Displaying 1 result from an estimated 1 matches for "island2".
Did you mean:
island
2007 Mar 30
0
problem using mcmcsamp() with glmer models containing interaction terms in fixed effects
...it <- c(rep(0, 50), rep(1,50)) # sample visit number, 0 or 1 to indicate
pre and post no-fishing reserve establishment
zone <- as.factor(rep(c(rep("blue", 25), rep("green",25)),2)) #
inside/outside no-fishing reserve
region <- as.factor(rep(c("island1", "island2"),50)) # spatial grouping
factor
data1 <- data.frame(counts, visit, zone, region)
rm (counts, visit, zone, region)
# a quick look
bwplot(I(sqrt(counts))~as.factor(visit)|region, data=data1)
# two lmer models
lmer1 <- lmer (counts~visit+visit:zone + (1|region), data=data1,
family=poiss...