在MCMCglmm中设置先验

如何解决在MCMCglmm中设置先验

我正在尝试通过使用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