问题描述
我正在尝试通过使用R中的“ MCMCglmm”包来拟合线性混合效果模型,但无法设置先验值。我是MCMCglmm软件包的新手,恐怕我既不了解该符号,也不了解如何正确设置先验。
这是我的代码:
mydata$XNEW<-(mydata$X-mean(mydata$X))/(2*sd(mydata$XNEW))
mydata$YNEW<-scale(mydata$Y)
prior1 <- list(G = list(G1 = list(V = 3,nu = 0.002)),R = list(V = 1,nu = 0.002))
model1 <-MCMCglmm(YNEW~XNEW,random=~us(1+XNEW):ID,data=mydata,prior = prior1)
Error in priorformat(if (nopriorG) { :
V is the wrong dimension for some prior$G/prior$R elements
我的印象是,对于G,V应该= 3,因为存在随机项拦截,ID和Temperaturez,尽管我可能完全误解了所读内容。请有人可以提供白痴指南来设定先验条件吗?
以下是一些数据:
ID Y X
A1.2.2 25 3
A1.2.2 7 3
A1.2.2 43 10
A1.2.2 20 10
A1.2.2 66 17
A1.2.2 54 17
A1.2.2 72 23
A1.2.2 31 23
A4.3.1 19 3
A4.3.1 13 3
A4.3.1 44 10
A4.3.1 33 10
A4.3.1 59 17
A4.3.1 26 17
A4.3.1 88 23
A4.3.1 49 23
A4.3.2 11 3
A4.3.2 16 3
A4.3.2 17 10
A4.3.2 11 10
A4.3.2 41 17
A4.3.2 37 17
A4.3.2 67 23
A4.3.2 2 23
A7.10.1 24 3
A7.10.1 15 3
A7.10.1 38 10
A7.10.1 23 10
A7.10.1 70 17
A7.10.1 49 17
A7.10.1 47 23
A7.10.1 47 23
A7.3.1 11 3
A7.3.1 4 3
A7.3.1 23 10
A7.3.1 6 10
A7.3.1 20 17
A7.3.1 27 17
A7.3.1 72 23
A7.3.1 47 23
A8.5.1 12 3
A8.5.1 12 3
A8.5.1 15 10
A8.5.1 19 10
A8.5.1 57 17
A8.5.1 25 17
A8.5.1 47 23
A8.5.1 37 23
A9.7.1 8 3
A9.7.1 23 3
A9.7.1 12 10
A9.7.1 20 10
A9.7.1 60 17
A9.7.1 42 17
A9.7.1 61 23
A9.7.1 15 23
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)