[R] summary.lme and anova question

Christoph Scherber Christoph.Scherber at agr.uni-goettingen.de
Thu Aug 21 10:02:41 CEST 2008


Dear all,

When analyzing data from a climate change experiment using linear mixed-effects models, I recently
came across a situation where:

- the summary(model) showed a significant difference between the levels of a two-level factor,
- while the anova(model) showed no significance for that factor (see below).

My question now is: Is the anova.lme() approach correct for that model? And why does the F-test for 
CO2 yield a non-significant P-value, while the t-test in the summary.lme() is significant?

Many thanks for your help!

Best wishes
Christoph

######################################################

mod11=lme(log(ind1+1) ~ CO2*DROUGHT*TEMP, random=~1|B/C,na.action=na.exclude)

summary(mod11)
Linear mixed-effects model fit by REML
Data: NULL
      AIC      BIC    logLik
  97.3077 115.6069 -37.65385

Random effects:
Formula: ~1 | B
         (Intercept)
StdDev: 1.303146e-05

Formula: ~1 | C %in% B
        (Intercept)  Residual
StdDev:   0.2466839 0.4846578

Fixed effects: log(ind1 + 1) ~ CO2 * DROUGHT * TEMP
                      Value Std.Error DF   t-value p-value
(Intercept)       1.9981490 0.2220158 29  9.000030  0.0000
CO2              -1.0308687 0.3139778  5 -3.283254  0.0219
DROUGHT          -0.9715216 0.2798173 29 -3.471986  0.0016
TEMP             -0.5592615 0.2954130 29 -1.893151  0.0684
CO2:DROUGHT       1.2196261 0.3957214 29  3.082032  0.0045
CO2:TEMP          0.9791044 0.4068987 29  2.406261  0.0227
DROUGHT:TEMP      0.6413038 0.4068987 29  1.576077  0.1259
CO2:DROUGHT:TEMP -1.1448624 0.5675932 29 -2.017047  0.0530
Correlation:
                 (Intr) CO2    DROUGHT TEMP   CO2:DROUGHT CO2:TE DROUGHT:
CO2              -0.707                                                 DROUGHT          -0.630 
0.446                                          TEMP             -0.597  0.422  0.474 
                    CO2:DROUGHT       0.446 -0.630 -0.707  -0.335                           CO2:TEMP 
          0.433 -0.613 -0.344  -0.726  0.486                    DROUGHT:TEMP      0.433 -0.306 
-0.688  -0.726  0.486       0.527        CO2:DROUGHT:TEMP -0.311  0.439  0.493   0.520 -0.697 
-0.717 -0.717
Standardized Within-Group Residuals:
       Min         Q1        Med         Q3        Max
-1.4631313 -0.5715171 -0.2024273  0.4592221  1.9568914

Number of Observations: 47
Number of Groups:
       B C %in% B
       6       12

######################################################

anova(mod11)
                 numDF denDF   F-value p-value
(Intercept)          1    29 162.95719  <.0001
CO2                  1     5   1.15108  0.3324
DROUGHT              1    29   5.53240  0.0257
TEMP                 1    29   0.04519  0.8331
CO2:DROUGHT          1    29   5.66686  0.0241
CO2:TEMP             1    29   1.88455  0.1803
DROUGHT:TEMP         1    29   0.03481  0.8533
CO2:DROUGHT:TEMP     1    29   4.06848  0.0530


######################################################



More information about the R-help mailing list