Boa tarde,
estou precisando fazer um gráfico chamado forest plot porem não posso
usar os já prontos no R. Já esta quase pronto utilizando este script.
require(gridExtra)
nomes <- sample(letters,5)
logOR <- abs(rnorm(5,0,.3))
sd_logOR <-rbeta(5,10,10)
OR <- round(exp(logOR),2)
lci_OR <- round(exp(logOR-1.96*sd_logOR),2)
uci_OR <- round(exp(logOR+1.96*sd_logOR),2)
OR <- paste(OR,' (',lci_OR,' - ',uci_OR,')',sep='')
df <- data.frame(nomes,OR,lci_OR,uci_OR)
p1 <- tableGrob(df[,1:2],cols=c('Estudos','OR (CI)'),row.just
='left',col.just ='center', core.just ='right')
require(ggplot2)
OR <- round(exp(logOR),2)
df <- data.frame(nomes,OR,lci_OR,uci_OR)
p2 <- p2 <- ggplot(df, aes(x=nomes, y=OR, ymin=lci_OR,
ymax=uci_OR))+geom_hline(yintercept=1,lty=2)+expand_limits(y =
c(0,max(uci_OR)*1.1))+geom_pointrange()+coord_flip()#+
geom_point(aes(size = sd_logOR))
grid.arrange(p1,p2,ncol=2)
Porém aqui no fim há uma desproporção entre os dados da tabela com os
dados do gráfico.
Porventura alguém sabe como resolver isto ?
[]s
Tura
_______________________________________________
R-br mailing list
[email protected]
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.
