Ei Éder, Hoje finalmente eu pude olhar e testar o exemplo que voce mandou e funcionou bem. Alguma pequenas modificações pra funcionar com os meus dados verdadeiros, Principalmente retirar a fixação do tamanho da janela e colocar essas timensões dentro da função tiff que salva o grafico. Eu percebi que a grande diferença entre o seu exemplo e as minhas tentativas frustradas eram o new=F dentro da forest que eu não estava usando, ao contrario eu estava tentando forçar o new = T, o que fez toda a diferença. Valeu mesmo. Abraço, Pedro Brasil 2015-04-16 9:10 GMT-03:00 Éder Comunello <comunello.eder@gmail.com>:
Pedro, bom dia!
A figura que segue (DoubleForest.png) é gerada pelo código colocado logo abaixo. A exibição on-screen só fica adequada no modo zoom do RStudio. Apesar de não aparecer corretamente na exibição da janela gráfica (windows/x11/quartz), quando direcionada pro dispositivo png() parece funcionar corretamente.
### <code r> if (!require("meta")) {install.packages("meta", dep=T); require("meta")}
b <- structure(list(ID2 = c("Lorca M-1992", "Carvalho MR-1993", "Hamerschlak N-1997","Oeleman WMR-1998", "Cabalero ZC-2007", "Otani MM-2009", "Otani MM-2009", "Remesar MC-2009", "Anes N-2010", "Flores Chaves M-2010"), VP = c(77, 224, 20, 515, 68, 167, 167, 340, 19, 66), VN = c(98, 182, 36, 481, 225, 257, 261, 1097, 15, 132), with = c(82, 224, 20, 539, 69, 168, 168, 356, 23, 66), without = c(98, 191, 40, 486, 238, 262, 262, 1099, 15, 157), group = c("ABBOTT-ABBOTT CHAGAS ELISA", "ABBOTT-ABBOTT CHAGAS ELISA", "ABBOTT-ABBOTT CHAGAS ELISA", "ABBOTT-ABBOTT CHAGAS ELISA", "BIOSCHILE-TEST ELISA CHAGAS", "ABBOTT-ABBOTT CHAGAS ELISA", "BIOSCHILE-TEST ELISA CHAGAS", "BIOSCHILE-TEST ELISA CHAGAS", "BIOSCHILE-TEST ELISA CHAGAS", "BIOSCHILE-TEST ELISA CHAGAS")), .Names = c("ID2","VP", "VN", "with", "without", "group"), row.names = c(50L, 42L,69L, 60L, 37L, 113L, 111L, 124L, 97L, 102L), class = "data.frame")
b$Q47="ELISA"
m.e.se <- metaprop(VP, with, studlab=ID2, data=b, subset=Q47=="ELISA", incr=.1, allincr=T, addincr=T, comb.fixed=F,hakn=T,title="ELISA for Chagas disease", complab='With Chagas', outclab="True positives", byvar=group, bylab='Trademark-name', print.byvar=T)
m.e.sp <- metaprop(VN, without, studlab=ID2, data=b, subset=Q47=="ELISA", incr=.1, allincr=T, addincr=T, comb.fixed=F,hakn=T,title="ELISA for Chagas disease", complab='With Chagas', outclab="True negatives", byvar=group, bylab='Trademark-name', print.byvar=T)
# parRS <- par(no.readonly=T)
require(grid) png(filename="DoubleForest.png", width=1200, height=460, units = "px", bg = "white") grid.newpage() grid.show.layout(l <- grid.layout(nrow=1, ncol=2, w=c(.6,.35), h=.9, default="npc"), vp.ex=.8)
# grid.show.layout(l <- grid.layout(nrow=1, ncol=2, w=c(1.6,1), h=1, default="npc", respect=T), vp.ex=.8)
vp0 <- viewport(layout=l)
pushViewport(vp0) grid.rect(gp=gpar(fill="white"))
pushViewport(viewport(layout.pos.col=1, layout.pos.row=1)) forest(new=F, m.e.se,overall=F,text.random="D&L random effects", pooled.events=T,xlab="Sensitivity", rightcols=F, leftcols=c("studlab", "event", "n","effect", "ci"), leftlabs=c("Study","TP","With","Sensitivity","95%-CI"), height=200) popViewport(1)
pushViewport(viewport(layout.pos.col=2, layout.pos.row=1)) forest(new=F,m.e.sp,overall=F,text.random="D&L random effects", pooled.events=T,xlab="Specificity", rightcols=c("effect", "ci"), leftcols=c("n","event"), leftlabs=c("TN","Without"), rightlabs=c("Specificity","95%-CI")) popViewport(1) dev.off()
file.info("DoubleForest.png")
### </code>
_______________________________________________ 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.