統計とか,研究とか,日常についての備忘録的なもの。この広告は、90日以上更新していないブログに表示しています。2016-06-11Mplusで潜在クラス分析を行うときのめもMplusこのノートの荒っぽいまとめ。https://www.statmodel.com/examples/webnotes/webnote14.pdf3ステップで潜在クラスの数を決める。14の挿入ブートストラップ尤度比検定を行うには結構時間がかかるので,いきなりやるより,手順を踏んだほうが良いようだ。クラス数kとk-1での最適な対数尤度を探す。120;などとすると,初期値を600個ランダムに生成して,初期の対数尤度がいけてるもののうち120個を推定を行い,その中でベストなものを最終的なアウトプットにするというもの。普通はそこまででいいのだけれど,潜在クラス数を決める場合には,誤ったseedだと検定がうまくできない。そこで,このステップではイケてる乱数のSEEDを探すということになる。推定が終わったら,outputの中身を"seed"などで検索して,対数尤度のランキング表を見つけること。そこにイケてるseedがある。つぎに,このseedを使って次のステップに入る。Rubin(LMR)検定が出来る。al.(2001)にあるらしい。ブートストラップ尤度比検定(BLRT)より計算量が少ないので,先にやってあたりをつける感じ。最後はBLRTの登場このseedを使って次のステップに入る。具体的には,先ほどとどうように。11はもはや入れなくて良い。これで推定しておしまい。ところが警告が出ることがある。20;などとして初期値の個数を増やす。前半2つはクラス数k-1について,後半2つはクラス数kについて。駄目だったらモデル自体を再考するところから始める。putについて月別アーカイブはてなブログをはじめよう!kazzstatさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。