www.cmm.bristol.ac.uk/forum

Hello everyone,

I want to estimate a multilevel mixed type model: y1 is Normal and y2 is Poisson.

I’m using the following code:

model= c(y1, log(y2)) ~ 1 + x1 + (1 | group) + (1[1] | idn )
mymodel3 <- runMLwiN( model, D = c("Mixed", "Normal", "Poisson"),
estoptions = list(EstM = 0, resi.store ...