Displaying 1 result from an estimated 1 matches for "expl1".
Did you mean:
exp1
2011 Feb 17
0
Multi-response MCMCglmm (gaussian and zapoisson)
...ables and consequently to specify completely parameterized covariance matrices.
My current model looks like this:
prior<-list(R=list(V=diag(4),nu=0.004),G=list(G1=list(V=diag(4)/4,n=4),G2=list(V=diag(4)/4,n=4),G3=list(V=diag(4)/4,n=4)))
m1<-MCMCglmm(fixed=cbind(resp1,resp2,resp3) ~
trait:(expl1 + expl2+ expl3+ expl4+ expl5...)+trait-1,
random=~us(trait):rand1+us(trait): rand1: rand2+us(trait): rand1: rand2: rand3,
rcov=~us(trait):units,family=c("gaussian","gaussian","zapoisson"),nitt=20000,burnin=5000,thin=10,prior=prior,data=dat)
However, for zero-altered m...