Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy

Abstract

Introduction. Cost-effectiveness analyses are vital in guiding decisions on treatment reimbursement. Natural history models are central to these, enabling the estimation of long-term costs and quality-adjusted life-years (QALYs) in the absence of lifetime trial data. Rare disease data are often scarce, resulting in disease progression being estimated through clinical assumptions. This study aims to evaluate how different modeling approaches influence cost-effectiveness estimates in rare disease health technology assessments (HTAs), using Duchenne muscular dystrophy (DMD) as a case study. Methods. A published economic model was used to compare 2 approaches for estimating disease progression: an assumption-based method relying on clinical plausibility and data-driven methods using data from 1,005 patients with DMD across 8 studies. Transition probabilities were estimated assuming increasing flexibility of study heterogeneity and compared with a simulated treatment cohort. Models were evaluated by comparing incremental cost-effectiveness ratios (ICERs) across approaches. No gold standard exists, so the plausibility of predictions was evaluated by comparing survival and disease progression estimates to published milestones. Results. Results showed that although the assumption-based model was clinically plausible, it predicted higher QALY gains (0.77) and lower ICERs (£1.96M per QALY) than data-driven methods did, which estimated QALY gains of 0.25, 0.26, 0.27, and 0.28 and ICERs of £6.2M, £6.2M, £5.8M, and £5.7M per QALY for the least to most flexible models, respectively. Limitations. No covariate effects or updated cost and utility data were incorporated, as the study purpose was a methodological comparison between approaches. Analyses were deterministic not probabilistic. Conclusions and Implications. This study emphasizes the critical role of model selection for HTA in rare diseases, showing that cost-effectiveness estimates from robust data-driven approaches can differ from clinically plausible assumption-based models.

Highlights

The choice of a natural history modeling method can drastically alter the cost-effectiveness results in rare disease evaluations.

A case study in Duchenne muscular dystrophy demonstrates how different modeling approaches yield divergent cost-effectiveness outcomes.

Assumption-based models, even when clinically plausible, may underestimate measures of cost-effectiveness and result in less reliable guidance for decision makers.

Data-driven models using real-world patient data provide more reliable estimates for health technology assessment (HTA).

This study offers practical guidance for analysts and HTA bodies on selecting robust modeling approaches in rare disease contexts.

Visual Abstract

This is a visual representation of the abstract.

Keywords

rare diseases economic evaluation Duchenne muscular dystrophy

Introduction

Cost-effectiveness analyses (CEAs) are a cornerstone of health technology assessment (HTA), informing decisions about the reimbursement of new treatments and therapies. Central to a CEA is the appropriate modeling of disease progression in the absence of lifetime clinical trial data. Natural history models (NHMs) enable the estimation of long-term costs and quality-adjusted life-years (QALYs) by predicting how patients move through health states over time. These are especially important when treatment effects extend beyond the duration of clinical trials. Data-driven approaches are likely to lead to more reliable estimation, but for rare disease CEAs, access to natural history data can be problematic,¹ meaning some models instead rely on clinical assumptions to estimate disease progression. Even with data access, standard statistical methods may not be applicable due to the scarcity and heterogeneity inherent in rare disease data.²

In the United Kingdom, the National Institute for Health and Care Excellence (NICE) employs specific thresholds, typically between £20,000 and £30,000 per QALY,³ to determine whether a new technology offers sufficient value for money. However, for rare diseases, these thresholds are often less stringent⁴ due to the higher costs associated with treatment development and acquisition and the limited evidence base inherent to such conditions.⁵ The highly specialized technology (HST) evaluation route allows for higher cost-effectiveness thresholds of £100,000 to £300,000 per QALY in exceptional cases.⁶ The upper limit is reserved for significant QALY gains of up to 30 QALYs or more,⁷ which is a very high threshold to meet. To be classed as an HST, the disease must be rare, with a very small number of patients eligible for treatment, no alternative approved treatment, and a shortened life expectancy.⁶ NICE are currently reviewing the HST criteria.⁸

Duchenne muscular dystrophy (DMD) is a rare neuromuscular disease that nearly exclusively affects males, with a global incidence of 1 in 3,600 to 6,000 live male births⁹ and a UK prevalence of about 2,500.¹⁰ Ataluren received NICE approval for the treatment of DMD in 2014 through HST3¹¹ and HST22.¹² Previous research has compared different statistical modeling approaches to a NHM for DMD.¹³ An economic model framework for DMD has also been published,¹⁴ although this is solely reliant on (biologically plausible) clinical assumptions to model disease progression. The original aim of the study was to describe different cost-effectiveness models based on varying measures of disease progression of patients.¹⁴ The DMD Functional Ability Self-Assessment Tool (DMDSAT)¹⁵ was compared with commonly used model structures based on ambulatory status and ventilation support.

Building on this work, this study focuses on DMD as a case study, using the previous economic model framework based on ambulatory status to show how data-driven models can improve on assumption-based models, even when the assumptions are based on clinically plausible evidence, and assesses the impact of this on the results of a CEA. The assumption-based and data-driven approaches can yield substantially different results due to the nature of their underlying inputs. Assumption-based models rely on expert opinion and clinical plausibility, which may oversimplify disease progression and omit variability observed in real-world data. In contrast, data-driven models incorporate observed patient-level transitions, capturing heterogeneity and potentially more accurate estimates of disease progression. These differences can significantly affect cost-effectiveness outcomes, particularly in rare diseases where small changes in assumptions or data can lead to large shifts in incremental cost-effectiveness ratios (ICERs). This distinction is central to the contribution of this article. While data-driven approaches are expected to provide more robust and credible estimates, consideration is required of the reasons behind deviations between differences (both between data-driven methods and when compared with assumption-based approaches). Unlike prior studies that focused on either model development or economic evaluation in isolation, this analysis integrates both elements to highlight methodological implications for HTA. In addition, the potential for testing and identifying cost-effective treatments targeted at particular stages of a disease (rather than across the whole NHM) is explored. Sensitivity analyses explore situations in which the model choice is most pertinent and how cost-effective thresholds might be reached in rare disease HTA.

