進行荟萃分析時,您很可能必須使用通用度量将效果大小計算或轉換為效果大小。有多種工具可以執行此操作 。
計算效果大小
R語言涵蓋了大多數效果尺寸計算和轉換選項,可以讓您更好地了解。例如,從t檢驗中擷取效果大小:
1. esc_t(t, p, totaln, grp1n, grp2n,
2.
3. es.type = c("d", "g", "or", "logit", "r", "cox.or", "cox.log"),
4.
5. study = NULL, ...)
然後,您可以根據可用參數來計算效果大小,如下所示:
1. # 不相等的樣本量
2.
3. esc_t(t = 3.3, grp1n = 100, grp2n = 150)
4.
5.
6.
7. # 樣本大小相等
8.
9. esc_t(t = 3.3, totaln = 200)
10.
轉換效果大小
軟體提供了多種功能,可将一種效應量轉換為另一種效應量大小:( 标準差 平均對數比),(标準差平均對數比), (标準差平均對數r), (奇數比)到标準差的平均值), (将相關系數r轉換為Fisher的z)和 (将Fisher的z轉換為相關系數的r)。
工作流程
效果大小計算函數的結果以清單形式傳回 。
1. e1 <- esc(grp1yes = 30, grp1no = 50, grp2yes = 40,
2.
3. grp2no = 45, study = "Study 1")
4.
5.
6. e4 <-mean_sd(grp1m = 7, grp1sd = 2, grp1n = 50, grp2m = 9, grp2sd = 3,
7.
8. grp2n = 60, es.type = "logit", study = "Study 4")
現在,mydat 包含一個資料幀, 其中包含幾種效果大小計算的結果:
> mydat
然後按如下方式計算荟萃分析(請注意,不同的效應量度量僅用于示範目的–通常,您應該隻有一個共同的效應量才能進入荟萃分析):
rm(yi = es, sei = se, method = "REML", data = mydat)