search for: trapfield

Dear helpers, I tried these models to run in the package nlme, but allways got the same error message... I have a correlation in 5 sessions within a field (n=12) with ten traps in one field. res2a <- lme(response~x+y+z+treatment),correlation = corARMA(form = ~ session|trapfield, p = 1, q = 0), random=~1|field, na.action=na.omit, data=plates, method="ML") 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 &...

...lysis. I have presence absence data of 5 sampling sessions within 3 years, of 12 fields. Each field had 12 traps. I would like to analyse the data with a Generalized Estimating Equations (GEE) Model in R. For the abundance data I used a gls with the function: correlation = corARMA(form = ~ session|trapfield, p = 1, q = 0) But now I want to use presence absence data?I know about the problem with correlated -binary data- but maybe somebody already created a solution? Something like that, which is not working?. rgee1 <- geeglm(Alus01 ~ treatment + scale(LAI) + scale(log(Tot_Sp)) + scale(AveH) , data...