Capítulo 50 Meta-análise

50.1 Características

50.1.1 O que é meta-análise?

• Meta-análise é um método estatístico que combina quantitativamente os resultados de múltiplos estudos independentes sobre uma mesma questão de pesquisa, aumentando o poder estatístico e a precisão das estimativas de efeito.³⁸⁹

• Meta-análise sintetiza evidências considerando o peso de cada estudo (geralmente inversamente proporcional à variância) e permite avaliar a consistência dos resultados, identificar fontes de heterogeneidade e estimar um efeito global.³⁸⁹

50.2 Modelos de meta-análise

• Modelo de efeitos fixos: assume que todos os estudos avaliam o mesmo efeito verdadeiro, e a variação observada é apenas devido ao erro de amostragem. É adequado quando os estudos são homogêneos e as diferenças entre eles são pequenas.³⁸⁹

• Modelo de efeitos aleatórios: assume que os estudos avaliam efeitos verdadeiros diferentes, com uma distribuição normal. É mais apropriado quando há heterogeneidade entre os estudos, pois considera a variação entre eles além do erro de amostragem.³⁸⁹

• Modelo de efeitos de rede: estende a meta-análise para comparar múltiplas intervenções simultaneamente, mesmo que não tenham sido comparadas diretamente em estudos. É útil para avaliar a eficácia relativa de várias intervenções.^REF?

O pacote metafor³⁹⁰ fornece a função forest para criar figuras tipo forest plot.

Figura 50.1: Comparação entre modelos de efeito fixo e aleatório com 10 ensaios clínicos simulados.

50.3 Conversão de Medidas em Meta-análises

50.3.1 O que fazer quando os estudos apresentam resultados com diferentes parâmetros?

• Quando os estudos reportam médias e desvios-padrão, os dados podem ser usados diretamente na metanálise.^REF?

• Quando apresentam mediana e intervalo interquartil (ou mínimo–máximo), existem métodos estatísticos para converter em média e DP.³⁹¹

• Hozo et al. (2005) propuseram fórmulas para estimar a média e o desvio-padrão a partir da mediana, amplitude e tamanho da amostra.³⁹¹

• Wan et al. (2014) aperfeiçoaram essas estimativas, oferecendo métodos mais precisos para converter mediana e IQR em média e DP.³⁹²

O pacote metafor³⁹⁰ fornece a função conv.fivenum para converter mínimo/mediana/máximo ou Q1/mediana/Q2 em média e desvio-padrão.

50.4 Interpretação de efeitos em meta-análise

50.4.1 Como avaliar a variação do tamanho do efeito?

• O intervalo de predição contém informação sobre a variação do tamanho do efeito.³⁹³

• Se o intervalo de predição não contém a hipótese nula (\(H_{0}\)), podemos concluir que (a) o tratamento funciona igualmente bem em todas as populações, ou que ele funciona melhor em algumas populações do que em outras.³⁹³

• Se o intervalo de predição contém a hipótese nula (\(H_{0}\)), podemos concluir que o tratamento pode ser benéfico em algumas populações, mas prejudicial em outras, de modo que a estimativa pontual (geralmente a média) torna-se amplamente irrelevante. Nesse caso, é recomendado investigar em que populações o tratamento seria benéfico e em quais causaria danos.³⁹³

50.4.2 Como avaliar a heterogeneidade entre os estudos?

• A heterogeneidade - variação não-aleatória - no efeito do tratamento entre os estudos incluídos em uma meta-análise pode ser avaliada pelo \(I^{2}\).^393,394

• \(I^{2}\) representa qual proporção da variância observada reflete a variância nos efeitos verdadeiros em vez do erro de amostragem.³⁹³

• \(I^{2}\) mede a proporção da variância total que pode ser atribuída à heterogeneidade entre os estudos incluídos.³⁹⁴

• \(I^{2}\) não depende da quantidade de estudos incluídos na meta-análise. Entretanto, \(I^{2}\) aumenta com a quantidade de participantes incluídos nos estudos meta-analisados.³⁹⁴

• A heterogeneidade entre estudos é explicada de modo mais confiável utilizando dados de pacientes individuais, uma vez que a direção verdadeira da modificação de efeito não pode ser observada a partir de dados agregados no estudo.³⁹⁵

O pacote psychmeta²⁷⁸ fornece a função ma_d para meta-analisar valores d.

