Good day,
I'm ran aov () and lme() on a split-plot using a 64-bit machine. For aov() I
don't see the values for ErrorA, F-value and p-value in the output. For lme(),
output is different from results from a 32-bit. Please see codes used and
corresponding output. Is my code wrong and/or not sufficient or is this a
compatibility issue?
**************
model1<-aov(Y~Main*Sub + Error(Block/Main),data=splitplot)
anova(model1)
Output:
Error: Block
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 1 4.9116 4.9116
Error: Block:Main
Df Sum Sq Mean Sq
Main 5 27.4793 5.4959
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Main 5 14.35 2.87 1.2704 0.2757
Sub 2 3.89 1.95 0.8620 0.4231
Main:Sub 10 14.77 1.48 0.6541 0.7669
Residuals 408 921.43 2.26
*****************
model1 <- lme(y~Main*Sub, random=~1|Block/Sub,data=splitplot)
anova(model1)
Output:
numDF denDF F-value p-value
(Intercept) 1 405 150.84639 <.0001
Main 5 405 3.74269 0.0025
Sub 2 6 0.55782 0.5995
Main:Sub 10 405 0.69772 0.7268
Regards,
Emily
Emily Deomano
Biometrician
BSES Limited
Postal Address: PO Box 86
Street Address: 50 Meiers Road
Indooroopilly QLD 4068 Australia
Ph: +61 (0) 7 3331 3306
Mobile: +61 (0) 408 656 452
Fax: +61 (0) 7 38710383
Email: edeom...@bses.org.au<blocked::mailto:jstrin...@bses.org.au>
Web: www.bses.org.au<http://www.bses.org.au>
BSES Limited Disclaimer
--------------------------------------------------------------------------------
This email and any files transmitted with it are confide...{{dropped:15}}
______________________________________________
R-help@r-project.org 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.
