stanとRでMCMCした結果をいじりたい(fitの結果の取り出し方)

確立プログラミング言語stanでMCMCやった後の結果を個別に取り出したいことがある。例えば、複数の仮説の同時確立を出したい時(連言命題の確立)など…
そのまま取り出そうとしても無理。ちょっと迷ったので、忘れないようにメモしておく。
以下のコードが大事。

rstan::extract(fit)

でOK。ちなみにfitはstan()で計算された中身のこと。つまり、stan()でMCMCやった結果のことをfitで表している。

fit_data = rstan::extract(fit)

あるいは

fit_data < - rstan::extract(fit)

みたいにすればよい。

コメント