This study aims to evaluate how different natural history modeling approaches influence cost-effectiveness estimates in the context of rare diseases, using DMD as a case study. Specifically, we compare assumption-based and data-driven methods for estimating disease progression and assess their impact on ICERs. The findings are intended to inform health economists, analysts, and decision makers involved in HTA, particularly in rare disease contexts in which data limitations often necessitate methodological compromises. By highlighting the implications of model choice, this article provides practical guidance for selecting robust modeling approaches in future economic evaluations.

Methods

This study evaluates the impact of different natural history modeling approaches on cost-effectiveness estimates. It focuses on a population of patients with DMD, beginning with diagnosis at age 5 y. The intervention is a hypothetical treatment assumed to reduce disease progression rates by 25% and cost £100,000 per patient per year, compared against a standard-of-care cohort based on either historical disease progression data or clinically plausible assumptions. The outcomes of interest are lifetime costs, QALYs and ICERs.

Economic Model Framework

To highlight differences in natural history modeling methods on decision models, the ambulatory status model from the published economic framework is used as a case study. This model defines 4 health states—early ambulatory (where patients are aged approximately 5–7 y), late ambulatory (aged approximately 8–11 y), early nonambulatory (aged approximately 12–15 y), and late nonambulatory (aged approximately 16 y and older)—with mortality as the absorbing state (Figure 1).

Figure 1

Natural history model structure based on ambulatory status, with hazard function h_k (t) on transition k.

Due to the rarity of DMD, disease progression data were unavailable to the authors publishing the economic model framework, in particular for the newly developed DMDSAT.¹⁴ Transition probabilities between states were instead estimated through a number of clinically justified assumptions identified from a targeted literature review. Patients in the economic model framework were assumed to enter the early ambulatory state at age 5 y, since this is the average age of diagnosis,¹⁶ and progress through subsequent states to reflect increasing disease severity. As such, no backwards transitions were allowed. Mortality was allowed from all disease states, and an increasing hazard of death as disease severity increases was assumed.¹⁷ Specifically, patients were assumed to spend on average 4 y in each of the intermediate states,¹⁸ corresponding to a constant annual transition probability of 0.159 for hazards $h_{1}$ to $h_{3}$ . A median survival of 25 was assumed,¹⁹ leading to constant annual transition probabilities of 0.094 for hazards $h_{4}$ to $h_{7}$ . There was no assumed mortality between the ages of 5 and 18 y, meaning the annual transition probabilities for these hazards were 0 for the first 13 y of the model. The mortality probabilities were also assumed to increase by 15% each year beyond the age of 35 y, so that no patients were observed to survive beyond 50 y in the standard-of-care cohort. These transitions were assumed to be consistent, with no variation across the study population.

In this study, the standard-of-care cohort simulation was replicated under these disease progression assumptions. To verify the accuracy of this replication, a comparison of the disease progression estimates (presented as transition probabilities) was conducted. For additional accuracy, the plots were uploaded to WebPlotDigitizer,²⁰ and transition probabilities at a range of ages were compared.

The economic model adopted a UK health care perspective, incorporating direct medical costs and quality-of-life estimates for each health state. Annual costs and utilities associated with being in the 4 intermediate health states that were calculated by Landfeldt et al.¹⁴ are shown in Table 1. The costs shown are direct medical costs, that is, from the health care perspective, and utilities are for patients (rather than caregivers).

Table 1

Annual Direct Medical Costs and Patient Utilities for the Duchenne Muscular Dystrophy Natural History Model Based on Ambulation

Health State	Annual Cost	Utility
Early ambulatory	£10,670	0.699
Late ambulatory	£11,190	0.607
Early nonambulatory	£16,490	0.224
Late nonambulatory	£27,590	0.146

The costs were obtained from 2012 data from the United States, exchanged to UK currency and inflated to 2015 values.¹⁴ The values in Table 1 have been taken directly from the original study and have not been altered or updated for this analysis, as the purpose is to compare different modeling techniques rather than to conduct a genuine CEA of a DMD health technology.

Simulation of Treatment Effects

A hypothetical treatment cohort was simulated in the economic model framework to represent a new health technology for DMD.¹⁴ The treatment was assumed to cost £100,000 per patient per year (in line with some rare disease treatments at the time of the economic framework publication²¹ and approximately half the current cost of ataluren¹²). The treatment was assumed to reduce all intermediate annual transition probabilities by 25% (consistent with the observed efficacy of glucocorticoid treatment in patients with DMD²²) and to not affect utility values. A treatment cohort was simulated in the same way in this study to enable a demonstrative CEA to be conducted. The hypothetical treatment was assumed not to affect utility values or introduce adverse effects. To verify the accuracy of the replication of treatment and standard-of-care cohorts, the estimated total costs and QALYs for each cohort were compared with the estimates of Landfeldt et al.¹⁴

Disease Progression Data

As a comparison to the assumption-based approach adopted in the economic model framework, transition probabilities were also estimated from a collection of datasets on disease progression in patients with DMD. The datasets were made available from the Cooperative International Neuromuscular Research Group (CINRG)²³ and C-PATH Duchenne Regulatory Science Consortium (D-RSC)²⁴ and are described in detail in previous studies.^13,25 To summarize, disease progression data were available from observational and natural history studies, placebo arms of clinical trials, and clinical test data. A total of 1,005 patients from 8 studies were included across a range of ages (4–34 y), study periods (2004–2018), and locations (North America, South America, Europe and, Asia) and contained patients on and not on steroids. The following variables were used to map patients to disease states: ability to stand from supine, ability to walk/run 10 m, North Star Ambulatory Assessment score,²⁶ Brooke score,²⁷ forced vital capacity percent predicted, and ventilation status (none/nighttime/full-time). Further details of the C-Path datasets are in the Supplementary Materials, and D-RSC data access applications can be made via https://c-path.org/tools-platforms/d-rsc-database/.

Methods for Estimating Disease Progression

Four previously evaluated methods for combining individual patient data from multiple sources¹³ were used to calculate disease progression through transition probabilities between the states described in Figure 1. These consisted of a model that does not adjust for study source, a 1-stage frailty model, a 2-stage model assuming proportional baselines between studies, and a 2-stage model assuming stratified baselines. All the models assumed a Weibull baseline hazard function, and no patient-level covariate information was included (covariate information could be incorporated into all of these models but was not consistently available from the datasets, a common rare disease issue¹). The 4 models are formulated as follows (where $h_{jk}$ corresponds to the hazard $h_{k}$ on transition $k$ for study $j$ in Figure 1):

