Ethische Zustimmung und Zustimmung zur Teilnahme
Diese Studie war eine Analyse bestehender, öffentlich zugänglicher und nicht identifizierter Zusammenfassungsdaten aus der Global Burden of Disease Study 2023 (GBD 2023). Da die Forschung keine direkte Interaktion mit menschlichen oder tierischen Probanden beinhaltete, keine neuen Primärdaten sammelte oder auf individuell identifizierbare Informationen zugreifen konnte, galt sie als von der Pflicht zur formellen Genehmigung durch ein institutionelles Überprüfungsgremium (IRB) oder Ethikkomitee ausgenommen. Die Studie wurde gemäß den in der Helsinkier Erklärung dargelegten Grundsätzen durchgeführt. Die Nutzung von GBD-Daten entspricht deren Nutzungsbedingungen und Datenzugriffsrichtlinien.
Datenquellen
Diese Studie ist eine sekundäre Analyse der deidentifizierten aggregierten Daten, die in der GBD-2023-Studie veröffentlicht wurden. Das Online-Abfragetool GHDx (https://ghdx.healthdata.org/gbd–results–tool) wurde von 1990 bis 2023 verwendet, um Daten zur Krankheitslast im Zusammenhang mit Gebärmutterprolaps in China und G20-Ländern zu erhalten. Fälle von Gebärmutterprolaps wurden anhand der GBD 2023 Ursachenhierarchie unter der Kategorie "Genitalprolaps" identifiziert, die auch den Gebärmutterprolaps umfasst. Die entsprechenden ICD-10-Codes, die vom Institute for Health Metrics and Evaluation (IHME) kartiert wurden, umfassen N81.2 (unvollständiger Uterovaginalprolaps), N81.3 (vollständiger Uterovaginalprolaps), N81.4 (nicht näher bezeichneter Uterovaginalprolaps) und N81.9 (nicht näher bezeichneter weiblicher Genitalprolaps). Die GBD-Ursachenerkennung für einen Gebärmutterprolaps wurde direkt aus dem GBD Results Tool gewonnen, und die Schätzungen der Krankheitslast wurden anhand der standardisierten GBD-Ursachendefinitionen und Codierungsverfahren extrahiert. Die Krankheitslastanalyse in GBD 2023 schätzte Inzidenz, Prävalenz, Tod, YLDs, YLLs und DALYs. Die Studie untersuchte epidemiologische Merkmale, indem geografische und altersgruppenspezifische Unterschiede in der Krankheitslast über Zeit und Raum analysiert wurden. Die Daten wurden mittels EAPC-Analyse (Estimated Annual Percent Change), dem Joinpoint-Modell, dem ARIMA-Modell, dem BAPC-Modell und der Dekompositionsanalyse verarbeitet. Die spezifischen Berechnungsmethoden stimmen mit denen überein, die in früheren Studien 27,33,34,35 berichtet wurden.
EAPC-Analyse
EAPC war ein weithin anerkanntes Maß zur Quantifizierung des Trends von ASRs über bestimmte Zeitintervalle und wurde auf Basis des Regressionsmodells berechnet, das dem natürlichen Logarithmus der Raten angepasst war. Das Regressionsmodell wurde definiert als: ln (Rate) = α+βx+ε, und EAPC wurde berechnet als 100×(exp(β)- 1). Das 95%-Konfidenzintervall (KI) wurde ebenfalls durch das lineare Regressionsmodell bestimmt.
wobei ln (Rate) der natürliche Logarithmus von ASR ist, x das Kalenderjahr bezeichnet, α der Intercept ist, β den Steigungskoeffizienten darstellt und ε der Fehlerterm ist. Wenn EAPC 0 > und sein 95%-KI ebenfalls > 0 ist, gibt es einen deutlichen Aufwärtstrend. Wenn der EAPC 0 < und das 95%-KI < 0 liegt, gibt es einen deutlichen Abwärtstrend. Wenn das 95%-KI 0 enthält, ist der Trend statistisch nicht signifikant; das heißt, die Veränderung war über die Zeit stabil.
Joinpoint-Analyse
Das Joinpoint-Regressionsmodell wurde verwendet, um die jährliche prozentuale Veränderung (APC) und die durchschnittliche jährliche prozentuale Veränderung (AAPC) zusammen mit ihren jeweiligen 95%-Konfidenzintervallen (KI) zu berechnen, um die langfristigen Trends signifikanter Veränderungen für UP in China und G20-Ländern von 1990 bis 2023 zu bestimmen. Dieses Modell etabliert eine segmentierte Regression basierend auf den zeitlichen Eigenschaften der Krankheitsverteilung und teilt den Zeitbereich in verschiedene Intervalle auf, von denen jedes Trendanpassung und Optimierung durchläuft, wodurch die Subjektivität typischer Trendanalysen auf Basis linearer Trends effektiv vermieden werden. Die Trendrichtung wurde auf Basis des AAPC bestimmt, das aus dem endgültigen Modell berechnet wurde; insbesondere bedeutet das, wenn der 95%-KI 0 nicht enthält, dass der Trend signifikant ist, AAPC > 0 einen steigenden Trend, und APPC < 0 einen abnehmenden Trend.
BAPC-Analyse
Das APC-Modell wurde mit bayesschen statistischen Methoden (BAPC) erweitert, die historische Datenmuster, Unsicherheitsfaktoren und Vorwissen integrieren, um die Vorhersagegenauigkeit und Robustheit zu verbessern und so Rauschen und Unklarheiten in den Daten effektiv zu bewältigen. BAPC-Prognosen für 2024 bis 2038 wurden mit einem Glättungsparameter von 5 für Alter, Periode und Kohorteneffekte erstellt, mit standardisierten Gewichtungen basierend auf der weltweiten Standardpopulation.
ARIMA-Analyse
Das ARIMA-Modell ist ein beliebtes Modell in der Ökonometrie, das das Verhalten stationärer und nicht-stationärer Zeitreihen sowie die Auswirkungen von Plänen und Richtlinien auf bestimmte Ergebnisse über die Zeit analysieren kann34. Im ARIMA-(p, d, q)-Modell steht 'p' für die Anzahl der autoregressiven Terme, 'd' für die Differenzierungsordnung und 'q' für die Anzahl der gleitenden Durchschnittsterme. Die spezifischen Berechnungsmethoden stimmen mit denen in früheren Studien überein. Für jede Vorhersage werden der vorhergesagte Wert, der obere Grenzwert, das untere Limit und der 95%-CI jedes Jahres sowie die Leistungsparameter des Vorhersagemodells ermittelt (_eval.csv Datei). Das ARIMA-Zeitreihenmodell wurde verwendet, um die nächsten 27 Jahre (2024–2050) vorherzusagen, und die Modellleistung wurde bewertet. Die Modellleistung wurde mittels standardisierter Prognosediagnostik bewertet, darunter das Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), Residual-Autokorrelationsanalyse und Goodness-of-Fit-Bewertung. Residualdiagramme zeigten keine wesentlichen Verstöße gegen die Modellannahmen, was die Eignung der ausgewählten Prognosemodelle für die Lastprojektion unterstützt. Auto.arima wurde verwendet, um automatisch die optimalen ARIMA(p,d,q)-Parameter auszuwählen, und AICc (modifizierte AIC) wurde verwendet, um das beste Modell auszuwählen. Basierend auf der Bootstrap-Vorhersagemethode wurden der Vorhersagewert und das Konfidenzintervall durch Simulation möglicher zukünftiger Pfade berechnet.
Zerlegungsanalysen
Wir führten eine Das-Gupta-Zerlegung (Kitagawa–Das Gupta-Faktorisierung) der Änderung der absoluten Zählungen zwischen 1990 und 2023 durch. Die Analyse wurde hauptsächlich in drei Komponenten unterteilt: den Alterungseffekt (Veränderungen der Altersstruktur), den Effekt der Bevölkerungsgröße (Wachstum der Gesamtbevölkerung) und den epidemiologischen Veränderungseffekt (Veränderungen altersspezifischer Raten). Die Datenspaltennamen lauten wie folgt: 'overall_difference' bezeichnet die Gesamtänderung; 'a_effect' steht für den Alterungseffekt; 'p_effect' steht für den Effekt der Populationsgröße; 'r_effect' steht für den epidemiologischen Wandeleffekt; 'a_percent', 'p_percent' und 'r_percent' bezeichnen den prozentualen Beitrag jedes jeweiligen Effekts; 'val_1990' und 'val_2023' stehen für die Anzahl der Fälle in den Jahren 1990 bzw. 2023; und 'Veränderung' bezieht sich auf die Differenz in 'Diff1'. Die Größe jedes Faktors im resultierenden Diagramm veranschaulicht seinen proportionalen Einfluss, wobei positive und negative Werte jeweils einen Steigerungs- oder Minderungseffekt anzeigen. Schwarze Punkte auf dem Plot dienen als Markierungen für die vollständige Veränderung.
Analysesoftware
Die datenstatistische Analyse und Visualisierung in dieser Studie wurden mit dem R-Softwarepaket und dem Softwareprogramm Joinpoint durchgeführt (siehe Materialtabelle für Details). Ein p-Wert < 0,05 galt als statistisch signifikant.