Displaying 1 result from an estimated 1 matches for "mcmc2".
Did you mean:
mcmc
2007 Mar 30
0
problem using mcmcsamp() with glmer models containing interaction terms in fixed effects
...data1)
# two lmer models
lmer1 <- lmer (counts~visit+visit:zone + (1|region), data=data1,
family=poisson(link="log"))
lmer2 <- lmer (counts~visit+visit:zone + (visit+visit:zone|region),
data=data1, family=poisson(link="log"))
mcmc1 <- mcmcsamp(lmer1, n=10000) # works
mcmc2 <- mcmcsamp(lmer2, n=10000) # doesn't work
[[alternative HTML version deleted]]