Rで統計(メモ)

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

ベイズ

段階反応モデルにおける特性値θのベイズ推定

困難度・識別力パラメータの読み込み 特性値θの推定 困難度・識別力パラメータの読み込み #Q7,9,10,11,12,13,14 雲財中村(2018)に基づくEAP推定 Extrmt1=c(-2.09,-1.55,-1.99,-1.86,-2.35,-1.48,-1.38) Extrmt2=c(-0.79,-0.55,-0.90,-0.78,-1.35,-0.49,-0.…

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

データの読み込み 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…

対応のある分散分析(1要因)forベイズ

対応のある分散分析(1要因) データの読み込み stanコード格納 連言命題が正しい確率 特定のパラメーターだけプロット 対応のある分散分析(1要因) データの読み込み #作業ディレクトリの設定(適宜変更) setwd("C:/Users/turid/OneDrive/ドキュメント"…

頻度論的な統計手法のオルタナティブ

相関分析 データの格納 ピアソンの積率相関 ピアソンの積率相関 ver.ベイズ stanコード格納 MCMC実行 2群の平均値の差の検討(性差の検討を例に) データの格納 対応のないt検定 対応のないt検定 ver.ベイズ stanコード格納 相関分析 データの格納 score_1<…