O pacote psychmeta²⁷⁸ fornece a função ma_r para meta-analisar correlações.

50.5 Interpretação do Forest Plot

50.5.1 O que é um forest plot?

• Um forest plot é uma representação gráfica dos achados de uma meta-análise. Ele resume os resultados de estudos individuais e apresenta uma estimativa combinada do efeito, permitindo interpretação visual da magnitude, direção e significância estatística dos resultados.³⁹⁶

Figura 50.2: Forest plot de uma meta-análise de efeito fixo com 10 ensaios clínicos simulados.

Figura 50.3: Forest plot de uma meta-análise de efeito aleatório com 10 ensaios clínicos simulados.

50.5.2 2. Quais são as seis colunas básicas que um forest plot geralmente apresenta?

• As seis colunas básicas incluem: estudos incluídos (e subgrupos, se analisados); dados do grupo de intervenção, dados do grupo controle; peso de cada estudo; medida numérica do efeito; representação gráfica do efeito.³⁹⁶

50.5.3 Como diferenciar um desfecho binário de um contínuo em um forest plot?

• Em desfechos binários, são mostrados número de eventos e total da amostra, sendo o efeito medido por risk ratio (\(RR\)) ou odds ratio (\(OR\)).³⁹⁶

• Em desfechos contínuos, apresentam-se médias, desvios-padrão e tamanhos amostrais, com o efeito medido pela diferença de médias.³⁹⁶

50.5.4 4. O que representa o ponto central da caixa e o tamanho desta no gráfico?

• O ponto central indica a estimativa pontual do efeito (melhor estmativa para o efeito real).³⁹⁶

• O tamanho da caixa é proporcional ao peso do estudo na meta-análise, geralmente maior para estudos com amostras maiores.³⁹⁶

50.5.5 5. Qual é o significado da linha vertical do “nenhum efeito”?

• É a linha de referência que indica efeito nulo.³⁹⁶

• Para desfechos binários, corresponde ao valor 1 (\(RR\) ou \(OR\) = 1).³⁹⁶

• Para desfechos contínuos, corresponde ao valor 0 (diferença de médias = 0).³⁹⁶

• Se o intervalo de confiança de um estudo ou do resultado combinado cruza essa linha, o resultado não é estatisticamente significativo.³⁹⁶

50.5.6 6. Como interpretar o diamante na parte inferior do forest plot?

• O diamante representa o efeito combinado dos estudos incluídos.³⁹⁶

• O ponto central do diamante é a estimativa global.³⁹⁶

• A largura do diamante representa o intervalo de confiança de 95% para o efeito combinado.³⁹⁶

50.5.7 Como a heterogeneidade pode ser avaliada no forest plot?

• A variabilidade nos resultados dos estudos incluídos é avaliada pela sobreposição dos intervalos de confiança dos estudos; pelo teste do qui-quadrado (\(chi^2\)) e pelo valor de \(I^{2}\).³⁹⁶

50.5.8 Quais são as interpretações usuais para os valores de \(I^{2}\)?

• 0% a 40%: pode não ser importante; 30% a 60%: heterogeneidade moderada; 50% a 90%: heterogeneidade substancial; 75% a 100%: heterogeneidade considerável.³⁹⁶

50.6 Viés de publicação em meta-análises

50.6.1 O que é viés de publicação?

• O viés de publicação ocorre quando estudos com resultados não significativos ou contrários à hipótese tendem a não ser publicados, afetando a estimativa final da meta-análise e podendo levar a conclusões incorretas.³⁹⁷

50.6.2 Quais métodos podem ser usados para identificar viés de publicação?

• O método mais simples é o funnel plot, que representa a estimativa de efeito de cada estudo em função de sua precisão (\(1/SE\)).³⁹⁸

• Na ausência de viés, espera-se uma distribuição simétrica (“forma de funil”). Assimetria pode indicar viés de publicação, heterogeneidade entre estudos ou efeitos de tamanho de estudo (small-study effects).³⁹⁸

• Para odds ratios (\(OR\)), a correlação entre \(ln(OR)\) e seu erro padrão pode gerar assimetria mesmo sem viés, por isso recomenda-se, nesses casos, plotar em função do tamanho amostral.³⁹⁹

50.6.3 O que é um gráfico de funil (funnel plot)?

