[R] an error message with 't.test' with R under Unix

Tao Shi shitao at hotmail.com
Thu Sep 14 00:34:35 CEST 2006


Hi Peter,

Thank you for the reply!

Really sorry for the formating (forgot to change to plain text).  Here is basically what I'm talking about:

> R.version
         _                       
platform x86_64-unknown-linux-gnu
arch     x86_64                  
os       linux-gnu               
system   x86_64, linux-gnu       
status                           
major    2                       
minor    1.1                     
year     2005                    
month    06                      
day      20                      
language R                       
> t.test(1:6~rep(1:2,each=3))
Error in terms.formula(formula[-3]) : invalid model formula in ExtractVars
> 

Again, I don't see the error in R under Windows (R 2.3.0) or Unix (R2.3.1).  Is this the bug you were talking about?

...Tao




----------------------------------------
> To: shitao at hotmail.com
> CC: r-help at stat.math.ethz.ch
> Subject: Re: [R] an error message with 't.test' with R under Unix
> From: p.dalgaard at biostat.ku.dk
> Date: Wed, 13 Sep 2006 23:52:04 +0200
> 
> "Tao Shi" <shitao at hotmail.com> writes:
> 
> > Hi list,Could you please help me to explain the following error
>  messages with 't.test' in R Unix 2.1.1? 
> 
> This is completely unreadable! However, yes, there was at some point a
> bug where the LHS of model formulas was checked more rigorously than
> need be, using the same rules that apply to the RHS. This fact can be
> found in a recent NEWS file.
> 
> 
> 
> > I don't see it in R under Windows (R 2.3.0) or Unix (R2.3.1).  Is it really due to the different R versions?Thanks,...TaoUnix session: (R.2.1.1)========================> R.version         _                       platform x86_64-unknown-linux-gnuarch     x86_64                  os       linux-gnu               system   x86_64, linux-gnu       status                           major    2                       minor    1.1                     year     2005                    month    06                      day      20                      language R                       > t.test(extra ~ group, data = sleep)        Welch Two Sample t-testdata:  extra by group t = -1.8608, df = 17.776, p-value = 0.0794alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -3.3654832  0.2054832 sample estimates:mean in group 1 mean in group 2        !
> >      0.75            2.33 > t.test((1:6)~rep(1:2,each=3))Error in terms.formula(formula[-3]) : invalid model formula in ExtractVars> t.test(2^(1:6)~rep(1:2,each=3))Error in terms.formula(formula[-3]) : invalid power in formula> t.test(2^extra ~ group, data = sleep)Error in terms.formula(formula[-3]) : invalid power in formulaUnix session: (R 2.3.1)============================> R.version               _                         platform       x86_64-unknown-linux-gnu  arch           x86_64                    os             linux-gnu                 system         x86_64, linux-gnu         status                                   major          2                         minor          3.1                       year           2006                      month          06                        day            01                        svn rev        38247                     language       R                         version.string Version 2.3.1 (2006-06-01)> t.test(1:6~rep(1:2,each
>  =!
> >  3))        Welch Two Sample t-testdata:  1:6 by rep(1:2, each = 3) t =
> >  -3.6742, df = 4, p-value = 0.02131alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -5.266958 -0.733042 sample estimates:mean in group 1 mean in group 2               2               5 Windows session:===========================> R.version               _                         platform       i386-pc-mingw32           arch           i386                      os             mingw32                   system         i386, mingw32             status                                   major          2                         minor          3.0                       year           2006                      month          04                        day            24                        svn rev        37909                     language       R                         version.string Version 2.3.0 (2006-04-24)> help(t.test)> t.test(extra ~ group, data = sleep)        Welch Two Sample t-testdata:  extra by group t = -1.8608, df = 17.
>  7!
> >  76, p-value = 0.0794alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -3.3654832  0.2054832 sample estimates:mean in group 1 mean in group 2            0.75            2.33 > t.test(2^extra ~ group, data = sleep)        Welch Two Sample t-testdata:  2^extra by group t = -1.6362, df = 10.718, p-value = 0.1308alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -18.641557   2.773344 sample estimates:mean in group 1 mean in group 2        3.448644       11.382751 > t.test(1:6~rep(1:2,each=3))        Welch Two Sample t-testdata:  1:6 by rep(1:2, each = 3) t = -3.6742, df = 4, p-value = 0.02131alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -5.266958 -0.733042 sample estimates:mean in group 1 mean in group 2               2               5 
> > _________________________________________________________________
> > Use Messenger to talk to your IM friends, even those on Yahoo!
> > 
> > 685ee3e858fe
> > 	[[alternative HTML version deleted]]
> > 
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> > 
> 
> -- 
>    O__  ---- Peter Dalgaard             Øster Farimagsgade 5, Entr.B
>   c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
>  (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45) 35327918
> ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)                  FAX: (+45) 35327907

_________________________________________________________________
Use Messenger to talk to your IM friends, even those on Yahoo!



More information about the R-help mailing list