[R-sig-ME] prior error in mcmcglmm multi-response model

[Victor Oribamise]

Hello,

I tried running this multi-response model in mcmcglmm with this prior:

priorsd <- list(R=list(V=diag(8), nu=7.001), G = list(G1 =
list(V=diag(8),nu=7.001), G2 = list(V=diag(8),nu=7.001)))

and this model:

modelTrait<-MCMCglmm(cbind(Trait1,Trait2, Trait3,Trait4,Trait5, Trait6,
Trait7, Trait8)~year+pit,random= ~ID+Dam ,ginverse=list(ID=Ainv1,
Dam=M1),family=rep("gaussian",8),prior=priorsd, , rcov=~units,
data=data,nitt=10000,burnin=1000,thin=10)


but I ran into this error:


"Error in priorformat(if (NOpriorG) { :V is the wrong dimension for some
prior$G/prior$R elements."


What am I doing wrong? How do I correct this?


And, this is a subset of my data:


Trait1 Trait2 Trait3 Trait4 Trait5 Trait6 Trait7 Trait8 year pit
1 1 1 1 1 1 1 1 1954 1
1 1 1 1 1 1 1 1 1956 2
1 1 1 1 1 1 1 1 1957 2
1 1 1 1 1 1 1 1 1957 1
1 1 1 1 1 1 1 0 1978 2
1 1 1 1 1 1 1 1 1978 2

Thanks.


Victor

	[[alternative HTML version deleted]]