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

尊敬するHoxo_mさんらが翻訳された”みんなのR”にあるGelman's secret weaponを真似して自分の武器を作成してみようかと思う。 d.hatena.ne.jp https://www.amazon.co.jp/みんなのR-データ分析と統計解析の新しい教科書-Jared-P-Lander-ebook/dp/B012Z86Q6Qや…

Rから今流行のスクリプト言語Juliaを使いたい。 普通にJuliaにパスを通してRのsystem関数を使えばいいのだけれど，パスを通さずにRで下のような関数を書いて作業をする。jlfileにはJuliコードを書いた拡張子.jlであるファイル名を入れる。jlfilepathは作業フ…

このノートの荒っぽいまとめ。 https://www.statmodel.com/examples/webnotes/webnote14.pdf３ステップで潜在クラスの数を決める。 最適な対数尤度探し optseedとtech 11の挿入 optseedとtech 14の挿入 ブートストラップ尤度比検定を行うには結構時間がかか…

Mplusはミクスチャーモデルの分析もらくらくできるし，計算も結構速い。しかし，人間業が深いもので，もっと早く分析して欲しい，なんていう欲望が出てくる。 そんな時はamalysisコマンドの中にprocessors = 4;などと記述すると，複数コアで計算を回してくれ…

MCMCサンプラーのうちの1つであるJAGSをRから利用する新しいパッケージができたようだ。基本はrjagsのラッパーのようだが，並列計算の指定などがやりやすそうだ。今後，これを中心に使っていくかもしれないのでヘルプにあった例をコピペしておく。初期値の指…

学会に行く意味とはなんだろうか。今回はこのことに関して自分にとって何が良いことなのか考えてみる。 まず，第1に締め切りがあることによる効果が挙げられる。締め切りがあることによって，適当なところで研究をまとめなければならない。コレによって，1つ…

いきなり新しい趣味を始めたくなった。昔から見るのが好きだったジャグリングをやろうと思う。一番基本的なボールというやつから手をつける。3つのボールを小気味良くポコポコと投げ続けられるようになりたいものだ。 ネットでやり方を調べてみたり，youtube…

自分で文章を書くときにやりがちなことをまとめておく。文章を書くのは苦手なので，注意点をキチンと書いてみる。この文章の目的は，自分にとっての現在把握している気をつけるべき点を明確にし，書き残すことで自分の作文能力を向上させるための資料を作成…

Mplusは潜在変数モデリングに適した解析ソフトで，心理学，社会学，疫学などの分野で近年人気が出ているソフトウェアである。ホームページは次のものがある。 Muthén & Muthén, Mplus Home Pageこのソフトは非常に便利で，文法も簡潔なのであるが，文法が特…

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

今日やったこと 某ポスター作成 発表資料の作成 MCMCの観察 なす購入 1.2.については直前に再度確認を行うこと。3. について，バーンインは結構早い。サンプリングがえらく時間がかかる。いらないパラメータはサンプリングしないでいたほうがいいか？ 欲しい…

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

前回から随分間が空いているが，気にしない。 思ったらやってみるのが大事な気がする。最近の困り 研究のオリジナリティが良くわからない。何を研究したら面白いのか全くわからない。 査読のコメントへの対応が困った。オリジナリティの主張の仕方，伝え方。…