[R] lme outer
yclough at gwdg.de
Tue Jul 19 15:15:31 CEST 2005
Dear R users
a question about "outer" explanatory variables in lme:
I have measured size of a population of insects in fields.
These fields were spread out over a large region.
The fields are grouped (spatially) in pairs: one with fertiliser "high",
the other one "low".
I want to test effect of mean temperature and fertiliser on popsize.
Meantemp was measured for each field, but measurements are correlated
and this should be taken into account to avoid pseudoreplication (in
other words, I d like meantemp to be considered an "outer" variable).
Do I need to replace the temperature measurements by the means for each
Or can I leave in the measurements per field pair?
this is my data and model (with meantemp values for each field):
I now create a vector with the values of meantemp averaged per pair
rerun a model with that explanatory variable:
summary.lme and anova.lme suggest minute differences in the estimated
parameters and DF (!) between model1 and model2.
How do I explain these differences, especially in the DF?
Is there a model to prefer?
email: yclough at gwdg.de
More information about the R-help