Skip to main content

Expected lifetime numbers and costs of fractures in postmenopausal women with and without osteoporosis in Germany: a discrete event simulation model

Abstract

Background

Osteoporotic fractures cause a large health burden and substantial costs. This study estimated the expected fracture numbers and costs for the remaining lifetime of postmenopausal women in Germany.

Methods

A discrete event simulation (DES) model which tracks changes in fracture risk due to osteoporosis, a previous fracture or institutionalization in a nursing home was developed. Expected lifetime fracture numbers and costs per capita were estimated for postmenopausal women (aged 50 and older) at average osteoporosis risk (AOR) and for those never suffering from osteoporosis. Direct and indirect costs were modeled. Deterministic univariate and probabilistic sensitivity analyses were conducted.

Results

The expected fracture numbers over the remaining lifetime of a 50 year old woman with AOR for each fracture type (% attributable to osteoporosis) were: hip 0.282 (57.9%), wrist 0.229 (18.2%), clinical vertebral 0.206 (39.2%), humerus 0.147 (43.5%), pelvis 0.105 (47.5%), and other femur 0.033 (52.1%). Expected discounted fracture lifetime costs (excess cost attributable to osteoporosis) per 50 year old woman with AOR amounted to €4,479 (€1,995). Most costs were accrued in the hospital €1,743 (€751) and long-term care sectors €1,210 (€620). Univariate sensitivity analysis resulted in percentage changes between -48.4% (if fracture rates decreased by 2% per year) and +83.5% (if fracture rates increased by 2% per year) compared to base case excess costs. Costs for women with osteoporosis were about 3.3 times of those never getting osteoporosis (€7,463 vs. €2,247), and were markedly increased for women with a previous fracture.

Conclusion

The results of this study indicate that osteoporosis causes a substantial share of fracture costs in postmenopausal women, which strongly increase with age and previous fractures.

Peer Review reports

Background

With a lifetime risk of 40% to 50%, an osteoporotic fracture is one of the most likely negative health events in the remaining lifetime of a 50 year old woman [1]. The individual probability for a 50 year old woman of experiencing an osteoporotic fracture is higher than the probability of developing an atherosclerotic cardiovascular disease (4.1% to 30.7% [2] (depending on risk profile)) or breast cancer (around 11% [3]).

Studies have shown that a decrease in bone mineral density (BMD) is a strong predictor for fractures on different sites of the human body [46]. According to the World Health Organization (WHO) a person suffers from osteoporosis when the person’s BMD is at least 2.5 standard deviations lower (T≤-2.5) than the average BMD of healthy adults [7]. Clinically, osteoporosis is characterized ”by low bone mass and microarchitectural deterioration of bone tissue, leading to enhanced bone fragility and consequent increase in fracture risk” [8]. While women are more often affected than men, a strong increase in osteoporosis prevalence rates beyond the age of 50 occurs in both genders [9, 10].

Osteoporotic fractures have serious consequences on the individual and the societal level. Fracture patients often have a decreased functional mobility [11], health related quality of life [12] and are faced with an increased mortality [13]. Especially for advanced elderly persons, fracture related functional impairment often necessitates external help from professional home care services, relatives or friends [11]. In some cases institutionalization in a nursing home (NH) is unavoidable [14]. On the societal level fractures have repeatedly been shown to cause a high economic burden [1521]. As a consequence of the (expected) demographic change in Germany [22], the economic burden due to osteoporotic fractures is likely to increase strongly in the next decades [17].

Internationally, the economic burden of osteoporosis and osteoporotic fractures has frequently been investigated in the last decade [16, 1821, 23]. In contrast to the large number of cost studies worldwide, the evidence for Germany is not abundant. Four of the available cost studies on osteoporotic fractures were prevalence based, estimating fracture costs for a specific year [2427]. One study used an incidence based approach to estimate lifetime costs caused by hip fractures [28]. Another study used a Markov-cohort model to estimate the cumulative (2010–2050) and yearly (2010, 2030 and 2050) osteoporosis-attributable costs of different fracture types [17]. None of these studies estimated expected numbers and costs for all relevant fractures occurring over the remaining lifetime. Moreover, expected costs of specific subgroups, e.g. persons already suffering from osteoporosis, having a prevalent fracture or living in a NH, were not analyzed.

Our study had different aims. First, a discrete event simulation model (DES) was developed to estimate the expected lifetime fracture numbers and costs of six osteoporotic fracture types for 50 year old women in Germany. Second, the same estimates for hypothetical women who never suffer from osteoporosis over their remaining lifetime were made. Third, lifetime fracture costs attributable to the risk factor osteoporosis (excess cost) were calculated which indicate the economic potential of preventing osteoporosis as the main risk factor for fractures. Finally, lifetime fracture costs were estimated for women with different risk profiles (with previous fractures, with osteoporosis, or living in a NH) and at different starting ages. In summary, our model should give a deeper insight into the economics of six osteoporotic fractures in Germany, and should provide new and valuable information to the international literature.

Methods

Modeling approach

A DES model [2931] was developed to estimate the expected lifetime fracture numbers and costs of postmenopausal women (aged 50 years and older) in the German general population (at average risk getting osteoporosis) and in those who never suffer from osteoporosis. The simulation starts in the year 2009 and simulates 50 year old women until the age of 100 years or death (lifetime horizon). The included costs are considered from a societal perspective.

A DES is an individual based simulation technique where all individuals are simulated one by one [30]. The main components in DES are entities, attributes, events, time, resources (optional) and queues (optional) [29]. In a disease simulation context entities are individuals (e.g. patients). Individuals have attributes which reflect different personal characteristics, for example age, gender, health status, risk factors or event history. Depending on the number and combinations of attributes the model determines the individual’s probability to experience an event. The individual’s attributes profile is dynamic and can change at each time point in the model. Changes can be triggered by time alone (e.g. due to natural aging) or by events (e.g. occurrence of osteoporosis). Events are dichotomous (discrete) and can occur simultaneously at each time point (interval) in the model [29]. Event initialization is done by comparing a uniform (pseudo) random number between 0 and 1 with an event probability or a value of a hazard function. An event occurs when the random number is smaller (or equal) than the event probability value [30]. Time can be handled in two different ways in DES: In an “event driven” simulation time jumps from event to event, whereas in a “time driven” simulation time progresses by constant intervals (e.g. years, months, days) [31]. Optionally, DES can also handle resource constraints using queuing-systems and allow for interaction between individuals [29].

In order to estimate the expected lifetime fracture numbers and costs of a 50 year old woman, we conceptualized a “time driven” DES with time intervals of one year, no resource constraints and no interaction between individuals. The patient flow diagram in Figure 1 shows the mechanism of the DES-model: “Circles” represent start- and end-point of each simulation run (simulation of one individual n); in “rectangles” events are initialized whereas “rhombuses” represent decision points which determine the individual’s way through the model based on the individual’s attributes. The main events in this DES are the occurrences of hip, other femur, clinical vertebral, humerus, pelvis and wrist fractures. The probability of a fracture event is influenced by four risk factors (attributes) known from the literature: Age [32, 33], osteoporosis (BMD) [46], prevalent fractures (prev. fx) [34], and living in a NH [35, 36]. The occurrence of a risk factor itself is modeled as an event, e.g. institutionalization in a NH, and can arise in each time interval. At the beginning of each simulation run for an individual n the model determines if the individual directly starts with osteoporosis or a prevalent fracture (preload) [29], based on the age specific prevalence rates of these risk factors. All women start with the age of 50 in the base case.

Figure 1
figure 1

Model process and structure (Flow Chart). Legend of Figure 1: prev. fx = previous fracture, NH = nursing home, av. Inst. = average NH institutionalization probability, acute fx = acute fracture, fx status = acute or previous fracture, t = time interval, n = individual.

There is evidence that fractures increase the probability of institutionalization in a NH [14]. Therefore, the probability of the event “admission to a NH” (NH entry) is dependent on acute fracture events (acute fx) and the age-dependent average institutionalization probability (av. Inst.). Fractures increase mortality in the year of the fracture and in the subsequent years [13]; furthermore it is known that individuals living in a NH have higher mortality rates compared to those who do not at the same age [37]. Therefore, whether an individual survives the actual time interval (“Survival”) depends on age, the presence of fractures (fx status) in the actual and previous time interval and the living situation (living in NH or not). Each individual is simulated until one of the two endpoint conditions, age of 100 years or death, is fulfilled. The simulation was programmed in Matlab R2012a (The MathWorks, Inc., Natick, MA, USA) in combination with Excel 2010 (Microsoft Cooperation, Redmond, WA, USA).

Cost tracking

