[R] test for equal distributions with small numbes of observations

2003-08-19 Thread Christoph Lehmann
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

Re: [R] test for equal distributions with small numbes of observations

2003-08-19 Thread Carlos J. Gil Bellosta
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

[R] WARNING. You tried to send a potential virus or unauthorisedcode

2003-08-19 Thread alert
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

[R] Re: Thank you!

2003-08-19 Thread The MathWorks Inc. Documentation
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

Re: [R] Command line R / PHP?

2003-08-19 Thread Thomas W Blackwell
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

[R] R-1.7.1 gets installed without default packages withoutreadline

2003-08-19 Thread Laurent Faisnel
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

Re: [R] R-1.7.1 gets installed without default packages withoutreadline

2003-08-19 Thread Peter Dalgaard BSA
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

[R] Re: glmmPQL() and memory limitations

2003-08-19 Thread Elliott Moreton
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]

Re: [R] Command line R / PHP?

2003-08-19 Thread Martin Maechler
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

RE: [R] R Install on Solaris 9

2003-08-19 Thread Akpodigha Filatei
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]

Re: [R] R-1.7.1 gets installed without default packages withoutreadline

2003-08-19 Thread Thomas W Blackwell
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

Re: [R] R-1.7.1 gets installed without default packages withoutreadline

2003-08-19 Thread Laurent Faisnel
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

[R] On the Use of the nnet Library

2003-08-19 Thread yukihiro ishii
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

[R] Mgw: Blocked mail Re: Your application fromr-help@lists.r-project.org

2003-08-19 Thread mailmarshal
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

[R] Return Message: Thank you!

2003-08-19 Thread ORAPOST
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

[R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Padmanabhan, Sudharsha
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.

[R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Padmanabhan, Sudharsha
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.

Re: [R] Variance Computing - HELP

2003-08-19 Thread Thomas W Blackwell
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.

Re: [R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Jonathan Baron
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

RE: [R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Liaw, Andy
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

Re: [R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread James MacDonald
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

Re: [R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Tony Plate
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:

RE: [R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Wiener, Matthew
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

[R] for those of you who want a patent...

2003-08-19 Thread A.J. Rossini
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

[R] On the Use of the nnet Library

2003-08-19 Thread yukihiro ishii
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

[R] logistic regression without intercept

2003-08-19 Thread Ross Boylan
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 -

Re: [R] Variance Computing-

2003-08-19 Thread Frank Roberts - SOTL
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

Re: [R] logistic regression without intercept

2003-08-19 Thread Thomas Lumley
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

Re: [R] logistic regression without intercept

2003-08-19 Thread Spencer Graves
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.

Re: [R] Variance Computing-

2003-08-19 Thread A.J. Rossini
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

Re: [R] Variance Computing-

2003-08-19 Thread Frank Roberts - SOTL
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

Re: [R] Variance Computing-

2003-08-19 Thread A.J. Rossini
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,

Re: [R] compare two linear regression slopes

2003-08-19 Thread John Christie
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

[R] question for loop on matrix row level.

2003-08-19 Thread Lily
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

Re: [R] question for loop on matrix row level.

2003-08-19 Thread Spencer Graves
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

Re: [R] Variance Computing- - HELP!!!!!!!!!!!!!!!!!!

2003-08-19 Thread Richard A. O'Keefe
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

Re: [R] test for equal distributions with small numbes of observations

2003-08-19 Thread kjetil brinchmann halvorsen
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

Re: [R] for those of you who want a patent...

2003-08-19 Thread Patrick Connolly
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