Pessoal, bom dia! estou com problemas de colocar priori para os efeitos fixos do modelo II. Alguém saberia onde estaria o erro?
rm(list=ls(all=TRUE)) # Modelo I data(PlodiaPO) model1<-MCMCglmm(PO ~ plate, random=~FSfamily, data=PlodiaPO, verbose=FALSE, nitt=1300, burnin=300, thin=1) summary(model1) # Modelo II prioriC0 <- list(B = list(B1=list(mu = 0, V = 1e+8), B2=list(mu = 0, V = 1e+8), R = list(V=1, nu=1), G = list(G1=list(V=1, nu=1, alpha.mu = 0, alpha.V = 25^2))))model2<-MCMCglmm(PO ~ plate, random=~FSfamily, data=PlodiaPO, verbose=FALSE, nitt=13000, burnin=300, thin=2, prior=prioriC0)Error in MCMCglmm(PO ~ plate, random = ~FSfamily, data = PlodiaPO, verbose = FALSE, : either both or neither R and G structures need a priorsummary(model2) att,.André
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
