[Stan] multinomialをつかった回帰(3)

前回のsoftmaxをつかったモデルでは信用区間がひろかったので、何か制約をいれた方がよさそうということで、パラメータの事前分布をせまくしてみました。



結果は信用区間がせまくなりました。まあそれはそうだ。

          mean se_mean   sd  2.5%   25%   50%   75% 97.5% n_eff Rhat
alpha[1]  1.30    0.01 0.59  0.11  0.90  1.30  1.70  2.43  2514    1
alpha[2]  0.13    0.01 0.59 -1.05 -0.27  0.13  0.53  1.28  2540    1
alpha[3] -1.50    0.01 0.60 -2.70 -1.90 -1.49 -1.08 -0.36  2541    1
beta[1]  -2.66    0.01 0.62 -3.84 -3.09 -2.66 -2.25 -1.44  2514    1
beta[2]   0.01    0.01 0.61 -1.17 -0.41  0.01  0.44  1.21  2436    1
beta[3]   2.63    0.01 0.62  1.42  2.21  2.64  3.04  3.87  2434    1

タグ：STAn