2017-10-03 7 views
0

rstanarmに基本的なGLMコードを書きました。私はそれが行われたすべての反復の予測変数の平均値を与えている次のコードを実行すると:rstanarmの反復ごとに異なる値を見つけるにはどうすればよいですか?

summary(TestGLM) 
    print(TestGLM) 

を私は私の解決策のために一人ひとりの反復で予測変数の値を望んでいました。誰も私がこの問題を解決するのを助けることができますか?

+0

は(https://cran.r-project.org/web/packages/rstanarm [マニュアル]で素早く見ていたが/rstanarm.pdf)。 'prior_summary'と' summary.stanreg'に運がありましたか?もしそうでなければ、あなたは 'str(TestGLM)'をして突き進むことができると思います。 –

+0

@EricFailお返事ありがとうございます。私はstr(TestGLM)を使って値を得ることができ、値を得るためにggmcmc(TestGLM)も見つけました。私はposterior_predict(TestGLM)が私に値を与えてくれるのは何だろうと思っていました。 posterior_predict()について考えているかどうか教えてください。 –

+0

上記リンク先のrstanarmマニュアルの23ページをご覧ください。 –

答えて

2

あなたは次のように係数のそれぞれについて、描画抽出することができます。

require(rstanarm) 

TestGLM <- stan_glm(mpg ~ wt, data = mtcars) 
draws <- as.data.frame(TestGLM) 
print(colnames(draws)) 
str(draws) 

# See http://mc-stan.org/rstanarm/reference/as.matrix.stanreg.html 
関連する問題