RのJAGSを使用して離散値にワイブルモデルをフィッティングしています。連続データにワイブルを当てはめるのに問題はありませんが、離散値。ここで は、いくつかのデータであり、ぎざぎざでワイブルモデルに合うようにコード: #draw data from a weibull distribution
y <- rweibull(200, shape = 1, scale = 0.9)
#y <- r
R2WinBUGSを使用していくつかのシミュレーションスタディを実行する方法を理解する上で問題があります。目的は、n個のデータセットをシミュレートすることです(1000を目標としますが、10から始まります)。それらをすべてR2WinBUGSコードに入れて、WinBUGSに移植すると、n個のデータセットの見積もりを実行します。 モデル::ここに私は現在持っているものです model{
a
以前はOpenBUGS/WinBUGSを使用してベイジアン統計を行ってきましたが、PythonでPYMC3パッケージを使用することにしました。だから私はかなり新しいpacakageですが、まだ完全に使い方を学んでいます。私はバグコードをPYMC3に変換するのにいくつかの困難を抱えています。次のようにオリジナルのバグコードは次のとおりです。 model {
for (i in 1 : N) {