Plot Forest

1. Inputs

index=c(1:5)
label=c("LDL-c", "LDL-p", "Large LDL-p", "Small LDL-p", "IDL")
study=c("5 LDL-c", "4 LDL-p", "3 Large LDL-p", "2 Small LDL-p", "1 IDL")

OR= c(1.44, 1.46, 1.08, 1.34, 1.06)
LL= c(1.20, 1.22, 0.91, 1.12, 0.90)
UL= c(1.72, 1.75, 1.29, 1.60, 1.25)
dat=data.frame(index, label, study, OR, LL, UL)

2. Code by layer

library(ggplot2)
p = ggplot(data=dat, aes(y=index, x=OR, label=index, col=index))
p = p + geom_point(size=4, shape=19) + geom_errorbarh(aes(xmin=LL, xmax=UL), height=0.3, size=2.0)
p = p + geom_vline(xintercept=1, linetype="dashed")
p = p + xlab("Odds Ratio") + ylab("") + theme_bw() + theme(legend.position="none")
p + scale_y_continuous(name = "", breaks=1:5, trans = "reverse", labels = dat$label)   

3. Simple code

ggplot(data=dat, aes(y=index, x=OR, label=index, col=index))+geom_point(size=2, shape=19) + geom_errorbarh(aes(xmin=LL, xmax=UL), height=0.2, size=1.0)+ geom_vline(xintercept=1, linetype="dashed", col="red")+scale_y_continuous(name = "", breaks=1:5, labels = dat$label, trans = "reverse")+ xlab("Odds Ratio") + ylab("") + theme_bw() + theme(legend.position="none")

4.

# And photoshop
p = ggplot(data=dat, aes(y=study, x=OR, label=study, col=study))
p = p + geom_point(size=4, shape=19) + 
   geom_errorbarh(aes(xmin=LL, xmax=UL), height=0.3, size=2.0)
p = p + geom_vline(xintercept=1, linetype="dashed", color = "red")
p + xlab("ORs for the presence of ICAS per 1SD higher of lipid indices") + ylab("") + theme_bw() + 
   theme(legend.position="none")