Aprovação Ética e Consentimento para Participar
Este estudo foi uma análise de dados resumidos, disponíveis publicamente e desidentificados do Estudo Global Burden of Disease 2023 (GBD 2023). Como a pesquisa não envolveu interação direta com humanos ou animais, coleta de novos dados primários ou acesso a qualquer informação individualmente identificável, foi considerada isenta da exigência de aprovação formal por um conselho de revisão institucional (IRB) ou comitê de ética. O estudo foi conduzido de acordo com os princípios estabelecidos na Declaração de Helsinque. O uso dos dados GBD está em conformidade com seus termos de uso e políticas de acesso a dados.
Fontes de Dados
Este estudo é uma análise secundária de dados agregados desidentificados publicados no estudo GBD 2023. A ferramenta de consulta online GHDx (https://ghdx.healthdata.org/gbd–results–tool) foi empregada para obter dados sobre a carga da doença relacionados ao prolapso uterino na China e em países do G20 de 1990 a 2023. Casos de prolapso uterino foram identificados usando a hierarquia de causas GBD 2023 sob a categoria "Prolapso genital", que abrange o prolapso uterino. Os códigos ICD-10 correspondentes mapeados pelo Instituto de Métricas e Avaliação de Saúde (IHME) incluem N81.2 (prolapso uterovaginal incompleto), N81.3 (prolapso uterovaginal completo), N81.4 (prolapso uterovaginal não especificado) e N81.9 (prolapso genital feminino não especificado). A identificação de causa do GBD para prolapso uterino foi obtida diretamente da Ferramenta de Resultados do GBD, e as estimativas de carga da doença foram extraídas seguindo as definições padronizadas de causa do GBD e os procedimentos de codificação. A análise da carga da doença no GBD 2023 estimou a incidência, prevalência, morte, YLDs, YLLs e DALYs. O estudo examinou características epidemiológicas analisando diferenças geográficas e etárias na carga da doença ao longo do tempo e do espaço. Os dados foram processados usando análise EAPC (Estimada de Variação Percentual Anual), o modelo de ponto de junção, o modelo ARIMA, o modelo BAPC e análise de decomposição. Os métodos de cálculo específicos são consistentes com os relatados em estudosanteriores 27,33,34,35.
Análise EAPC
O EAPC foi uma medida amplamente aceita para quantificar a tendência dos ASRs em intervalos de tempo específicos, e foi calculado com base no modelo de regressão ajustado ao logaritmo natural das taxas. O modelo de regressão foi definido como: ln (taxa)=α+βx+ε, e o EAPC foi calculado como 100×(exp(β)- 1). O intervalo de confiança (IC) de 95% também foi determinado pelo modelo de regressão linear.
onde ln (taxa) é o logaritmo natural do ASR, x denota ano civil, α é o intercepto, β representa o coeficiente de inclinação e ε é o termo de erro. Se o EAPC > 0 e seu IC 95% também for > 0, há uma tendência significativa de alta. Se o EAPC < 0 e o IC 95% for < 0, há uma tendência significativa de queda. Se o IC de 95% inclui 0, a tendência não é estatisticamente significativa; ou seja, a mudança foi estável ao longo do tempo.
Análise de Ponto de Junção
O modelo de regressão Joinpoint foi empregado para calcular a variação percentual anual (APC) e a variação percentual anual média (AAPC), juntamente com seus respectivos intervalos de confiança (IC) de 95%, a fim de determinar as tendências de longo prazo das mudanças significativas para UP na China e nos países do G20 de 1990 a 2023. Esse modelo estabelece regressão segmentada baseada nas características temporais da distribuição da doença, dividindo a faixa de tempo em diferentes intervalos, cada um dos quais passa por ajuste e otimização de tendência, evitando efetivamente a subjetividade das análises típicas de tendências baseadas em tendências lineares. A direção da tendência foi determinada com base no AAPC calculado a partir do modelo final; especificamente, quando o IC de 95% não inclui 0, significa que a tendência é significativa, AAPC > 0 indica uma tendência de aumento, e APPC < 0 indica uma tendência decrescente.
Análise do BAPC
O modelo APC foi ampliado usando métodos estatísticos bayesianos (BAPC), que integram padrões históricos de dados, fatores de incerteza e conhecimento prévio para melhorar a precisão e robustez das predições, lidando assim de forma eficaz com ruído e ambiguidades nos dados. As projeções do BAPC para 2024 a 2038 foram criadas usando um parâmetro de suavização de 5 para efeitos de idade, período e coorte, com pesos padronizados baseados na população padrão mundial.
Análise ARIMA
O modelo ARIMA é um modelo popular em econometria que pode analisar o comportamento de séries temporais estacionárias e não estacionárias, bem como o impacto de planos e políticas em resultados específicos ao longo dotempo 34. No modelo ARIMA (p, d, q), 'p' representa a contagem de termos autoregressivos, 'd' indica a ordem de diferença, e 'q' denota a contagem de termos da média móvel. Os métodos de cálculo específicos são consistentes com os relatados em estudos anteriores. Para cada previsão, serão obtidos o Valor previsto, limite superior, limite inferior e IC de 95% de cada ano, assim como os parâmetros de desempenho do modelo de previsão (arquivo _eval.csv). O modelo de série temporal ARIMA foi usado para prever os próximos 27 anos (2024-2050), e o desempenho do modelo foi avaliado. O desempenho do modelo foi avaliado usando diagnósticos padrão de previsão, incluindo o Critério de Informação de Akaike (AIC), o Critério de Informação Bayesiano (BIC), análise de autocorrelação residual e avaliação de boa do ajuste. Os gráficos residuais não demonstraram violações significativas das suposições do modelo, apoiando a adequação dos modelos de previsão selecionados para a projeção de carga. O Auto.arima era usado para selecionar automaticamente os parâmetros ótimos do ARIMA(p,d,q) e o AICc (AIC modificado) era usado para selecionar o melhor modelo. Com base no método de predição bootstrap, o valor da predição e o intervalo de confiança foram calculados simulando os possíveis caminhos futuros.
Análises de decomposição
Realizamos uma decomposição de Das Gupta (fatoração de Kitagawa–Das Gupta) da variação nas contagens absolutas entre 1990 e 2023. A análise foi principalmente dividida em três componentes: o efeito do envelhecimento (mudanças na estrutura etária), o efeito do tamanho populacional (crescimento da população total) e o efeito de mudança epidemiológica (mudanças nas taxas específicas por idade). Os nomes das colunas de dados são os seguintes: 'overall_difference' indica a mudança total; 'a_effect' indica o efeito do envelhecimento; 'p_effect' representa o efeito do tamanho da população; 'r_effect' significa o efeito de mudança epidemiológica; 'a_percent', 'p_percent' e 'r_percent' denotam a contribuição percentual de cada efeito respectivo; 'val_1990' e 'val_2023' representam o número de casos em 1990 e 2023, respectivamente; e 'mudança' refere-se à diferença em 'diff1'. A magnitude de cada fator no gráfico resultante ilustra sua influência proporcional, com valores positivos e negativos indicando um efeito de aumento ou mitigação, respectivamente. Pontos pretos no gráfico servem como marcadores para a mudança total.
Software de análise
A análise estatística e visualização de dados neste estudo foram realizadas utilizando o pacote de software R e o programa Joinpoint (Consulte a Tabela de Materiais para detalhes). Um valor p < 0,05 foi considerado estatisticamente significativo.