[R] Forestplot

[francogrex]

Answering my own query here's a solution (if any reader was curious):

forestplot=function(mean,std,conf,threshold){
z=-qnorm((1-conf)/2)
CI.H <- mean+z*sqrt(std) # Calculate upper CI 
CI.L <- mean-z*sqrt(std) # Calculate lower CI 
plot(mean,1:length(mean),xlim=c(min(CI.L),max(CI.H)),pch=15,cex=1,
ylab="",xlab="mean & CIs",main="Forest plot")
arrows(mean,1:length(mean),CI.H,1:length(mean), angle=90,length=0.05)
arrows(mean,1:length(mean),CI.L,1:length(mean), angle=90,length=0.05)
abline(v=threshold,lty=3)
}


francogrex wrote:
> 
> I know there is a function forestplot from rmeta package and also the
> plot.meta from the meta package and maybe others, but they are rather
> complicated with extra plot parameters that I do not need and also they
> process only objects created with other package functions. 
> But I wonder if anyone has a much simpler function using the basic plot to
> make a forestplot with only a median (or mean) and just the confidence
> intervals...
> 

-- 
View this message in context: http://www.nabble.com/Forestplot-tp14404133p14461981.html
Sent from the R help mailing list archive at Nabble.com.