Predicting the intervention threshold for initiating osteoporosis treatment among postmenopausal women in China: a cost-effectiveness analysis based on real-world data.

This study built a micro-simulation Markov model to determine the treatment threshold of osteoporosis in postmenopausal women in Mainland China. Treatment with zoledronate is cost-effective when FRAX-based (Fracture risk assessment tool) fracture probability is over 7%.
INTRODUCTION: The purpose of this study is to estimate FRAX-based fracture probabilities in Mainland China using real-world data, at which intervention could be cost-effective.
METHODS: We developed a micro-simulation Markov model to capture osteoporosis states and relevant morbidities including hip fracture, vertebral fracture, and wrist fracture. Baseline characteristics including incidences of osteoporosis and distribution of risk factors were derived from the Peking Vertebral Fracture study, the largest prospective cohort study of postmenopausal women in Mainland China. We projected incidences of fractures and deaths by age groups under two treatment scenarios: 1) no treatment, and 2) zoledronate. We also projected total quality-adjusted life-years (QALY) and total costs including fracture management and osteoporosis drugs for cost-effectiveness analysis. Cost-effective intervention thresholds were calculated based on the Chinese FRAX model.
RESULTS: Treatment with zoledronate was cost-effective when the 10-year probability of major osteoporotic fracture based on FRAX was above 7%. The FRAX threshold increased by age from 51 to 65 years old, and decreased in elder age groups, ranging from 4% to 9%.
CONCLUSIONS: Using real-world data, our model indicated that widespread use of zoledronate was of both clinical and economic benefit among Chinese postmenopausal women. Using a FRAX-based intervention threshold of 7% with zoledronate should permit cost-effective access to therapy to patients and contribute to reducing the disease burden of osteoporosis in Mainland China.

Background

Osteoporosis is a major public health concern in Mainland China. It was estimated that there were up to 70 million osteoporotic patients, and over 200 million patients with osteopenia in China in 2006 [1]. Osteoporotic fracture is the most severe complication of osteoporosis and results in hospitalization, disability, and mortality. The prevalence of vertebral fracture in postmenopausal Chinese women increased from 13.4% at ages 50–59 years, to 58.1% at age 80 years or older [2]. The overall major osteoporotic fracture in Chinese population was estimated to be 2.69 million cases in 2015 [3], indicating an enormous disease burden of osteoporotic fractures in China.

Osteoporosis treatment beforehand exhibits efficacy in fracture risk reduction. Of all treatment regimens, bisphosphonate is the most commonly used drug and was predicted to be cost-effective in other models [4]. In one of our previous studies (unpublished), alendronate was dominated by zoledronate in cost-effectiveness analysis, suggesting zoledronate a better option for treating Chinese postmenopausal women. In Mainland China, there is no universally accepted policy for treating patients at high risk, and only a minority of osteoporotic people receive proper treatment due to lack of evidence-based guidelines and recommendations [5]. Fracture risk assessment tool (FRAX) is a popular tool worldwide, used to predict 10-year probability of osteoporotic fracture based on clinical risk factors and femoral neck bone mineral density (BMD) [6]. Several cost-effectiveness analyses based on Western settings suggest that osteoporotic treatment should be initiated when the 10-year major fracture probability, based on FRAX, is over 7-15% [7-10]. Previous studies on estimating the FRAX threshold for the Chinese population were only limited to clinical observations or expert consensus [11-14].

We introduced real-world data from the Peking Vertebral Fracture (PK-VF) study into modeling to improve the robustness of cost-effectiveness analysis. The PK-VF study is currently the largest prospective cohort study of postmenopausal osteoporotic women in Mainland China, which involved 1100 women and follow-up from 2008 to 2013, to capture baseline epidemiological data of osteoporosis and osteoporotic vertebra fracture [2, 15–17]. The epidemiological data of the cross-sectional survey in 2008 has been published, including prevalence of osteoporosis and osteoporotic vertebra fracture and the distribution of risk factors of osteoporosis in the population [2].

