The present study is a full economic evaluation that compared the MS rehabilitation strategy with the non-intervention strategy in Iran from a social perspective.
Modeling
We developed a Markov model to simulate the natural history of the disease and assess the cost-effectiveness of MS rehabilitation services. To design the model, different sessions were formed with the presence of clinical specialists and economic team, and after a discussion of different models, the best model was finalized.
To evaluate the clinical efficacy of rehabilitation in this study, based on the available evidence as well as other economic evaluation studies on MS, the Expanded Disability Status Scale (EDSS) was considered.
The designed decision-analytic model can be seen in Figure 1. The Markov model had five mutually exclusive health states including severe disease health state (EDSS> = 8), moderate disease health state (6 <EDSS <7.5), mild disease health state (3 <EDSS <5.5), very-mild health state (EDSS <3) and death. Patients in the first cycle were distributed between health states according to the initial distribution. The Markov process was such that individuals in each group and health state remain in the same state at the end of each cycle, transfer to another health state or die. The model time horizon was 5-year and the cycle duration was one month. Outcomes were Quality Adjusted Life Years (QALYs), and costs.
According to the available evidence on the rehabilitation interventions, one month was considered as cycle length. Also, the effectiveness of rehabilitation interventions (18) was considered in different transition probabilities of compared strategies. The time horizon in the base case analysis as mentioned was considered to be 5 years equivalent to 60 cycles.
The final outcome of the evaluation was the quality-adjusted life years (QALYs) and QALY values are defined in each health state according to the EDSS score and health-related quality of life.
To extract evidence on mortality rate and other evidence, according to the average age of MS patients in Iran, the average age of 30 years was considered for patients at the beginning of Markov cycles. The period and frequency of the rehabilitation in each group have been done in accordance with Iran’s Ministry of Health standards for MS rehabilitation services in each subgroup.
Model Assumptions and Parameters
According to the natural course of the disease and the rehabilitation protocol for MS patients in Iran, the efficacy of MS rehabilitation (as a transition to better health status) was considered in only one cycle, and in the following cycles, rehabilitation was considered with the aim of stabilizing the health condition.
In each health state, in addition to the costs of rehabilitation interventions, the costs of patients’ care and maintenance and the costs of absenteeism were also considered.
The costs of inpatient and outpatient rehabilitation costs were calculated based on the unit cost used in each intervention. To calculate the average costs per patient and to calculate the total cost of rehabilitation services, a bottom-up approach was used. For this purpose, the extracted data of relevant evidence and Iran’s Ministry of Health standards for MS rehabilitation services were used and the list of rehabilitation services and their frequency divided different levels of the disease severity was extracted based on the EDSS index and also the type of rehabilitation services, then by multiplying the number of rehabilitation services per patient in the tariff of services, the average cost of rehabilitation services per patient was calculated.
In this study, all tariffs, including rehabilitation services, were extracted using the public sector official tariffs of Iran’s Ministry of Health in 2019. Accordingly, the average cost of MS rehabilitation services per patient was calculated based on public sector tariffs. Rehabilitation cost items included hospitalization costs, physiotherapy, speech therapy, occupational therapy, monitoring, and maintenance costs.
Home care and maintenance costs were calculated based on the average salary of a full-time and part-time nurse at a home care center. Absenteeism costs were also calculated based on the monthly number of days off work of different groups and also based on the minimum wage in 2019. For the severe and moderate disease group, due to the relatively high level of disability and lack of relative independence, the number of days off work was considered 30 days per month, and for the mild group, 15 days were considered.
In general, according to Iran’s Ministry of Health standards for MS rehabilitation services, as well as the opinions of clinical specialists, the length of the rehabilitation period was considered 2 years, so the cost of rehabilitation is only included in the model for 2 years (24 cycles).
Other related parameters and variables include the transition probabilities between health states in compared strategies, mortality risk in each health state, initial distributions, the effectiveness of rehabilitation interventions in each health state and other parameters extracted from international evidence. In this regard, each parameter was searched separately based on keywords and specific strategies in scientific databases, studies that had appropriate evidence were classified and studied, and finally, the best available evidence was extracted. The values of parameters and variables of the model and their sources can be seen in Table 1.
Table 1: Model Inputs and Sources
Statistic variable
|
Base case
|
SD/(CI)
|
Distribution
|
Source
|
Initial Prevalence (Severe)
|
26.4%
|
|
|
(19)
|
Initial Prevalence (Moderate)
|
58.6%
|
|
|
Initial Prevalence (Mild/Very Mild)
|
15%
|
|
|
Annual discount rate
|
0.05
|
(0.03-0.12)
|
Beta
|
|
Time Horizon(years)
|
5 years (60 cycles)
|
(+5 years)
|
|
|
DMTs Efficacy
|
51%
|
±13%
|
Beta
|
(20)
|
Transition Probabilities
(Non-rehabilitation strategy)
|
|
|
|
|
Mild to Moderate
|
0.008
|
|
|
(21)and Calibration
|
Mild to Severe
|
0.00042
|
|
|
Mild to Very Mild
|
0.0126
|
|
|
Moderate to Mild
|
0.0069
|
|
|
Moderate to Severe
|
0.0045
|
|
|
Moderate to Very Mild
|
0.00041
|
|
|
Severe to Mild
|
0.000089
|
|
|
Severe to Very Mild
|
0.00000314
|
|
|
Severe to Moderate
|
0.0051
|
|
|
Very Mild to Mild
|
0.0046
|
|
|
Very Mild to Moderate
|
0.00043
|
|
|
Very Mild to Severe
|
0.0000099
|
|
|
Transition Probabilities
(MS Rehabilitation strategy)
|
|
|
|
|
Mild to Moderate
|
0.008
|
|
|
(21)and Calibration
|
Mild to Severe
|
0.00042
|
|
|
Mild to Very Mild
|
0.047
|
±0.0094
|
Beta
|
(19)
|
Moderate to Mild
|
0.047
|
±0.0094
|
Beta
|
(19)
|
Moderate to Severe
|
0.0045
|
|
|
(21) and Calibration
|
Moderate to Very Mild
|
0.00041
|
|
|
Severe to Mild
|
0.000089
|
|
|
Severe to Very Mild
|
0.00000314
|
|
|
Severe to Moderate
|
0.04
|
±0.008
|
Beta
|
(19)
|
Very Mild to Mild
|
0.0046
|
|
|
(21) and Calibration
|
Very Mild to Moderate
|
0.00043
|
|
|
Very Mild to Severe
|
0.0000099
|
|
|
Mortality
|
|
|
|
|
Probability of Death (Severe State)
|
0.00040875
|
|
Beta
|
(22, 23)
|
Probability of Death (Moderate State)
|
0.000201375
|
|
Beta
|
Probability of Death (Mild State)
|
0.00012875
|
|
Beta
|
Probability of Death (Very Mild State)
|
0.000113625
|
|
Beta
|
Costs($)
|
|
|
|
|
Monthly cost of care and maintenance and patient absenteeism (Mild State)
|
11.164
|
±1.116
|
Gamma
|
Survey and Calibration
|
Monthly cost of care and maintenance and patient absenteeism (Moderate State)
|
515.364
|
±51.536
|
Gamma
|
Survey and Calibration
|
Monthly cost of care and maintenance and patient absenteeism (Severe State)
|
807.454
|
±80.745
|
Gamma
|
Survey and Calibration
|
Monthly Outpatient Rehabilitation Cost (Mild State)
|
10.191
|
±1.019
|
Gamma
|
Survey and Calibration
|
Monthly Outpatient Rehabilitation Cost (Moderate State)
|
13.604
|
±1.36
|
Gamma
|
Survey and Calibration
|
Monthly Outpatient Rehabilitation Cost (Severe State)
|
13.604
|
±1.36
|
Gamma
|
Survey and Calibration
|
Monthly Inpatient Rehabilitation Cost (Mild State)
|
462.408
|
±46.241
|
Gamma
|
Survey and Calibration
|
Monthly Inpatient Rehabilitation Cost (Moderate State)
|
520.249
|
±52.025
|
Gamma
|
Survey and Calibration
|
Monthly Inpatient Rehabilitation Cost (Severe State)
|
520.249
|
±52.025
|
Gamma
|
Survey and Calibration
|
Annual Utilities
|
|
|
|
|
Severe Health State
|
0.155
|
±0.14
|
Beta
|
(24)
|
Moderate Health State
|
0.385
|
±0.09
|
Beta
|
Mild Health State
|
0.606
|
±0.081
|
Beta
|
Very Mild Health State
|
0.81
|
±0.065
|
Beta
|
Economic Evaluation Analysis
The incremental cost-effectiveness ratio (ICER) index was used to analyze and determine the most cost-effective strategy according to the cost and outcome of each strategy. The formula of this index is as follows:
ICER=C1-C2/E1-E2
In this regard, C represents the cost of strategies 1, 2 and E represents the effectiveness of the strategies. The calculated ICER was compared with the cost-effectiveness threshold and the most cost-effective strategy was determined. It should be noted that in this study, the cost-effectiveness threshold was considered according to the World Health Organization proposal for developing countries as 1-3 times the GDP per capita. According to the World Bank report in 2019, the threshold was considered equivalent to $3114.623 (Iran’s GDP per capita in 2019) (25).
Sensitivity Analysis
Due to the uncertainty about some parameters used in the model, Deterministic Sensitivity Analysis (DSA) and Probabilistic Sensitivity Analysis (PSA) were performed. DSA was performed by changing the Markov time horizon from 5 years to 10 years as scenario analysis. Accordingly, in the new scenario, 120 monthly cycles were considered for the implementation of the Markov model.
The distributions of the uncertain parameters used in the PSA are given in Table 2. In cases where no evidence was found regarding the range of the variable amount, 10 to 20% of the mean of the parameters was considered as the standard deviation, and the appropriate distribution was selected according to the type of variable.
Preparation of preliminary data including cost data and cost calculations of different stages of rehabilitation services in each health state was conducted using Excel 2016 software. Modeling, analysis of base case cost-effectiveness results, as well as all steps of sensitivity analysis was performed using TreeAge software version 2020.