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

チラシの裏の落書き日記

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

list操作めも

R

Rのリスト操作について改めてまとめる。 特に,purrrなどを用いた操作を理解すること。purrrの使い方についてのページ purrr: ループ処理やapply系関数の決定版 - Heavy Watal Hadley神のListについてのページ Lists · R for Data Scienceforループ的なlappl…

Rの関数のソースコードなど見たいときの方法

R

一番有名な,関数名をタイプするものでは見れないものについて。 あまり日本語の記事が無いようなきがしたので,メモを残してておく。 library("Matrix") Matrix # Matrix関数の中にある謎の関数。このままだとエラーが返ってくる。 spV2M # コロン3つなの…

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

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

Lavaanによるパラメータ推定結果を見やすくグラフにする。(Gelman's secret weaponのパクリ)

R

尊敬する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を使う

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

jagsUIなるパッケージ

R

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

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} による非テーブルデータ…