In this study, we built a micro-simulation Markov model based on the PK-VF study. We aim to estimate the intervention threshold of FRAX in Mainland China, at which treatment with zoledronate could be cost-effective.

Methods

Model approach

We developed an individual-based micro-simulation Markov model in NetLogo 5.3.1 (Evanston, IL) to evaluate cost-effectiveness of osteoporotic treatment with zoledronate under different initiation thresholds. Subjects in the model were simulated on the following aspects: age, BMI, osteoporosis status, fracture status, risk factors of fracture (including previous fracture, parent fracture, smoking, alcohol, use of glucocorticoid, and the history of rheumatoid arthritis and secondary osteoporosis), and medical costs (including drug costs and fracture management costs) (Fig. 1a). Subjects in the model could transit between different states, including states of non-fracture, fracture, and postfracture, as shown by arrows (Fig. 1b). We simulated 20,000 people representing Chinese postmenopausal women as a closed cohort with a lifelong time horizon. A discount rate was set as 3% per year for both costs and effectiveness. Cost-effectiveness was considered as ICER less than 3 times of per capita GNP, which was $20,000/QALY for 2018 China.

Fig. 1

Schematic of the model structure. a Subjects in the model were simulated on the aspects of demographic feature, health feature, risk factor feature, and treatment feature. b Subjects could transit between different states, including states of non-fracture, fracture, and postfracture, as shown by arrows. Multiple fractures were allowed to happen in the same subjects during the lifetime. All states would also be absorbed into death state. c The anti-osteoporosis treatment would be initiated to subjects when FRAX score was higher than the threshold, or when fracture happened in the subject. Once initiated, we treated those subjects for five years under each of the two treatment scenarios, including (1) no treatment and (2) zoledronate. Subjects might drop-off during the 5-year treatment, and the remaining effect would decrease linearly in the next 5-year posttreatment

Baseline population from real-world data

We incorporated characteristics from the PK-VF study into the model with a stochastic sampling method for baseline population setup. By initial setting up the population, individuals’ characteristics were distributed according to the PK-VF study, including eight risk factors related to osteoporotic fracture in FRAX that is body mass index (BMI), previous fracture, parent fracture hip, current smoking, glucocorticoid use, rheumatoid arthritis, secondary osteoporosis, and alcohol use. For each age group, we randomly selected a number of individuals of the same age, together with their osteoporotic status, accompanied risk factors and treatment histories, from PK-VF participants. The percentage of each clinical risk factor based on PK-VF cohort was listed in Supplementary Table 1. As most of the risk factors (e.g., parent fracture hip, current smoking, glucocorticoid use, rheumatoid arthritis, BMI, secondary osteoporosis, and alcohol use) may have sustained impact on the osteoporosis-related fractures, we assumed that the risk factors of individuals simulated in the model remain unchanged except for previous fractures, which is updated if osteoporotic fractures occurred. For example, we assumed the effect of smoking on the increased risk of osteoporosis remains sustained no matter the individual stopped smoking or not.

Fracture risks and mortality

Our model began with women in the no-fracture state and simulated yearly transition till death. Osteoporosis developed with a yearly probability by age derived from the age-specific incidence of osteoporosis from PK-VF study (Table 1). The incidence of vertebral fracture, hip fracture, and wrist fracture were from pooled data of the Hong Kong Osteoporosis Study, the Hefei osteoporosis project, and a Norwegian study (multiplied by 0.72 to adjust for the Chinese population) [18-20] (Table 1). Osteoporotic subjects were more likely to experience fractures, with the relative risks published by Melton et al. [21]. The risk of a subsequent fracture increased when a prior fracture occurred [22]. Age distribution and natural mortality rate were from the 2012 National Sample Survey on Population Changes [33]. Standardized mortality ratios (SMR) of vertebral fracture, hip fracture, and wrist fracture were set as 1.82, 2.43, and 1.42 respectively [23] (Table 1).