• É um gráfico de dispersão que relaciona a estimativa de efeito de cada estudo com uma medida de seu tamanho ou precisão (por exemplo, erro-padrão no eixo vertical, em escala invertida). Em condições ideais (ausência de viés e heterogeneidade), os estudos se distribuem de forma simétrica, formando um “funil invertido”.⁴⁰⁰

Figura 50.4: Gráfico de funil de meta-análise com 10 ensaios clínicos simulados.

50.6.4 A assimetria do funnel plot indica sempre viés de publicação?

• Viéses de relato (reporting biases), como viés de publicação, viés de linguagem ou de citação.⁴⁰⁰

• Diferenças metodológicas entre estudos pequenos e grandes.⁴⁰⁰

• Heterogeneidade verdadeira (diferença real no efeito conforme o tamanho ou o contexto do estudo).⁴⁰⁰

• Artefatos estatísticos ou mero acaso.⁴⁰⁰

50.6.5 O que é trim and fill?

• O método trim and fill “apara” (trim) os estudos que causam assimetria no funnel plot, estima o número de estudos possivelmente ausentes (não publicados) e “preenche” (fill) o gráfico com esses estudos simulados, recalculando o efeito combinado.⁴⁰¹

• O método assume que a assimetria é causada unicamente por viés de publicação, podendo levar a conclusões equivocadas quando há outras causas, como heterogeneidade.⁴⁰¹

50.6.6 O que é o teste de Egger?

• É um teste estatístico amplamente utilizado que avalia a relação entre o efeito padronizado (\(efeito/SE\)) e a precisão (\(1/SE\)).³⁹⁸

• No entanto, para meta-análises com \(OR\), apresenta taxas de erro tipo I excessivas, especialmente quando o efeito é grande ou há alta heterogeneidade.³⁹⁹

50.6.7 O que é o teste de Peters?

• Uma regressão linear ponderada com \(ln(OR)\) como variável dependente e o inverso do tamanho total da amostra como variável independente (modificação do teste de Macaskill).³⁹⁹

• Essa abordagem reduz a correlação entre \(ln(OR)\) e seu \(SE\), resultando em taxas de erro tipo I mais adequadas (~10%) independentemente do tamanho do efeito, número de estudos ou heterogeneidade.³⁹⁹

• O teste de Peters é preferível ao teste de Egger quando o desfecho é expresso como OR, pois mantém taxas de erro tipo I adequadas e ainda apresenta poder comparável para detectar viés em condições de baixa heterogeneidade.³⁹⁹

• Em casos de alta heterogeneidade, o teste de Egger pode apresentar maior poder, mas sua alta taxa de falsos positivos compromete a interpretação.³⁹⁹

50.6.8 Quais são as recomendações para testar a assimetria?

• Evitar testes quando há menos de 10 estudos, devido ao baixo poder estatístico.⁴⁰⁰

• Inspecionar visualmente o funnel plot junto com os resultados dos testes.⁴⁰⁰

• Para desfechos contínuos (diferença de médias), o teste de Egger pode ser usado.⁴⁰⁰

• Para desfechos dicotômicos expressos como odds ratio (\(OR\)) com baixa heterogeneidade (\(\tau^2 < 0,1\)), considerar os testes de Harbord, Peters ou Rücker.⁴⁰⁰

• Para desfechos dicotômicos expressos como odds ratio(\(OR\)) com alta heterogeneidade (\(\tau^2 > 0,1\)), o teste de Rücker com transformação \(arcsine\) é mais indicado.⁴⁰⁰

50.6.9 Como interpretar os resultados de testes de viés de publicação?

• Um resultado não significativo não garante ausência de viés.³⁹⁹

• Recomenda-se complementar com inspeção visual do funnel plot e considerar métodos adicionais como trim and fill.^399,401

• Quando há suspeita de viés, discutir as implicações e interpretar o efeito combinado com cautela.³⁹⁹

O pacote psychmeta²⁷⁸ fornece a função plot_funnel para criar figuras tipo funnel plot.

50.7 Diretrizes para redação

50.7.1 Quais são as diretrizes para redação de meta-análises?

• Visite a rede Enhancing the QUAlity and Transparency Of health Research (EQUATOR Network) para encontrar diretrizes específicas para cada tipo de meta-análises.

  • The PRISMA 2020 statement: An updated guideline for reporting systematic reviews:⁴⁰² https://www.equator-network.org/reporting-guidelines/prisma/

