小弟最近学习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

哪位大侠指点一下，感激不尽，着急呀！