Table 1

Key parameters derived from Peking-Vertebral Fracture (PK-VF) study and from literatures

Parameters		Distribution
Prevalence of osteoporosis (%) (PK-VF study)	5.26 (50–54 years), 13.1 (55–59 years), 19.6 (60–64 years), 26.4 (65–69 years), 24.5 (70–74 years), 27.2 (75–79 years), 30.1 (80–84 years), 35.0 (85 + years)
Incidence of osteoporosis (%) (PK-VF study)	1.91 (50–54 years), 2.80 (55–59 years), 1.80 (60–64 years), 1.04 (65–69 years), 2.10 (70–74 years), 2.22 (75–79 years), 3.96 (80 + years)
Fracture incidence (‰)
Vertebral fracture [18]	2.19 (50–54 years), 3.13 (55–59 years), 5.16 (60–64 years), 5.64 (65–69 years), 8.74 (70–74 years), 12.05 (75–79 years), 21.19 (80–84 years), 26.89 (85–89 years), 27.10 (90 + years)	-
Hip fracture [19]	0.33 (50–54 years), 0.46 (55–59 years), 0.54 (60–64 years), 0.96 (65–69 years), 2.33 (70–74 years), 4.08 (75–79 years), 6.44 (80–84 years), 6.59 (85–89 years), 8.67 (90 + years)	-
Wrist fracture [20]a	4.76 (50–54 years), 7.32 (55–59 years), 11.16 (60–64 years), 12.95 (65–69 years), 13.17 (70–74 years), 13.87 (75–79 years), 15.01 (80–84 years), 15.10 (85–89 years), 13.97 (90 + years)	-
Osteoporosis attribution probabilities by fracture type [21]
Vertebral fracture	0.75 (range 0.40–0.80) (45–64 years), 0.75 (range 0.50–0.90) (65–84 years), 0.95 (range 0.60–0.95) (85 + years)	Triangular
Hip fracture	0.75 (range 0.20–0.85) (45–64 years), 0.85 (range 0.50–0.95) (65–84 years), 0.95 (range 0.60–0.95) (85 + years)	Triangular
Wrist fracture	0.60 (range 0.10–0.70) (45–64 years), 0.70 (range 0.35–0.80) (65–84 years), 0.70 (range 0.55–0.90) (85 + years)	Triangular
Relative risk of subsequent fracture following a prior fracture [22]
Vertebral fracture	2.0 (range 1.6–2.4)	Gamma
Hip fracture	2.0 (range 1.9–2.2)	Gamma
Wrist fracture	1.9 (range 1.6–2.2)	Gamma
Standardized mortality ratios (SMR) [23]
Vertebral fracture	1.82 (95% CI 1.52–2.17)	Gamma
Hip fracture	2.43 (95% CI 2.02–2.93)	Gamma
Wrist fracture	1.42 (95% CI 1.19–1.70)	Gamma
Quality-adjusted life year (QALY)
By age [24]	0.772 (55–59 years), 0.728 (range 0.582–0.874) (60–64 years), 0.702 (range 0.562–0.842) (65–69 years), 0.685 (range 0.548–0.822) (70–74 years), 0.669 (range 0.535–0.803) (75–79 years), 0.655 (range 0.524–0.786) (80–84 years), 0.643 (range 0.514–0.77) (85 + years)	Triangular
Vertebral fracture [25]	0.724 (range 0.667–0.779) (1st year), 0.868 (range 0.827–0.922) (subsequent years)	Beta
Hip fracture [25]	0.776 (range 0.720–0.844) (1st year), 0.855 (range 0.800–0.909) (subsequent years)	Beta
Wrist fracture [25]	1.000 (range 0.960–1.000) (1st year), 1.000 (range 0.930–1.000) (subsequent years)	Triangular
Treatment of zoledronate related parameters
Relative risk of fracture of zoledronate [26]	0.23 (range 0.14–0.37) (vertebral fracture), 0.59 (range 0.42–0.83) (hip fracture), 0.75 (range 0.64–0.87) (wrist fracture)	Lognormal
Adherence of zoledronateb	1
Persistence of zoledronatec [27]	1 (1st year), 0.748 (subsequent years)
COST (2019 US dollar)
Vertebral fracture 1st year [32]d	3411
Hip fracture 1st year [32]d	4514
Wrist fracture 1st year [32]d	1383
Zoledronate + Ca + Vit D [28]	492
Annual discount rate [29]	0.03 (range 0–0.05)	-

