Displaying 1 result from an estimated 1 matches for "lmedatab".
Did you mean:
lmedata
2004 Nov 16
1
lme, two random effects, poisson distribution
...batch to be another random-effect on intercept. Thus I
fit model:
************pumps_tij = state_t + slugs_i + batch_ij + error_tij
Slugs, batch and error are normally distributed, but slugs and batch
are not nested factors.
I had fit following (however I'm not sure if that is right):
slugs.lmedataB <- groupedData(pumps ~ state | slugs/batch,
data=as.data.frame(data.slugs))
csB <- corAR1(form= ~ x|slugs/batch)
res1B <- lme(pumps ~ state, random = ~1 | slugs/batch,
data=slugs.lmedataB, cor=csB)
----------------------------------------
QUESTIONS:
1) Are my models right? Do I model the...