R 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

我试图通过在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)
           
我的印象是,对于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