Skip to content
Prev 13337 / 20628 Next
Original Raw

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

Thread (2 messages)

Szymek Drobniak genetic effects and multiple membership in MCMCglmm (Alexandre Martin) May 14 Alexandre Martin genetic effects and multiple membership in MCMCglmm (Alexandre Martin) May 19