search for: multinomial2

...here the intercept, var1 and var2 are jointly distributed according to study. My question is now how I would translate this formula to the fixed and random terms in MCMCglmm. For the fixed part, I figured that I should make a variable nooutcome=abs(1-outcome) because it can then be modeled with a multinomial2 family as there is no binomial(logit) option available. Then, the fixed part would look as follows: cbind(outcome,nooutcome)~1+var1+var2 However, I am unsure how to specify the random effects over the intercept, var1 and var2 jointly. So far, I was able to generate the following: random=~us(var...

Hi, I have fitted a generalized linear mixed effects model using lmer (library lme4), and the family = quasibinomial. I have tried to obtain a MCMC sample, but on calling mcmcsamp(model1, 1000) I get the following error which I don't understand at all: Error in .local(object, n, verbose, ...) : Update not yet written traceback() delivers: 4: .Call(mer_MCMCsamp, ans, object) 3: