天天看點

拓端tecdat|R語言Meta分析效應量

進行荟萃分析時,您很可能必須使用通用度量将效果大小計算或轉換為效果大小。有多種工具可以執行此操作 。

計算效果大小

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      
拓端tecdat|R語言Meta分析效應量

然後按如下方式計算荟萃分析(請注意,不同的效應量度量僅用于示範目的–通常,您應該隻有一個共同的效應量才能進入荟萃分析):

rm(yi = es, sei = se, method = "REML", data = mydat)