search for: lmedata

Displaying 1 result from an estimated 1 matches for "lmedata".

Did you mean: medata
2004 Nov 16
1
lme, two random effects, poisson distribution
...ime slots (t=1,...,10) state_t = 1 for excited time slots (t=11,12) slugs_i - ith animal, where i = 1,...,25 I would like to find out if the # of pumps depends on the variable state, assuming the correlation AR(1) between x_t and slugs being a random-effect on intercept. slugs.lmedata <- groupedData(pumps ~ state | slugs, data=as.data.frame(data.slugs)) cs <- corAR1(form= ~ x|slugs) res1 <- lme(pumps ~ state, random = ~1 | slugs, data=slugs.lmedata, cor=cs) --------------------------------- Now, I would like to add a complication to the model: The slugs were observed...