Rで統計(メモ)

難しいことはわからないけど、とにかくRでやってみようというページ

2018-02-18から1日間の記事一覧

項目反応理論(段階反応モデル,ベイズ推定)

データの読み込み IRTによる尺度構成 局所独立の検討 単調性の確認 段階反応モデルの実行 基本情報 特性反応曲線 テスト情報関数 項目情報関数 個人の能力値 因子得点と能力値θのプロット 水平テストの構築 段階反応モデル verベイズ 呼び出し stanコード格…

2群の平均値の差の検討(t検定,ベイズ推定)

2群の平均値の差の検討(性差の検討を例に) データの格納 x<-SNFC_1$theta[SNFC_1$sex=="2"] #第1群の指定 y<-SNFC_1$theta[SNFC_1$sex=="1"] #第2群の指定 #平均値の高い方をyにしておくと後で解釈しやすい nx <- length(x) ny <- length(y) data <- list…

ピアソンの積立相関係数(ベイズ推定)

相関分析 データの格納 ピアソンの積率相関 ピアソンの積率相関 ver.ベイズ stanコード格納 MCMC実行 相関分析 データの格納 score_1<-SNFC_1$btheta score_2<-SNFC_1$kakusan score<-scale(cbind(score_1,score_2)) ni<-1317 nt<-2 data<-list(y=score,ni=n…