Rで統計(メモ)

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

ロジスティック回帰分析

ロジスティック回帰分析(2値) ロジスティック回帰分析(2値) #データの読み込み library(readxl) dat1 <- read_excel("C:/Users/2ken/Desktop/Umeda3.xlsx") View(dat1) dat1=na.omit(dat1) #基礎集計 library(psych) describe(dat1) #データ型の確認と変…

メタ分析(変量効果モデル)

データの読み込み 変量効果モデル 調整変数ごとの平均効果量 funnel plot fail-safe N trim and fill p.rep データの読み込み library(metafor) library (MAc) library (MAd) dat=read.csv("csv_all2.csv") summary(dat) 変量効果モデル res=rma(yi,vi,data=…

多重比較(ボンフェローニ)

群間比較(多重比較,ボンフェローニ) 群間比較(多重比較,ボンフェローニ) #学校ごとの平均算出 aggregate(X1.all~school,data=dat4,FUN = mean) #多重比較 pairwise.t.test(dat4$X1.all, dat4$school, p.adj = "bonf")

変数関係の視覚化(ペアプロット)

GGallyパッケージ psychパッケージ GGallyパッケージ library(ggplot2) library(GGally) dat6=subset(dat4,select = c(school,think1_pre,think1_post,eap_pre,eap_post,X4_8,Q2.8)) colnames(dat6) <- c("学校","思考プレ","思考ポスト","態度プレ","態度ポ…

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

困難度・識別力パラメータの読み込み 特性値θの推定 困難度・識別力パラメータの読み込み #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…

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/ドキュメント"…

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

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

通常の因子分析

KMOとMSAの確認 サンプリングの適切性を検討する 因子数の検討 平行分析(通常版/ポリコリック相関行列版) スクリープロット 固有値 探索的因子分析 因子数1の検討(通常版/ポリコリック相関行列版) 信頼性係数の算出(アルファとオメガ) 確認的因子分析 …

データの読み込みと整形

データの読み込みと整形~標準編~ データの読み込み データセットの作成と確認 データの読み込みと整形~大規模データ編~ データの読み込み データセットの作成 データの読み込みと整形~標準編~ データの読み込み #作業ディレクトリの設定(適宜変更) s…