메타분석 파이프라인

library(metafor)

# C6 내보내기에서 코드북 로드
data <- read.csv("diverga_codebook_v2.2.csv")

# 랜덤 효과 메타분석
res <- rma(yi = hedges_g,
           vi = variance,
           method = "REML",
           data = data)

# Forest plot
forest(res,
       header = "Study",
       xlab = "Hedges' g")

# Funnel plot (출판 편향)
funnel(res)
regtest(res)  # Egger 검정

import pandas as pd
from metaanalysis import MetaAnalysis

# 코드북 로드
df = pd.read_csv("diverga_codebook_v2.2.csv")

# 메타분석 초기화
ma = MetaAnalysis(df,
                  effect_size="hedges_g",
                  variance="variance")

# 랜덤 효과 모델 실행
results = ma.fit(method="REML")

# Forest plot 생성
ma.plot_forest()

# 이질성 확인
print(f"I²: {results.I2:.1f}%")
print(f"Q: {results.Q:.2f}, p={results.Q_pval:.3f}")