O pacote metagear⁴⁰³ fornece a função plot_PRISMA para gerar o fluxograma de uma revisão sistemática de acordo com o Preferred Reporting Items for Systematic Reviews and Meta-Analyses⁴⁰⁴.

O pacote PRISMA2020⁴⁰⁵ fornece a função PRISMA_flowdiagram para elaboração do fluxograma de revisões sistemáticas no formato padrão.

Citar como:
Ferreira, Arthur de Sá. Ciência com R: Perguntas e respostas para pesquisadores e analistas de dados. Rio de Janeiro: 1a edição,

Referências

278.

Dahlke JA, Wiernik BM. psychmeta: An R Package for Psychometric Meta-Analysis. Applied Psychological Measurement. 2018;43(3):415–416. doi:10.1177/0146621618795933

389.

Borenstein M, Hedges LV, Higgins JPT, Rothstein HR. A basic introduction to fixed-effect and random-effects models for meta-analysis. Research Synthesis Methods. 2010;1(2):97–111. doi:10.1002/jrsm.12

390.

Viechtbauer W. Conducting meta-analyses in R with the metafor package. Vol 36.; 2010. doi:10.18637/jss.v036.i03

391.

Hozo SP, Djulbegovic B, Hozo I. Estimating the mean and variance from the median, range, and the size of a sample. BMC Medical Research Methodology. 2005;5(1). doi:10.1186/1471-2288-5-13

392.

Wan X, Wang W, Liu J, Tong T. Estimating the sample mean and standard deviation from the sample size, median, range and/or interquartile range. BMC Medical Research Methodology. 2014;14(1). doi:10.1186/1471-2288-14-135

393.

Borenstein M. In a meta-analysis, the I-squared statistic does not tell us how much the effect size varies. Journal of Clinical Epidemiology. outubro 2022. doi:10.1016/j.jclinepi.2022.10.003

395.

Grooth HJ de, Parienti JJ. Heterogeneity between studies can be explained more reliably with individual patient data. Intensive Care Medicine. julho 2023. doi:10.1007/s00134-023-07163-z

396.

Dettori JR, Norvell DC, Chapman JR. Seeing the Forest by Looking at the Trees: How to Interpret a Meta-Analysis Forest Plot. Global Spine Journal. 2021;11(4):614–616. doi:10.1177/21925682211003889

397.

Song, Eastwood, Gilbody, Duley, Sutton. Publication and related biases. Health Technology Assessment. 2000;4(10). doi:10.3310/hta4100

398.

Egger M, Smith GD, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315(7109):629–634. doi:10.1136/bmj.315.7109.629

399.

Peters JL. Comparison of Two Methods to Detect Publication Bias in Meta-analysis. JAMA. 2006;295(6):676. doi:10.1001/jama.295.6.676

400.

Sterne JAC, Sutton AJ, Ioannidis JPA, et al. Recommendations for examining and interpreting funnel plot asymmetry in meta-analyses of randomised controlled trials. BMJ. 2011;343(jul22 1):d4002–d4002. doi:10.1136/bmj.d4002

401.

Duval S, Tweedie R. Trim and Fill: A Simple Funnel-PlotBased Method of Testing and Adjusting for Publication Bias in Meta-Analysis. Biometrics. 2000;56(2):455–463. doi:10.1111/j.0006-341x.2000.00455.x

402.

Page MJ, McKenzie JE, Bossuyt PM, et al. The PRISMA 2020 statement: An updated guideline for reporting systematic reviews. PLOS Medicine. 2021;18(3):e1003583. doi:10.1371/journal.pmed.1003583

403.

Lajeunesse MJ. Facilitating systematic reviews, data extraction, and meta-analysis with the metagear package for R. Methods in Ecology and Evolution. 2016;7(3):323–330. doi:10.1111/2041-210X.12472

404.

Moher D, Shamseer L, Clarke M, et al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015 statement. Systematic Reviews. 2015;4(1). doi:10.1186/2046-4053-4-1

405.

Haddaway NR, Page MJ, Pritchard CC, McGuinness LA. PRISMA2020: An R package and Shiny app for producing PRISMA 2020-compliant flow diagrams, with interactivity for optimised digital transparency and Open Synthesis. Campbell Systematic Reviews. 2022;18:e1230. doi:10.1002/cl2.1230