2016-10-27 12 views
1

JAGSモデル内の変数に確率p [i]の値を割り当て、別の値を確率1-p [i]に割り当てる必要があります。 p [i]はモデル内の変数です。JAGSの確率的ifelseステートメント

私は以下を試しましたが、それが適切であるかどうかはわかりません。私は、1つの値(0または1)が被験者iに帰され、この値は常に同じままであることを恐れている:与えられた被験者iについて、ベクトルa [i]の値は常に20または30に等しい私が代わりにしたいのは20と30を含むベクトルで、20がp [i]に等しい確率です。

for (i in 1:n) { 
    m[i] ~ dbin(p[i],1) 
    a[i] <- ifelse(m[i]==1, 20, 30)  
    } 
+0

BUGS言語なる - ぎざぎざが読み取るものである - モデル宣言型言語でありますそれはRに触発され、実際にR言語のサブセットではなく、Rのように見え、感じることを意味しました。これは手続き型言語です。したがって、R関数 'ifelse'はバグ言語の一部ではありません。 –

答えて

3

あなたのようなことができます。

for(i in 1:n){ 
m[i] ~ dbern(p[i]) 
a[i] <- (m[i] * 20) + ((1 - m[i]) * 30) 
} 

それは0をとる場合ベルヌーイ証跡が確率pで1をとる場合、aは20となるが30

関連する問題