0

我对 JAGS/OpenBUGS 完全陌生,所以在指定我的模型时,我真的很感激朝着正确的方向前进。我正在使用由 103 个国家在 15 年内编制的不平衡纵向数据,在这种情况下选择 12 年。DV 是基尼系数,它不应该被建模为 log-Normal,而可能是 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 个“技巧”可以应用于不平衡数据,但我仍然对所有这些工作原理感到有些困惑,即 JAGS 如何读取数据输入。

干杯

Ĵ

4

1 回答 1

0

您的数据框y只有 2 列。但是Time是 12。你在哪里

y[i,t]~dlnorm(mu[i,t],tau.1)

在一个循环内

for(t in 1:Time){

想想当t上升到 3 时会发生什么(在上升到Time=12 的路上)。

您要求 JAGS 查看不存在的 y[i,3]。因此“索引超出范围”。

于 2016-08-14T01:44:40.953 回答