1

以下の WinBUGS コードを動作させることができませんでした。通常の事前分布では機能しますが、一様事前分布では機能しません。クリックした後に表示されるエラー メッセージcompileはarray index is greater than array upper bound for age. どういう意味ですか？以下のコードの作業を手伝ってもらえますか？

model
{
for (i in 1:n) {
# Linear regression on logit
logit(p[i]) <- alpha + b.sex*sex[i] + b.age*age[i]
# Likelihood function for each data point
frac[i] ~ dbern(p[i])
}
alpha ~ dunif(0, 1) # Prior for intercept
b.sex ~ dunif(0, 1) # Prior for slope of sex
b.age ~ dunif(0, 1) # Prior for slope of age
}
Data
list(sex=c(1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0,     1,
1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1,     1, 0,
0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1,      1, 1,
0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1),
age= c(69, 57, 61, 60, 69, 74, 63, 68, 64, 53, 60, 58, 79, 56, 53, 74, 56, 76, 72,
56, 66, 52, 77, 70, 69, 76, 72, 53, 69, 59, 73, 77, 55, 77, 68, 62, 56, 68, 70, 60,
57, 51, 51, 63, 57, 80, 52, 65, 72, 80, 73, 76, 79, 66, 51, 76, 75, 66, 75, 78, 70,
67, 51, 70, 71, 71, 74, 74, 60, 58, 55, 61, 65, 52, 68, 75, 52, 53, 70),
frac=c(1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1,        0,
1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1,        1, 1,
1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1,      1, 1,
1, 0, 1, 1, 0, 0, 1, 0, 0, 1),
n=100)
Initial Values
list(alpha=0.5, b.sex=0.5, b.age=0.5)
4

1 に答える 1

5

ああ、それは明らかです。WinBUGS は言いarray index is greater than array upper bound for ageます。それは明らかにエラーを示唆しています->あなたが持っていることがわかりn = 100、ageリストは十分に長くありません:

> your_list <- list(...)
> str(your_list)
List of 4
 $ sex : num [1:100] 1 1 1 0 1 1 0 0 0 0 ...
 $ age : num [1:79] 69 57 61 60 69 74 63 68 64 53 ...
 $ frac: num [1:100] 1 1 1 0 1 1 0 1 1 0 ...
 $ n   : num 100

とにかく、ここではユニフォームを使用しません。自分が何をしているのかを実際に理解していない限り、フラットノーマルなどをお勧めしdnorm(0, 1.0E-10)ます。係数には負の値も許可する必要があります。「帰無仮説」は通常、係数がゼロであるということです。そのため、係数の事後分布の平均値がゼロになるには、「両側からある程度のスペースを確保する」必要があります(直感的に言います)。

于 2013-01-22T22:32:18.237 に答える