[R-sig-ME] MCMCglmm Heterogenous Residual Variance in a Categorial Model

FERNANDEZ Ana maria s1049258 at sms.ed.ac.uk
Fri Nov 7 12:18:04 CET 2014

Dear List,

I was wondering what the recommendation is for structuring residual variances in the context of a categorial model when you have a sex interaction in the random term of the model (where animal is the pedigree structure).

MCMCglmm course notes have a brief recommendation about sex interaction and allowing the sexes to have different residual variances: "3.4 Heterogenous Residual Variance To be started... In short - if you’ve fitted a sex by dam interaction I would always allow the sexes to have different residual variances. Use rcov=∼idh(sex):units.”

But the conflict that I am coming across is another recommendation about how the residual term should be- or needs to be- fixed in a categorical model.

I have a case-control binary response variable and I want to look at the sex effect interaction.

Let me know if you need more detail.


Ana Maria Fernández Pujals
A.Fernandez-2 at sms.ed.ac.uk<mailto:A.Fernandez-2 at sms.ed.ac.uk>

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20141107/30b4525e/attachment.pl>

More information about the R-sig-mixed-models mailing list