meta分析亚组分析R实现

library(meta)
bin.metagen<-read.table("clipboard",header=T,sep="\t")
bin.metagen

bin.metagen$RR <- log(bin.metagen$RR)
bin.metagen$lower <- log(bin.metagen$lower)
bin.metagen$upper <- log(bin.metagen$upper)
bin.metagen$seTE <- (bin.metagen$upper - bin.metagen$lower)/3.92

fit<-metagen(RR, 
        seTE, 
        studlab = Author,
        method.tau = "SJ",
        sm = "RR",
        data = bin.metagen)
forest(fit)

 

fit<-metagen(RR, 
             seTE, 
             studlab = Author,
             method.tau = "SJ",
             sm = "RR",
             data = bin.metagen,byvar = Group)
forest(fit)

 

 ref:https://bookdown.org/MathiasHarrer/Doing_Meta_Analysis_in_R/binary.html