全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3707 3
2014-03-30
悬赏 50 个论坛币 未解决
小弟最近学习CMCglmm, 在设置随机变量时遇到问题请打大侠帮忙,数据在附件中
11.xls
大小:(33.5 KB)

 马上下载


程序如下:
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

哪位大侠指点一下,感激不尽,着急呀!

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2014-3-31 12:35:29
有人指点一下吗?
急求帮助
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-3-31 22:48:01
设定一个先验分布参考https://stat.ethz.ch/pipermail/r-sig-mixed-models/2011q1/005208.html

prior.1<-list(R = list(V = diag(2), n = 2, fix = 2), G = list(G1 =list(V = diag(c(1, 1e-6)), n = 2, fix = 2)))
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-4-1 07:57:02
dongshengjay 发表于 2014-3-31 22:48
设定一个先验分布参考https://stat.ethz.ch/pipermail/r-sig-mixed-models/2011q1/005208.html

prior.1
董哥,怎么才能把这个做先验信息
mu.beta<-c(-2.037,2.415)
V.betas<-matrix(c(0.35716566, -0.08127909,-0.08127909 , 0.03360510),nrow=2)

n.iter<-25000
算出来,我还是没做出来呀
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群