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.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