No-adjustment model:

h_{jk} (t) = λ_{k} γ_{k} t^{γ_{k} - 1}

One-stage frailty model:

h_{jk} (t) = λ_{k} γ_{k} t^{γ_{k} - 1} \exp (α_{jk}), α_{jk} N (0, σ_{k}^{2})

Two-stage proportional model:

h_{jk} (t) = λ_{jk} γ_{k} t^{γ_{k} - 1}

\log ({\hat{λ}}_{jk}) N (\log ({\hat{λ}}_{k}), σ_{k}^{2})

Two-stage stratified model:

h_{jk} (t) = λ_{jk} γ_{jk} t^{γ_{jk} - 1}

(\begin{matrix} \log ({\hat{λ}}_{jk}) \\ \log ({\hat{γ}}_{jk}) \end{matrix}) ~ N ((\begin{matrix} 0 \\ 0 \end{matrix}), (\begin{matrix} σ_{λ_{k}}^{2} & Cov (\log ({\hat{λ}}_{k}), \log ({\hat{γ}}_{k})) \\ Cov (\log ({\hat{λ}}_{k}), \log ({\hat{γ}}_{k})) & σ_{γ_{k}}^{2} \end{matrix}))

The no-adjustment model assumes that the transition probabilities of moving from one state to the next are the same across all studies, meaning all studies have the same Weibull baseline hazard function for transition $k$ . The 1-stage frailty model assumes that studies vary by an unobserved random effect $α_{jk}$ , meaning their baseline hazards differ by a proportional term $\exp (α_{jk})$ . Between-study heterogeneity is estimated in this model through $σ_{k}^{2}$ . The 2-stage proportional model also assumes that baseline hazards in each study are proportional but estimates $λ_{jk}$ directly in each study before combining these estimates ${\hat{λ}}_{jk}$ (on the log-scale) in a second stage to obtain a pooled estimate of $λ_{k}$ (with between-study heterogeneity again estimated through $σ_{k}^{2}$ ). The 2-stage stratified model allows the baseline hazards to vary beyond proportionality, estimating a separate shape and scale parameter for each study on each transition, and combines these estimates ${\hat{λ}}_{jk}$ and ${\hat{γ}}_{jk}$ (again on the log-scale) in a second stage through a multivariate meta-analysis²⁸ to obtain pooled estimates of $λ_{k}$ and $γ_{k}$ . Between-study heterogeneity is now estimated through $σ_{λ_{k}}^{2}$ and $σ_{γ_{k}}^{2}$ .

A previous simulation study identified the 2-stage models (particularly the 2-stage stratified model) as the least biased estimators of transition probabilities among the rare disease scenarios considered, while the 1-stage frailty model showed low absolute bias conditional on convergence.¹³ The no-adjustment model was the most biased estimator of transition probabilities.

These 4 models were applied to the disease progression datasets. The datasets had previously been assigned to health states for DMD from a more granular NHM with 8 intermediate states²⁹ (according to clinical and test records), published as part of Project HERCULES.³⁰ Details of how patients were mapped to these states based on the available variables have been published previously.²⁹ Disease progression could be estimated more precisely using the HERCULES model; however, this would have prevented a direct comparison with the assumption-based estimates of disease progression from the economic model framework and would have required additional estimates of costs and utilities for the extra health states in the HERCULES model. The 8 intermediate health states from the HERCULES model were therefore collapsed into the 4 intermediate states in the ambulatory model of Figure 1 (with clinical input) as follows. States 1 (early ambulatory) and 2 (late ambulatory) of the ambulatory model are equivalent to states 1 and 2 of the HERCULES model. Patients in states 3 to 6 of the HERCULES model were assigned to state 3 (early nonambulatory) of the ambulatory model, since these health states correspond to patients who cannot walk but have better health than patients in states 7 or 8 of the HERCULES model, who were as such assigned to state 4 (late nonambulatory). Observations from the datasets were then mapped from the HERCULES model health states to the ambulatory model health states accordingly.

The 4 methods produced 4 models estimating disease progression under standard of care based on patient data. No treatment data were available, and so disease progression under a treatment were simulated in the same way as the original economic framework to allow the CEA to be conducted.¹⁴

To evaluate the impact of different modeling approaches on disease progression, predicted health state occupancies and lengths of stay (LOS) were calculated for each method. These visualizations allow for direct comparison of how patients are expected to transition through disease stages under each model. The LOS estimates feed directly into the calculation of costs and QALYs, meaning differences in LOS estimation between methods will alter the corresponding ICERs.

Cost-Effectiveness, Sensitivity, and Threshold Analyses

Once annual transition probabilities for the standard-of-care and treatment cohort had been estimated from the C-PATH data using the 4 models and from the previous assumption-based approach, the annual costs and utilities from Table 1 were applied to determine lifetime costs and QALYs. An annual discount rate of 3.5% was assumed. From these, ICERs were calculated for each method.

Sensitivity analyses were conducted to vary assumptions around treatment duration and efficacy, the discount rate, and annual costs/utilities. Scenarios reported in the original economic framework were replicated, including increasing or decreasing the discount rate, state utility values and baseline (direct [medical] costs), as well as treatment efficacy on mortality.¹⁴ In addition, treatment effects were assumed to affect only transitions out of one state at a time, representing a targeted treatment. The sensitivity analyses were conducted deterministically.

Finally, threshold analyses investigated what total treatment cost, annual treatment cost, and total QALY gain would be required to obtain a cost-effective ICER for each of the methods in the base-case scenario.

The primary aim of the analysis is to explore how different natural history modeling approaches influence cost-effectiveness outcomes. Rather than identifying a single superior model, the study compares 5 approaches, 1 assumption-based and 4 increasingly flexible data-based, within a consistent economic framework. While no formal gold standard exists, the plausibility of each model’s predictions is assessed against published literature and clinical expectations. This comparative approach highlights the sensitivity of cost-effectiveness results to the flexibility of modeling choices and aims to inform best practices in rare disease HTA.

