メタ分析（変量効果モデル）

データの読み込み

library(metafor)
library (MAc)
library (MAd)


dat=read.csv("csv_all2.csv")
summary(dat)

変量効果モデル

res=rma(yi,vi,data=dat)
res
forest(res)

調整変数ごとの平均効果量

macat(yi,vi,mod = subject,method ="random",data = dat)
macat(yi,vi,mod = grade,method ="random",data = dat)
macat(yi,vi,mod = century,method ="random",data = dat)
macat(yi,vi,mod = Measure,method ="random",data = dat)
macat(yi,vi,mod = purpose,method ="random",data = dat)
macat(yi,vi,mod = technology,method ="random",data = dat)

funnel plot

#funnel plot
funnel(res)
regtest(res,model = "rma")

fail-safe N

#fail-safe N
FSN=(res$zval*11)^2/2.706-11
FSN
fsn (y = res$yi, v = res$vi)

plot(x=dat$year,y=dat$yi)

trim and fill

#trim and fill
taf2=trimfill(res)
funnel(taf2)

p.rep

#p.rep
p.rep.old=pnorm(pnorm(1-0.0008)/sqrt(2))
p.rep.old #Lilleen(2005)
t=res$b/res$se
p.rep=pt(abs(t)/sqrt(2),df=9)
p.rep #Lecoutre(2010)