Message-ID: <14461981.post@talk.nabble.com>
Date: 2007-12-21T20:35:29Z
From: francogrex
Subject: Forestplot
In-Reply-To: <14404133.post@talk.nabble.com>
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.