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: edeomano@bses.org.au<blocked::mailto:jstringer@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}}


