0

私は JAGS/OpenBUGS をまったく初めて使用するので、モデルの指定に関して正しい方向に進んでいただければ幸いです。私は、103 か国が 15 年間にわたってまとめたアンバランスな縦断データを使用しています。この場合は 12 年が選択されています。DV はジニ係数であり、対数正規ではなく、おそらくベータでモデル化する必要がありますが、現在は 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 がデータ入力を読み取る方法についてはまだ少し混乱しています。

乾杯

J

4

1 に答える 1

0

データフレームyには 2 つの列しかありません。しかしTime、12です。

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

ループ内

for(t in 1:Time){

が 3になるとどうなるか考えてみてください ( =12tに向かう途中)。Time

存在しない y[i,3] を調べるように JAGS に要求しています。したがって、「範囲外のインデックス」。

于 2016-08-14T01:44:40.953 に答える