Not all fracture patients are treated in an inpatient setting (hospital). Therefore, the model determines on the basis of fracture specific hospitalization probabilities, whether a fracture patient is treated in a hospital or exclusively treated in an outpatient setting. If a patient is treated in a hospital, the model tracks acute hospital and post-hospital outpatient treatment costs; otherwise only outpatient costs are tracked. Also, not all patients receive inpatient rehabilitation after a hospital stay. For this reason, the same cost tracking approach, based on fracture specific post-hospital rehabilitation probabilities, was applied for rehabilitation costs. For tracking fracture-attributable long term care costs, a parallel background simulation approach was implemented [20]. The approach simulates the probability of institutionalization in a NH for any reason and the fracture specific probability in parallel. If a woman is institutionalized due to a fracture, long term care costs (i.e. yearly NH costs weighted by level of care as well as capital costs) are tracked as fracture-attributable only as long as no institutionalization occurs for any other reason [20]. A half-cycle correction was applied to long term care costs [38].

All other cost categories in the model are tracked as unit costs (rewards) per fracture event, similar to the common cohort model approach [38]. A detailed description of the epidemiological and cost specific input data used can be found in the electronic supplementary material (Additional file 1).

Determination of life time costs and fracture events

In total, two risk groups of 200,000 hypothetical women (simulation runs) each were simulated through the model. For the first risk group it was assumed that these women were at average osteoporosis risk over their remaining lifetime (based on osteoporosis incidence and prevalence), whereas the second risk group was assumed to be free of osteoporosis over their remaining lifetime. To determine osteoporosis-attributable fracture costs (excess costs), the difference in costs between the two considered risk groups was calculated [17]. To eliminate the stochastic noise between the two groups (with regard to simulated costs and fracture numbers), that is the randomly generated difference between the two groups because of using different random streams (sequence of random numbers) for each group, the same random stream (common random numbers) was used for both simulated groups [39].

Epidemiological input data

Fracture probabilities – general population

Age-dependent fracture probabilities for all considered fracture types were calculated based on official German population [22] and hospital discharge data [40] of the year 2009. The number of hospitalized fracture cases for women between 50 and 100 years of age (five-year age classes) were identified using the international classification of disease (ICD-10) with the following codes: hip S72.0-2, other femur S72.3-9, clinical vertebral S12.0-2/S12.7/S22.0-1/S32.0, humerus S42.2-4, pelvis S32.1-8 and wrist S52.5-6. Fracture probabilities based on hospital cases were estimated by dividing the identified fracture cases from the discharge statistic with the number of persons in the general population in the corresponding five-year age class (population at risk). While 100% of hip and other femoral fractures are treated in an inpatient setting (hospital) [41], some of clinical vertebral, humerus, pelvis and wrist fractures are treated exclusively in an outpatient setting [4144]. Therefore, the estimated fracture probabilities calculated from hospital discharge statistics and population data underestimate the “total” fracture probabilities, at least for clinical vertebral, humerus, pelvis and wrist fractures. To include these outpatient fracture cases in the total fracture probability, the fracture probabilities estimated based on hospital cases were divided by fracture specific and age-independent hospitalization probabilities from the literature (clinical vertebral [41, 44, 45], humerus [42], pelvis [43], wrist [42]) (see Additional file 1: A.3.a-b).

The calculated “total” fracture probabilities pertain to the German female general population (see Table 1). However, as women were assumed to be in specific living situations (women with or without osteoporosis, women with or without a previous fracture, community-dwelling women or women living in a NH), these general population fracture probabilities had to be adjusted.

Table 1 Total yearly fracture probability of the female general population by age and fracture type

Fracture probabilities for women with and without a previous fracture

Clinical studies showed that a previous fracture is an important risk factor for a subsequent fracture, even after adjustment for BMD [34]. The age-adjusted risk ratios for further hip (RRHipPreviousFx vs. NoPreviousFx) and other osteoporotic fractures (RROsteoPreviousFx vs. NoPreviousFx), comparing women with a previous fracture to those without a previous fracture, were taken from a meta-analysis [34]. However, the risk ratios reported in this meta-analysis are only applicable to fracture probabilities of women without a previous fracture. Therefore, in a first step, age-dependent relative risks for women without a previous fracture compared to the general population (RRNoPreviousFx vs. GpFx) were estimated. The calculation method was adopted from Schousboe et al. [46] using the following formula:

RR NoPreviousFx vs . GpFx = 1 / 1 + RR AnyPreviousFx vs NoPreviousFx 1 * Prev PreviousFx

The age-adjusted risk ratios for any fracture comparing women with previous fractures to those without a previous fracture adjusted for BMD (RRAnyPreviousFx vs. NoPreviousFx) and age-dependent prevalence rates for previous fractures (PrevPreviousFx) were both taken from the previously mentioned meta-analysis [34]. In the second step the estimated RRNoPreviousFx vs. GpFx were multiplied with the fracture probabilities of the general population to obtain fracture probabilities for women without a previous fracture. If a woman suffered from a fracture in the model, further fracture probabilities were calculated by combining fracture probabilities for women without a previous fracture with RRHipPreviousFx vs. NoPreviousFx and RROsteoPreviousFx vs. NoPreviousFx depending on the fracture type (see Additional file 1: A.3.c).

Fracture probabilities for women with and without osteoporosis

Fracture probabilities for women with osteoporosis were derived by multiplying the age-dependent female general fracture probabilities with age and fracture specific relative risks for women with osteoporosis compared to the female general population (RROST). The RROST were estimated based on a method described by Kanis et al. [10]. Necessary input data such as fracture specific relative risks by a decrease of one standard deviation in BMD measured at the femoral neck (RRfx) [46], the osteoporosis BMD threshold [10], and age-dependent population BMD values (reference values from NHANES III) [47] were taken from different international studies. To reflect lower fracture probabilities for women without osteoporosis compared to the female general population, probabilities were estimated based on a method described by Bleibler et al. [17] (see Additional file 1: A.3.d-e).

Fracture probabilities for women living and not living in a NH

There is evidence that women living in a NH have a higher fracture risk compared to community-dwelling women [35, 36]. To include this in our model, relative fracture risks for women residing in a NH and for those who do not, each compared to the female general population, were calculated. Basis for the calculation was a dataset (claims data) from a large German mandatory sickness fund (Allgemeine Ortskrankenkasse Bayern (AOK Bavaria)). This dataset contains individual data about level of long-term care (e.g. begin of care, admission to NH), hospital cases by ICD-10 and date of death for persons aged 65 years and older insured between January 2004 and June 2009. Relative fracture risks were determined on pooled cumulative fracture incidence rates from 2004–2008. Incident hospital fracture cases by fracture type were identified using the ICD-10 codes mentioned above. A detailed description of the dataset can be found in the electronic supplementary material (see Additional file 1: A.3.f-g).

Probability of developing osteoporosis

The age-dependent probability of developing osteoporosis was calculated based on osteoporosis prevalence rates originally estimated from the NHANES III study [9, 48]. The osteoporosis prevalence rate was transformed into a applicable transition probability (incidence) using a method described by Podgor and Leske [49] (see Additional file 1: A.2).

Probability of institutionalization in a NH

Especially hip fracture patients in higher age groups often have persistent reductions in their functional abilities and mobility after a fracture event [11]. Moreover, they frequently depend on external care or need to be institutionalized as a result of a fracture [11]. Similar to hip fractures, other types of osteoporotic fractures have been shown to increase the probability of institutionalization [14]. To model the association of a fracture event and institutionalization, crude age and fracture specific probabilities of admission to a NH (within 3 months after hospitalization) were estimated based on claims data from the AOK Bavaria. Besides a fracture related transition into a NH, the model allows institutionalization due to other reasons. To include this, official long-term care prevalence rates from 2009 [50] were transformed into NH institutionalization probabilities [49]. A transition between the status community-dwelling and status NH was possible for women 65 years and older. Similar to other simulation studies, it was assumed that persons in a NH will remain there for the rest of their lives [51] (see Additional file 1: A.4.a-b).

Mortality

As main source of mortality, an official generation life table (all-cause mortality) of women born in 1959 was used to reflect a realistic remaining life expectancy for a 50 year old woman living in the year 2009 (start year of the model) [52]. Depending on the actual model state (community dwelling, NH or fracture), age specific all-cause mortality was multiplied with state specific relative risks of mortality. Relative risks of mortality due to fractures were taken form a Canadian study [13]. The authors of this study determined the increase in relative mortality risk due to hip, clinical vertebral, proximal humerus, wrist and other fractures for five age classes and a follow up period of 10 years separated in 3 categories (1; 2–5; 6–10 years). They adjusted the calculated relative risks for co-morbidities and location of residence. Relative mortality risk for persons living in a NH or in a community-dwelling setting compared to the female general population were calculated based on the AOK Bavaria claims data (see Additional file 1: A.1.a-d).

