QGglmm multivariate with fixed effects

Hello everyone,

I am working with the QGglmm package by Pierre de Villemereuil to 
calculate heritability and additive variance on the observable scale 
since my response variable is binomial and I am using MCMCglmm. I read 
his 'How to' pdf but it does not precisely say if the output represents 
the posterior mode or the mean for each parameter. Is it possible to 
compute the mode of the additive genetic variance and heritability 
instead of the mean if that is the case.

In addition, in his tutorial, he explains how to calculate the 
parameters if there are fixed effects in the model by using predict 
instead of calculating the mean of the trait. However that is only in 
the case of a univariate model, I was not successful in applying it for 
a multivariate model with fixed effects. I would appreciate it if I 
could get some help on how to proceed in this case.

Thank you

Walid Mawass

Ph.D. candidate in Cellular and Molecular Biology

Population Genetics Laboratory

University of Qu?bec at Trois-Rivi?res
3351, boul. des Forges, C.P. 500
Trois-Rivi?res (Qu?bec) G9A 5H7
Telephone: 819-376-5011 poste 3384