[R-sig-ME] genetic effects and multiple membership in MCMCglmm (Alexandre Martin)

[Szymek Drobniak]

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]]