草薙の研究ログ

英語の先生をやってます。

metacorパッケージで相関係数のメタ分析

Rのメタ分析パッケージ多すぎて何がなんだかなんだけど,相関係数関係はmetacorってのがシンプル。rmeta風ね。


https://cran.r-project.org/web/packages/metacor/metacor.pdf


#パッケージの準備
install.packages("metacor")
library(metacor)

#データの準備
r<-c(.2,.3,.4,.2,.1,.7,.6,.3)
n<-c(30,28,100,20,300,20,30,100)
l<-c("Taro (1991)","Ziro (1993)","Saburo (2000)","Shiro (2011)","Goro (2012)","Rokuro (2014)","Nanaro (2015)","Hachiro (2016)")

#変量効果モデルでやる
DSL<-metacor.DSL(r,n,l)
DSL

#固定効果モデルでやる
OP<-metacor.OP(r,n,l)
OP


f:id:kusanagik:20160427160219p:plain


2 by 2 by xの層別データの分析のとき,マンテルヘンツェル法によるオッズ比のメタ分析は,テトラコリック相関係数のメタ分析と結構一緒になるだろうね。