Results

Disease Progression Estimates

Disease progression from the 4 models is shown in Figure 2, alongside the disease progression estimated from the assumptions-based approach. Patients were assumed to start in the early ambulatory state at age 5 y. Equivalent plots for the treatment cohort are available in Figure 1 of the Supplementary Materials. Table 2 gives the total LOS in each of the intermediate states (calculated as the shaded areas in each of the stacked plots in Figure 2) alongside the median survival from each of the 5 methods. Also tabulated are estimates from the simulated treatment cohort.

Figure 2

Health state occupancies after age 5 y across the ambulatory disease states from the assumption based and 4 model-based approaches for the standard of care cohort: (a) assumption-based approach,¹⁴ (b) no-adjustment model, (c) 1-stage frailty model, (d) 2-stage proportional model, and (e) 2-stage stratified model.

Table 2

Mean Lengths of Stay from the Age of 5 y in the 4 Intermediate Health States and Mean/Median Survival for the 5 Methods, across the Standard-of-Care and Simulated Treatment Cohorts

Method	Assumption Based		No Adjustment		One-Stage Frailty		Two-Stage Proportional		Two-Stage Stratified
Cohort	SoC	Trt	SoC	Trt	SoC	Trt	SoC	Trt	SoC	Trt
Early ambulatory	5.5	7.2	10.7	11.4	12.3	13.1	11.4	12.2	10.7	11.4
Late ambulatory	5.3	6.0	3.4	3.4	3.2	3.2	3.2	3.3	2.5	2.5
Early nonambulatory	4.1	4.0	3.1	3.2	3.4	3.5	3.2	3.3	3.9	4.1
Late nonambulatory	7.2	4.9	7.6	7.1	5.2	4.3	6.0	5.2	6.5	5.6
Mean survival	22.1	22.1	24.8	25.1	24.1	24.1	23.8	24.0	23.6	23.6
Median survival	25.0	25.0	25.9	26.1	28.5	28.5	28.5	28.5	28.2	28.3

SoC, standard-of-care cohort; Trt, treatment cohort.

One clear distinction between the probabilities derived from the assumption-based approach and those from the data-driven methods is the lack of mortality until age 18 y from the former; in Figure 2a, patients aged 18 y are distributed fairly evenly across the 4 intermediate health states, but no patients are predicted to die, whereas the other 4 methods predict a nonnegligible possibility of mortality from as early as age 7 y and probabilities of between 12% and 23% of mortality by age 18 y. These predictions are more closely aligned with literature estimates, which have estimated a probability of death by age 18 y of approximately 30%.¹⁷ The no-adjustment model estimates a 10% probability of still being alive at 50 y, which does not match clinical practice and is likely a reflection of the simpler modeling assumptions in this method. Transition probabilities from the original and reproduced assumption-based model were in agreement across a patient’s lifetime.

Patients spend longer in the early ambulatory state under the data-based models than in the assumption-based model, which could be due to different state definitions (clinical plausibility from published guidelines¹⁸ versus clinical test data). Patients are also distributed much more evenly across intermediate states using the assumption-based approach, which is expected since the transition probabilities were determined by assuming patients move through each state on average every 4 y. Predictions from the other models suggest patients in the C-Path datasets spend much longer in the early ambulatory state across their lifetime. There is also some variation between the predictions from the 4 models, with the no-adjustment model having much longer tails than the other 3 models. The median survival estimates from the more complex statistical models are more in line with literature estimates of median survival than from the assumption-based and no-adjustment models.^17,31

The overall mean and median survival across the standard-of-care and treatment cohorts differs notably for only the no-adjustment methods, with a mean and median difference of 0.3 y and 0.2 y, respectively. The lack of noticeable difference in median survival between the 2 cohorts is because the treatment cohort is assumed to have no impact on mortality; in general, patients in the treatment cohort spend longer in the ambulatory states and (as a result) shorter in the nonambulatory states than patients in the standard-of-care cohort.

CEA Results

Table 3 shows the lifetime costs and QALYs estimated for the standard-of-care and treatment cohorts, lifetime differences in costs and QALYs between the standard-of-care and (simulated) treatment cohort, and the ICER obtained from these differences for each of the 5 methods.

Table 3

Estimated Lifetime Costs and QALYs in the Standard-of-Care and Treatment Cohorts, Differences in Costs and QALYs between the 2 Cohorts, and ICERs Obtained from the 5 Methods

Method	Total SoC Costs	Total Trt Costs	Total SoC QALYs	Total Trt QALYs	Δ Costs (£)	Δ QALYs	ICER (£)
Assumption based	£360,000	£1,880,000	6.76	7.53	£1,520,000	0.77	£1,960,000
No adjustment	£280,000	£1,810,000	8.29	8.54	£1,530,000	0.25	£6,190,000
One-stage frailty	£310,000	£1,890,000	8.91	9.17	£1,580,000	0.26	£6,190,000
Two-stage proportional	£320,000	£1,900,000	8.89	9.16	£1,580,000	0.27	£5,820,000
Two-stage stratified	£340,000	£1,900,000	8.43	8.71	£1,560,000	0.28	£5,650,000

SoC, standard-of-care cohort; Trt, treatment cohort.

The estimates from the assumption-based model accurately replicate the estimates from the original economic framework model, which estimated a difference in costs of £1.52 million, a difference in QALYs of 0.79, and an ICER of £1.94 million. The small disparities are likely due to rounding differences in the software (Excel versus Stata) used to calculate the values.

The total difference in costs between the treatment and standard-of-care cohorts was fairly consistently estimated across all models (approximately £1.5–£1.6 million over the lifetime of a patient). This is due to the methods predicting that patients will live for approximately the same length of time; the median survival is slightly lower for the assumption-based and no-adjustment methods, leading to a slightly lower difference in costs. The differences are largely driven by the high cost of the treatment relative to the standard-of-care costs, and so the treatment cost was varied in a threshold analysis. The assumption-based approach also predicted a much higher QALY difference between the 2 cohorts (0.77 QALYs) than the other methods, all of which predicted a difference of 0.25 to 0.28 QALYs. This is mainly because the difference between the LOS for the treatment cohort and standard-of-care cohort was greater for the assumption-based method than for the data-based methods. In particular, under the assumption-based method, patients in the treatment cohort were estimated to spend 1.7 more years in the early ambulatory state and 2.3 fewer years in the late nonambulatory state than patients in the standard-of-care cohort were.

