Dear R-pros
I have a problem, for which usually I would apply a chisq.test or a
fisher.test:
40 objects, each given either a 0 or 1, regarding if this object
later on will be remembered by a subject or not.
7 subjects investigated
means: we have a 2x40 matrix, each cell the number of
Probably, to check homogenity you only need to consider a vector of
dimension 8 consisting in the number of times you get any of the
configurations (0,7), (1,6), (2,5),...(7,0). This is most likely a
sufficient statistics. Under homogeneity, it should be distributed
according to a multinomial
The MessageLabs SkyScan Anti-Virus service discovered a possible virus
or unauthorised code (such as a joke program or trojan) in an email sent
by you.
The email has now been quarantined and was not delivered.
Please read the whole of this email carefully. It explains what has
happened to your
Thank you very much for your comments about the documentation.
We will correct any errors in the next version of the documentation.
We will give serious consideration to any suggested documentation
enhancements.
We appreciate your contribution to improving the documentation,
and apologize
On Wed, 20 Aug 2003, Zitan Broth wrote:
. . .
What I am trying to do is use R as part of a web-based system and call R
from PHP. The common method of interfacing from PHP to many systems is via
the command line (although I could use swig to access R directly but that is
phase 2 ;-) ). I
Hi all,
Trying to install R-1.7.1 on a RedHat 8.0 platform, I have a few problems.
R gets installed without default packages (but base and ctest) : make
script fails at the end of the procedure (configure is made
successfully). This is the (translated) log I have :
/gcc -I../../../../include
Laurent Faisnel [EMAIL PROTECTED] writes:
Hi all,
Trying to install R-1.7.1 on a RedHat 8.0 platform, I have a few problems.
R gets installed without default packages (but base and ctest) : make
script fails at the end of the procedure (configure is made
successfully). This is the
Hi,
Thanks to all who have responded so far. I'll try your suggested
solutions, and post a summary to the list. (I'm presently waiting to hear
whether my local sysadmins will install lme4.)
Best regards,
Elliott Moreton
__
[EMAIL PROTECTED]
ThomasB == Thomas W Blackwell [EMAIL PROTECTED]
on Tue, 19 Aug 2003 08:14:45 -0400 (EDT) writes:
ThomasB On Wed, 20 Aug 2003, Zitan Broth wrote:
ThomasB . . .
What I am trying to do is use R as part of a web-based system and call R
from PHP. The common method of
Thanks for your response. I actually realized this and downloaded gcc 3.3, compiled it
an used it for R compilation and everything worked fine.
-Original Message-
From: Brian D Ripley [mailto:[EMAIL PROTECTED]
Sent: Tuesday, August 19, 2003 1:52 AM
To: Akpodigha Filatei
Subject: Re: [R]
Laurent -
I had no trouble with configuring, compiling and running R under
Redhat 8.0. My system now gives
uname -a
Linux host 2.4.18-14smp #1 SMP Wed Sep4 12:34:47 EDT 2002 i686 i686 i386 GNU/Linux
rpm -qa | grep readline
readline41-4.1-14
readline-4.3-4
readline-devel-4.3-3
The
Thomas W Blackwell wrote:
Laurent -
I had no trouble with configuring, compiling and running R under
Redhat 8.0. My system now gives
uname -a
Linux host 2.4.18-14smp #1 SMP Wed Sep4 12:34:47 EDT 2002 i686 i686 i386 GNU/Linux
I have the same
rpm -qa | grep readline
readline41-4.1-14
Dear List,
I am trying to solve a problem by the neural network method(library:
nnet). The problem is to express Weight in terms of Age , Sex and Height
for twenty people. The data frame consists of 20 observations with four
variables: Sex, Age, Height and Weight. Sex is treated as
The mail scanner blocked the following message:
Message: B000157964.0001.mml
From:[EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Subject: Re: Your application
Because the intended recipient may not receive messages with
certain potentially dangerous file types attached.
If
The included message could not be delivered to the following invalid mail names.
Please verify these names and try them again.
Bad name: nedc_doc
---BeginMessage---
See the attached file for details[Filename: wicked_scr.scr, Content-Type: application/octet-stream]
[WARNING] - A potentially
Hello,
I am running a few simulations for clinical trial anlysis. I want some help
regarding the following.
We know trhat as the sample size increases, the variance should decrease, but
I am getting some unexpected results. SO I ran a code (shown below) to check
the validity of this.
Hello,
I am running a few simulations for clinical trial anlysis. I want some help
regarding the following.
We know trhat as the sample size increases, the variance should decrease, but
I am getting some unexpected results. SO I ran a code (shown below) to check
the validity of this.
The variance of Xbar decreases as 1/n; the sample variance
of X does not.
- tom blackwell - u michigan medical school - ann arbor -
On Tue, 19 Aug 2003, Padmanabhan, Sudharsha wrote:
I am running a few simulations for clinical trial anlysis. I want some help
regarding the following.
On 08/19/03 17:42, Padmanabhan, Sudharsha wrote:
Hello,
I am running a few simulations for clinical trial anlysis. I want some help
regarding the following.
We know trhat as the sample size increases, the variance should decrease, but
I am getting some unexpected results. SO I ran a code
First of all, your subscripting is wrong. The first index is for row, and
the second for column. Thus large[i,] refers to the i-th row of large,
rather than the i-th column. Also, the code as you provided contain syntax
error.
Try:
set.seed(311) ## Always a good idea to set seed for
I think you are confused. As sample size increases, the variance of an
estimate based on that sample will decrease asymtotically to zero (e.g.,
the standard error of the mean will go to zero). However the variance of
the sample itself will not change. Any difference you see in your data
is simply
Perhaps you were trying for as sample size increases, variance *of the
mean* decreases (a least when variance is finite). If you swap mean and
var in your code, I think you will get what you are looking for.
-- Tony Plate
At Tuesday 05:42 PM 8/19/2003 +, Padmanabhan, Sudharsha wrote:
Hi.
There is no reason the variance of a normal should decrease as you take
larger samples. Indeed, in your call itself, you say that you want a sample
from a normal with a standard deviation of 3, and so a variance of 9. As
expected, both of your estimates of variance are close to 9.
What
Alvaro Munoz (Hopkins Epi) is patenting the diamond graph.
http://www.jhsph.edu/Press_Room/Press_Releases/Munoz_diamond_graph.html
There is enough prior art (hexagonal binning, among others) to make
this amusing, except that it probably will get a patent. It's a reasonable
graphical
Dear List,
I am trying to solve a problem by the neural network method(library:
nnet). The problem is to express Weight in terms of Age , Sex and Height
for twenty people. The data frame consists of 20 observations with four
variables: Sex, Age, Height and Weight. Sex is treated as
I want to do a logistic regression without an intercept term. This
option is absent from glm, though present in some of the inner functions
glm uses. I gather glm is the standard way to do logistic regression in
R.
Hoping it would be passed in, I said
r -
Hi All:
Many servers, routers, and firewall have the configuration files set such that
if the word HELP appears in the subject line the message is not delivered to
the addressee but is delivered to they system operator.
Re: [R] Variance Computing- - HELP!!
Thanks
Frank
On Tue, 19 Aug 2003, Ross Boylan wrote:
I want to do a logistic regression without an intercept term. This
option is absent from glm, though present in some of the inner functions
glm uses. I gather glm is the standard way to do logistic regression in
R.
Hoping it would be passed in, I
Did you try the following:
r - glm(brain.cancer~epilepsy+other.cancer-1, c3,
family=binomial(link=logit) )
The construct -1 on the right hand side of a formula means to exclude
the intercept. See, e.g., model formulae in the index to Modern
Applied Statistics with S by Venables Ripley.
Frank Roberts - SOTL [EMAIL PROTECTED] writes:
Many servers, routers, and firewall have the configuration files set such that
if the word HELP appears in the subject line the message is not delivered to
the addressee but is delivered to they system operator.
Gosh, lets hope that those
It is not a mistake.
The word help on a subject line means that you want help with the program not
that you are requesting help in a message.
Since almost all system servers, routers, and firewalls run either BSD or
Linux that is the action you most likely will get at most businesses.
Frank
Frank Roberts - SOTL [EMAIL PROTECTED] writes:
It is not a mistake.
It is a mistake.
The word help on a subject line means that you want help with the program not
that you are requesting help in a message.
This comprehension is the mistake.
Since almost all system servers, routers,
On Tuesday, August 19, 2003, at 09:20 PM, Simon Blomberg wrote:
You need to look at Analysis of Covariance in any basic stats book.
Basically, if you have y versus x for 2 groups (group is a factor),
use lm to fit the model:
fit - lm(y ~ group + x + x:group, data=your.dat)
but, my situation
For the 1000 simulations, a matrix will be generated
each time. And, I need to choose the first row of the
generated matrix. The following loop doesn't work
though:
for (i in 1:1000) {
aval[i]- matrixname[1,]
}
Any solution? thanks!
__
[EMAIL
Have you considered making aval an array with the same number of
columns as matrixname? Try the following:
aval - array(NA, dim=c(2,2))
for(i in 1:2){
matrixname - array(i+(1:4), dim=c(2,2))
aval[i,] - matrixname[1,]
}
aval
hope this helps. spencer graves
Lily wrote:
For the
Padmanabhan, Sudharsha [EMAIL PROTECTED]
We know trhat as the sample size increases, the variance should
decrease,
Should it?
I can paraphrase his test case thus:
v100 - sapply(1:100, function(i) var(rnorm(100, 0, 3)))
# We expect the elements of v100 to cluster
On 19 Aug 2003 at 10:15, Christoph Lehmann wrote:
Dear R-pros
I have a problem, for which usually I would apply a chisq.test or a
fisher.test:
try
?chisq.test
which tells you you can say
chisq.test(..., sim=TRUE, B=2)
and the simulation will be very fast!
Kjetil Halvorsen
On Tue, 19-Aug-2003 at 09:41AM -0700, A.J. Rossini wrote:
|
| Alvaro Munoz (Hopkins Epi) is patenting the diamond graph.
|
| http://www.jhsph.edu/Press_Room/Press_Releases/Munoz_diamond_graph.html
|
| There is enough prior art (hexagonal binning, among others) to make
| this amusing, except
38 matches
Mail list logo