aIncidence of wrist fracture was from Norwegian population, and multiplied by 0.72 to adjust for Chinese population [20]

bAdherence referred to the extent to which the subject acted in accordance with the prescribed drug, with respect to dosage and timing [30]. The adherence of zoledronate was hypothesized to be 1, as zoledronate was injected once a year

cPersistence referred to the duration from initiation until discontinuation of the treatment with a certain drug [30]

dCosts of vertebral, hip, and wrist fracture referred to all direct medical costs and direct non-medical costs reported by Qu et al. in 2012 [31]. Direct medical costs included costs of outpatient and inpatient care, examination, medication, rehabilitation and physical therapy after fracture; and direct non-medical costs included traffic fee, preventive care foods, and specific equipment. All costs (collected in 2012) were converted to 2019 RMB by an inflation rate of 1.1466 [32]

Treatment strategies

We incorporated the treatment of zoledronate into the model, dosed as 5 mg a year, combined with calcium and vitamin D3 supplement (calcium 600 mg + vitamin D3 125IU per tablet), dosed once a day. The effect of zoledronate on reducing fracture risks was taken from clinical trials [26, 34, 35].

We calculated the 10-year major fracture probability (FRAX score) for each subject based on FRAX equation [6], which incorporated covariates including age and risk factors. Those with an estimated FRAX score higher than the threshold set were eligible for treatment. The antiosteoporosis treatment was initiated for subjects when FRAX score was higher than the threshold, or when fracture occurred in the subject (Fig. 1c). Once initiated, those subjects would be treated for five years under each of the two treatment scenarios: (1) no treatment or (2) zoledronate. The remaining effect was assumed to decrease linearly in the next 5 years after the treatment was over [36].

Imperfect adherence and persistence may affect the cost-effectiveness of the treatment and were thus incorporated in the model (Table 1) [27, 37]. Adherence was defined as the extent to which the subject acted in accordance with the prescribed drug, with respect to dosage and timing [30]. We set the adherence of zoledronate to be 1, as zoledronate was injected once a year. Subjects in the model may also drop out during the 5-year treatment, and persistence was defined as the duration from initiation until discontinuation of the treatment with a certain drug [30].

Cost assumptions

Costs of vertebral, hip, and wrist fracture referred to all direct medical costs and direct nonmedical costs reported by Qu et al. in 2012 [31]. Direct medical costs included costs of outpatient and inpatient care, examination, medication, rehabilitation, and physical therapy after fracture; direct non-medical costs included traffic fee, preventive care foods, and specific equipment. All fracture costs (data collected in 2012) were converted to 2018 RMB by a total inflation rate of 1.1466 reported by inflation.eu [32]. Costs of zoledronate (100 ml: 5 mg per year) and calcium + vitamin D supplement (calcium 600 mg + vitamin D3 125IU per tablet per day) was obtained from the Chinese market drug price database (access date Feb 07, 2019) [28]. RMB was converted to US dollars at the exchange rate of 6.7426 (Feb 07, 2019).

Utilities

Quality-adjusted life year (QALY) is a measurement of both quality and quantity of life lived, varying from 0 to 1, where 0 referring to death and 1 referring to perfect health. The age-specific QALY for the general population (Table 1) was based on the National Health Services Survey 2008 in China [24]. The QALY of different fracture types in the general population and osteoporotic population (Table 1) were obtained from a meta-analysis [25].