The 2-stage methods (which were previously observed to be the least biased estimators of transition probabilities in a simulation study focused on DMD disease progression¹³) gave similar estimates for the differences in costs and QALYs. The 1-stage frailty model estimated a slightly higher difference in costs, whereas the assumption-based and no-adjustment approaches gave a lower difference. Similar ICERs were obtained for the 2-stage methods (£5.7–£5.8 million per QALY), whereas the no-adjustment and 1-stage frailty models estimated higher ICERs of £6.2 million per QALY. Due to the difference in QALYs being so much higher for the assumption-based approach, the ICER from this method was vastly lower, at £2.0 million per QALY.

Sensitivity Analysis Results

Given the disparity between the assumption-based and data-based ICER estimates, the sensitivity analyses highlighted under which scenarios these differences might be more or less apparent in this case study. None of the methods were close to achieving a cost-effective ICER, even at the increased threshold of £100,000 per QALY, due to the very high treatment cost relative to quite a low efficacy. The resulting differences in costs, QALYs, and ICERs from sensitivity analysis scenarios relating to treatment duration and efficacy are displayed in Table A1 of the Supplementary Materials.

Figure A1 in the Supplementary Materials shows the resulting differences in lifetime costs and QALYs, and corresponding ICERs, from the remaining scenarios.

Assuming that the treatment also improved survival increased QALYs and costs. In the scenarios considered, treatments that lead to higher QALY gains were less affected by modeling differences, although mortality assumptions drove the greatest variation in ICERs. Varying discount rates and baseline costs had little impact, whereas varying utility values had a high impact on ICERs, with higher utilities reducing ICERs by up to 30%. These findings again underline the finding that greater disparities between restrictive and more flexible approaches are likely to be observed when QALY gains (or treatment costs) from the standard-of-care cohort to the treatment cohort are lower.

Threshold Analysis Results

Figure A2 in the Supplementary Materials shows the results of threshold analyses to determine the required difference in total costs, annual costs, and total QALYs, respectively, between the standard-of-care and treatment cohorts to obtain various ICERs, including a cost-effective ICER threshold of £100,000 per QALY. Threshold analyses relating to treatment cost (Figure A2a and A2b) showed clear differences between assumption-based and data-based approaches, while threshold analyses on QALYs were consistent across these (as expected, since lifetime costs were consistently estimated across models). In order to be cost-effective, the different models would require a QALY gain of between 8.2 and 8.8, a long way from the observed difference in QALYs in Table 3.

Discussion

This study demonstrated the impact that natural history modeling approaches can have on the results of CEAs for rare diseases, using DMD as a case study. By comparing assumption-based and data-driven methods, the choice of methodology for estimating disease progression was shown to significantly influences predictions of QALYs, costs, and ICERs. A model based on ambulatory status was used due to similarities between the model structure and the available disease progression data and because of its common use in clinical practice.¹⁸ The disease stages are specified in published DMD clinical care guidelines to assist with diagnosis and symptom management, including neuromuscular, cardiac, psychosocial management, and many others.^18,32,33

The assumption-based approach yielded substantially different results from the data-driven methods. It predicted longer durations in less severe health states for the treatment cohort relative to the standard-of-care cohort and perhaps most significantly assumed no mortality until the age of 18 y, leading to higher QALY differences and much lower ICERs (£1.96M per QALY versus £5.65–£6.19M per QALY for data-driven methods). This difference was driven by a much higher incremental QALY gain estimated by the assumption-based approach, with an approximately 3-fold increase in the QALY gain from this approach compared with the data-driven methods. Incremental costs were estimated relatively consistently between methods due to much higher lifetime treatment costs compared with standard-of-care costs. This finding underscores the importance of measuring both costs and QALYs alongside overall ICERs to determine what is driving differences in results between methodologies. The difference in QALYs between the assumption-based and data-based methods could be driven by assumptions around mortality, specifically, that the assumption-based model assumes no mortality until the age of 18 y. The large difference in LOS in the early ambulatory stay between the treatment and standard-of-care cohorts for this method is likely to be driven by this assumption. This may be compounded by the assumption of patients all entering the model at age 5 y, meaning differences in earlier health states of the model (corresponding to younger patients) have a greater impact on the CEA than for later health states, due to the economic discounting.

However, the optimistic outlook of the assumption-based approach is inconsistent with literature estimates and demonstrates the risks of relying on simplified, assumption-driven models. In particular, excess mortality has been observed in patients with DMD well before the age of 18 y.¹⁷ This study highlights that assumptions that appear clinically plausible and lead to believable transition probabilities may well be inaccurate, but this may not be obvious until a data-driven approach is applied. More thorough elicitation approaches to populate assumption-based models may lead to more similar ICER estimates with data-driven models (although a targeted literature review was conducted in the original economic model framework¹⁴).

While the primary conclusion that model choice significantly influences outcome estimates remains robust, we acknowledge some inconsistencies between the assumption-based and data-based models and observations in clinical practice regarding other clinical milestones beyond mortality. For instance, early ambulatory status was observed to be prolonged in the data-based models, but the median loss of ambulation was no lower than 20 y across all 5 models, much higher than has been observed in previous literature.³⁴ Although the data-driven approach may still be limited by the quality and completeness of available data, it arguably offers a closer reflection of real-world disease progression than the assumption-based alternative. These limitations can be partially addressed through the application of a 2-stage stratified model, which helps mitigate bias and improve the robustness of the findings.

Among the data-driven approaches, the 2-stage stratified model offered the most flexible and robust approach to accounting for heterogeneity across data sources. The results underline the importance of accurately adjusting for study-specific differences when pooling individual patient data, particularly for rare diseases for which data are often sparse and heterogeneous. Predictions from the frailty model are taken as conditional on zero frailty, in other words, an “average” study population relative to the studies that were used to populate the NHM. A new health technology would likely necessitate a new study, separate from the populations used to estimate disease progression in the NHM. Therefore, while it is possible that this average study population represents a reasonable baseline against which to compare a new study population, it is important to evaluate where in these populations the new study sits (especially given the heterogeneity that is often present in these study populations¹). If it is similar in nature to a more extreme population from one of the original studies, then study-specific predictions conditional on the frailty from this study would be more representative as a baseline NHM.

