統計とか,研究とか,日常についての備忘録的なもの。高速なコードにはなっていないが,これで基本的にはOK。詳しくはヘルプを見たり,exampleを走らせればすぐにわかる。簡単に使えるので,良さそう。統計IRTモデルで予測分布や対数尤度の計算などをおこなってみる。対数密度の記述は以下の頁の質問に少し書かれている。統計やり方は簡単で,予測分布を未知の変数としてモデルに書けばいい。sigma)のように乱数を生成してくれる...JAGSでの予測分布の作成2017-12-21JAGSで対数尤度を得る方法(Stanも)RWAICなどを計算するには,対数尤度のMCMCサンプルが必要。個々のMCMCサンプルからがんばって計算する方法もあるが,より簡単にはMCMCサンプラーのなかで対数尤度を得ると楽。quantityが容易されているが,JAGSでは少し工夫が必要。R気が向いたらStanやJAGSなどとの対応関係をきちんと取りならが検算を行っていきたい。MCMCサンプルは「.gh5」ファイルに入っているので,これを読み込む。以下の資料のp.65あたりの「mplus.get.bayesian.」を利用すると,ギブスサンプラーのMCMCサンプルが...MplusのMCMCサンプルをGetする方法(+Mplusの結果をRでプロットする方法)その12017-11-21よく必要になるけど,よく忘れるR(とExcel)の小技メモR文字(数字)の桁数を揃えたい。不要な空白を消したい。全角半角の指定...よく必要になるけど,よく忘れるR(とExcel)の小技メモ2017-10-20平均値の信頼区間のシミュレーションプロットRRで,信頼区間を書いていいプロットするコードのメモ。layout関数で,上手くプット画面を分割して,さも1枚のグラフかのように見せかけるのがポイント。可愛い感じのグラフを描画できる。縦向きの,信頼区間の図。ポイントは,レイアウト...平均値の信頼区間のシミュレーションプロット2017-03-25list操作めもRRのリスト操作について改めてまとめる。purrrなどを用いた操作を理解すること。Scienceforループ的なlappl...list操作めも2017-03-11Rの関数のソースコードなど見たいときの方法R一番有名な,関数名をタイプするものでは見れないものについて。日本語の記事が無いようなきがしたので,メモを残してておく。Matrix関数の中にある謎の関数。このままだとエラーが返ってくる。統計ちょっと気になったので,近年一部界隈で話題になっているベイズ情報仮説の評価をやってみる。普通に計算すると大変なベイズファクターも,情報仮説の元では簡単に計算できる。http://www.tandfonline.com/doi/abs/10.1080/17405629.201...MplusでMCMCした結果を使ったSEMでの情報仮説の評価2016-06-20Lavaanによるパラメータ推定結果を見やすくグラフにする。weaponを真似して自分の武器を作成してみようかと思う。https://www.amazon.co.jp/みんなのR-データ分析と統計解析の新しい教科書-Jared-P-Lander-ebook/dp/B012Z86Q6Qや...Lavaanによるパラメータ推定結果を見やすくグラフにする。JuliaRから今流行のスクリプト言語Juliaを使いたい。普通にJuliaにパスを通してRのsystem関数を使えばいいのだけれど,パスを通さずにRで下のような関数を書いて作業をする。jlfileにはJuliコードを書いた拡張子.jlであるファイル名を入れる。jlfilepathは作業フ...RからJuliaを使う2016-06-07jagsUIなるパッケージRMCMCサンプラーのうちの1つであるJAGSをRから利用する新しいパッケージができたようだ。基本はrjagsのラッパーのようだが,並列計算の指定などがやりやすそうだ。今後,これを中心に使っていくかもしれないのでヘルプにあった例をコピペしておく。統計3レベルの階層線形を行う時のlmerのメモ学校のなかに,個人がネストされていて,個人の中に測定時点がネストされている場合。最近のRの界隈で流行っている(気がする)purrrパッケージの使いかたについて。参照させてもらうので,とりあえずまとめておく。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。