Vinicius, eu pesquisei os pacotes que fazem o forest plot. Parece que existem vários. No entanto, vendo a dinâmica das sintaxes, dei preferência para o metafor.
O metafor é um pacote bastante robusto e permite a construação manual do forest plot. Sendo assim, a primeira coisa que deve ser feita é ativar o pacote metafor, utilizando a sintaxe abaixo:
## Comando para ativar o pacote metafor
library(metafor)
## Loading required package: Matrix
## Loading 'metafor' package (version 2.1-0). For an overview
## and introduction to the package please type: help(metafor).
Agora que baixamos o pacote, o próximo passo é avisar para ele quais os dados que devem ser apresentados no forest plot. Nós podemos fazer isso de duas formas: a primeira é importar o banco de dados no mesmo formato que utilizamos no pacote meta e, a segunda forma é criar manualmente um banco de dados no R.
Para este exemplo, vou utilizar a segunda forma:
## Criando o data.frame (banco de dados)
covid.prev <- c(18.2, 14.9, 23.0, 18.7)
lower <- c(12.9, 9.9, 21, 14.7)
upper <- c(19.6, 17.2, 25.2, 22.7)
autor <- c("Revisão 1", "Revisão 2", "Revisão 3",
"Revisão 4")
Veja que no exemplo, temos os resultados de 4 estudos. Veja que nós não estamos preocupados em retirar os estudos que se repetem. Nós apenas estamos apresentando os resultados dos estudos no gráfico!
Vamos agora criar o Forest Plot desses dados que eu simulei com a sintaxe, a seguir:
forest(x = covid.prev,
ci.lb = lower,
ci.ub = upper,
slab = autor,
xlab = "Prevalence of...",
refline = 16, psize = 1.2, cex = .8) ## Eu coloquei a linha de referência em 16
text(35.5, 5.5, "Prevalence [95%-IC]", cex = 0.8)
text(-49, 7.5, "Outcomes")
O único trabalho que teremos é o de colocar a unidade do efeito em cada forest. Isso vai mudar de forest plot para forest plot.