[R-br] Definição de Prior do modelo

[Andre Oliveira]

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é 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200306/c7255695/attachment.html>