Rで統計(メモ)

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

2019-03-01から1ヶ月間の記事一覧

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

データの読み込み 変量効果モデル 調整変数ごとの平均効果量 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.…