The costs were not altered from the conversion from 2012 US currency to 2015 UK currency, as originally published,¹⁴ meaning they are not representative of current medical costs in the United Kingdom. Similarly, a DMD-specific quality-of-life questionnaire, the DMD-QoL,³⁵ has since been developed that may allow more accurate state utility values to be determined.

For HTA in rare diseases, where randomized controlled trial data are often unavailable,^2,36 NHMs play a critical role in estimating disease progression and informing decision making. This study highlighted the dangers of assumption-driven methods, which can misrepresent disease progression and treatment effects, potentially compromising the robustness of cost-effectiveness estimates. Conversely, data-driven methods offer greater reliability but require careful selection of the most appropriate modeling approach based on biological plausibility, model checking, and evaluating robustness through sensitivity analysis. This is especially pertinent when it is not clear which approach is most favorable and if results from different methods vary. The methods in this study have been rigorously assessed in a rare disease setting, minimizing this risk. We recognize that many rare diseases lack the data infrastructure available for DMD, such as that provided by Project HERCULES. In such contexts, assumption-based models may be the only feasible option, and efforts should focus on improving transparency and elicitation methods to strengthen their validity.

A key limitation of this analysis is the absence of covariate effects in the NHMs. Just as the assumption-based approach oversimplified disease progression, omitting covariates may have the same impact. Incorporating patient-specific factors, such as genotype, baseline functional status, or corticosteroid use, could improve the precision of transition probability estimates and better capture variability in disease progression. This would be particularly important if costs and utilities varied substantially across patient subgroups. However, data availability often constrains the complexity of model structures (as in this case study), necessitating tradeoffs between realism and feasibility. Rare disease models, in particular, are frequently criticized for oversimplifications, such as omitting key health states. For instance, in previous NICE appraisals of DMD treatments, the exclusion of scoliosis as a health state was highlighted as a limitation.¹¹ Further stratification of the ambulatory states could enhance the granularity of future models, although this would increase data requirements. More flexible model specifications for between-state transitions beyond Weibull distributions could also have been considered, although the adopted approach did offer more flexibility than previous multistate modeling of DMD, where constant transition rates between states were assumed.^14,29

Another limitation is that analyses were deterministic rather than probabilistic. A probabilistic sensitivity analysis (PSA) could quantify uncertainty in transition probabilities, costs, and utilities and assess whether differences in ICERs between modeling approaches remain meaningful when uncertainty is considered. In this study, we did not conduct a PSA because the primary aim was methodological, comparing deterministic outputs under identical assumptions rather than providing a comprehensive cost-effectiveness evaluation and maintaining comparability with previous analysis.¹⁴ Future methodological work in this area would be of interest, and applications of the methods to a CEAs should incorporate a PSA to better understand the robustness of conclusions and the potential overlap in ICER distributions across methods.

The results also emphasize the value of targeting treatments to specific disease stages. Sensitivity analyses revealed that interventions focused on states with high(er) utilities, such as the early ambulatory phase of DMD, reduced ICERs more effectively than treatments applied across all health states. This targeted approach may help novel therapies achieve cost-effectiveness thresholds, particularly in the context of high costs and modest QALY gains. An important further consideration highlighted by this study is the need for robust utility data, as utility values were shown to significantly influence the outcomes in sensitivity analyses. This underscores the value of early modeling, which can help identify critical data gaps such as the need for high-quality utility estimates. Addressing these gaps may require the proactive planning of targeted studies to generate the necessary evidence.

This study was not intended as a comprehensive cost-effectiveness evaluation but rather as a methodological investigation into the effects of natural history modeling approaches on economic outcomes for a particular case study. The cost and utility estimates, derived from older data, may not reflect current clinical and economic conditions. In addition, the analyses were deterministic rather than probabilistic, limiting the evaluation of uncertainty in model inputs. Future work could incorporate probabilistic sensitivity analyses and value-of-information analyses to better quantify the uncertainty in cost-effectiveness estimates. Moreover, extending the analysis to include Bayesian frameworks may allow for more robust modeling of heterogeneity and uncertainty.^37,38 Alternative diseases should also be explored to widen the context of the research.

Considerations for Model Selection

Selecting an appropriate modeling approach for estimating disease progression in CEAs, particularly in rare diseases, is critical as it can substantially influence estimates. When estimating disease progression, researchers should consider the following key factors. First, what is the level of data availability and quality—are individual patient data available, and are the data heterogeneous across data sources? Second, balance clinical plausibility with statistical robustness by ensuring assumptions are evaluated and the model is externally validated where possible. Third, consider model flexibility, such as adjustments for between-study variation or covariate inclusion/stratification. Fourth, ensure the methods are reported transparently to ensure reproducibility. Finally, evaluate how sensitive the cost-effectiveness results are to model selection and assumptions. For instance, would different models lead to different reimbursement decisions?

Conclusion

In the context of rare diseases, where robust evidence is often limited, the selection of natural history modeling methods is critical. This study emphasizes the importance of adopting advanced, data-driven approaches to accurately estimate disease progression and inform HTA decision making, highlighting the need for timely preparation of data and evidence generation plans. By demonstrating the impact of methodology on cost-effectiveness outcomes, this research provides a framework for improving the validity and reliability of economic evaluations in rare diseases, ultimately aiding in the development and reimbursement of new treatments.

Supplemental Material

sj-docx-1-mpp-10.1177_23814683261447231 – Supplemental material for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy

Supplemental material, sj-docx-1-mpp-10.1177_23814683261447231 for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy by Jonathan Broomfield, Keith R. Abrams, Michael J. Crowther, Michela Guglieri, Nicholas R Latimer and Mark J. Rutherford in MDM Policy & Practice

Supplemental Material

sj-docx-2-mpp-10.1177_23814683261447231 – Supplemental material for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy

