Displaying 1 result from an estimated 1 matches for "batch_ij".
2004 Nov 16
1
lme, two random effects, poisson distribution
..._9, slugs_10}
Batch_4 = {slugs_11}
.
.
.
Batch_12 = {slugs_24, slugs_25}
Notice that there are 12 batches, and the number of slugs in each
batch differ, from 1 slug to 4 slugs.
I consider 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...