Raw Message

Message-ID: <CANXb-o5zZBsrPa_6OmQhJcgLQpSxggydTeMhEcGKMFYhsRSWNQ@mail.gmail.com>
Date: 2015-05-14T21:38:45Z
From: Szymek Drobniak
Subject: genetic effects and multiple membership in MCMCglmm (Alexandre Martin)

Hello Alexandre,

using ped assigns a given correlation structure only to the "animal" random
effect. You have to create an inverse of A matrix:

my_inverse <- inverseA(ped)$Ainv

and assign it in MCMCglmm to specific random effects:

MCMCglmm(Y~1,
random=~animal+mm(m1+m2+m3),
ginverse=list(animal=my_inverse, m1=my_inverse,
m2=my_inverse,m3=my_inverse), data=dat, pr=T).

Cheers
Szymek

	[[alternative HTML version deleted]]