Presentation of results

We projected QALY, as a measurement of health and total costs of both fracture management and osteoporosis drugs for the cost-effectiveness analysis. The incremental cost-effectiveness ratio (ICER) was calculated as (cost A–cost B)/(QALY A–QALY B). The willingness to pay (WTP) was assumed to be 3 times of per capita GNP, which was $20,000/QALY for 2018 China.

Sensitivity analysis

To test the robustness of the model to uncertainties of input data, one-way sensitivity analysis comparing the treatment of zoledronate to no treatment was performed on key model parameters, including fracture incidence, incidence of osteoporosis, standard mortality ratio, relative risk of subsequent fracture following a prior fracture, QALY by age, QALY by fracture types, relative risk of fracture with treatment, persistence, treatment year, adherence, drug cost, and annual discount. The upper and lower limits of each key parameter are listed in Table 1 and Appendix.

Results

Validation

Face validation was done by consulting experts on osteoporosis prevention and treatment. Internal validation was performed on comparison of annual fracture incidences from model prediction to those from input data (Supplementary Table 2). The projected incidences were all within ± 5% of the reference value. External validation was performed on comparison of life expectancy, prevalence of osteoporosis, and lifetime fracture risk predicted from model to those from literature (Supplementary Table 3) [38-40]. Cross-model validation was not performed, as no modeling work has been done on estimating the intervention threshold of Chinese mainland population.

Intervention FRAX threshold estimation

Comparing the no-treatment scenario to treatment with zoledronate, the total costs increased from $1407.4 to $1822.9 per person, and the corresponding QALYs increased from 9.232 to 9.269. The ICER was $26637.3, $22129.7, $20338.4, $19285.0, $18181.0, $16680.2, $15047.7, and $14447.7 per QALY at FRAX threshold 0.02, 0.06, 0.07, 0.08, 0.09, 0.1, 0.5, and 1, respectively. As Fig. 2 indicated, the overall treatment threshold was 7% for treatment with zoledronate to be cost-effective in comparison to no treatment.

Fig. 2

FRAX threshold for zoledronate. Treatment with zoledronate for patients with 10-year major fracture probability by FRAX over 0.07 was cost-effective in comparison with no treatment

Furthermore, we conducted cost-effectiveness analysis by age groups. As Fig. 3 and Supplementary Table 4 indicated, FRAX thresholds in different age groups ranged from 0.040 to 0.090. The threshold increased by age from 51 to 65 years old, and decreased in elder age groups.

Fig. 3

Heatmap of incremental cost-effectiveness ratio (ICER) for zoledronate by age-group and by FRAX threshold. Red lines indicated the willingness to pay (WTP) as $20,000/QALY

We also performed a subgroup cost-effectiveness analysis on subjects all with previous fracture history (Supplementary Table 5, Fig. 4). By comparing Supplementary Table 5 to Supplementary Table 4, treatment with zoledronate was found to be cost-saving in the youngest and oldest age groups with previous fracture, while in other age groups, the treatment threshold based on FRAX was similar between people with and without previous fracture.

Fig. 4

Heatmap of incremental cost-effectiveness ratio (ICER) for zoledronate in subjects with previous fracture by age-group and by FRAX threshold. Red lines indicated the willingness to pay (WTP) as $20,000/QALY

We performed sensitivity analysis on key model parameters, comparing the treatment of zoledronate to no treatment (Supplementary Table 6). As shown in the tornado plot (Fig. 5), the drug cost of zoledronate, incidence of vertebral fracture in osteoporosis, and QALY by age and persistence were the most sensitive parameters, but ICERs for all parameters remain below $23,000/QALY.

Fig. 5

Tornado plot of one-way sensitivity analysis on key model parameters, comparing treatment with zoledronate to no-treatment scenario. The upper and lower limits of all key model parameters run for the sensitivity analysis were listed in Table 1. ICER, incremental cost-effectiveness ratio; RR, relative risk; QALY, quality-adjusted life year

