install.packages("rmeta")
library(rmeta)
data(cochrane)
svg("cochrane_logo_data.svg",width=9,height=4.35) 
	steroid <- meta.MH(n.trt, n.ctrl, ev.trt, ev.ctrl,
		names=name, data=cochrane)
	tabletext<-cbind(
	c("","Study",steroid$names,NA,"Summary"), #Just tags, not single-author papers
#	c("","Year",c(1111,1111,1111,1111,1111,1111,1111),NA,NA), #Don't have this data
	c("Deaths","(corticosteroid)",cochrane$ev.trt,NA,NA),
	c("Deaths","(placebo)", cochrane$ev.ctrl, NA,NA),
	c("","Odds ratio",format(exp(steroid$logOR),digits=2),NA,format(exp(steroid$logMH),digits=2))
	)
	m<- c(NA,NA,steroid$logOR,NA,steroid$logMH)		#Mean
	l<- m-c(NA,NA,steroid$selogOR,NA,steroid$selogMH)*2	#Lower
	u<- m+c(NA,NA,steroid$selogOR,NA,steroid$selogMH)*2	#Upper
	forestplot(tabletext,m,l,u,zero=0,is.summary=c(TRUE,TRUE,rep(FALSE,8),TRUE),
	clip=c(log(0.015),log(3)), xlog=TRUE,
	xticks=c(0.01,0.1,1,2,4),
	col=meta.colors(box="royalblue",line="darkblue", summary="royalblue"))
	forestplot(tabletext,m,l,u,zero=0,is.summary=c(TRUE,TRUE,rep(FALSE,8),TRUE),
	clip=c(log(0.015),log(3)), xlog=TRUE, 
	#boxsize=0.75,
	xticks=c(0.01,0.1,1,2,4),
	xlab="\nOdds ratio with 95% confidence interval\n(1=no effect, <1=treatment has fewer deaths)",
	col=meta.colors(box="royalblue",line="darkblue", summary="royalblue"))
dev.off()