Cost input data

Fracture related unit costs were estimated from a societal perspective for Germany. The base year of cost evaluation was 2009 and a discount rate of 3% was applied [53]. Included fracture related direct cost were separated in three main categories: inpatient costs include costs due to acute hospital care, rehabilitation after hospitalization and long term care, outpatient costs include costs due to physician and physiotherapist visits, analgesics (medication) as well as home care, and family costs include informal care costs. An overview over the main direct unit costs is shown in Table 2. Productivity costs were included in the model and estimated based on the human capital [54] and friction cost approach [55].

Table 2 Overview of direct unit costs in € by cost category and fracture type

Inpatient costs

In Germany, hospitals are reimbursed according to a dual system: Operating costs due to direct resource uses are paid on the basis of a German version of diagnosis related groups (G-DRG) by health insurance funds (private or mandatory), whereas capital costs are compensated by federal states [53]. Therefore, applied hospital costs consist of fracture related DRGs and capital costs. To determine the DRG part of total hospital costs per fracture type, the G-DRG Browser V2010 [56] was used. This dataset includes information on a large representative sample of German hospitals in 2009, including number of DRG-cases, age, gender, disease diagnosis (ICD-10), DRG relative cost weights and mean length of stay. On the basis of this information an average relative cost weight and mean length of stay was calculated for each considered fracture type and multiplied with a state weighted base rate of 2009 [57] to obtain DRG costs. Capital costs were calculated by combining the fracture specific length of stay with a daily capital cost rate [58] (inflated to 2009 [59]) (see Additional file 1: B.1.a). Rehabilitation costs after a hospital stay were estimated by combining the fracture type related duration of an inpatient rehabilitation treatment [60] with a daily cost rate [61] (see Additional file 1: B.1.b). Long term care costs are based on the official care statistic 2009 [50]. This statistic provides information on daily care costs and number of persons living in long term care by level of care in 2009. Therefore, level of care weighted yearly inpatient long term care costs were calculated. Additionally, capital costs for long term care [62] were added to the yearly unit costs (see Additional file 1: B.1.c).

Outpatient costs

Fracture related costs for physician visits (including outpatient surgeons), physiotherapeutic treatments and analgesics for patients exclusively treated in the outpatient sector, as well as outpatient cost for post-hospital treatment were estimated using German unit costs [58], inflated to 2009 [59]. Necessary resource use data in the three outpatient categories were taken from a German cost-effectiveness study [63]. The authors of this study collected information on resource use associated with hip, vertebral and wrist fractures. It was assumed that outpatient costs for other femur as well as pelvis fractures are similar to hip fractures and outpatient costs for humerus fractures are similar to wrist fractures (see Additional file 1: B.1.d).

Costs for home care per fracture type were determined based on information about fracture type specific hours of home care needed [19] and hourly unit costs [64]. As no information on fracture type specific hours of home care was available for Germany, a study from Austria was used [19]. It was assumed that only women older than 65 years utilize home care. No home care costs were applied to women living in a NH (see Additional file 1: B.1.e).

Informal care costs

To reflect the fracture related costs of informal care [65], fracture type specific amount of hours spend by relatives of a fracture patient were evaluated with the market cost approach (proxy good) [65]. Each hour of care was monetarily valued by the hourly gross salary of an employee in the field of care for elderly and disabled persons [66], corrected by the employer share of social contribution [54, 67]. As no information about number of hours spend by relatives of the fracture patients was available for Germany, Austrian data [19] were used. Similar to home care costs, it was assumed that women aged 65 or younger and women living in NH receive no informal care (see Additional file 1: B.1.f).

Productivity costs

The model considers two categories of productivity costs. First, loss in productivity due to the fracture related inability to perform paid work, and second, the productivity loss due to osteoporosis related premature deaths [68]. Productivity loss was evaluated using the human capital [54] and the friction cost approach [55]. Necessary information, such as employment rate [69], yearly average gross earnings [66] corrected by employer share of social contribution [67], time away from work (due to disability and rehabilitation [60]), length of friction period [70] for the year 2009, were taken from official statistics. A yearly net wage increase of 2% was assumed [54, 71]. If necessary, a half-cycle correction was applied to determine costs [38] (see Additional file 1: B.2.a-b).

Model assumptions

In order to improve model transparency, all model assumptions are presented in Table 3, distinguishing between assumptions regarding event probabilities and costs. Model assumptions were evaluated subjectively with respect to their expected impact on modeling results and scored from 1 asterisk (*, low impact) to 3 asterisks (***, high impact).

Table 3 Overview of the model assumptions

Scenario analysis

A number of scenarios of direct lifetime fracture costs were analyzed for women with different start characteristics. Each considered scenario consists of a combination of four of the following (model) start characteristics: Age (50 or 75), residential status (community dwelling (C) or NH (N)), disease status (already suffering from osteoporosis (O), at average risk of getting osteoporosis (avO) or not at risk of getting osteoporosis (nO)) and fracture history (with a previous fracture (P) or with no previous fracture (nP)). In total we modeled 18 scenarios, and each scenario was labeled based on the combination of the four start characteristics e.g. a 50 year old woman, community dwelling with average risk of getting osteoporosis and no previous fracture was labeled as 50_C_avO_nP in the results section. All scenarios were simulated separately, direct costs are presented undiscounted and discounted by 3%.

Sensitivity analysis

Deterministic one way and probabilistic sensitivity analyses were conducted. We performed 20 deterministic sensitivity analyses (S1-S20) to evaluate the impact on excess costs due to osteoporosis-attributable fractures:

S1: We valued informal care based on the opportunity cost approach [65] with an average German 2009 hourly wage rate (corrected by employer share) [66, 67]. S2: We valued informal care with a minimum hourly wage rate for employees in the area of health care and nursing introduced in 2010 in Germany [72]. S3: We assumed no excess mortality due to fractures. S4: We assumed baseline excess mortality only in the first year but not thereafter. S5: We assumed 30% of baseline excess mortality only in the first year, but none thereafter. S6: We did not consider a previous fracture as a risk factor for further fractures. S7: We applied the lower relative risk value of the 95% confidence interval to model the risk factor “previous fracture”. S8: We applied the upper relative risk value of the 95% confidence interval to model the risk factor “previous fracture”. S9: We applied the lower value for the fracture specific relative risks by a decrease of one standard deviation in BMD measured at the femoral neck. S10: We applied the upper value for the fracture specific relative risks by a decrease of one standard deviation in BMD measured at the femoral neck. S11: We increased the osteoporosis prevalence rate by +30%. S12: We decreased the osteoporosis prevalence rate by -30%. S13: We increased all input fracture probabilities by +30%. S14: We decreased all input fracture probabilities by -30%. S15: We assumed that only hip fracture patients have an increased risk to be institutionalized in a NH. S16: We assumed a yearly increase of modeled input fracture incidence rates by +1%. S17: We assumed a yearly decrease of modeled input fracture incidence rates by -1%. S18: We assumed a yearly increase of modeled input fracture incidence rates by +2%. S19: We assumed a yearly decrease of modeled input fracture incidence rates by -2%. S20: We applied a discount rate of 5%.

In order to perform a probabilistic sensitivity analysis, we simulated 10,000 hypothetical women for each risk group (simulation runs (n)) 1,500 times (simulations (m)) by varying parameters simultaneously for each of the 1,500 simulations (second-order uncertainty) [73]. Common distributional assumptions were used for costs, probabilities and relative risks [73] (see Additional file 1: C). The results of the probabilistic sensitivity analyses are reported as uncertainty intervals, estimated based on the percentile method (2.5 and 97.5 percentile) [73].

Results

Base case

Figure 2 shows the expected lifetime numbers of the six analyzed fracture types of a 50 year old woman for the two risk groups. Considering a 50 year old woman with an average risk of getting osteoporosis, there are 0.282 hip fractures expected for her remaining life, whereas 0.119 hip fractures are expected for a woman never suffering from osteoporosis. Thus, around 57.9% of all expected hip fractures are attributable to osteoporosis. Wrist fractures are the second most frequent fracture type, with expected 0.229 fractures for women at average risk of osteoporosis and 0.188 for those not at risk. With around 18.2% of wrist fractures attributable to osteoporosis, this fracture type had the lowest osteoporosis risk attribution. The lifetime numbers of other femur, clinical verterbral, humerus and pelvis fractures for a woman with average osteoporosis risk are expected at 0.033, 0.206, 0.147 and 0.105 respectively, which corresponds to an osteoporosis risk attribution of around 52.1%, 39.2%, 43.5% and 47.5% for these fracture types.

Figure 2
figure 2

