小弟最近学习CMCglmm, 在设置随机变量时遇到问题请打大侠帮忙,数据在附件中
程序如下:
library(coda)
library(MCMCpack)
library(MCMCglmm)
rd<-read.csv("11.csv")
head(rd)
attach(rd)
mu.beta<-c(-2.037,2.415)
V.betas<-matrix(c(0.35716566, -0.08127909,-0.08127909 , 0.03360510),nrow=2)
n.iter<-25000
prior.betas <- list(B = list(mu = mu.beta, V = V.betas),
R = list(V = 1, nu = 0.001))
model1.sim <- MCMCglmm(TOTAL~D, random=NULL, rcov=~units, family="gaussian",
data=rd, prior=prior.betas ,nitt=n.iter,burnin=floor(n.iter/2),
thin=3,verbose=F)
#在不设置随机变量时,model1.simk可以做出来
model2.sim <- MCMCglmm(TOTAL~D, random = ~ region,, rcov=~units, family="gaussian",
data=rd, prior=prior.betas ,nitt=n.iter,burnin=floor(n.iter/2),
thin=3,verbose=F)
#在设置随机变量时,出现:
Error in MCMCglmm(TOTAL ~ D, random = ~D | region, , rcov = ~units, :
either both or neither R and G structures need a prior
哪位大侠指点一下,感激不尽,着急呀!