統計とか,研究とか,日常についての備忘録的なもの。この広告は、90日以上更新していないブログに表示しています。R気が向いたらStanやJAGSなどとの対応関係をきちんと取りならが検算を行っていきたい。MCMCサンプルは「.gh5」ファイルに入っているので,これを読み込む。以下の資料のp.65あたりの「mplus.get.bayesian.」を利用すると,ギブスサンプラーのMCMCサンプルが得られるはず。http://www.statmodel.com/mplus-R/Mplus%20R%20tutorial.pdf以下のinpファイルを作成して,Mplusで推定を行う。データはhttps://www.statmodel.com/usersguide/chap4/ex4.1a.datにあるものを,ダウンロードしてくる。MacではRを使ってプロットをする必要がある。以下のsourceを読み込むと必要な関数がいろいろと読み込まれる。Rstudioのコンソールでmplusと入力して,タブを押すと色々な関数が候補として出てくる。プロットできるもの一覧を表示する。データの取得:1変数だけ。mplus.get.data("ex4.1bayes.gh5","Y1")ベイズ推定の結果をプロットするためのもの。利用しきれていない関数もあるので,今後追記する。ついでに,パラメタ名が必要な関数も表示してくれる。特定のパラメタのMCMCサンプルを取得する。列がチェインに対応している。par_list[1])Mplusは便利なので,もう少し使いこなしていきたいところ。modelコマン...2017-10-20平均値の信頼区間のシミュレーションプロットRで,信頼区間を書いていいプロットするコードのメモ。layout...2017-08-15Mplusによるモンテカルロ・シミュレーション覚書2今回は前回紹介しなかったコマンドなどについて,紹介する。putについて月別アーカイブはてなブログをはじめよう!kazzstatさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。