R MCMCglmm中的先验设置
我试图通过在R中使用“MCMCglmm”包来拟合线性混合效果模型,但设置优先级时遇到问题。我是MCMCglmm软件包的新手,恐怕我不懂符号,也不知道如何正确设置优先级 这是我的密码:R MCMCglmm中的先验设置,r,mcmc,R,Mcmc,我试图通过在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
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)
我的印象是,对于G,V应该=3,因为有随机项intercept,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
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