Discussion

In this study, we built an individual-based Markov model of osteoporosis simulating postmenopausal women in Mainland China, based on real-world data from the PK-VF study. We predicted the intervention threshold of FRAX in Mainland China at which treatment with zoledronate could be cost-effective and calculated the variation of the threshold in different age groups.

In Mainland China, there is no accepted criterion for initiating the treatment of osteoporosis. Current expert consensus in China is that anti-osteoporosis treatment should begin in any of the following conditions: (1) the occurrence of vertebral fragile fracture or hip fragile fracture; (2) BMD of femoral neck, whole hip, lumbar vertebrae, or distal radius < − 2.5; (3) − 2.5 < BMD < − 1 plus (1) occurrence of the proximal humerus, distal forearm, or pelvis fragile fracture or (2) 10-year hip fracture probability > 3% or 10-year major fracture probability > 20% predicted by FRAX [41]. However, there has been no cost-effectiveness analysis to determine a Chinese-specific threshold to initiate osteoporosis treatment based on FRAX. Thus, building a cost-effectiveness model that predicts FRAX thresholds for the Chinese population is both important and urgent.

Several cost-effectiveness analyses have been performed on predicting the disease burden of osteoporosis in postmenopausal women in Mainland China [3, 42–45], but none of them were based on real-world data of osteoporotic cohorts in the Chinese population. This study, however, was based on our previous Peking Vertebral Fracture (PK-VF) [2, 15–17], in which 1100 Chinese postmenopausal women were followed up from 2008 to 2013, and data were collected on the incidence and prevalence of osteoporosis and osteoporotic vertebra fracture, and the distribution of risk factors of osteoporosis in the population [2]. Introducing real-world data into modeling would significantly improve the simulation.

The determination of cost-effectiveness FRAX intervention threshold should be country-specific, considering the large variation of epidemiology and economic conditions across countries. For example, osteoporosis treatment was predicted to be cost-effective when 10-year major fracture exceeded 7% in UK [7], 15% in Switzerland [8], 10% in Greece [9], 8.8% in Portugal [10], and when 10-year hip fracture exceeded 3% in USA [46]. There have been a few studies estimating the intervening threshold of osteoporosis for Chinese population based on FRAX, but they were only limited to clinical observations or expert consensus without the evaluation of the economic and disease burden of osteoporosis [11-14].

No cost-effectiveness analysis has been done on estimating FRAX threshold in China so far. One of the reasons was that real-world data of the distribution of clinical risk factors were not available in China before [45]. There were eight clinical risk factors in FRAX model to determine a 10-year probability of major fracture and hip fracture of an individual, including age, sex, BMI, previous fracture, parent fracture hip, current smoking, glucocorticoid use, rheumatoid arthritis, secondary osteoporosis, and alcohol use (with or without the addition of femoral neck BMD). As risk factors of osteoporosis were always not independent of each other, it would lead to potential bias if we were to bring each risk factor separately into the model. For example, current smokers are more likely to be a regular alcohol drinker. To solve this issue, we generated a sampling table of risk factors based on our PK-VF cohort, and randomly sampled a group of risk factors as a whole from the sampling table for each subject simulated in the model.

There are various anti-osteoporosis drugs among which bisphosphonate is the most widely used one. Bisphosphonate functions in inhibiting the bone reabsorption process of osteoclasts. Alendronate, an oral bisphosphonate, is dosed once per week, and zoledronate is dosed once a year intravenously. Alendronate and zoledronate reduced the risk of vertebral fracture by 48% and 70% respectively in clinical trials [26, 47], and in one of our previous studies (unpublished), alendronate was dominated by zoledronate in cost-effectiveness analysis, suggesting zoledronate a better option for treating Chinese postmenopausal women.

As indicated in Fig. 2, FRAX threshold for the whole Chinese population was predicted to be 7%. FRAX thresholds differed in different age groups (Fig. 3 and Supplementary Table 4), as the number first increased in age groups younger than 65 years, and decreased in elder age groups. The reason behind this phenomenon is the bidirectional effect of aging on the resultant gain of QALYs, and thus on the ICERs, in a treatment-as-prevention mode, (1) early intervention with low FRAX threshold among younger patients, who have longer expected lifetime, resulted in more cumulative QALYs gained. (2) Intensive intervention with low FRAX threshold among older patients, who have higher fracture risk, resulted in a marked reduction of osteoporotic fractures and related deaths and thus in more QALYs gained. Therefore, we suggested additional resources and attentions shift to the middle aged or advanced aged.

In the external validation table (Supplementary Table 3), a gap was found between the projected life expectancy (16.49) and the reference life expectancy (16.9) in the age group of 65–69 years. We think the reason for the gap was that we used the overall death rate by age from a national survey [33], and applied osteoporosis, fractures, and deaths due to fractures to the simulated population additionally. Therefore, the simulated population presented a slightly higher death rate and thus a shorter life expectancy, especially among the elders. However, such an overestimated death rate would result in an underestimated benefit of anti-osteoporosis treatment, because it relies on long-term follow-up to observe the benefit of treatment of osteoporosis in preventing fractures, and overestimation of the death rate tends to underestimate the preventive gain of QALYs and the cost-effectiveness of treatment. Currently, treatment of osteoporosis in clinical practice was rather insufficient. With such a conserved estimation, we put forward a FRAX-based intervention threshold as 7%.

There are certain limitations to this study. First, not all data needed for modeling were available for the Chinese population. For example, the incidence of wrist fracture was adopted from the Norwegian population, and multiplied by 0.72 to adjust to the Chinese population [20]. And we generally stratified transition-probabilities to osteoporotic status from healthy status by age, regardless of those listed risk factors such as BMI, parent fracture, current smoking. Since only 67 new cases of osteoporosis were observed in the 5-year follow-up of PK-VF study, the impact of these covariates on the incidences was not available. And for the same reason, we took only age, previous fractures, treatment, and osteoporotic status as impact factors on transition probabilities of fractures but did not considered the impact of all listed risk factors above. Second, our real-world data were collected in Beijing population, which may not represent the complete picture of the Chinese population. Our group is now performing an ongoing national-wide epidemiological study on osteoporosis with over 20,000 people recruited, and will attempt to incorporate these data into modeling in the future. Third, side effects of zoledronate, such as osteonecrosis of the jaw and atypical femoral fracture, were not considered in the model. The reason is that the health state utility of these side effects was unclear for Chinese population, and these effects occur rarely at the dose of treating osteoporosis. Fourth, one-way sensitivity analysis was performed in the study, which may not be as robust as probabilistic sensitivity analysis. However, we found ICERs tested on all variables in the current sensitivity analysis fell below $23,000, suggesting the results were not sensitive to the uncertainty of those input data. Fifth, we used visual analogue scale (VAS)-based QALY by fracture types as the input data of the model. Though a trade time-off (TTO)-based QALY model was available for Chinese population [48], we were unable to use it to calculate QALY by fracture types, as we did not collect the information in 5 dimensions (i.e., mobility, self-care, usual activities, pain/discomfort, and anxiety/depression) and 3 levels (i.e., 1 = no problems, 2 = some/moderate problems, and 3 = extreme problems) for each individual in PK-VF study. As an alternative, we used the VAS-based QALY data by fracture types in the model [25], which has also been used in other previous cost-effectiveness analysis [43, 45].

Conclusions

This study is the first cost-effectiveness analysis simulating osteoporosis and osteoporotic fractures in Chinese postmenopausal women based on real-world data. The study indicated that widespread use of zoledronate was of clinical and economic benefit to Chinese postmenopausal women with a 10-year major fracture probability by FRAX over 7%.

Electronic supplementary material

(DOCX 46 kb)

(PDF 1308 kb)