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...