data(anticoagulation_df)
m1 <- meta::metagen(yi, sqrt(vi), sm = "OR", data = anticoagulation_df,
studlab = LETTERS[1:18])
set.seed(1090) # Make sampled effect sizes reproducible
sample1 <- getsamples(m1)
dis1 <- inc(sample1, dt1 = 16, dt2 = 31, dt3 = 60, br = 0.504, sm = "OR")
dis1
if(requireNamespace("metafor")){
# Same sampled effect sizes using R package 'metafor' (must be installed)
m2 <- metafor::rma(anticoagulation_df, measure = "OR", slab = LETTERS[1:18])
set.seed(1090) # Make sampled effect sizes reproducible
sample2 <- getsamples(m2)
dis2 <- inc(sample2, dt1 = 16, dt2 = 31, dt3 = 60, br = 0.504, sm = "OR")
dis2
all.equal(sample1, sample2) # Only difference: package name
}
Run the code above in your browser using DataLab