[R-sig-ME] lme:correlations tructure AR1 and random factor
Yvonne Fabian
yvo.fabian at gmx.de
Wed Mar 30 17:26:52 CEST 2011
Dear mixed models specialists,
I tried these model to run in the package nlme, but allways got the same
error message...
I have a timeseries correlation in 5 sessions within 12 fields with 12 traps per field.
res2a <- lme(response~x+y+z+treatment,correlation =
corCAR1(form = ~ session|trapfield), random=~1|field, na.action=na.omit,
data=plates, method="ML")
Error: Incompatible formulas for groups in "random" and "correlation"
I do not want to use trapfield as a random factor and using only session|field in the correlation term leaves me with a errormessage that there have to be unique values.
Do you have any suggestions what to do?
thanks a lot!
Yvonne
--
Yvonne Fabian
Département de biologie
Unité d'écologie & évolution
Université de Fribourg
Chemin du Musée 10
Ch-1700 Fribourg
Switzerland
e-mail: yvo.fabian at gmx.de
Mobile Phone: 0041-76-2744664
Home 0041-26-4361312
Work: 0041-26-3008849
More information about the R-sig-mixed-models
mailing list