タグ付けされた質問 「winbugs」

2
二項分布の
この質問は、技術的なフォローアップでこの質問。 Raftery(1988)でNNN提示されたモデルの理解と複製に問題があります:二項Nパラメーターの推論: WinBUGS / OpenBUGS / JAGSの階層ベイズアプローチ。ただし、コードだけではないので、ここでトピックを取り上げる必要があります。 バックグラウンド ましょうの未知との二項分布から、成功回数の集合Nとθ。さらに、Nはパラメーターμのポアソン分布に従うと仮定します(論文で説明)。その後、各xは、私は平均のポアソン分布持っλ = μ θを。λとθの観点から事前確率を指定したい。x=(x1,…,xn)x=(x1,…,xn)x=(x_{1},\ldots,x_{n})NNNθθ\thetaNNNμμ\muxixix_{i}λ=μθλ=μθ\lambda = \mu \thetaλλ\lambdaθθ\theta またはθに関する十分な事前知識がないと仮定して、λとθの両方に情報のない事前分布を割り当てたいと思います。言う、私の事前確率であるλ 〜G A M M A(0.001 、0.001 )およびθ 〜U N I F O R M(0 、1 )。NNNθθ\thetaλλ\lambdaθθ\thetaλ∼Gamma(0.001,0.001)λ∼Gamma(0.001,0.001)\lambda\sim \mathrm{Gamma}(0.001, 0.001)θ∼Uniform(0,1)θ∼Uniform(0,1)\theta\sim \mathrm{Uniform}(0, 1) 著者は不適切な事前分布を使用していますが、WinBUGSは不適切な事前分布を受け入れません。p(N,θ)∝N−1p(N,θ)∝N−1p(N,\theta)\propto N^{-1} 例 紙(226ページ)では、観測されたウォーターバックの次の成功数が用意されています:。母集団のサイズであるNを推定したい。53,57,66,67,7253,57,66,67,7253, 57, 66, 67, 72NNN WinBUGS(@StéphaneLaurentのコメントの後に更新)の例を解決しようとした方法は次のとおりです。 model { # Likelihood for …

1
ネットワークのメタ分析に最適な方法はどれですか?
現在、ネットワークのメタ分析または混合治療の比較を実行するためのいくつかの異なるアプローチがあります。 最も一般的に使用され、アクセス可能なものは、おそらく次のとおりです。 ベイジアンフレームワークで: WinBUGSの処理ごとの相互作用アプローチ(例、Jackson et al); WinBUGSでの階層的な腕ベースのベイズモデリング(たとえば、Zhao et al); 階層的なコントラストに基づく(すなわち、ノード分割)WinBUGS伴うまたは介してのいずれかでベイジアンモデリング、gemtc及びrjagsR(例えばディアスら若しくはバンValkenhoefら)。 WinBUGSに統合されたネストされたラプラス近似(INLA)(例えば、Sauter et al); 常連主義の枠組みで: SASの要因分散分析(例:Piepho); SASでのマルチレベルネットワークメタ分析(例:Greco et al); mvmetaStataまたはRでの多変量メタ回帰(例、White et al); R lmeとのネットワークメタ分析netmeta(例:Lumley、ただし2群試験に限定、またはRuckerら)。 私の質問は、単純です:それらはほぼ同等ですか、またはほとんどの場合に一次分析に好ましいものがありますか(したがって、他のものを補助分析に予約します)? 更新 時間の経過とともに、ネットワークメタ分析の方法に関するいくつかの比較分析が行われてきました。 カーリンBP、ホンH、シャムヤンTA、セイントF、ケインRL。複数の治療を比較するためのベイジアンアプローチとフリークエンティストアプローチの比較に関する事例研究。Healthcare Research and Quality(米国)の代理店。2013。

3
Winbugsおよびその他のMCMC(事前配布の情報なし)
パラメータの分布がわからない場合はどうなりますか?どのアプローチを使用する必要がありますか? ほとんどの場合、特定の変数が特定の種の存在/不在に何らかの影響を及ぼし、その変数が変数の重要度に従って受け入れられるかどうかについては、十分に検討することを目指しています。これは、ほとんどの場合、パラメーターが持つはずの推定分布については考えていません。 b1、b2、b3およびb4が-2と2の間で変化し、b0が-5と5の間で変化する可能性があることを私が知っているすべてのパラメーターが正規分布に従うと仮定することは正しいですか? model { # N observations for (i in 1:N) { species[i] ~ dbern(p[i]) logit(p[i]) <- b0 + b1*var1[i] + b2*var2[i] + b3*var3[i] + b4*var4[i] } # Priors b0 ~ dnorm(0,10) b1 ~ dnorm(0,10) b2 ~ dnorm(0,10) b3 ~ dnorm(0,10) b4 ~ dnorm(0,10) }
10 r  bayesian  mcmc  bugs  winbugs 

1
WinBUGSで特定の変数の予測を取得する方法は?
私はWinBUGSの新しいユーザーであり、あなたの助けに1つの質問があります。次のコードを実行した後、私はのパラメータを持っbeta0てbeta4(統計情報、密度)が、私は最後の値の予測を取得する方法がわからないh私がするように設定、NAコード内でモデル化することを。 誰かが私にヒントを与えることができますか?何かアドバイスをいただければ幸いです。 model { for(i in 1: N) { CF01[i] ~ dnorm(0, 20) CF02[i] ~ dnorm(0, 1) h[i] ~ dpois (lambda [i]) log(lambda [i]) <- beta0 + beta1*CF03[i] + beta2*CF02[i] + beta3*CF01[i] + beta4*IND[i] } beta0 ~ dnorm(0.0, 1.0E-6) beta1 ~ dnorm(0.0, 1.0E-6) beta2 ~ dnorm(0.0, 1.0E-6) beta3 ~ dnorm(0.0, …
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.