[R] mixed-effects model with two fixed effects: interaction
Ilona Leyer
ileyer at yahoo.de
Tue Jun 29 11:09:33 CEST 2010
Dear all,
In a greenhouse experiment we tested performance of 4 different species (B,H,P,R) under 3 different water levels in 10 replications. As response variable e.g. the number of emerging sprouts were measured on three dates. A simple Anova considering every measurement date separately shows a higly significant effect of species and moisture (and partly the interaction of both). The mixed-effects model with species and moisture shows a highly significant effect of species and moisture as well. However, when I included the interaction the t-values of the species dropped strongly and the SE increase and the results for the species are not significant anymore. For me this does not seem plausible. Has anybody an idea, how this can be interpreted and if I have done a mistake in calculating the data?
Thanks in advance for any help!
Ilona
model1<-lme(sprouts~species+moisture,random=~time|ID)
model2<-lme(sprouts~species*moisture,random=~time|ID)
Fixed effects: sprouts ~ species + moisture
Value Std.Error DF t-value p-value
(Intercept) 7.971267 1.330500 240 5.991180 0.0000
speciesH -6.459344 1.536329 114 -4.204400 0.0001
speciesP -10.063604 1.536329 114 -6.550421 0.0000
speciesR -5.051894 1.536329 114 -3.288288 0.0013
moisturemoist 2.228835 1.330500 114 1.675185 0.0966
moisturewaterlogged 17.111149 1.330500 114 12.860688 0.0000
Fixed effects: sprouts ~ species * moisture
Value Std.Error DF t-value p-value
(Intercept) 4.831965 1.750970 240 2.759594 0.0062
speciesH -4.464197 2.476245 108 -1.802809 0.0742
speciesP -3.986787 2.476245 108 -1.610013 0.1103
speciesR -0.809376 2.476245 108 -0.326856 0.7444
moisturemoist 3.505506 2.476245 108 1.415654 0.1598
moisturewaterlogged 24.766934 2.476245 108 10.001811 0.0000
speciesH:moisturemoist -0.457291 3.501939 108 -0.130582 0.8963
speciesP:moisturemoist -2.458125 3.501939 108 -0.701932 0.4842
speciesR:moisturemoist -2.555356 3.501939 108 -0.729697 0.4672
speciesH:moisturewaterlogged -5.597498 3.501939 108 -1.598400 0.1129
speciesP:moisturewaterlogged -15.538272 3.501939 108 -4.437048 0.0000
speciesR:moisturewaterlogged -10.206874 3.501939 108 -2.914635 0.0043
More information about the R-help
mailing list