我使用ggplot制作森林情节。 这是我的样本数据
D<–data.frame(metabolite,expb,lower95,upper95)
$ metabolite: Factor w/ 70 levels "IDL.C_proc","IDL.CE_proc",..: 69 66 67 68 70 59 56 57 58 60 ...
$ expbeta : num -1.355 0.545 0.3 1.224 -0.861 ...
$ lower95 : num -2.589 1.034 0.809 -0.346 -0.666 ...
$ upper95 : num -0.2617 0.0916 -0.0561 1.3656 -0.5855 ...
此处metabolite是包含变量名称的列，expb是expcoef，lower95，upper95是conf.ivals。
我使用ggplot制作森林情节
forestplot <- function(d, xlab="HR", ylab=""){
require(ggplot2)
p <- ggplot(d, aes(x=d$metabolite, y=d$expb, ymin=d$lower95, ymax=d$upper95)) +
geom_pointrange() +
coord_flip() +
geom_hline(aes(x=0,yintercept=1), lty=2) +
ylab(xlab) +
xlab(ylab) #switch because of the coord_flip() above
return(p)}
forestplot(d)
如何修改森林图以具有对等和pvalues的边表？
感谢您的任何提示。