Definição de Prior do modelo

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é
participantes (1)
-
Andre Oliveira