Supplemental material, sj-docx-2-mpp-10.1177_23814683261447231 for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy by Jonathan Broomfield, Keith R. Abrams, Michael J. Crowther, Michela Guglieri, Nicholas R Latimer and Mark J. Rutherford in MDM Policy & Practice

Supplemental Material

sj-docx-3-mpp-10.1177_23814683261447231 – Supplemental material for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy

Supplemental material, sj-docx-3-mpp-10.1177_23814683261447231 for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy by Jonathan Broomfield, Keith R. Abrams, Michael J. Crowther, Michela Guglieri, Nicholas R Latimer and Mark J. Rutherford in MDM Policy & Practice

Supplemental Material

sj-docx-4-mpp-10.1177_23814683261447231 – Supplemental material for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy

Supplemental material, sj-docx-4-mpp-10.1177_23814683261447231 for Evaluating the Impact of Different Natural History Modeling Methods on Cost-Effectiveness Decisions: A Case Study in Duchenne Muscular Dystrophy by Jonathan Broomfield, Keith R. Abrams, Michael J. Crowther, Michela Guglieri, Nicholas R Latimer and Mark J. Rutherford in MDM Policy & Practice

Footnotes

Acknowledgements

The authors would like to thank C-Path Duchenne Regulatory Science Consortium members and study team members from participating CINRG sites, details of whom can be found in the Supplementary Materials. This study was supported by the National Institute for Health and Care Research (NIHR) Applied Research Collaboration East Midlands (ARC EM) and Leicester NIHR Biomedical Research Centre (BRC). The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care.

The authors declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: JB, KRA, and MJC have received consultancy fees from Duchenne UK. MG has received research funding from Sarepta, PTC, Muscular Dystrophy UK, NIH, and H2020; is (or has been) chief investigator/principal investigator for clinical trials of Pfizer, Italfarmaco, Santhera, Roche, ReveraGen, Dynacure, and Dyne; is a member of advisory boards for Pfizer, NS Pharma, and Dyne (honoraria through Newcastle University); and has received speaker honoraria from Sarepta, Italfarmaco, Novartis, and Roche. NRL received funding from Yorkshire Cancer Research (award reference number: S406NL). NRL reports receiving grants and consultancy fees paid to his institution from Bristol-Myers Squibb, MSD, Pfizer, Daiichi Sankyo, GSK, Sanofi, Abbvie, and Novartis outside the submitted work. The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Financial support for this study was provided by a fellowship from the National Institute for Health and Care Research (NIHR; Doctoral Research Fellowship NIHR300984) and an NIHR Development and Skills Enhancement Award (NIHR305039). The funding agreement ensured the authors’ independence in designing the study, interpreting the data, writing, and publishing the report.

Ethical Considerations and Consent to Participate

All studies within the D-RSC dataset are fully anonymized and compliant with the terms and conditions specified by the data owners, with more information at .

Consent for Publication

Not applicable.

ORCID iDs

Jonathan Broomfield

Michael J. Crowther

Data Availability Statement

Further details of the C-Path D-RSC datasets used are in the Supplementary Materials. D-RSC data access applications can be made via .

References

Griggs

Batshaw

Dunkle

, et al. Clinical research for rare disease: opportunities, challenges, and solutions. Mol Genet Metab. 2009;96(1):20–6.

Schieppati

Henter

Daina

Aperia

Why rare diseases are an important medical and social issue. Lancet. 2008;371(9629):2039–41.

McCabe

Claxton

Culyer

AJ.

The NICE cost-effectiveness threshold: what it is and what that means. Pharmacoeconomics. 2008;26(9):733–44.

Griffiths

Hendrich

Stoddart

Walsh

SC.

Acceptance of health technology assessment submissions with incremental cost-effectiveness ratios above the cost-effectiveness threshold. Clinicoecon Outcomes Res. 2015;31(7):463–76.

Clarke

Ellis

Brownrigg

The impact of rarity in NICE’s health technology appraisals. Orphanet J Rare Dis. 2021;16:218.

NICE. Interim process and methods of the highly specialised technologies programme updated to reflect 2017 changes. Available from: https://www.nice.org.uk/media/default/about/what-we-do/nice-guidance/nice-highly-specialised-technologies-guidance/hst-interim-methods-process-guide-may-17.pdf. [Updated 2017; accessed 1 February, 2025].

Zhang

Garau

International cost-effectiveness thresholds and modifiers for HTA decision making. 2020. Available from: https://www.ohe.org/publications/international-cost-effectiveness-thresholds-and-modifiers-hta-decision-making/. [Accessed 8 January 2025].

Perez-Kempner

Chowdhury

Lubojemska

Unpacking NICE’s review of the HST routing criteria: implications for manufacturers. 2024. Available from: https://www.parexel.com/insights/blog/unpacking-nices-review-of-the-hst-routing-criteria-implications-for-manufacturers. [Accessed 14 January, 2025].

Fox

Millington

Mahabeer

Ruiten

Duchenne muscular dystrophy. BMJ. 2020;368:1–9. DOI:10.1136/bmj.l7012

10.

National Health Service (NHS). Muscular dystrophy. 2021. Available from: https://www.nhs.uk/conditions/muscular-dystrophy/. [Accessed 14 January, 2025].

11.

NICE. Ataluren for treating Duchenne muscular dystrophy with a nonsense mutation in the dystrophin gene (NICE guideline). Available from: https://www.nice.org.uk/guidance/hst3. [Updated 2016; Accessed 28 January, 2025].

12.

NICE. Routine funding recommended for Duchenne muscular dystrophy gene therapy ataluren. Available from: https://www.nice.org.uk/news/articles/routine-funding-recommended-for-Duchenne-muscular-dystrophy-gene-therapy-ataluren. [Updated 2022; accessed 14 January, 2025].

13.

Broomfield

Abrams

Freeman

, et al. Modelling the multi-state natural history of rare diseases with heterogeneous individual patient data: a simulation study. Stat Med. 2024;43(1):184–200.

14.

Landfeldt

Alfredsson

Straub

Lochmuller

Bushby

Lindgren

Economic evaulation in Duchenne muscular dystrophy: model frameworks for cost-effectiveness analysis. Pharmacoeconomics. 2017;35(2):249–58.

