1 Abstract
Background: Bladder cancer (BC) is a common malignancy with substantial morbidity, mortality and health-care cost, but comparable, state-level estimates of its burden and of sex differences in the U.S. are limited.
Aim: To describe Global Burden of Disease (GBD) 2021 estimates of BC incidence, prevalence, mortality, disability-adjusted life years (DALYs), years lived with disability (YLDs) and years of life lost (YLLs) in the U.S. and its states, their change over time, and male-to-female differences. This is a descriptive analysis of modelled estimates; it does not test causes or policy effects.
Methods: We extracted age-standardized rates (ASRs; per 100,000) with 95% uncertainty intervals (UIs) from GBD 2021 for the U.S., the 50 states and the District of Columbia, and for the European Union, high-income and high-SDI aggregates (which include the U.S.). Series cover 1990–2021 (deaths and YLLs also 1980–2021). Percent changes (1990–2021) and male:female rate ratios (2021) were derived from published point estimates and UIs under approximating assumptions.The final version and code of this study are also available from www.bladder-cancer.ir.
Results: In 2021, the U.S. ASR of BC incidence was 14.7 [13.6, 15.4] per 100,000 (men 24.0 [22.4, 25.2]; women 7.1 [6.3, 7.4]) and the ASR of mortality was 3.4 [3.0, 3.6] per 100,000 (men 5.8 [5.2, 6.1]; women 1.7 [1.4, 1.8]); DALYs were 68.1 [62.7, 72.5] per 100,000. Compared with 1990, the 2021 ASR changed by +3.6% [−3.7%, +11.4%] for incidence, −3.0% [−12.8%, +7.9%] for mortality and −8.6% [−16.4%, −0.1%] for DALYs (percent change [approximate 95% UI]). State-level male:female incidence rate ratios in 2021 ranged from 2.5 to 4.1 (point estimates); the approximate UIs of the lowest and highest jurisdictions overlapped, and the published UIs of the 2021 incidence ASR overlapped for 97% of all pairs of jurisdictions, so rankings should be interpreted cautiously.
Conclusion: GBD 2021 estimates show a male excess in BC burden and variation across U.S. states and over time. Because these are modelled, descriptive estimates, they do not identify the causes of the differences; hypotheses about exposures, health-care access and policy require study designs that measure them.
Keywords: bladder cancer; global burden of disease; United States of America; incidence; mortality; prevalence
2 Introduction
Bladder cancer (BC) is among the most commonly diagnosed cancers in the U.S. The American Cancer Society projects approximately 84,530 new cases and 17,870 deaths from BC in 2026 (1), and BC is the tenth leading cause of cancer death (1,2). It is considerably more common in men than in women (3). Globally, the absolute burden of BC is increasing despite declining age-standardized rates (4).
Beyond its epidemiological impact, BC imposes a substantial economic burden. Its management commonly requires lifelong surveillance and repeated therapeutic interventions, and BC has been described as the most costly cancer to treat on a per-patient lifetime basis (5–7). In addition to direct healthcare expenditures, BC affects patients’ quality of life and the emotional, social, and practical well-being of caregivers and families (8).
Substantial geographic heterogeneity in BC incidence and outcomes has been documented across the U.S. (9–12). For instance, Maine reported an age-standardized incidence rate of 26.7 per 100,000 during 2017–2021, nearly double that observed in Hawaii (13.4) (13). Registry-based reports describe the highest incidence in northeastern states (11) and disproportionately high mortality in Southern states (14); proposed explanations include differences in environmental exposures, socioeconomic conditions, age and racial composition, rurality, and access to care (15).
Despite these observations, most prior BC research has focused on national trends or single-institution cohorts, and contemporary, age-standardized, state-level comparisons remain limited. Although the GBD project produces estimates of disability-adjusted life years (DALYs), years lived with disability (YLDs), and years of life lost (YLLs) for BC, these measures have rarely been used in disease-specific descriptions of the U.S. National studies indicate that men have approximately three to four times higher BC incidence than women (3), whereas women often present with more advanced disease and have poorer survival (16); the geographic distribution of sex differences has not, to our knowledge, been described systematically at the state level.
We therefore used GBD 2021 estimates, which provide internally consistent, modelled estimates of key BC burden measures for all U.S. states, to describe (i) the estimated burden of BC in the U.S. and its states over time, benchmarked against comparator regions; (ii) its age and sex patterns; and (iii) male-to-female differences across states. The analysis is descriptive. It does not evaluate health policies, exposures, or their effects.
3 Methodology
3.1 Study design and data sources
This is a descriptive, secondary analysis of publicly available GBD 2021 estimates. GBD produces modelled estimates for diseases and injuries by age, sex, location and year, which allows standardized comparison across U.S. states and over time (17,18). Estimates were downloaded from the Global Health Data Exchange (GHDx) GBD Results Tool. The download provides, for each combination of measure, location, sex, age group and year, a point estimate and the lower and upper bounds of the 95% uncertainty interval (UI); posterior draws are not part of this public download (see Section 3.3.2).
We extracted incidence, prevalence, mortality (deaths), YLLs, YLDs and DALYs for males, females and both sexes, for all age groups, five-year age groups and the age-standardized rate, for the U.S., the 50 states and the District of Columbia (DC; not a state, hence “jurisdictions” when pooled), and for three comparator regions: the European Union, the high-income super-region and the high-SDI quintile. Rates are expressed per 100,000 population. Years covered differ by measure: mortality and YLLs are available for 1980–2021, whereas the non-fatal measures (incidence, prevalence, YLDs) and therefore DALYs start in 1990. All quantified changes use 1990–2021 (the first year for which all six measures exist and the conventional GBD baseline) as the primary contrast, and 2010–2021 as a secondary contrast; the 1980–2021 change for deaths and YLLs is provided in Supporting Info 1. Figures show all available years.
GBD 2021 was used because the analysis pipeline and all state-level extracts were built on this release. GBD 2023 results, including cancer estimates, have since been published (19); this is discussed as a limitation (Section 5.6). Given the value codes and data structure remains the same in GBD 2023 data, this study will be easily reproducible by replacing the data. Cancer estimation in GBD is anchored on mortality: incidence is derived from mortality using mortality-to-incidence ratios, and prevalence and YLDs are derived from incidence and survival (20). The six measures are therefore not independent sources of information.
3.2 Definitions
- Bladder cancer
-
As defined in GBD 2021 using the International Classification of Diseases (ICD): ICD-10 codes C67–C67.9 and ICD-9 codes 188–188.9, V10.51, V16.52 and V76.3. These codes are reproduced from the GBD 2021 cause list and cause-mapping appendices (17,18) and were not modified by the authors.
- Disability-adjusted life years (DALYs)
-
A summary measure of population health that combines YLLs and YLDs.
- Years of life lost (YLLs)
-
The fatal component of DALYs: years lost due to death occurring before the reference life expectancy.
- Years lived with disability (YLDs)
-
The non-fatal component of DALYs: time lived in health states less than full health, weighted by severity.
- Socio-Demographic Index (SDI)
-
A composite measure developed within the GBD framework reflecting income per capita, average educational attainment among people aged 15 years or older, and fertility among those younger than 25 years (21). It ranges from 0 to 1 and is used to classify locations into five quintiles (low to high SDI).
- Comparator regions
-
The European Union, the high-income super-region and the high-SDI quintile, used as reference points for the U.S. The U.S. is itself a constituent of the high-income and high-SDI aggregates and contributes a substantial share of their population; contrasts between the U.S. and these aggregates are therefore not independent country-to-country comparisons and are expected to be attenuated. The European Union does not contain the U.S.
- Age-standardized rate (ASR)
-
A rate adjusted to a common age structure so that populations with different age distributions can be compared. By construction, an ASR is not affected by changes in the age structure of the population. ASRs were obtained directly from the GBD Results Tool (they were not recalculated by the authors). GBD uses direct standardization to the GBD world standard population:
\[\mathrm{ASR} = \frac{\sum_{i=1}^{A} a_i\, w_i}{\sum_{i=1}^{A} w_i} \times 100{,}000\]
where \(a_i\) is the age-specific rate (events per person-year) in age group \(i\), \(w_i\) is the corresponding standard-population weight, and \(A\) is the number of age groups.
3.3 Statistical analysis
3.3.1 Summary measures
All statements about trends refer either to the plotted series or to the following fixed summary measures; no trend model (e.g., Joinpoint regression, annual percent change, or an age–period–cohort model) was fitted, and no test of differing slopes between the U.S. and comparators was performed.
- Percent change in the ASR between two years, \(100\times(a_{\text{end}}/a_{\text{start}}-1)\), for 1990–2021 and 2010–2021 (and 1980–2021 for deaths and YLLs), for the U.S., each jurisdiction and each comparator region, by sex and measure.
- Male:female rate ratio of ASRs in 2021 for each jurisdiction and measure.
- Rank of jurisdictions by point estimate of the 2021 ASR (both sexes).
The 95% UIs reported for GBD estimates are Bayesian uncertainty intervals (2.5th and 97.5th percentiles of 1,000 posterior draws) that reflect data sparsity, measurement error and model specification (18). They are not frequentist confidence intervals. We therefore do not use the word “significant”, and describe whether an approximate UI includes zero only as a description of that interval.
3.3.2 Approximate uncertainty intervals for derived quantities
Percent changes and rate ratios are derived quantities for which the public GBD download does not provide UIs. The posterior-draw approach, in which the quantity is computed in each draw and then summarised, is the preferred way of obtaining UIs for derived quantities in GBD (22). Because draws were not part of our download, we approximated these UIs from the published bounds using standard methods for ratios of independent estimates (23,24).
For two ASRs \(a_1\) and \(a_2\) with published 95% UIs \([l_1,u_1]\) and \([l_2,u_2]\), the standard error of each log-rate was back-calculated from the bounds, assuming that the UI is symmetric on the log scale (this is a direct back-calculation, not a delta-method approximation):
\[SE_{\log(a_k)} \approx \frac{\log(u_k) - \log(l_k)}{2\, z_{0.975}}, \quad k = 1, 2, \qquad z_{0.975}\approx 1.96 .\]
The ratio \(R=a_1/a_2\) (for percent change, \(a_1\) is the later and \(a_2\) the earlier year; percent change is \(100(R-1)\); for sex differences, \(a_1\) is the male and \(a_2\) the female rate) has \(\log R=\log a_1-\log a_2\) and
\[SE_{\log R} = \sqrt{SE_{\log a_1}^2 + SE_{\log a_2}^2}, \qquad \text{UI}_{95\%} = \exp\!\big(\log R \pm z_{0.975}\, SE_{\log R}\big).\]
The point estimate is the ratio of the published point estimates, not the median of draw-wise ratios.
Limitations of the approximation. (i) The estimates being compared (two years, or the two sexes, of the same location) come from shared models, covariates and input data and are therefore not independent. The sign of the dependence cannot be determined from published summaries: positive correlation would make the approximate UIs too wide, negative correlation too narrow. (ii) Published UIs are often asymmetric on the log scale, and the GBD point estimate is a mean of draws rather than a median. The resulting intervals are therefore labelled approximate, are intended to convey the order of magnitude of uncertainty, and are not used for formal tests. In particular, we do not interpret non-overlap (or overlap) of two jurisdictions’ intervals as a test of a difference between them.
3.4 Reporting and reproducibility
This study follows the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) (25); the completed checklist, with manuscript locations, is provided as Supporting Info 10. All numbers quoted in the Abstract and Results are generated by the analysis code from the downloaded data, rather than typed. Analyses used R version 4.6.1 (2026-06-24); package versions are listed in Supporting Info 9. Complete code, cleaned data (in RDS format) and figure files are in the study’s GitHub repository; the submitted figures and tables were generated from commit 16f52d8, archived at 10.5281/zenodo.23814418 (26). The study’s website (bladder-cancer.ir) hosts the manuscript, its source code, and authors’ correspondence with the reviewers.
4 Results
4.1 U.S. compared with comparator regions
In 2021, the U.S. age-standardized incidence rate of BC was 24.0 [22.4, 25.2] per 100,000 in men and 7.1 [6.3, 7.4] in women, and the age-standardized mortality rate was 5.8 [5.2, 6.1] in men and 1.7 [1.4, 1.8] in women (point estimate [95% UI]). For both sexes combined, the 2021 mortality rate was 3.4 [3.0, 3.6] in the U.S., compared with 4.6 [4.1, 4.9] in the European Union, 3.5 [3.1, 3.7] in the high-income aggregate and 3.4 [3.0, 3.6] in the high-SDI aggregate (all per 100,000).
The U.S. ASR of incidence (both sexes) peaked in 2011 (16.5 [15.4, 17.0] per 100,000) and the ASR of prevalence in 2011 (113.9 [108.3, 117.2]). The U.S. point estimate of the incidence ASR exceeded that of each comparator in the following numbers of years: European Union, 32 of 32 years; High-income, 32 of 32 years; High SDI, 32 of 32 years. For mortality, the U.S. point estimate exceeded that of each comparator in: European Union, 0 of 42 years; High-income, 0 of 42 years; High SDI, 7 of 42 years. Figure 1 shows the series; Table 1 gives, for each region and measure, the 1990 and 2021 rates and the percent changes, and is the source for the figures quoted below.
Abbreviations: DALYs: disability-adjusted life years; YLDs: years lived with disability; YLLs: years of life lost; EU: European Union; SDI: Socio-Demographic Index.
Measure | Region | 1990 | 2021 | % change 1990–2021 | % change 2010–2021 |
|---|---|---|---|---|---|
Incidence | United States | 14.2 [13.5, 14.7] | 14.7 [13.6, 15.4] | +3.6% [−3.7%, +11.4%] | −10.5% [−17.4%, −3.1%] |
European Union | 13.4 [12.8, 13.8] | 13.1 [12.1, 13.9] | −2.3% [−9.6%, +5.7%] | −8.7% [−15.9%, −0.9%] | |
High-income | 12.6 [12.0, 13.0] | 11.5 [10.6, 12.0] | −9.1% [−15.7%, −1.9%] | −10.1% [−17.1%, −2.4%] | |
High SDI | 11.2 [10.7, 11.6] | 10.6 [9.8, 11.1] | −5.3% [−12.2%, +2.1%] | −7.7% [−14.9%, 0.0%] | |
Prevalence | United States | 94.0 [90.6, 96.8] | 100.6 [95.0, 104.6] | +7.0% [+1.0%, +13.5%] | −11.4% [−16.8%, −5.7%] |
European Union | 72.2 [69.5, 74.5] | 80.6 [75.8, 85.2] | +11.6% [+4.3%, +19.5%] | −7.6% [−13.7%, −1.1%] | |
High-income | 73.0 [70.5, 75.2] | 73.8 [69.1, 77.0] | +1.1% [−5.0%, +7.6%] | −9.4% [−15.1%, −3.3%] | |
High SDI | 63.6 [61.2, 65.4] | 66.8 [62.8, 69.6] | +5.1% [−1.1%, +11.6%] | −7.1% [−12.7%, −1.2%] | |
Mortality | United States | 3.5 [3.3, 3.7] | 3.4 [3.0, 3.6] | −3.0% [−12.8%, +7.9%] | −6.3% [−16.6%, +5.2%] |
European Union | 5.8 [5.5, 6.0] | 4.6 [4.1, 4.9] | −21.2% [−28.4%, −13.2%] | −8.9% [−17.9%, +1.1%] | |
High-income | 4.7 [4.4, 4.9] | 3.5 [3.1, 3.7] | −25.7% [−33.0%, −17.8%] | −10.2% [−19.8%, +0.6%] | |
High SDI | 4.4 [4.1, 4.5] | 3.4 [3.0, 3.6] | −21.9% [−29.5%, −13.6%] | −8.0% [−17.5%, +2.7%] | |
YLDs | United States | 8.0 [5.8, 10.4] | 8.4 [6.1, 10.8] | +4.1% [−30.6%, +56.4%] | −11.0% [−40.4%, +33.0%] |
European Union | 6.7 [4.9, 8.7] | 7.1 [5.2, 9.4] | +6.2% [−29.7%, +60.3%] | −7.7% [−39.0%, +39.7%] | |
High-income | 6.6 [4.8, 8.5] | 6.4 [4.6, 8.2] | −3.4% [−35.5%, +44.7%] | −9.5% [−39.6%, +35.7%] | |
High SDI | 5.8 [4.2, 7.5] | 5.8 [4.3, 7.5] | +0.1% [−33.1%, +49.7%] | −7.2% [−37.8%, +38.4%] | |
YLLs | United States | 66.6 [63.2, 68.8] | 59.8 [55.1, 62.6] | −10.1% [−16.8%, −2.9%] | −9.9% [−17.1%, −2.1%] |
European Union | 112.7 [108.1, 115.8] | 82.9 [76.9, 87.8] | −26.4% [−31.7%, −20.7%] | −11.8% [−18.4%, −4.6%] | |
High-income | 88.5 [84.4, 91.1] | 60.4 [55.6, 63.6] | −31.7% [−36.7%, −26.2%] | −13.5% [−20.4%, −5.9%] | |
High SDI | 82.1 [78.3, 84.6] | 59.3 [54.7, 62.4] | −27.8% [−33.1%, −22.1%] | −11.4% [−18.3%, −4.0%] | |
DALYs | United States | 74.6 [70.4, 78.1] | 68.1 [62.7, 72.5] | −8.6% [−16.4%, −0.1%] | −10.1% [−18.1%, −1.2%] |
European Union | 119.4 [114.7, 123.4] | 90.1 [83.3, 96.0] | −24.6% [−30.4%, −18.3%] | −11.5% [−18.7%, −3.6%] | |
High-income | 95.1 [90.8, 98.6] | 66.8 [61.6, 70.9] | −29.7% [−35.2%, −23.8%] | −13.1% [−20.4%, −5.1%] | |
High SDI | 87.9 [84.0, 91.2] | 65.1 [60.1, 68.7] | −26.0% [−31.6%, −20.0%] | −11.1% [−18.3%, −3.2%] |
Abbreviations: UI, uncertainty interval; SDI, Socio-Demographic Index; YLDs, years lived with disability; YLLs, years of life lost; DALYs, disability-adjusted life years.
For both sexes combined, the percent change in the ASR of mortality between 1990 and 2021 was, by region (percent change [approximate 95% UI]): United States −3.0% [−12.8%, +7.9%]; European Union −21.2% [−28.4%, −13.2%]; High-income −25.7% [−33.0%, −17.8%]; High SDI −21.9% [−29.5%, −13.6%]. For DALYs it was: United States −8.6% [−16.4%, −0.1%]; European Union −24.6% [−30.4%, −18.3%]; High-income −29.7% [−35.2%, −23.8%]; High SDI −26.0% [−31.6%, −20.0%]. For incidence it was: United States +3.6% [−3.7%, +11.4%]; European Union −2.3% [−9.6%, +5.7%]; High-income −9.1% [−15.7%, −1.9%]; High SDI −5.3% [−12.2%, +2.1%]; and between 2010 and 2021: United States −10.5% [−17.4%, −3.1%]; European Union −8.7% [−15.9%, −0.9%]; High-income −10.1% [−17.1%, −2.4%]; High SDI −7.7% [−14.9%, 0.0%]. Sex-specific values, and values for all other measures and for every jurisdiction, are in Supporting Info 1 and, as a graphical table, Supporting Info 3. Because the percent-change intervals are approximate and no trend model was fitted, we describe these as differences between years, not as tests of trend. Because the high-income and high-SDI aggregates contain the U.S., differences between the U.S. and these aggregates are attenuated by construction.
4.2 Age and sex patterns in the U.S.
As illustrated in Figure 2, in 2021 the male rate exceeded the female rate in 84 of 96 age-group–by–measure combinations. The highest 2021 age-specific rates (per 100,000) were: DALYs, women (80+ years): 462.3 [354.7, 521.8]; men (80+ years): 1,617.1 [1,365.6, 1,752.9]. mortality, women (80+ years): 42.1 [31.9, 47.5]; men (80+ years): 142.2 [118.3, 154.5]. incidence, women (80+ years): 84.1 [63.8, 94.5]; men (80+ years): 309.3 [258.4, 335.0]. prevalence, women (75-79 years): 411.0 [358.2, 441.0]; men (75-79 years): 1,481.3 [1,366.4, 1,568.7]. YLDs, women (80+ years): 37.8 [26.7, 48.8]; men (80+ years): 136.1 [99.5, 174.7]. YLLs, women (80+ years): 424.6 [323.9, 477.4]; men (80+ years): 1,481.0 [1,241.4, 1,606.1].
Abbreviations: DALYs: disability-adjusted life years; YLDs: years lived with disability; YLLs: years of life lost.
State-level population pyramids for 2021 are provided in Supporting Info 4, and animated pyramids showing change by year for each state and the U.S. in Supporting Info 5. Figure 3 displays age-specific rates in 1990, 2000, 2010 and 2021 by sex. It is a descriptive age-by-period display; no age–period–cohort model was fitted. Corresponding plots for each state are in Supporting Info 6.
Abbreviations: DALYs: disability-adjusted life years; YLDs: years lived with disability; YLLs: years of life lost.
Absolute case and death counts were not part of the extracted dataset; we therefore cannot report changes in the number of cases.
4.3 Bladder cancer burden across U.S. states
State-specific, sex-stratified temporal trends in age-standardized incidence rates are displayed in Figure 4; geofaceted plots for the other five measures are provided in Supporting Info 7.
Abbreviations: Alabama: AL; Alaska: AK; Arizona: AZ; Arkansas: AR; California: CA; Colorado: CO; Connecticut: CT; Delaware: DE; District of Columbia: DC; Florida: FL; Georgia: GA; Hawaii: HI; Idaho: ID; Illinois: IL; Indiana: IN; Iowa: IA; Kansas: KS; Kentucky: KY; Louisiana: LA; Maine: ME; Maryland: MD; Massachusetts: MA; Michigan: MI; Minnesota: MN; Mississippi: MS; Missouri: MO; Montana: MT; Nebraska: NE; Nevada: NV; New Hampshire: NH; New Jersey: NJ; New Mexico: NM; New York: NY; North Carolina: NC; North Dakota: ND; Ohio: OH; Oklahoma: OK; Oregon: OR; Pennsylvania: PA; Rhode Island: RI; South Carolina: SC; South Dakota: SD; Tennessee: TN; Texas: TX; Utah: UT; Vermont: VT; Virginia: VA; Washington: WA; West Virginia: WV; Wisconsin: WI; Wyoming: WY.
Table 2 lists the 2021 age-standardized incidence and mortality rates (both sexes), the rank of each jurisdiction by point estimate (1 = highest), and the percent change since 1990. By point estimate, the highest 2021 incidence ASRs (both sexes, per 100,000) were in Maine (19.9), Nevada (18.9), Delaware (18.6), Kentucky (18.0), and West Virginia (17.5) and the lowest in Hawaii (11.0), Utah (12.0), North Dakota (12.3), California (12.8), and South Dakota (12.9); the highest 2021 mortality ASRs were in Maine (4.5), Nevada (4.3), Kentucky (4.2), Delaware (4.2), and West Virginia (4.1) and the lowest in Hawaii (2.3), North Dakota (2.8), California (2.9), New York (2.9), and Utah (3.0). These are rankings of point estimates; the published 95% UIs of the 2021 ASR overlapped for 97% of all pairs of jurisdictions for incidence and 95% for mortality (Table 2), so the ranking itself is uncertain. Between 1990 and 2021, the largest point-estimate decreases in incidence ASR were in Massachusetts (−21.2%), New York (−14.1%), New Jersey (−11.9%), Rhode Island (−11.3%), and Wisconsin (−9.9%), and the smallest decreases or largest increases in West Virginia (+25.4%), Louisiana (+24.6%), Oklahoma (+24.5%), Alabama (+23.5%), and New Mexico (+22.4%). The approximate UI of the 1990–2021 percent change in incidence lay entirely below zero in 1 of 51 jurisdictions, entirely above zero in 6, and included zero in 44; for mortality the corresponding numbers were 2, 1 and 48. Rates by sex for all six measures are in Supporting Info 1.
Incidence | Mortality | |||||
|---|---|---|---|---|---|---|
Jurisdiction | 2021 rate | Rank | % change 1990–2021 | 2021 rate | Rank | % change 1990–2021 |
United States of America | 14.7 [13.6, 15.4] | – | +3.6% [−3.7%, +11.4%] | 3.4 [3.0, 3.6] | – | −3.0% [−12.8%, +7.9%] |
Alabama | 15.0 [12.6, 17.8] | 25 | +23.5% [+2.2%, +49.2%] | 3.8 [3.2, 4.4] | 9 | +15.5% [−3.8%, +38.8%] |
Alaska | 14.6 [12.5, 16.9] | 34 | +2.8% [−13.3%, +21.9%] | 3.3 [2.8, 3.8] | 42 | −13.6% [−28.1%, +3.9%] |
Arizona | 15.3 [12.8, 18.3] | 19 | +6.6% [−12.3%, +29.5%] | 3.4 [2.8, 4.0] | 34 | −2.7% [−19.8%, +18.0%] |
Arkansas | 15.0 [12.5, 17.7] | 26 | +15.7% [−4.1%, +39.6%] | 3.7 [3.1, 4.3] | 13 | +7.7% [−10.6%, +29.8%] |
California | 12.8 [10.7, 14.9] | 48 | +3.6% [−13.5%, +24.1%] | 2.9 [2.4, 3.4] | 49 | −8.8% [−23.7%, +8.9%] |
Colorado | 13.0 [10.5, 15.6] | 44 | +1.4% [−18.1%, +25.6%] | 3.0 [2.4, 3.6] | 45 | −1.0% [−20.8%, +23.7%] |
Connecticut | 14.9 [12.1, 17.8] | 28 | −8.3% [−25.5%, +12.9%] | 3.5 [2.8, 4.1] | 29 | −9.5% [−26.3%, +11.2%] |
Delaware | 18.6 [15.8, 21.5] | 3 | +9.5% [−7.5%, +29.6%] | 4.2 [3.6, 4.9] | 4 | −2.8% [−18.1%, +15.3%] |
District of Columbia | 13.7 [11.3, 16.4] | 40 | +11.5% [−8.4%, +35.7%] | 3.3 [2.7, 3.9] | 41 | −25.1% [−38.1%, −9.4%] |
Florida | 16.1 [13.3, 19.0] | 9 | +9.1% [−9.9%, +32.1%] | 3.5 [2.9, 4.2] | 20 | +0.2% [−17.7%, +21.9%] |
Georgia | 14.7 [12.2, 17.6] | 31 | +21.0% [−0.3%, +46.9%] | 3.6 [3.0, 4.2] | 18 | +7.2% [−10.8%, +29.0%] |
Hawaii | 11.0 [9.1, 13.2] | 51 | +12.2% [−8.8%, +38.1%] | 2.3 [1.9, 2.8] | 51 | −2.9% [−22.0%, +20.8%] |
Idaho | 15.0 [12.3, 17.8] | 24 | +16.1% [−5.0%, +41.9%] | 3.7 [3.0, 4.3] | 15 | +13.7% [−6.8%, +38.8%] |
Illinois | 15.3 [12.8, 18.1] | 21 | +7.3% [−10.9%, +29.3%] | 3.5 [2.9, 4.2] | 19 | −1.9% [−18.8%, +18.5%] |
Indiana | 16.3 [13.5, 19.3] | 8 | +10.9% [−8.6%, +34.5%] | 3.9 [3.2, 4.6] | 7 | +9.8% [−9.3%, +32.9%] |
Iowa | 14.5 [12.1, 17.4] | 36 | +11.1% [−8.9%, +35.5%] | 3.3 [2.8, 4.0] | 37 | +10.7% [−8.1%, +33.4%] |
Kansas | 15.0 [12.4, 18.1] | 23 | +19.6% [−2.6%, +46.9%] | 3.4 [2.9, 4.0] | 33 | +20.0% [−0.9%, +45.2%] |
Kentucky | 18.0 [15.1, 21.4] | 4 | +16.3% [−3.5%, +40.1%] | 4.2 [3.5, 5.0] | 3 | +15.5% [−4.0%, +39.0%] |
Louisiana | 15.8 [13.3, 18.6] | 12 | +24.6% [+4.0%, +49.2%] | 3.8 [3.2, 4.4] | 10 | +12.7% [−5.3%, +34.1%] |
Maine | 19.9 [16.5, 24.0] | 1 | +5.8% [−13.7%, +29.9%] | 4.5 [3.8, 5.4] | 1 | +3.6% [−14.7%, +25.9%] |
Maryland | 15.3 [12.6, 18.1] | 20 | −4.9% [−21.8%, +15.7%] | 3.5 [2.9, 4.1] | 23 | −15.3% [−30.3%, +2.8%] |
Massachusetts | 14.6 [12.0, 17.5] | 33 | −21.2% [−35.7%, −3.4%] | 3.4 [2.8, 4.1] | 31 | −15.2% [−30.2%, +3.0%] |
Michigan | 14.5 [12.3, 17.1] | 35 | +1.6% [−15.0%, +21.6%] | 3.5 [2.9, 4.0] | 27 | −6.1% [−21.3%, +12.0%] |
Minnesota | 13.7 [11.2, 16.4] | 41 | −4.4% [−22.0%, +17.1%] | 3.3 [2.6, 3.8] | 40 | +4.2% [−15.0%, +27.7%] |
Mississippi | 13.0 [11.0, 15.4] | 45 | +16.9% [−2.7%, +40.3%] | 3.4 [2.9, 3.9] | 36 | +8.6% [−8.8%, +29.3%] |
Missouri | 15.7 [13.1, 18.6] | 14 | +14.9% [−4.9%, +38.9%] | 3.6 [3.1, 4.2] | 17 | +8.2% [−8.9%, +28.6%] |
Montana | 14.9 [12.3, 17.6] | 27 | +6.7% [−12.2%, +29.6%] | 3.5 [2.9, 4.2] | 24 | +0.5% [−17.6%, +22.4%] |
Nebraska | 14.2 [11.8, 16.8] | 39 | +3.4% [−15.0%, +25.9%] | 3.3 [2.7, 3.9] | 39 | +6.0% [−13.4%, +29.9%] |
Nevada | 18.9 [15.4, 22.5] | 2 | +14.0% [−7.0%, +39.8%] | 4.3 [3.5, 5.0] | 2 | +1.2% [−17.1%, +23.4%] |
New Hampshire | 16.5 [13.4, 19.8] | 7 | −9.1% [−26.3%, +12.1%] | 3.9 [3.2, 4.6] | 8 | −6.8% [−24.4%, +15.1%] |
New Jersey | 15.4 [12.7, 18.4] | 15 | −11.9% [−27.8%, +7.5%] | 3.5 [2.9, 4.2] | 26 | −16.6% [−31.5%, +1.6%] |
New Mexico | 13.2 [11.0, 15.6] | 43 | +22.4% [+0.8%, +48.6%] | 3.0 [2.5, 3.6] | 44 | +2.7% [−15.8%, +25.4%] |
New York | 13.0 [10.6, 15.4] | 46 | −14.1% [−29.5%, +4.7%] | 2.9 [2.4, 3.5] | 48 | −24.5% [−38.4%, −7.5%] |
North Carolina | 14.8 [12.4, 17.6] | 30 | +19.5% [−0.9%, +44.2%] | 3.5 [2.9, 4.1] | 28 | +8.7% [−10.2%, +31.7%] |
North Dakota | 12.3 [10.3, 14.3] | 49 | −3.4% [−20.0%, +16.6%] | 2.8 [2.3, 3.3] | 50 | −5.0% [−21.3%, +14.8%] |
Ohio | 16.8 [13.9, 19.8] | 6 | +7.9% [−10.7%, +30.4%] | 4.0 [3.4, 4.7] | 6 | +6.0% [−11.4%, +26.8%] |
Oklahoma | 15.4 [12.9, 18.2] | 16 | +24.5% [+3.2%, +50.2%] | 3.7 [3.1, 4.3] | 11 | +23.7% [+3.1%, +48.4%] |
Oregon | 16.0 [13.2, 19.3] | 11 | +7.7% [−12.3%, +32.2%] | 3.5 [3.0, 4.2] | 22 | +1.5% [−16.9%, +24.1%] |
Pennsylvania | 15.7 [13.0, 18.7] | 13 | −2.6% [−20.0%, +18.5%] | 3.7 [3.0, 4.4] | 14 | −6.0% [−23.0%, +14.9%] |
Rhode Island | 14.9 [12.2, 17.6] | 29 | −11.3% [−27.6%, +8.7%] | 3.5 [2.9, 4.2] | 25 | −11.5% [−27.2%, +7.6%] |
South Carolina | 14.7 [12.1, 17.1] | 32 | +20.8% [+0.5%, +45.2%] | 3.5 [3.0, 4.1] | 21 | +4.2% [−12.5%, +24.0%] |
South Dakota | 12.9 [10.9, 15.2] | 47 | +2.8% [−14.7%, +24.0%] | 3.0 [2.5, 3.6] | 46 | −1.2% [−18.9%, +20.4%] |
Tennessee | 15.3 [12.8, 18.1] | 18 | +19.6% [−0.7%, +44.1%] | 3.7 [3.1, 4.2] | 16 | +11.0% [−7.2%, +32.6%] |
Texas | 13.3 [11.3, 15.5] | 42 | +9.7% [−7.4%, +29.9%] | 3.1 [2.6, 3.6] | 43 | +0.8% [−15.4%, +20.1%] |
Utah | 12.0 [10.0, 14.3] | 50 | +11.2% [−8.5%, +35.3%] | 3.0 [2.5, 3.5] | 47 | +6.2% [−12.7%, +29.3%] |
Vermont | 16.1 [13.6, 18.6] | 10 | −3.1% [−18.9%, +15.7%] | 3.7 [3.0, 4.3] | 12 | −9.8% [−25.3%, +9.0%] |
Virginia | 15.1 [12.5, 18.0] | 22 | +3.8% [−14.5%, +26.1%] | 3.4 [2.8, 4.1] | 32 | −3.9% [−21.1%, +17.1%] |
Washington | 14.3 [11.8, 16.9] | 37 | −0.4% [−17.9%, +20.9%] | 3.3 [2.7, 3.9] | 38 | −1.7% [−19.1%, +19.3%] |
West Virginia | 17.5 [14.7, 20.3] | 5 | +25.4% [+5.0%, +49.7%] | 4.1 [3.5, 4.8] | 5 | +15.3% [−3.1%, +37.0%] |
Wisconsin | 14.2 [11.8, 17.2] | 38 | −9.9% [−26.3%, +10.2%] | 3.4 [2.8, 4.1] | 35 | −6.5% [−23.9%, +14.8%] |
Wyoming | 15.4 [13.0, 17.7] | 17 | +18.7% [−0.1%, +41.1%] | 3.5 [2.9, 3.9] | 30 | +5.9% [−10.9%, +25.9%] |
Figure 5 displays male-to-female rate ratios of ASRs in 2021 for each jurisdiction and measure. Because incidence is derived from mortality and prevalence and YLDs are derived from incidence and survival in the GBD cancer models (20), the six ratios are strongly related by construction and are not independent evidence; we therefore describe incidence and mortality. Across the 50 states and DC, male:female incidence rate ratios ranged from District of Columbia, 2.5 [1.8, 3.4] (lowest point estimate) to North Dakota, 4.1 [3.1, 5.4] (highest point estimate), and mortality ratios from District of Columbia, 2.4 [1.8, 3.2] to South Dakota, 4.2 [3.2, 5.6] (ratio [approximate 95% UI]). The approximate UIs of the lowest and highest jurisdictions overlapped for incidence and overlapped for mortality. These UIs assume independence between the male and female estimates (Section 3.3.2) and must not be read as a test of a difference between jurisdictions. For the U.S. overall, the 2021 male:female incidence ratio was 3.4 [3.1, 3.8]. All ratios are in Supporting Info 2.
Abbreviations: Alabama: AL; Alaska: AK; Arizona: AZ; Arkansas: AR; California: CA; Colorado: CO; Connecticut: CT; Delaware: DE; District of Columbia: DC; Florida: FL; Georgia: GA; Hawaii: HI; Idaho: ID; Illinois: IL; Indiana: IN; Iowa: IA; Kansas: KS; Kentucky: KY; Louisiana: LA; Maine: ME; Maryland: MD; Massachusetts: MA; Michigan: MI; Minnesota: MN; Mississippi: MS; Missouri: MO; Montana: MT; Nebraska: NE; Nevada: NV; New Hampshire: NH; New Jersey: NJ; New Mexico: NM; New York: NY; North Carolina: NC; North Dakota: ND; Ohio: OH; Oklahoma: OK; Oregon: OR; Pennsylvania: PA; Rhode Island: RI; South Carolina: SC; South Dakota: SD; Tennessee: TN; Texas: TX; Utah: UT; Vermont: VT; Virginia: VA; Washington: WA; West Virginia: WV; Wisconsin: WI; Wyoming: WY; DALYs: disability-adjusted life years; YLDs: years lived with disability; YLLs: years of life lost.
Bump plots comparing the U.S., all jurisdictions and the comparator regions across measures in 1990, 2000, 2010 and 2021 are in Supporting Info 8.
5 Discussion
5.1 Principal findings
In this descriptive analysis of GBD 2021 estimates, the U.S. age-standardized incidence rate of BC in 2021 was 14.7 [13.6, 15.4] per 100,000 and the age-standardized mortality rate 3.4 [3.0, 3.6] per 100,000 (both sexes). The incidence ASR peaked in 2011; between 2010 and 2021 it changed by −10.5% [−17.4%, −3.1%], and between 1990 and 2021 by +3.6% [−3.7%, +11.4%]. The mortality ASR changed by −3.0% [−12.8%, +7.9%] between 1990 and 2021 and by −6.3% [−16.6%, +5.2%] between 2010 and 2021 (percent change [approximate 95% UI]). We did not fit a trend model, so these are differences between years and not statements about the statistical significance of a trend. Across jurisdictions, estimated changes since 1990 varied (Table 2); for 44 of 51 jurisdictions the approximate UI of the incidence change included zero. All of these results are modelled estimates, and no exposure, policy or health-system variable was analysed; explanations offered below are hypotheses drawn from the literature, not findings of this study.
5.2 Comparison with other regions
The numbers of years in which the U.S. point estimate of incidence exceeded that of each comparator are given in Section 4.1. The high-income and high-SDI aggregates contain the U.S., so those contrasts are attenuated by construction and are not independent country-to-country comparisons; the European Union is the only comparator that excludes the U.S. A higher estimated incidence in the U.S. may partly reflect more intensive diagnostic activity (cystoscopy, imaging, and evaluation of hematuria) and case ascertainment (27–29), which raises recorded incidence without necessarily changing mortality, and/or differences in underlying risk, such as smoking history. This study cannot separate these possibilities, which are not mutually exclusive. Registry-based age-adjusted rates (for example, the SEER rate of about 18 per 100,000 for 2018–2022 (2)) are not directly comparable with GBD ASRs, which use a different standard population and are modelled.
5.3 Age
Age-standardized and age-specific rates are by design unaffected by population ageing; the rates reported here therefore cannot be attributed to population ageing. Ageing influences the number of cases: the U.S. population aged 65 years and older increased by 38.6% between 2010 and 2020 (30). Absolute case and death counts were not part of the extracted dataset; we therefore cannot report changes in the number of cases. The age groups with the highest age-specific rates are reported in Section 4.2; a rise in BC burden with age is in line with prior literature (4). An earlier analysis reported that the oldest adults, especially men aged 85 years and older, had not shared in the improvements in cancer mortality seen in younger seniors (31). We did not examine whether differences between states in age-specific rates reflect true differences, and we do not interpret them.
5.4 Sex differences
Male rates exceeded female rates in 84 of 96 age-group–by–measure combinations in 2021 (Section 4.2). The 2021 U.S. male:female incidence ratio was 3.4 [3.1, 3.8], and prior studies consistently report three- to four-fold higher incidence in men (3). For mortality, the 1990–2021 percent change was −6.7% [−15.3%, +2.8%] in men and −10.5% [−23.0%, +4.0%] in women. State-level ratios varied, with extremes of District of Columbia, 2.5 [1.8, 3.4] and North Dakota, 4.1 [3.1, 5.4] for incidence, but the approximate UIs are wide, assume independence of the male and female estimates, and the ratios are partly related across measures by construction (Section 3.1). GBD estimates contain no information on stage at diagnosis, treatment or survival. Studies outside GBD report that women more often present with advanced disease and have poorer survival (16,32,33), and that hematuria in women is more often attributed to urinary tract infection or gynecologic conditions (32,34). Differences in immune surveillance, hormonal milieu and genetic or epigenetic alterations have also been proposed (35). These are hypotheses that this analysis can neither support nor refute.
5.5 Geographic variation
By point estimate, the highest 2021 incidence ASRs were in Maine (19.9), Nevada (18.9), Delaware (18.6), Kentucky (18.0), and West Virginia (17.5) and the lowest in Hawaii (11.0), Utah (12.0), North Dakota (12.3), California (12.8), and South Dakota (12.9); the highest mortality ASRs were in Maine (4.5), Nevada (4.3), Kentucky (4.2), Delaware (4.2), and West Virginia (4.1) and the lowest in Hawaii (2.3), North Dakota (2.8), California (2.9), New York (2.9), and Utah (3.0) (Table 2). Because the published UIs overlap for 97% (incidence) and 95% (mortality) of all pairs of jurisdictions, these rankings are uncertain, and our data do not establish that any particular jurisdiction differs from another. Earlier registry-based reports describe high incidence in northeastern states (11) and elevated mortality in parts of the South (14,36). Explanations that have been proposed, none of which we tested, include arsenic in private-well drinking water in New England (37,38), smoking-attributable burden (39), rurality (40,41), and reduced access to timely diagnosis and treatment (41,42); studies of other cancers, such as colorectal and cervical cancer, have linked longer travel distance to worse care (43). GBD state estimates are not stratified by race or ethnicity, so we cannot assess whether differences in population composition (for example, lower BC risk reported in Asian and Pacific Islander populations (44)) contribute to differences between states. DC is a federal district with a small, predominantly urban population, and is not directly comparable with the states.
5.6 Strengths and limitations
Strengths include the use of a single, internally consistent modelling framework across all states, six burden measures, and fully scripted, openly available, reproducible analyses. Limitations are as follows. First, GBD estimates are modelled; for many states the underlying data are sparse, and incidence, prevalence and YLDs depend on mortality-based modelling (20), so the six measures are not independent. Second, derived UIs (percent change and rate ratios) are approximations from published bounds that assume independence and log-symmetry, and the direction of any resulting bias is unknown (Section 3.3.2); draw-based UIs would be preferable. Third, the U.S. is part of two of the three comparator aggregates, and we did not include comparators that exclude the U.S. (for example, individual peer countries or Western Europe). Fourth, no trend model was fitted, and no exposure, insurance, socioeconomic or policy variable was analysed, so the study cannot address causes or policy effects. Fifth, the analysis is limited to rates; population ageing affects case counts, which are reported only where available. Sixth, estimates are not stratified by race or ethnicity. Seventh, GBD estimates are revised between rounds, and GBD 2023 estimates, including for cancer (19), are now available; values here may differ from GBD 2023, and repeating the scripted pipeline on the newer release is a natural next step. Finally, discrepancies with surveillance systems that report crude or differently standardized rates likely reflect methodological differences, including the choice of standard population.
6 Conclusions
GBD 2021 modelled estimates describe a male excess in BC burden in the U.S. in 2021, changes in age-standardized rates between 1990 and 2021 that are summarised in Table 1, and variation across states for which uncertainty is substantial. These findings are descriptive. Whether differences between states, sexes or regions are related to exposures, health-care access or policy cannot be determined from these data, and such questions require analyses that measure those factors directly. Since the submission of this manuscript, GBD 2023 data has been released. Should the data structure and dictionary remain the same as GBD 2021, all visualizations of this study are reproducible for GBD 2023 as well. We encourage future researchers to further the current study by forking our GitHub repository and applying its source code to GBD 2023.
Supporting Info
Supporting Info are provided as appendices to the manuscript and are also publicly accessible at the links below (repository: GitHub; archived version: 10.5281/zenodo.23814418).
Percent change (numeric table)
Percent change in age-standardized rates, with approximate 95% UIs, for 1990–2021 and 2010–2021 (all measures) and 1980–2021 (deaths and YLLs), for the U.S., all jurisdictions and comparator regions, by sex, together with the start- and end-year rates. CSV file.Male:female rate ratios (2021, numeric table)
Male and female age-standardized rates and their ratio, with approximate 95% UIs (independence assumed), for every jurisdiction and measure. CSV file.Percent change 1990–2021 (graphical table)
Percent change in age-standardized rates between 1990 and 2021, by sex, for the U.S., jurisdictions and comparator regions. Bars show the point estimate and lines the approximate 95% UI. Green, yellow and red indicate an interval entirely below zero, including zero, and entirely above zero, respectively (a description of the interval, not a significance test). Available from the project’s GitHub repository.State-level population pyramids (2021)
Age- and sex-stratified bladder cancer burden for each jurisdiction in 2021, available from the project’s GitHub repository.Animated population pyramids
Temporal change in age- and sex-stratified burden for each jurisdiction, available from the project’s GitHub repository.Age-by-period plots
Line plots of age- and sex-specific rates in 1990, 2000, 2010 and 2021 for each jurisdiction (descriptive; no age–period–cohort model). The repository folder retains its original name, “APC plot 1”, available from the project’s GitHub repository.Geographically faceted state-level plots (all six measures)
Available from the project’s GitHub repository.
Abbreviations: DALYs: disability-adjusted life years; YLDs: years lived with disability; YLLs: years of life lost; EU: European Union; SDI: Socio-Demographic Index; Alabama: AL; Alaska: AK; Arizona: AZ; Arkansas: AR; California: CA; Colorado: CO; Connecticut: CT; Delaware: DE; District of Columbia: DC; Florida: FL; Georgia: GA; Hawaii: HI; Idaho: ID; Illinois: IL; Indiana: IN; Iowa: IA; Kansas: KS; Kentucky: KY; Louisiana: LA; Maine: ME; Maryland: MD; Massachusetts: MA; Michigan: MI; Minnesota: MN; Mississippi: MS; Missouri: MO; Montana: MT; Nebraska: NE; Nevada: NV; New Hampshire: NH; New Jersey: NJ; New Mexico: NM; New York: NY; North Carolina: NC; North Dakota: ND; Ohio: OH; Oklahoma: OK; Oregon: OR; Pennsylvania: PA; Rhode Island: RI; South Carolina: SC; South Dakota: SD; Tennessee: TN; Texas: TX; Utah: UT; Vermont: VT; Virginia: VA; Washington: WA; West Virginia: WV; Wisconsin: WI; Wyoming: WY.State ranking bump plots
Ranks of jurisdictions by point estimate in 1990, 2000, 2010 and 2021, available from the project’s GitHub repository.Software
R version and attached package versions used for the analysis (table below).Table 3: Software used for the analysis (Supporting Info 9). Package
Version
R
4.6.1
flextable
0.9.12
paletteer
1.7.0
gganimate
1.0.11
furrr
0.4.0
future
1.70.0
ggnewscale
0.5.2
geofacet
0.2.4
glue
1.8.1
scales
1.4.0
lubridate
1.9.5
forcats
1.0.1
stringr
1.6.0
dplyr
1.2.1
purrr
1.2.2
readr
2.2.0
tidyr
1.3.2
tibble
3.3.1
ggplot2
4.0.3
tidyverse
2.0.0
GATHER checklist
Completed GATHER checklist with manuscript locations: filegather_checklist.md, also available from the project’s GitHub repository.
7 List of abbreviations
- ACS – American Cancer Society
- ASR – age-standardized rate
- BC – bladder cancer
- DALYs – disability-adjusted life years
- GATHER – Guidelines for Accurate and Transparent Health Estimates Reporting
- GBD – Global Burden of Disease
- GHDx – Global Health Data Exchange
- SDI – Socio-Demographic Index
- SEER – Surveillance, Epidemiology, and End Results
- UI – uncertainty interval
- YLDs – years lived with disability
- YLLs – years of life lost
- AL – Alabama
- AK – Alaska
- AZ – Arizona
- AR – Arkansas
- CA – California
- CO – Colorado
- CT – Connecticut
- DE – Delaware
- DC – District of Columbia
- FL – Florida
- GA – Georgia
- HI – Hawaii
- ID – Idaho
- IL – Illinois
- IN – Indiana
- IA – Iowa
- KS – Kansas
- KY – Kentucky
- LA – Louisiana
- ME – Maine
- MD – Maryland
- MA – Massachusetts
- MI – Michigan
- MN – Minnesota
- MS – Mississippi
- MO – Missouri
- MT – Montana
- NE – Nebraska
- NV – Nevada
- NH – New Hampshire
- NJ – New Jersey
- NM – New Mexico
- NY – New York
- NC – North Carolina
- ND – North Dakota
- OH – Ohio
- OK – Oklahoma
- OR – Oregon
- PA – Pennsylvania
- RI – Rhode Island
- SC – South Carolina
- SD – South Dakota
- TN – Tennessee
- TX – Texas
- UT – Utah
- VT – Vermont
- VA – Virginia
- WA – Washington
- WV – West Virginia
- WI – Wisconsin
- WY – Wyoming
References
8 Declarations
8.1 Ethics approval and consent to participate
This study utilized publicly available, de-identified data from the Global Burden of Disease (GBD) study. As the data do not involve direct human participation and are anonymized, institutional ethics approval and informed consent were not required. The study was conducted in accordance with the principles of the Declaration of Helsinki.
8.2 Consent for publication
Not applicable.
8.4 Acknowledgments
None.
8.5 AI usage
No AI tools were used in the preparation of the first version of this manuscript’s code. All code was written manually by the first author, Alireza Sadeghi, with no AI involvement in the underlying logic. During the first revision, agentic AI was used solely for linting, styling, formatting, and commenting. Every code change was reviewed and tested by the first author. Readers can inspect the changes between the commit at submission and the commit at revision via the project’s GitHub repository.
8.6 Transparency Statement
The lead author, Alireza Sadeghi, affirms that this manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.
8.7 Availability of data and materials
All data are publicly available through the Global Health Data Exchange (GHDx). The complete analytic workflow, including cleaned datasets and all code required to reproduce the analyses, figures, and tables presented in this manuscript, is available in the study’s GitHub repository (commit 16f52d8; archived at 10.5281/zenodo.23814418). Data are provided in R’s native binary format (RDS). Full reproducibility can be achieved by cloning the repository and rendering the project using the supplied scripts and documentation. The manuscript and its source code are available via the manuscript’s website: www.bladder-cancer.ir
8.8 Competing interests
The authors declare no conflicts of interest.
8.9 Funding
None.