[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