Expected lifetime numbers of fractures (95%UI) of a 50 year old woman by fracture type and risk class.

Table 4 displays the expected direct costs for each of the six fracture types accrued in all considered healthcare sectors, excluding costs due to long term care. In women with average osteoporosis risk, hip fractures cause €1,277 (39.1%) of all discounted costs. Considering discounted excess cost due to osteoporosis, hip fractures cause €692 (50.4%) of total excess costs. The fracture types with the second and third largest contribution to discounted excess costs are clinical vertebral (15.9%) and humerus (14.9%) fractures. The lowest share of discounted excess costs was found for wrist fractures with 4.4%.

Table 4 Undiscounted and discounted fracture lifetime costs (€) of a 50 year old woman, by fracture type for two risk classes of osteoporosis and excess

Table 5 displays the expected direct lifetime fracture costs of a 50 year old woman by healthcare sector. Considering all six fracture types, discounted total costs of a 50 year old woman at average risk getting osteoporosis amounted to €4,479, of which around €1,995 (44.5%) were excess costs. With regard to the different healthcare sectors, we found that around 71.3% of the discounted total excess costs resulted from inpatient treatment, mainly caused by acute hospital treatments (37.6%) and fracture-attributable long term care (31.1%). Treatments in the outpatient sector were responsible for 17.2% and informal care for 11.5% of these costs.Figure 3 depicts the expected annual undiscounted total direct fracture costs per capita for postmenopausal women by age. It shows that women with average osteoporosis risk aged 50–70 years have relatively low annual costs per capita of €14 to €142, with low osteoporosis cost attributions of 9.7% to 25.5%. However, a strong increase in these costs and corresponding cost attributions can be observed beyond the age of 75. The annual costs per capita increase from €260 for women aged 75 years to around €1,535 for those aged 95, with osteoporosis cost attributions of 31.3% to 68.0%. For women aged 82 and older, excess costs exceed costs for women never suffering from osteoporosis.

Table 5 Undiscounted and discounted fracture lifetime costs (€) of a 50 year old woman, by healthcare sector for two risk classes of osteoporosis and excess
Figure 3
figure 3

Annual fracture costs per capita by age for two risk classes and excess.

Considering discounted productivity costs due to the inability to perform paid work, we found only low excess lifetime costs per capita of around €63 applying the human capital approach and €22 applying the friction cost approach.

Scenario analysis

The scenario analysis in Figure 4 shows that a 50 year old community-dwelling woman with osteoporosis and no previous fracture (50_C_O_nP) causes approximately 3.3 fold (€7,463 vs. €2,247) the discounted direct fractures lifetime costs of a 50 year old community-dwelling woman not at risk getting osteoporosis without a previous fracture (50_C_nO_nP). If a 50 year old community-dwelling woman with osteoporosis has additionally experienced a previous fracture (50_C_O_P), the fracture costs are 4.2 times (€9,452 vs. €2,247) higher. If a 50 year old community-dwelling woman never suffering from osteoporosis had a previous fracture (50_C_nO_P), discounted direct lifetime costs are increased by 46% compared to the same woman with no previous fracture (50_C_nO_nP). In general, we found that the expected discounted direct lifetime costs of a woman living in a NH are markedly lower than those of a woman living in a community dwelling setting. This is mainly explained by the high mortality in women institutionalized in a NH. Moreover, no costs for facture attributable long term care, professional home care and informal care were applied to institutionalized women.

Figure 4
figure 4

Undiscounted and discounted direct fracture lifetime costs (€) of different scenarios (women with different start characteristics).

Sensitivity analysis

The results of the deterministic sensitivity analyses are presented in Figure 5, whereas the results of the probabilistic sensitivity analysis (PSA) can be found in Tables 4, 5 and Figure 2 reported as 95% uncertainty intervals. Figure 5 shows the percentage difference (-%/+%) between each univariate sensitivity analysis and the discounted total direct fracture excess costs from the base case analysis. The largest percentage difference between costs was found when considering a yearly increase and decrease of fracture rates by -/+2% (S18/S19) (-48.4%/+83.5%), varying osteoporosis prevalence by -/+30% (S11/S12) (-44.1%/+64.4%), and by increasing the discount factor to 5% (S20) (-47.9%). The sensitivity analysis shows that neglecting the risk factor “previous fractures” (S6) reduces costs by 26.6%. Assuming that only hip fractures induce a fracture related institutionalization (S15) decreased costs by -13.9%. Applying minimal wage rates to value informal care (S2) led to a cost decrease of -4.8%, whereas applying the opportunity cost approach (S1) increased costs by +4.6%. Moreover, assuming no fracture excess mortality (S3) increased costs by 15.4%, while ignoring fracture related mortality effects beyond the year of fracture occurrence (S4) lead to a cost increase of 1.9%.

Figure 5
figure 5

Deterministic sensitivity analyses, difference (%) between direct fracture excess costs due to osteoporosis from base case excess costs. Legend of Figure 5: CL = confidence limit; RR = relative risk; SD = standard deviation; BMD = bone mineral density.

Validation

For an internal validation of our model we compared the modeled fracture incidence rates with the expected fractures rates (input data) of the female general population. We also reported the distribution of the number of fractures per woman to support internal validity and transparency of the model estimates. Table 6 shows the expected and modeled hip fracture rates as well as the proportion of women with 0 to 6 hip fractures (the same data for all other fracture types is available in the electronic supplementary material (see Additional file 1: D.1.a-b)).

Table 6 Modeled and expected hip fracture rates and proportions of women with 0 to 6 hip fractures

The expected and modeled fracture incidence rates are very similar (Table 6), which suggests that the proportions of women with different fracture probabilities (e.g. with and without a previous fracture, with and without osteoporosis, living in NH or not) are modeled appropriately.

As an external model validation, the modeled hip fracture incidence rates were compared to epidemiological studies from Sweden, USA and UK (Figure 6). Additionally, the fracture lifetime risk of a 50 year old woman at average risk of getting osteoporosis was compared to other modeling and epidemiological studies. Table 7 shows the (first ever) hip fracture lifetime risk of our modeling study in comparison with other epidemiological and modeling studies from different countries. Our estimated (first ever) fracture lifetime risk of 19.8% for hip fracture is in the range of other simulated lifetime hip fracture risk estimates. Furthermore, in comparison with epidemiological studies our estimated hip fracture lifetime risk and fracture rates (Table 7 and Figure 6) show a good fit. From a systematic review, comparing worldwide hip fracture risks, it is known that the hip fracture risk level in Germany should be between that of Sweden and USA [74], which is confirmed by our estimates. (The same data for all other fracture types is available in the electronic supplementary material (see Additional file 1: D.1.a-b)).

Figure 6
figure 6

Comparison of modeled hip fracture rates with hip fracture rates from Sweden, UK, and USA.

Table 7 Hip fracture lifetime risk from modeling and epidemiological studies

Discussion

Summary and comparison to the literature

The goal of this study was to estimate the expected lifetime fracture numbers and costs per capita for two risk groups of German postmenopausal women, that is, women with average risk getting osteoporosis and women never suffering from osteoporosis. As key results we found that discounted lifetime fracture costs of €4,479 can be expected for a 50 year old woman with average osteoporosis and previous fracture risk, whereof €1,995 or 44.5% are attributable to the risk factor osteoporosis (excess costs). For both risk groups, around 70% of costs were due to inpatient treatment, like acute hospital treatments and long term care. The fracture type with the largest proportion of direct lifetime costs were hip fractures.

In a previous study [17], we estimated the impact of demographic change on osteoporosis-attributable fractures costs in Germany from 2010–2050. Since aggregated costs were calculated, a direct comparison in terms of lifetime costs per capita was not possible. The cost distribution between fracture types and healthcare sectors are quite similar in both studies, though. In contrast, in the present study it was possible to implement the event “NH entry” directly in the model, which allowed a more precise estimate of fracture attributable long term care costs. Internationally, few studies estimated and reported the average lifetime costs per capita of osteoporotic fractures. A study from Switzerland [20] found discounted (3%) inpatient lifetime costs for a person aged 50 (both genders) of CHF5,400 in the year 2000. In 2009 US$ purchasing power parities (US$PPP) [79] this equals US$PPP4,007, whereas we found average lifetime costs of US$PPP5,599 (2009) for a 50 year old woman. The authors of the Swiss model [20] used a first-order lifetime Markov-model (individual based simulation) considering three fracture types (distal forearm, (all) vertebral and hip). They modeled all-cause mortality and excess fracture mortality for at least 5 years following hip and vertebral fractures. They considered three fracture related cost categories from a societal perspective: acute hospital care, inpatient rehabilitation and hip fracture attributable NH costs. The difference to our average fracture lifetime cost may be due to the difference in the number of modeled fracture types (3 vs. 6) and in included cost categories. We additionally included costs for fracture related informal care, professional home care and outpatient care. Furthermore we assumed that all hospitalized fractures have an influence on NH institutionalization, whereas the authors of the Swiss model assumed only hip fractures to have an impact. Additionally we modeled NH mortality, which may have had a decreasing effect on fracture attributable NH costs. Further differences in average lifetime costs may be due to different life expectancy, fracture incidence and unit costs assumptions.

A modeling study from Hiligsmann et al. [80] conducted in Belgium found discounted (3%) direct lifetime fracture costs of €10,288 (US$PPP12,102 (2009)) for women aged 55 already suffering from osteoporosis, whereas we found €7,463 (US$PPP9,329 (2009)) for a woman aged 50 years. The study’s main purpose was to evaluate the cost-effectiveness of osteoporosis screening followed by anti-osteoporotic treatment. However, the authors also reported simulated lifetime fracture costs of a 55 year old woman with osteoporosis (no intervention group). The authors used a first-order lifetime Markov model (Markov-microsimulation). They included hip, clinical vertebral, forearm and other fractures in their model. Excess mortality for hip and clinical vertebral fractures was assumed for the first and subsequent year [51]. All direct fracture costs from a health-care payer perspective were included. For hip fractures first year and long term care costs (due to a NH admission) were considered. For other fracture types only first year costs were modeled [51]. The difference in fracture lifetime cost estimates between our model and the model of Hiligsmann et al. may have various reasons. We used a parallel simulation to examine all fracture attributable NH costs with “real world” NH mortality. Hilligsmann et al. pre-calculated NH costs due to hip fractures using average life expectancy assumptions [51]. Also it seems that Belgian unit costs for hip fracture in the first year (without NH costs) are higher than our unit costs. Furthermore, we modeled excess mortality for all fracture types (except wrist) and assumed fracture excess mortality beyond the second year after a fracture, whereas Hilligsmann et al. considered excess mortality only for hip and vertebral fractures in the first two years after fracture. The difference in the cost perspective between the two modeling studies (health payer perspective vs. societal perspective) may also explain the divergence in lifetime costs.

Strengths of the model

In comparison to other models in the field of osteoporotic fractures, our model has some strengths. Compared to a Markov cohort model our model allows calculating event probabilities dependent on the history of a modeled individual, which overcomes the classical Markov assumption of “memorylessness” [38]. Another model type which can handle patient history is the first-order Markov model (Monte Carlo simulation) where individuals are simulated separately through a classical Markov structure (health states) [38, 73]. In comparison to this model type our model can handle multiple events at the same time which is not possible in a classical first-order Markov model because of the mutual exclusivity of health states [38]. Furthermore, a first-order Markov model would require a very large number of health states to reflect the same patient pathways (6 fracture types, previous fracture (yes/no), osteoporosis (yes/no), NH (yes/no)) that we considered. Applying the event based structure overcomes this “flood of health states” [30]. From an epidemiological point of view our model has strengths in terms of the number of modeled fracture types. In many models only 3 fracture types (hip, vertebra and wrist) are considered [81], whereas we modeled 6 fracture types. In addition, our model allows separating the effect of osteoporosis, previous fractures, and NH institutionalization on costs. Finally, the model validation showed that our model reflects the real world epidemiological data very well and therefore has high validity.

In terms of cost modeling, a main strength of our model is the large number of cost categories included and the way how fracture attributable long term costs are modeled. Using the parallel simulation of NH entries under real world NH mortality allowed us to estimate fracture attributable long term care costs more precisely.

Limitations of the model

The results of model-based studies should be interpreted against the background of the model assumptions, e.g. model type, time interval, input data. In the model, we applied age-dependent fracture specific RRfx only to hip fractures, but not to non-hip fractures. If the age related patterns of RRfx for hip fractures are similar to other fracture types, the fracture risk for non-hip fractures in women with osteoporosis may be overestimated in older, and underestimated in younger age classes. In the model we applied osteoporosis prevalence rates and BMD data from a US population (NHANES III study). Applying these data was necessary, because no German study determined osteoporosis prevalence rates in a comparable design to the NHANES III study (in terms of osteoporosis definition based on BMD (T ≤ -2.5), sample size and study quality [41]). Furthermore, the prediction of osteoporosis risk attributions was exclusively based on BMD values measured at the femoral neck, but not at other sides, which may have led to an underestimation of osteoporosis risk attribution on other fracture sites. Fracture probabilities were mainly estimated based on the German hospital discharge statistic. In order to include fracture cases exclusively treated in an outpatient setting, fracture probabilities from the hospital discharge statistic were divided by fracture specific hospitalization probabilities. In our model these hospitalization probabilities were assumed to be age-independent. However, in a real world setting these probabilities may be age dependent and increasing with age, which would lead to an overestimation of total fracture probabilities in younger and an underestimation in older age classes. Additionally, we neglected vertebral fractures not coming to clinical attention. However, these fractures likely have a small influence on costs, but have been shown to increase the risk of subsequent fractures [82]. In terms of NH probabilities our model also underlies limitations. We assumed that only fractures treated in the hospital would increase the risk of an institutionalization in a NH. However, fractures treated in the outpatient setting may also increase the probability of an institutionalization into a NH. Hence, we may have underestimated fracture attributable long term care costs as a result. Also, we assumed that only institutionalizations within 3 months were fracture related, which could have similarly led to an underestimation of fracture attributable long term care costs.

In terms of cost data, there are further limitations. We did not consider costs due to transport, medical aids and early retirement. Cost for osteoporosis specific medication was also not considered. The main reason for neglecting medication costs was based on the aim of the study to determine fracture costs attributable to the risk factor osteoporosis, and not the cost of treating osteoporosis. A study which analyzed German insurance claims data found that the share of osteoporosis medication costs was low in comparison to other fracture attributable cost categories [83]. However, this may be due to the low osteoporosis specific treatment prevalence in Germany [83]. A further limitation is the assumption that women already residing in a NH at the time of a fracture event were assumed to not accrue fracture related long term care costs. In a real-world setting this assumption may not hold, as fracture events in a NH can go along with an increase in care needs and related costs [84]. Fracture specific long term care costs could not be presented, because all of these costs were technically attached to a single ”NH” event. Because of a lack of data we assumed age constant rehabilitation probabilities after a fracture. However, in a real world setting these probabilities may differ by age. Also we assumed that rehabilitation probabilities after a hospitalization are identical for institutionalized and community-dwelling women, which may have had an effect on rehabilitation costs. In reality institutionalized women may have lower rehabilitation probabilities. As data on outpatient resource use was only available for three fracture types (hip, vertebral and wrist), we assumed the same resource use for fractures in near bone regions (other femur, pelvis, humerus).

The estimation of expected lifetime fracture numbers and costs is highly uncertain in terms of changes in age-adjusted fracture incidence rates (secular trends). Studies from Europe showed an inconsistent picture of changes in age-adjusted fracture incidence rates, with increasing and decreasing trends reported in the last years. Potential reasons for an increase in age-adjusted hip fracture rates may be higher survival rates in very old and fragile persons, a decrease in physical activity and deficiency in vitamin D [85]. However, there are also potential reasons for a decrease in age-adjusted fracture rates, which may be due to improved medical treatment of osteoporosis (e.g. introduction of bisphosphonates), changes in habits related to fracture increasing risk factors (e.g. smoking or heavy drinking) or due to higher prevalence rates in obesity in the past years [85]. Considering the current situation of treating osteoporosis in Germany, it seems that there is room for improvements in terms of detecting and treating osteoporosis. Currently only around 18% of osteoporosis patients receive an osteoporosis specific treatment [83]. Furthermore, only 14% of fracture patients receive an osteodensitometry, which may indicate an under-detection of osteoporosis [83]. Moreover, the adherence to oral osteoporosis medications is suboptimal [86]. Improving detection and medical treatment of osteoporosis in the future could lead to a decreasing secular trend in fracture incidence rates in the future in Germany, which should also lead to a decrease in expected fracture lifetime numbers and costs. On the other hand, the predicted demographic change in Germany [22] may have an increasing effect on future age-adjusted incidence rates and therefore may increase lifetime fracture numbers and costs.

Implications for decision makers

The results of this study show that fractures cause relevant costs in postmenopausal women, which strongly increase with age, previous fractures and osteoporosis status. The high costs in age classes over 70, in combination with demographic change, could in addition intensify the economic burden due to fractures in Germany [17]. Thus, tackling fractures in a cost-saving or cost-effective way should be an important goal for German health care decision makers.

In the future, our model can be used to perform cost-effectiveness analyses of fracture prevention interventions in Germany. Beside the described model functions, we optionally implemented an osteodensitometry module. This feature makes the modeling of different screening to treat scenarios for Germany possible. Moreover, the real-world adherence of medical anti-osteoporotic therapy can be realistically implemented in the model. As prior studies found cost-effectiveness results to be highly sensitive to medication adherence assumptions, a realistic implementation of the respective empirical evidence seems particularly important in health economic modeling of osteoporotic fractures [87].

Conclusion

This is the first analysis modeling the epidemiology and costs of fractures for Germany. The analysis shows that 44.5% of all discounted direct fracture lifetime costs of a 50 year old woman are attributable to osteoporosis. For a 50 year old woman already suffering from osteoporosis, fracture lifetime costs are expected to be 3.3 fold the fracture costs of a healthy woman at the same age. The main cost drivers are hospital and long term care costs, causing around 70% of total fracture related lifetime costs. Beyond cost of illness, the model can be used to evaluate the cost-effectiveness of fracture prevention.

References

  1. Johnell O, Kanis J: Epidemiology of osteoporotic fractures. Osteoporos Int. 2005, 16 (Suppl 2): S3-S7.

    Article  PubMed  Google Scholar 

  2. Berry JD, Dyer A, Cai X, Garside DB, Ning H, Thomas A, Greenland P, Van Horn L, Tracy RP, Lloyd-Jones DM: Lifetime risks of cardiovascular disease. N Engl J Med. 2012, 366: 321-329.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Armstrong K, Eisen A, Weber B: Assessing the risk of breast cancer. N Engl J Med. 2000, 342: 564-571.

    Article  CAS  PubMed  Google Scholar 

  4. Marshall D, Johnell O, Wedel H: Meta-analysis of how well measures of bone mineral density predict occurrence of osteoporotic fractures. BMJ. 1996, 312: 1254-1259.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Stone KL, Seeley DG, Lui LY, Cauley JA, Ensrud K, Browner WS, Nevitt MC, Cummings SR, Osteoporotic Fractures Research G: BMD at multiple sites and risk of fracture of multiple types: long-term results from the Study of Osteoporotic Fractures. J Bone Miner Res. 2003, 18: 1947-1954.

    Article  PubMed  Google Scholar 

  6. Johnell O, Kanis JA, Oden A, Johansson H, De Laet C, Delmas P, Eisman JA, Fujiwara S, Kroger H, Mellstrom D, Meunier PJ, Melton LJ, O’Neil T, Pols H, Reeve J, Silman A, Tenenhouse A: Predictive value of BMD for hip and other fractures. J Bone Miner Res. 2005, 20: 1185-1194.

    Article  PubMed  Google Scholar 

  7. World Health Organization: Assessment of fracture risk and its application to screening for postmenopausal osteoporosis. Report of a WHO Study Group. World Health Organ Tech Rep Ser. 1994, 843: 1-129.

    Google Scholar 

  8. Anonymous: Consensus development conference: prophylaxis and treatment of osteoporosis. Am J Med. 1991, 90: 107-110.

    Article  Google Scholar 

  9. Wildner M: Osteoporose [Osteoporosis]. Dtsch Med Wochenschr. 2001, 126: A1170-A1172.

    Google Scholar 

  10. Kanis JA, Johnell O, Oden A, Jönsson B, De Laet C, Dawson A: Risk of hip fracture according to the World Health Organization criteria for osteopenia and osteoporosis. Bone. 2000, 27: 585-590.

    Article  CAS  PubMed  Google Scholar 

  11. Osnes EK, Lofthus CM, Meyer HE, Falch JA, Nordsletten L, Cappelen I, Kristiansen IS: Consequences of hip fracture on activities of daily life and residential needs. Osteoporos Int. 2004, 15: 567-574.

    Article  CAS  PubMed  Google Scholar 

  12. Hiligsmann M, Ethgen O, Richy F, Reginster JY: Utility values associated with osteoporotic fracture: a systematic review of the literature. Calcif Tissue Int. 2008, 82: 288-292.

    Article  CAS  PubMed  Google Scholar 

  13. Morin S, Lix LM, Azimaee M, Metge C, Caetano P, Leslie WD: Mortality rates after incident non-traumatic fractures in older men and women. Osteoporos Int. 2011, 22: 2439-2448.

    Article  CAS  PubMed  Google Scholar 

  14. Morin S, Lix LM, Azimaee M, Metge C, Majumdar SR, Leslie WD: Institutionalization following incident non-traumatic fractures in community-dwelling men and women. Osteoporos Int. 2012, 23: 2381-2386.

    Article  CAS  PubMed  Google Scholar 

  15. Burge R, Dawson-Hughes B, Solomon DH, Wong JB, King A, Tosteson A: Incidence and economic burden of osteoporosis-related fractures in the United States, 2005–2025. J Bone Miner Res. 2007, 22: 465-475.

    Article  PubMed  Google Scholar 

  16. Burge RT, Worley D, Johansen A, Bhattacharyya S, Bose U: The cost of osteoporotic fractures in the UK: Projections for 2000–2020. J Med Econ. 2001, 4: 51-62.

    Article  Google Scholar 

  17. Bleibler F, Konnopka A, Benzinger P, Rapp K, König HH: The health burden and costs of incident fractures attributable to osteoporosis from 2010 to 2050 in Germany–a demographic simulation model. Osteoporos Int. 2013, 24: 835-847.

    Article  CAS  PubMed  Google Scholar 

  18. Börgstrom F, Sobocki P, Ström O, Jönsson B: The societal burden of osteoporosis in Sweden. Bone. 2007, 40: 1602-1609.

    Article  PubMed  Google Scholar 

  19. Dimai HP, Redlich K, Schneider H, Siebert U, Viernstein H, Mahlich J: Direkte und indirekte Kosten von osteoporotisch bedingten Frakturen in Osterreich [Direct and indirect costs of fractures due to osteoporosis in Austria]. Gesundheitswesen. 2012, 74: e90-e98.

    PubMed  Google Scholar 

  20. Schwenkglenks M, Lippuner K, Häuselmann HJ, Szucs TD: A model of osteoporosis impact in Switzerland 2000–2020. Osteoporos Int. 2005, 16 (6): 659-671.

    Article  PubMed  Google Scholar 

  21. Tarride JE, Hopkins RB, Leslie WD, Morin S, Adachi JD, Papaioannou A, Bessette L, Brown JP, Goeree R: The burden of illness of osteoporosis in Canada. Osteoporos Int. 2011, 23: 2591-2600.

    Article  Google Scholar 

  22. Statistisches Bundesamt [Federal Statistical Office]: Bevölkerung Deutschland bis 2060: Ergebnisse der 12. Koordinierten Bevölkerungsvorausbrechnung [Population in Germany up to the Year 2060: Results of the 12th Coordinated Projection of Population]. 2009, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  23. Budhia S, Mikyas Y, Tang M, Badamgarav E: Osteoporotic fractures: a systematic review of U.S. healthcare costs and resource utilization. Pharmacoeconomics. 2012, 30: 147-170.

    Article  PubMed  Google Scholar 

  24. Brecht JG, Schädlich PK: Burden of illness imposed by osteoporosis in Germany. HEPAC Health Economics in Prevention and Care. 2000, 1: 26-32.

    Article  Google Scholar 

  25. Häussler B, Gothe H, Gol D, Glaeske G, Pientka L, Felsenberg D: Epidemiology, treatment and costs of osteoporosis in Germany–the BoneEVA Study. Osteoporos Int. 2007, 18: 77-84.

    Article  PubMed  Google Scholar 

  26. Konnopka A, Jerusel N, König HH: The health and economic consequences of osteopenia- and osteoporosis-attributable hip fractures in Germany: estimation for 2002 and projection until 2050. Osteoporos Int. 2009, 20: 1117-1129.

    Article  CAS  PubMed  Google Scholar 

  27. Bleibler F, Benzinger P, Lehnert T, Becker C, König HH: Cost of fractures in German hospitals - What role does osteoporosis play?. Gesundheitswesen. 2014, 76 (3): 163-168.

    CAS  PubMed  Google Scholar 

  28. Weyler EJ, Gandjour A: Sozioökonomische bedeutung von hüftfrakturen in deutschland [socioeconomic burden of Hip fractures in germany]. Gesundheitswesen. 2007, 69: 601-606.

    Article  PubMed  Google Scholar 

  29. Karnon J, Stahl J, Brennan A, Caro JJ, Mar J, Möller J: Modeling using discrete event simulation a report of the ISPOR-SMDM modeling good research practices task force–4. Med Decis Mak. 2012, 32: 701-711.

    Article  Google Scholar 

  30. Caro JJ, Moller J, Getsios D: Discrete event simulation: the preferred technique for health economic evaluations?. Value Health. 2010, 13: 1056-1060.

    Article  PubMed  Google Scholar 

  31. Sobolev B, Sanchez V, Kuramoto L: Health Care Evaluation Using Computer Simulation: Concepts, Methods, and Applications. 2012, New York, Heidelberg, Dordrecht, London: Springer

    Book  Google Scholar 

  32. Icks A, Haastert B, Wildner M, Becker C, Meyer G: Trend of hip fracture incidence in Germany 1995–2004: a population-based study. Osteoporos Int. 2008, 19: 1139-1145.

    Article  CAS  PubMed  Google Scholar 

  33. Kanis J, Johnell O, Oden A, Sernbo I, Redlund-Johnell I, Dawson A, De Laet C, Jonsson B: Long-term risk of osteoporotic fracture in Malmö. Osteoporos Int. 2000, 11: 669-674.

    Article  CAS  PubMed  Google Scholar 

  34. Kanis JA, Johnell O, De Laet C, Johansson H, Oden A, Delmas P, Eisman J, Fujiwara S, Garnero P, Kroger H, McCloskey EV, Mellstrom D, Melton LJ, Pols H, Reeve J, Silman A, Tenenhouse A: A meta-analysis of previous fracture and subsequent fracture risk. Bone. 2004, 35: 375-382.

    Article  CAS  PubMed  Google Scholar 

  35. Rapp K, Becker C, Cameron ID, Klenk J, Kleiner A, Bleibler F, König HH, Büchele G: Femoral fracture rates in people with and without disability. Age Ageing. 2012, 41: 653-658.

    Article  PubMed  Google Scholar 

  36. Benzinger P, Becker C, Kerse N, Bleibler F, Büchele G, Icks A, Rapp K: Pelvic fracture rates in community-living people with and without disability and in residents of nursing homes. J Am Med Dir Assoc. 2013, 14: 673-678.

    Article  PubMed  Google Scholar 

  37. Lubitz J, Cai L, Kramarow E, Lentzner H: Health, life expectancy, and health care spending among the elderly. N Engl J Med. 2003, 349: 1048-1055.

    Article  CAS  PubMed  Google Scholar 

  38. Sonnenberg FA, Beck JR: Markov models in medical decision making: a practical guide. Med Decis Making. 1993, 13: 322-338.

    Article  CAS  PubMed  Google Scholar 

  39. Stout NK, Goldie SJ: Keeping the noise down: common random numbers for disease simulation modeling. Health Care Manag Sci. 2008, 11: 399-406.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Statistisches Bundesamt [Federal Statistical Office]: Tiefgegliederte Diagnosedaten der Krankenhauspatentinnen und -Patienten 2009 [Deep Stratified Diagnosis-Data from Female and Male Inpatients 2009]. 2011, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  41. Kanis JA, on behalf of the World Health Organization Scientific Group: Assessment of Osteoporosis at the Primary Health-Care Level. 2007, Technical Report. World Health Organization Collaborating Centre for Metabolic Bone Diseases

    Google Scholar 

  42. Einsiedel T, Becker C, Stengel D, Schmelz A, Kramer M, Daxle M, Lechner F, Kinzl L, Gebhard F: Do injuries of the upper extremity in geriatric patients end up in helplessness? A prospective study for the outcome of distal radius and proximal humerus fractures in individuals over 65. Z Gerontol Geriatr. 2006, 39: 451-461.

    Article  CAS  PubMed  Google Scholar 

  43. Boufous S, Finch C, Close J, Day L, Lord S: Hospital admissions following presentations to emergency departments for a fracture in older people. Inj Prev. 2007, 13: 211-214.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Finnern HW, Sykes DP: The hospital cost of vertebral fractures in the EU: estimates using national datasets. Osteoporos Int. 2003, 14: 429-436.

    Article  PubMed  Google Scholar 

  45. Cooper C, O’neill T, Silman A: The epidemiology of vertebral fractures. Bone. 1993, 14: 89-97.

    Article  Google Scholar 

  46. Schousboe J, Gourlay M, Fink H, Taylor B, Orwoll E, Barrett-Connor E, Melton L, Cummings S, Ensrud K: Cost-effectiveness of bone densitometry among Caucasian women and men without a prior fracture according to age and body weight. Osteoporos Int. 2013, 24: 163-177.

    Article  CAS  PubMed  Google Scholar 

  47. Looker AC, Wahner HW, Dunn WL, Calvo MS, Harris TB, Heyse SP, Johnston CC, Lindsay R: Updated data on proximal femur bone mineral levels of US adults. Osteoporos Int. 1998, 8: 468-489.

    Article  CAS  PubMed  Google Scholar 

  48. Jerusel N: Analyse und Prognose der Krankheitslast Osteoporose-Attributabler Proximaler Femurfrakturen in Deutschland [Analysis and Prognosis of the Health Burden of Osteoporosis Attributable Proximal Femur Fractures in Germany]. 2009, PhD thesis. University Leipzig (Germany): Health Economics Research Unit

    Google Scholar 

  49. Podgor MJ, Leske MC: Estimating incidence from age-specific prevalence for irreversible diseases with differential mortality. Stat Med. 1986, 5: 573-578.

    Article  CAS  PubMed  Google Scholar 

  50. Statistisches Bundesamt [Federal Statistical Office]: Pflegestatistik 2009 Pflege im Rahmen der Pflegeversicherung Deutschlandergebnisse [Statistics on care 2009. Care in the Context of the Mandatory Care Insurance]. Results from Germany. 2011, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  51. Hiligsmann M, Ethgen O, Bruyere O, Richy F, Gathon HJ, Reginster JY: Development and validation of a Markov microsimulation model for the economic evaluation of treatments in osteoporosis. Value Health. 2009, 12: 687-696.

    Article  PubMed  Google Scholar 

  52. Statistisches Bundesamt [Federal Statistical Office]: Generationensterbetafeln für Deutschland: Modellrechnung für die Geburtenjahrgänge 1871–2004 [Generation Lifetables for Germany: Model Calculation for the Birth Cohorts 1871–2004]. 2006, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  53. Institut für Qualität und Wirtschaftlichkeit im Gesundheitswesen [Institute for Quality and Efficiency in Health Care]: Arbeitspapier Kostenbestimmung Version 1.0 [Working paper cost determination version 1.0]. 2009, Köln (Germany): Institute for Quality and Efficiency in Health Care

    Google Scholar 

  54. Rice DP, Cooper BS: The economic value of human life. Am J Public Health Nations Health. 1967, 57: 1954-1966.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Koopmanschap MA, Rutten FF, Van Ineveld BM, Van Roijen L: The friction cost method for measuring indirect costs of disease. J Health Econ. 1995, 14: 171-189.

    Article  CAS  PubMed  Google Scholar 

  56. Institut für das Entgeltsystem im Krankenhaus [Institute for the Hospital Renumeration System]: G-DRG V2010 Browser (2009 § 21 KHEntgG). 2011, Siegburg (Germany): Institute for the Hospital Renumeration System, http://www.g-drg.de/cms/Archiv/Systemjahr_2011_bzw._Datenjahr_2009#sm15.

    Google Scholar 

  57. Deutsche Krankenhaus Gesellschaft [German Hospital Association]: Landesbasisfallwerte der Bundesländer [Baserates of the German Federal States]. 2010, Berlin (Germany): German Hospital Association, http://www.dkgev.de/media/file/7827.LBFW_2005_2010_Stand_160610.pdf.

    Google Scholar 

  58. Krauth C, Hessel F, Hansmeier T, Wasem J, Seitz R, Schweikert B: Empirical standard costs for health economic evaluation in Germany – a proposal by the working group methods in health economic evaluation. Gesundheitswesen. 2005, 67: 736-746.

    Article  CAS  PubMed  Google Scholar 

  59. Statistisches Bundesamt [Federal Statistical Office]: Preise - Verbraucherpreisindizes für Deutschland (Lange Reihen ab 1948) [ Prices - Consumer Price Indices for Germany (Long Series from 1948)]. 2013, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  60. AOK Bundesverband [AOK Federal Association]: Krankheitsartenstatistik 2008 [Disease statistic 2008]. 2009, Berlin (Germany): AOK Federal Association

    Google Scholar 

  61. Deutsche Rentenversicherung [German statutory pension insurance]: Rehabilitation 2009. 2010, Berlin (Germany): German statutory pension insurance

    Google Scholar 

  62. Bundesministerium für Gesundheit [Federal Ministry of Health]: Drucksache 17/8332 - Fünfter Bericht Über die Entwicklung der Pflegeversicherung und den Stand der Pflegerischen Versorgung in der Bundesrepublik Deutschland [Printed Matter 17/8332- Fifth Report on the Development of the Mandaratory Care Incurance and the Current Status of Nursing Care in the Federal Republic of Germany]. 2012, Berlin (Germany): Federal Ministry of Health

    Google Scholar 

  63. Kreck S, Klaus J, Leidl R, Von Tirpitz C, Konnopka A, Matschinger H, König HH: Cost effectiveness of ibandronate for the prevention of fractures in inflammatory bowel disease-related osteoporosis: cost-utility analysis using a Markov model. Pharmacoeconomics. 2008, 26: 311-328.

    Article  PubMed  Google Scholar 

  64. Rothgang H: Vergütung von Leistungen der Häuslichen Krankenpflege nach § 37 SGB V. Reimbursement of Home Care Services in Accordance With § 37 SGB V] In Oral Presentation at: 4Bayrischer Tag der Ambulanten Pflege des bpa [Forth Bavarian day of Home Care by bpa] Munich (Germany). 2011, http://www.sfb597.uni-bremen.de/homepages/rothgang/downloads/110328_Rothgang_Verguetung_von_Leistungen_der_Haeuslichen_Krankenpflege%20_bpa.pdf.

    Google Scholar 

  65. van den Berg B, Brouwer WB, Koopmanschap MA: Economic valuation of informal care. An overview of methods and applications. Eur J Health Econ. 2004, 5: 36-45.

    Article  PubMed  Google Scholar 

  66. Statistisches Bundesamt [Federal Statistical Office]: Verdienst und Arbeitskosten 2009 [Earning and labour costs 2009]. 2010, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  67. Eurostat: Social security and other labour costs paid by employer % of total labour costs (10 employees or more). http://epp.eurostat.ec.europa.eu/tgm/refreshTableAction.do;jsessionid=9ea7d07d30da009798d90cda41d8b89fadddd0393cc1.e34MbxeSaxaSc40LbNiMbxeNaxuRe0?tab=table&plugin=1&pcode=tps00114&language=en.

  68. Luce B, Manning W, Siegel J, Lipscomb J: Estimating costs in cost-effectiveness analysis. Cost-Effectiveness in Health and Medicine. Edited by: Gold M, Siegel J, Russel L, Weinstein M. 1996, New York: Oxford University Press, 176-209.

    Google Scholar 

  69. Statistisches Bundesamt [Federal Statistical Office]: Mikrozensus Bevölkerung und Erwerbstätigkeit Stand und Entwicklung der Erwertstätigkeit Deutschland 2009 [Microcensus Population and Employment Status quo and developement of employment Germany 2009]. 2010, Wiesbaden (Germany): Federal Statistical Office

    Google Scholar 

  70. Heckmann M, Kettner A, Rebien K, Vogler-Ludwig K: Unternehmensbefragung im IV. Quartal 2009: Stellenbesetzung in Zeiten der Krise (IAB Kurzbericht 23/2010). 2010, Nürnberg (Germany): [The German Job Vacancy Survey in the IVth quarter of 2009: Staffing in times of crisis (IAB short report 23/2010)]

    Google Scholar 

  71. Eurostat: Annual Net Earnings. http://appsso.eurostat.ec.europa.eu/nui/show.do?dataset=earn_nt_net&lang=eng.

  72. Leicht H, Heinrich S, Heider D, Bachmann C, Bickel H, van den Bussche H, Fuchs A, Luppa M, Maier W, Mösch E, Pentzek M, Rieder-Heller SG, Tebarth F, Werle J, Weyerer S, Wiese B, Zimmermann T, König HH, AgeCoDe study group: Net costs of dementia by disease stage. Acta Psychiatr Scand. 2011, 124: 384-395.

    Article  CAS  PubMed  Google Scholar 

  73. Briggs A, Claxton K, Sculpher M: Decision Modelling for Health Economic Evaluation. 2006, Oxford: Oxford University Press, 1

    Google Scholar 

  74. Kanis JA, Odén A, McCloskey E, Johansson H, Wahl DA, Cooper C: A systematic review of hip fracture incidence and probability of fracture worldwide. Osteoporos Int. 2012, 23: 2239-2256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Hiligsmann M, Bruyere O, Ethgen O, Gathon HJ, Reginster JY: Lifetime absolute risk of hip and other osteoporotic fracture in Belgian women. Bone. 2008, 43: 991-994.

    Article  PubMed  Google Scholar 

  76. Doherty DA, Sanders KM, Kotowicz MA, Prince RL: Lifetime and five-year age-specific risks of first and subsequent osteoporotic fractures in postmenopausal women. Osteoporos Int. 2001, 12: 16-23.

    Article  CAS  PubMed  Google Scholar 

  77. Nayak S, Roberts MS, Greenspan SL: Impact of generic alendronate cost on the cost-effectiveness of osteoporosis screening and treatment. PLoS One. 2012, 7: e32879.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Melton LJ, Chrischilles EA, Cooper C, Lane AW, Riggs BL: How many women have osteoporosis?. J Bone Miner Res. 2005, 20: 886-892.

    Article  PubMed  Google Scholar 

  79. OECD: OECD Health Data: Economic References. http://dx.doi.org/10.1787/data-00548-en.

  80. Hiligsmann M, Gathon HJ, Bruyere O, Ethgen O, Rabenda V, Reginster JY: Cost-effectiveness of osteoporosis screening followed by treatment: the impact of medication adherence. Value Health. 2010, 13: 394-401.

    Article  PubMed  Google Scholar 

  81. Müller D, Pulm J, Gandjour A: Cost-effectiveness of different strategies for selecting and treating individuals at increased risk of osteoporosis or osteopenia: a systematic review. Value Health. 2012, 15: 284-298.

    Article  PubMed  Google Scholar 

  82. Kanis JA, Johnell O, Oden A, Börgstrom F, Zethraeus N, De Laet C, Jönsson B: The risk and burden of vertebral fractures in Sweden. Osteoporos Int. 2004, 15: 20-26.

    Article  CAS  PubMed  Google Scholar 

  83. Klein S, Gothe H, Häussler B, Kless T, Linder R, Schmidt T, Steinle T, Verheyen F, Hadji P: Versorgungsinanspruchnahme und Behandlungskosten bei Osteoporose. Analyse von Krankenkassen-Routinedaten im Rahmen der Bone Evaluation Study (BEST). Osteologie. 2013, 22: 277-282.

    Google Scholar 

  84. Heinrich S, Rapp K, Rissmann U, Becker C, König H-H: Service use and costs of incident femoral fractures in nursing home residents in Germany: the Bavarian Fall and Fracture Prevention Project. J Am Med Dir Assoc. 2011, 12: 459-466.

    Article  PubMed  Google Scholar 

  85. Cooper C, Cole Z, Holroyd C, Earl S, Harvey N, Dennison E, Melton L, Cummings S, Kanis J: Secular trends in the incidence of hip and other osteoporotic fractures. Osteoporos Int. 2011, 22: 1277-1288.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Hadji P, Claus V, Ziller V, Intorcia M, Kostev K, Steinle T: GRAND: the German retrospective cohort analysis on compliance and persistence and the associated risk of fractures in osteoporotic women treated with oral bisphosphonates. Osteoporos Int. 2012, 23: 223-231.

    Article  CAS  PubMed  Google Scholar 

  87. Hiligsmann M, Kanis JA, Compston J, Cooper C, Flamion B, Bergmann P, Body JJ, Boonen S, Bruyere O, Devogelaer JP, Goemaere S, Kaufman JM, Rozenberg S, Reginster JY: Health technology assessment in osteoporosis. Calcif Tissue Int. 2013, 93 (1): 1-14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgement

This work was funded by the German Federal Ministry of Education and Research (BMBF), Germany, FKZ:01EC1007C. The BMBF had no further role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; and preparation, review, or approval of the manuscript. Also we will thank Thomas Lehnert for proofreading the final version of our manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Florian Bleibler.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FB and HHK conceived the analysis. FB developed the simulation model, conducted the analyses and drafted the manuscript. FB and AJ prepared data as input for the model. KR, CB, HHK, FB, AJ participated in study design and implementation. HHK, AJ, KR and CB critically revised the manuscript for important intellectual content. All authors read and approved the final manuscript.

Electronic supplementary material

12913_2013_3421_MOESM1_ESM.pdf

Additional file 1: File name: Bleibler_BMCHSR_ESM.pdf (Electronic supplementary material (ESM). Includes all used input data of the model). (PDF 1 MB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bleibler, F., Rapp, K., Jaensch, A. et al. Expected lifetime numbers and costs of fractures in postmenopausal women with and without osteoporosis in Germany: a discrete event simulation model. BMC Health Serv Res 14, 284 (2014). https://doi.org/10.1186/1472-6963-14-284

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1472-6963-14-284

Keywords