Rで統計(メモ)

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

2018-01-01から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…

3人以上の評価者間の一致度(Fleissのκ係数)

Fleissのκ係数 #データ読み込み library(readr) dat1 <- read_csv("csv1.csv") #欠損値除去 dat1=na.omit(dat1) #変数型の変更 dat1$Mscore=as.numeric(dat1$Mscore) dat1$Msitu=as.numeric(dat1$Msitu) dat1$Nscore=as.numeric(dat1$Nscore) dat1$Nsitu=as.…

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

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