Ethical Approval and Consent to Participate
This study was an analysis of existing, publicly available, and de-identified summary data from the Global Burden of Disease Study 2023 (GBD 2023). As the research did not involve direct interaction with human or animal subjects, the collection of new primary data, or access to any individually identifiable information, it was deemed exempt from the requirement for formal approval by an institutional review board (IRB) or ethics committee. The study was conducted in accordance with the principles outlined in the Declaration of Helsinki. The use of GBD data complies with its terms of use and data access policies.
Data Sources
This study is a secondary analysis of de-identified aggregated data published in the GBD 2023 study. The GHDx online query tool (https://ghdx.healthdata.org/gbd–results–tool) was employed to obtain disease burden data related to uterine prolapse in China and G20 countries from 1990 to 2023. Uterine prolapse cases were identified using the GBD 2023 cause hierarchy under the category “Genital prolapse,” which encompasses uterine prolapse. The corresponding ICD-10 codes mapped by the Institute for Health Metrics and Evaluation (IHME) include N81.2 (incomplete uterovaginal prolapse), N81.3 (complete uterovaginal prolapse), N81.4 (unspecified uterovaginal prolapse), and N81.9 (unspecified female genital prolapse). The GBD Cause ID for uterine prolapse were obtained directly from the GBD Results Tool, and disease burden estimates were extracted following the standardized GBD cause definitions and coding procedures. The disease burden analysis in GBD 2023 estimated the incidence, prevalence, death, YLDs, YLLs, and DALYs. The study examined epidemiological traits by analyzing geographic and age-group differences in disease burden over time and space. Data were processed using EAPC (Estimated Annual Percentage Change) analysis, the joinpoint model, the ARIMA model, the BAPC model, and decomposition analysis. The specific calculation methods are consistent with those reported in previous studies27,33,34,35.
EAPC Analysis
EAPC was a widely accepted measure to quantify the trend of ASRs over specific time intervals, and it was calculated based on the regression model fitted to the natural logarithm of the rates. The regression model was defined as: ln (rate)=α+βx+ε, and EAPC was calculated as 100×(exp(β)- 1). The 95% confidence interval (CI) was also determined by the linear regression model.
where ln (rate)is the natural logarithm of ASR, x denotes calendar year, α is the intercept, β represents the slope coefficient, and ε is the error term. If EAPC > 0 and its 95% CI is also > 0, there is a significant upward trend. If EAPC < 0 and 95% CI is < 0, there is a significant downward trend. If the 95% CI includes 0, the trend is not statistically significant; that is, the change was stable over time.
Joinpoint Analysis
The Joinpoint regression model was employed to calculate the annual percentage change (APC) and the average annual percentage change (AAPC), along with their respective 95% confidence intervals (CI), in order to determine the long-term trends of significant changes for UP in China and G20 countries from 1990 to 2023. This model establishes segmented regression based on the temporal characteristics of disease distribution, dividing the time range into different intervals, each of which undergoes trend fitting and optimization, effectively avoiding the subjectivity of typical trend analyses based on linear trends. The trend direction was determined based on the AAPC calculated from the final model; specifically, when the 95% CI does not include 0, it means that the trend is significant, AAPC > 0 indicates an increasing trend, and APPC < 0 indicates a decreasing trend.
BAPC Analysis
The APC model was extended using Bayesian statistical methods (BAPC), which integrate historical data patterns, uncertainty factors, and prior knowledge to improve prediction accuracy and robustness, thereby effectively handling noise and ambiguities in the data. BAPC projections for 2024 to 2038 were created using a smoothing parameter of 5 for age, period, and cohort effects, with standardized weights based on the world standard population.
ARIMA Analysis
The ARIMA model is a popular model in econometrics that can analyze the behavior of stationary and non-stationary time series as well as the impact of plans and policies on specific outcomes over time34. In the ARIMA (p, d, q) model, 'p' represents the count of autoregressive terms, 'd' signifies the order of differencing, and 'q' denotes the count of moving average terms. The specific calculation methods are consistent with those reported in previous studies. For each prediction, the predicted Value, upper limit, lower limit and 95%CI of each year and the performance parameters of the prediction model will be obtained (_eval.csv file). The ARIMA time series model was used to predict the next 27 years (2024-2050), and the model performance was evaluated. Model performance was assessed using standard forecasting diagnostics, including Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), residual autocorrelation analysis, and goodness-of-fit evaluation. Residual plots demonstrated no major violations of model assumptions, supporting the suitability of the selected forecasting models for burden projection. Auto.arima was used to automatically select the optimal ARIMA(p,d,q) parameters and AICc (modified AIC) was used to select the best model. Based on the bootstrap prediction method, the prediction value and confidence interval were calculated by simulating the possible future paths.
Decomposition analyses
We performed a Das Gupta decomposition (Kitagawa–Das Gupta factorization) of the change in absolute counts between 1990 and 2023. The analysis was primarily decomposed into three components: the aging effect (changes in age structure), the population size effect (growth in total population), and the epidemiological change effect (changes in age-specific rates). The data column names are as follows: 'overall_difference' denotes the total change; 'a_effect' indicates the aging effect; 'p_effect' represents the population size effect; 'r_effect' signifies the epidemiological change effect; 'a_percent', 'p_percent', and 'r_percent' denote the percentage contribution of each respective effect; 'val_1990' and 'val_2023' represent the number of cases in 1990 and 2023, respectively; and 'change' refers to the difference in 'diff1'. The magnitude of each factor on the resulting plot illustrates its proportional influence, with positive and negative values indicating an increase or mitigation effect, respectively. Black dots on the plot serve as markers for the total change.
Analysis software
Data statistical analysis and visualization in this study were performed using the R software package and the Joinpoint software program (Refer to Table of Materials for details). A p-value < 0.05 was considered statistically significant.