SSブログ

[Stan] multinomialをつかった回帰 [統計]

3種の生物が生息している地域のN地点でそれぞれM個体を採取する。採取した3種の個体数(の割合)と共変量Xとの関係を分析する。という設定をかんがえてみたテストです。

コードです。

データです。

Rplot01.png

結果です。

            mean se_mean   sd    2.5%     25%     50%     75%   97.5% n_eff Rhat
alpha[1]    3.35    0.04 1.15    1.17    2.64    3.32    4.09    5.62   881 1.00
alpha[2]   -0.08    0.02 0.48   -0.98   -0.34   -0.07    0.20    0.77   905 1.00
alpha[3]   -3.36    0.01 0.48   -4.27   -3.65   -3.35   -3.07   -2.47  1233 1.00
beta[1]    -6.55    0.04 1.34   -9.27   -7.42   -6.48   -5.68   -4.13  1162 1.00
beta[2]    -0.05    0.03 0.72   -1.42   -0.40   -0.02    0.38    1.12   621 1.00
beta[3]     6.14    0.06 1.59    2.88    5.27    6.21    7.09    9.00   697 1.00
lp__     -934.31    0.06 1.92 -938.95 -935.31 -933.93 -932.92 -931.71   900 1.01

Rplot02.png


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

nice! 1

コメント 0

コメントを書く

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

Facebook コメント