Caterina, Did you get an answer to this question? I'm trying to do something similar. Jason -- View this message in context: http://r.789695.n4.nabble.com/Multilevel-pseudo-maximum-likelihood-tp878413p3213583.html Sent from the R help mailing list archive at Nabble.com.