Hi all.

As far as I can tell newdata is still not implemented for this nice
package. Thus I wonder what would be the best way to get predictions "by
hand". My model is actually very simple. Still I need to marginalize the
random effects. Any hints? Thanks in advance, Antonio Pedro.


glm.MC.2 <- MCMCglmm(mortality.under.2 ~ maternal_age_c +
I(maternal_age_c^2)  +
                       as.factor(birth_year) + residence +
                       sex + wealth,
                     nitt=20000, thin=10, burnin=1000,
                     random= ~CASEID, prior=prior.2,data=egypt2,
family='categorical')

	[[alternative HTML version deleted]]

