Dear R users,
I'm trying figuring out with forest plot and rmeta. Here's my code...
library(rmeta)
tabletext <- structure(c("", NA, "Sex", "[F]",
"[M]", NA,
"Age class", "(0,60]",
"(60,80]", NA,
"Karnofsky class", "[70;90]",
"(90;100]", NA, "Tumor dimension",
"(0,5]", "(5,10]",
"(10,30]", NA,
"Primary tumor", "[T1]",
"[T2]",
"[T3a]", "[T3b-T3c]", NA,
"Histology",
"[ClearCell]", "[Papill-Chromo]",
"[Others]", NA, "Grading",
"[G1]",
"[G2]", "[G3]", "[G4]",
NA, "n", NA, NA,
"105", "198", NA, NA,
"161", "142", NA,
NA, "39", "264", NA, NA,
"121", "143",
"39", NA, NA, "131",
"51", "68", "53",
NA, NA, "242", "34",
"19", NA, NA, "53",
"136", "68", "16", NA,
"HR", NA, NA,
"0.67", "0.92", NA, NA,
"0.6", "0.98",
NA, NA, "1.19", "0.72", NA, NA,
"1.19",
"0.5", "1.11", NA, NA,
"0.82", "1.16",
"0.56", "1.21", NA, NA,
"0.82", "0.5",
"0.95", NA, NA, "0",
"0.52", "0.96",
"1.17", NA, "95 CI", NA, NA,
"[0.31-1.44]", "[0.53-1.62]", NA,
NA,
"[0.31-1.18]", "[0.53-1.82]", NA,
NA,
"[0.44-3.21]", "[0.43-1.21]", NA,
NA,
"[0.42-3.41]", "[0.27-0.93]",
"[0.44-2.76]", NA, NA,
"[0.28-2.46]",
"[0.39-3.46]", "[0.26-1.25]",
"[0.53-2.76]", NA, NA,
"[0.49-1.37]",
"[0.06-4.31]", "[0.27-3.39]", NA,
NA,
"[0-Inf]", "[0.25-1.1]",
"[0.44-2.07]",
"[0.31-4.38]", NA, "P-value", NA,
"0.46", NA, NA, NA, "0.3", NA, NA,
NA,
"0.39", NA, NA, NA, "0.21", NA, NA,
NA,
NA, "0.63", NA, NA, NA, NA, NA,
"0.85",
NA, NA, NA, NA, "0.22", NA, NA, NA, NA,
NA), .Dim = c(36L, 5L))
m <- c(NA, NA, NA, 0.67, 0.92, NA, NA, 0.6, 0.98, NA, NA, 1.19,
0.72, NA, NA, 1.19, 0.5, 1.11, NA, NA, 0.82, 1.16, 0.56, 1.21,
NA, NA, 0.82, 0.5, 0.95, NA, NA, 0, 0.52, 0.96, 1.17, NA)
l <- c(NA, NA, NA, 0.31, 0.53, NA, NA, 0.31, 0.53, NA, NA, 0.44,
0.43, NA, NA, 0.42, 0.27, 0.44, NA, NA, 0.28, 0.39, 0.26, 0.53,
NA, NA, 0.49, 0.06, 0.27, NA, NA, 0, 0.25, 0.44, 0.31, NA)
u <- c(NA, NA, NA, 1.44, 1.62, NA, NA, 1.18, 1.82, NA, NA, 3.21,
1.21, NA, NA, 3.41, 0.93, 2.76, NA, NA, 2.46, 3.46, 1.25,
2.76, NA, NA, 1.37, 4.31, 3.39, NA, NA, Inf, 1.1, 2.07,
4.38, NA)
forestplot(labeltext=tabletext,
m=m,
l=l,
u=u,
zero=1,
clip=c(0, 3.5),
align=c("r",rep("c",ncol(tabletext)-1)),
xlab="Hazard Ratio",
xlog=FALSE,
col=meta.colors(box="royalblue",line="darkblue",
summary="royalblue"),
boxsize=.4
)
Not considering for the moment a few problems with grading
section estimation, on my linux box labeltext columns overlap in the
graph produced;
is there a way to prevent this? What did i miss?
Thank you
Luca