15.

Landfeldt

Mayhew

Eagle

, et al. Development and psychometric analysis of the Duchenne muscular dystrophy functional ability self-assessment tool (DMDSAT). Neuromuscul Disord. 2015;25(12):937–44. Available from: https://www.clinicalkey.es/playcontent/1-s2.0-S0960896615007397; DOI:10.1016/j.nmd.2015.09.012

16.

Thomas

Conway

Fapo

, et al. Time to diagnosis of Duchenne muscular dystrophy remains unchanged: findings from the muscular dystrophy surveillance, tracking, and research network, 2000-2015. Muscle Nerve. 2022;66(2):193–7. Available from: https://onlinelibrary.wiley.com/doi/abs/10.1002%2Fmus.27532; DOI:10.1002/mus.27532

17.

Broomfield

Hill

Guglieri

Crowther

Abrams

KR.

Life expectancy in Duchenne muscular dystrophy: reproduced individual patient data meta-analysis. Neurology. 2021;97(23):e2304–14.

18.

Bushby

Finkel

Birnkrant

, et al. Diagnosis and management of Duchenne muscular dystrophy, part 1: diagnosis, and pharmacological and psychosocial management. Lancet Neurol. 2010;9:77–93.

19.

Eagle

Bourke

Bullock

, et al. Managing Duchenne muscular dystrophy—the additive effect of spinal surgery and home nocturnal ventilation in improving survival. Neuromuscul Disord. 2007;17(2007):470–5.

20.

Rohatgi

WebPlotDigitizer. Available from: https://automeris.io/. [Updated 2020; accessed 30 January, 2025].

21.

Picavet

Morel

Cassiman

Simoens

Shining a light in the black box of orphan drug pricing. Orphanet J Rare Dis. 2014;9(1):62. Available from: https://www.ncbi.nlm.nih.gov/pubmed/24767472; DOI:10.1186/1750-1172-9-62

22.

Wang

Silverstein Fadlon

Ulm

, et al. Online self-report data for Duchenne muscular dystrophy confirms natural history and can be used to assess for therapeutic benefits. PLoS Curr. 2014;1. DOI:10.1371/currents.md.e1e8f2be7c949f9ffe81ec6fca1cce6a

23.

The Cooperative International Neuromuscular Research Group. Duchenne natural history. Available from: https://cinrgresearch.org/Duchenne-natural-history/. [Accessed 8 July, 2025].

24.

Critical Path Institute. Duchenne regulatory science consortium (D-RSC). Available from: https://c-path.org/program/duchenne-regulatory-science-consortium/. [Updated 2017; accessed 8 July, 2025].

25.

Broomfield J. Modelling the natural history of rare diseases using disparate data sources: application to patients with Duchenne muscular dystrophy. 2024. Available from: https://figshare.com/articles/thesis/Modelling_the_Natural_History_of_Rare_Diseases_Using_Disparate_Data_Sources_pplication_to_Patients_with_Duchenne_Muscular_Dystrophy/27073477

26.

Ricotti

Ridout

Pane

, et al. The NorthStar ambulatory assessment in Duchenne muscular dystrophy: considerations for the design of clinical trials. J Neurol Neurosurg Psychiatry. 2016;87(2):149–55. Available from: https://jnnp.bmj.com/content/87/2/149.full; DOI:10.1136/jnnp-2014-309405

27.

Brooke

Griggs

Mendell

Fenichel

Shumate

Pellegrino

RJ.

Clinical trial in Duchenne dystrophy. I. The design of the protocol. Muscle Nerve. 1981;4(3):186–97. Available from: https://api.istex.fr/ark:/67375/WNG-42K8HSNJ-2/fulltext.pdf; DOI:10.1002/mus.880040304

28.

Jackson

Riley

White

Multivariate meta-analysis: potential and promise. Stat Med. 2011;30(20):2481–98.

29.

Broomfield

Hill

Chandler

, et al. Developing a natural history model for Duchenne muscular dystrophy. Pharmacoecon Open. 2024;8:79–89.

30.

Duchenne

UK.

Project HERCULES: HEalth research collaboration united in leading evidence synthesis. Available from: https://www.duchenneuk.org/project-hercules/. [Updated 2022; accessed 1 February, 2025].

31.

Landfeldt

Thompson

Sejersen

McMillan

Kirschner

Lochmuller

Life expectancy at birth in Duchenne muscular dystrophy: a systematic review and meta-analysis. Eur J Epidemiol. 2020;35:643–53.

32.

Bushby

Finkel

Birnkrant

, et al. Diagnosis and management of Duchenne muscular dystrophy, part 2: implementation of multidisciplinary care. Lancet Neurol. 2010;9:177–89.

33.

Birnkrant

Bushby

Bann

, et al. Diagnosis and management of Duchenne muscular dystrophy, part 1: diagnosis, and neuromuscular, rehabilitation, endocrine, and gastrointestinal and nutritional management. Lancet Neurol. 2018;17(3):251–67.

34.

Zambon

Ayyar Gupta

Ridout

, et al. Peak functional ability and age at loss of ambulation in Duchenne muscular dystrophy. Dev Med Child Neurol. 2022;64(8):979–88. Available from: https://onlinelibrary.wiley.com/doi/abs/10.1111%2Fdmcn.15176; DOI:10.1111/dmcn.15176

35.

Powell

Carlton

Rowen

Chandler

Guglieri

Brazier

JE.

Development of a new quality of life measure for Duchenne muscular dystrophy using mixed methods: the DMD-QoL. Neurology. 2021;96(19):e2438–50.

36.

Behera

Kumar

Soares

Sokol

Djulbegovic

Evidence-based medicine for rare diseases: implications for data interpretation and clinical trial design. Cancer Control. 2007;14(2):160–6.

37.

Biao

Bayesian Methods in Health Economics. Boca Raton (FL): Chapman & Hall; 2012.

38.

Fornacon-Wood

Mistry

Johnson-Hart

Faivre-Finn

O’Connor

Price

GJ.

Understanding the differences between Bayesian and frequentist statistics. Int J Radiat Oncol Biol Phys. 2022;112:1076–82.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.03 MB

0.81 MB

0.71 MB

0.03 MB