SSブログ

[Stan] beta_proportion 分布 [統計]

Stanにbeta_proportionという分布がはいっていたのでした。平均をつかってベータ分布のあてはめができます。

R notebookです。


データと
Rplot01.png

事後分布の要約と
> print(fit, pars = c("beta", "kappa"))
Inference for Stan model: stan-1af33ea9a9e9.
4 chains, each with iter=2000; warmup=1000; thin=1; 
post-warmup draws per chain=1000, total post-warmup draws=4000.

         mean se_mean   sd  2.5%   25%   50%   75% 97.5% n_eff Rhat
beta[1] -0.95    0.00 0.12 -1.19 -1.04 -0.96 -0.87 -0.72  2787    1
beta[2]  2.25    0.00 0.22  1.82  2.11  2.25  2.39  2.69  2385    1
kappa    1.89    0.01 0.27  1.41  1.70  1.88  2.07  2.46  2423    1


事後予測チェックの結果です。
Rplot02.png

タグ:STAn
nice!(1)  コメント(0) 
共通テーマ:日記・雑感

nice! 1

コメント 0

コメントを書く

お名前:
URL:
コメント:
画像認証:
下の画像に表示されている文字を入力してください。

Facebook コメント