確立プログラミング言語stanでMCMCやった後の結果を個別に取り出したいことがある。例えば、複数の仮説の同時確立を出したい時(連言命題の確立)など…
そのまま取り出そうとしても無理。ちょっと迷ったので、忘れないようにメモしておく。
以下のコードが大事。
rstan::extract(fit)
でOK。ちなみにfitはstan()で計算された中身のこと。つまり、stan()でMCMCやった結果のことをfitで表している。
fit_data = rstan::extract(fit)
あるいは
fit_data < - rstan::extract(fit)
みたいにすればよい。
そのまま取り出そうとしても無理。ちょっと迷ったので、忘れないようにメモしておく。
以下のコードが大事。
rstan::extract(fit)
でOK。ちなみにfitはstan()で計算された中身のこと。つまり、stan()でMCMCやった結果のことをfitで表している。
fit_data = rstan::extract(fit)
あるいは
fit_data < - rstan::extract(fit)
コメント