предположим, у меня есть df со следующей информацией, и я хочу создать лесной участок следующим образом:
df <- data.frame(Author = c("A", "B", "C", "D", "E"),
Year = c(2012, 2015, 2019, 2018, 2011),
Exposure = c("AAB", "GTD", "ABA", "CGH", "PLF"),
Outcome = c("XXX", "YYY", "ZZZ", "KKK", "QQQ"),
Case_Yes = c(20, 1, 3, 12, 5),
Case_No = c(15,3,5,7,100),
Control_Yes = c(140,405, 2048,45,10),
Control_No = c(120,234,563,56,776))
library (metafor)
dat <- escalc(ai = Case_Yes, bi = Case_No, ci = Control_Yes, di =Control_No,
data = df, measure = "OR")
res1 <- rma(yi = yi, vi = vi, data = dat, measure = "OR", add = 1/2, to = "only0",
vtype = "LS", slab=paste(Author, Year, sep=", "), method="REML")
Outcome_labels <- paste0(c(expression(italic("XXX")), "YYY",
expression(italic("ZZZ")), "KKK", expression(italic("QQQ"))))
forest(res1, xlim=c(-30,10), ylim = c(-2, 8), at=log(c(0.002, 0.25, 1, 100)), atransf = exp, cex=1, xlab = "Odds Ratio", mlab = "", psize = 1,
ilab=cbind.data.frame(df$Exposure, Outcome_labels, df$Case_Yes, df$Case_No, df$Control_Yes, df$Control_No),
ilab.xpos = c(-23, -18, -13, -10, -8, -6), font = 1)
По какой-то причине мой лесной участок не включает итоговые метки курсивом.
Любая помощь будет высоко ценится.Заранее спасибо.
Текущий лесной участок выглядит так:
