Question

我试图根据书＆＃34;做贝叶斯数据分析：R，JAGS和斯坦的教程（2015）＆＃34;来研究贝叶斯分析。

在本书中，有一些例子。所以，我试图在R中复制这个例子。但是，我在此示例中收到了错误消息。

具体而言，这是示例数据。

data
   y        s
1  1 Reginald
2  0 Reginald
3  1 Reginald
4  1 Reginald
5  1 Reginald
6  1 Reginald
7  1 Reginald
8  0 Reginald
9  0     Tony
10 0     Tony
11 1     Tony
12 0     Tony
13 0     Tony
14 1     Tony
15 0     Tony

y<-data$y
s<-as.numeric(data$s)
Ntotal=length(y)
Nsubj=length(unique(s))

dataList=list(y=y, s=s, Ntotal=Ntotal, Nsubj=Nsubj)

此外，这是我的模特。

modelString=" 
model{
  for(i in 1:Ntotal){
    y[i] ~ dbern(theta[s[i]])
  }
  for(s in 1:Nsubj){
    theta[s] ~ dbeta(2,2)
  }
}
"
writeLines(modelString, con="TEMPmodel.txt")

library(rjags)
library(runjags)
jagsModel=jags.model(file="TEMPmodel.txt",data=dataList)

在这种情况下，我收到了一条错误消息。

Error in jags.model(file = "TEMPmodel.txt", data = dataList) : 
  RUNTIME ERROR:
Cannot insert node into theta[1...2]. Dimension mismatch

我不知道我在这段代码中弄错了什么。请给我建议。

提前致谢。

Answer 1

正如@nicola所建议的那样，问题在于您将s作为数据传递给模型，而且还使用s作为迭代1:Nsubj的计数器。正如JAGS指出的那样，这会引起对theta维度的混淆......它的长度是15还是2？

以下作品：

model{
  for(i in 1:Ntotal){
    y[i] ~ dbern(theta[s[i]])
  }
  for(j in 1:Nsubj){
    theta[j] ~ dbeta(2,2)
  }
}

Rjags错误消息：尺寸不匹配

1 个答案: