Question

我对JAGS / OpenBUGS完全陌生，所以我真的很感激在指定我的模型方面正朝着正确的方向前进。我使用的是一个不平衡的纵向数据，这个数据由超过15年的103个国家编制，在这种情况下选择了12年。 DV是基尼系数，不应该建模为日志 - 正常但可能更确切地说是Beta，尽管现在的重点是仅仅了解如何在JAGS中编译模型。我暂时使用的是固定效果模型。

我正在运行的数据和代码：

> head(x)
          Year     II2       II3       II4     ..... II24
1          1       2.956233  40.90458 4.475183       16.443553
8          1       1.257794  85.47378 2.395186       19.333433
19         1       4.139706 141.07899 2.544640       25.555404
37         1       2.233664  98.51313 3.902835       42.533333
49         1       2.879734  61.39000 1.471334       18.884444
71         1       3.381762  60.23783 3.432614       16.334222


> head(y)
        Year       II1
1         1     0.3240000
8         1     0.2576667
19        1     0.3132500
37        1     0.2700000
49        1     0.2744286
71        1     0.3250000

dim(x)
1224   23

length(y)
1224  

Time <- 12, N <- length(y$II1)#No. of Obs.

dat <- list(x=x, y=y, N=N, Time=Time, p=dim(x)[2]), 
inits <- funtion(){list(tau.1=1, tau.2=1, eta=1, alpha=0, beta1=0, beta2=0, beta3=0)}

model6 <- "model{
for(i in 1:N){for(t in 1:Time){
y[i,t]~dlnorm(mu[i,t],tau.1)
mu[i,t] <- inprod(x[i,t],beta[])+alpha[i]}
alpha[i]~dnorm(eta, tau.2)}

for (j in 1:p) {
b[j]~dnorm(0,0.001)
}

eta~dnorm(0, 0.0001)
tau.2~dgamma(0.01,0.01)
tau.1~dgamma(0.01,0.01)


 }"

reg.jags <- jags.model(textConnection(model), data=dat, inits=inits, n.chains=1, n.adapt=1000)

我不断收到此运行时错误：

Error in jags.model(textConnection(model), data = dat, inits =   inits,  : 
RUNTIME ERROR:
Compilation error on line 3.
Index out of range taking subset of  y

任何关于我应该做些什么的建议都会非常感激！我知道有3个“技巧”你可以应用于不平衡的数据，但我仍然有点困惑所有这些如何工作，e.i。 JAGS如何读取数据输入。

干杯

Ĵ

用JAGS进行纵向层次贝叶斯回归

0 个答案: