Ola Ary, De uma olhada nesse código e veja se atende Abr Humberto ####################################################################################### ##Odds Ratio to bycatch (group) library(meta) library(mvmeta) library(metafor) dat.bcg<-read.csv("meta.csv" , header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) summary(dat.bcg) names(dat.bcg) res <- metabin(tpos,tneg,cpos,cneg, data=dat.bcg, sm="OR", method="I",studlab=paste(esp)) res summary(res) pdf(file='Figure OR3.pdf', height=120, width=120, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=10) forest(res, xlim="s",comb.fixed=FALSE,xlab="Odds of Cacthability",text.fixed="Pooled effect estimate") dev.off() metainf(res, pooled="random") ############################################################################# Enviada do meu iPad
Em 23/07/2014, às 21:29, Ary Serpa Neto <aryserpa@terra.com.br> escreveu:
Olá amigos,
Estou precisando construir um forest plot (similar a imagem em anexo) usando "pooled incidence rate".
Alguem pode me dar uma luz de como fazer isso no R?
Desde ja agradeço a atenção
Grande abraço
Ary Serpa Neto, MD, MSc
<Imagem1.jpg> _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.