読者です 読者をやめる 読者になる 読者になる

チラシの裏の落書き日記

統計とか,研究とか,日常についての備忘録的なもの。

統計

最高の初期値職人を目指して

潜在クラス分析(有限混合モデル)やら,k-meansとかいろいろと初期値に依存する分析するがある。また,Latent variableモデリングではEMアルゴリズムを使った最尤推定を行うことが多いが,複雑なモデルになるほど解が収束しないこともしばしばある。 とくに…

MplusでMCMCした結果を使ったSEMでの情報仮説の評価

ちょっと気になったので,近年一部界隈で話題になっているベイズ情報仮説の評価をやってみる。 普通に計算すると大変なベイズファクターも,情報仮説の元では簡単に計算できる。元ネタは以下の資料 http://www.tandfonline.com/doi/abs/10.1080/17405629.201…

lmerでの3レベルモデルの実行

3レベルの階層線形を行う時のlmerのメモ学校のなかに,個人がネストされていて,個人の中に測定時点がネストされている場合。学校IDをschool_id, 個人IDがind_id,あとは独立変数をx, 従属変数をyとする。簡単のため,中心化は行われているとして,レベル2以…

Rの使いかたメモ(purrrの使いかたについて)

最近読んだ良記事 最近のRの界隈で流行っている(気がする)purrrパッケージの使いかたについて。 よく参照させてもらうので,とりあえずまとめておく。sinhrks.hatenablog.com sinhrks.hatenablog.com uribo.hatenablog.com {purrr} による非テーブルデータ…