Expected lifetime numbers and costs of fractures in postmenopausal women with and without osteoporosis in Germany: a discrete event simulation model
© Bleibler et al.; licensee BioMed Central Ltd. 2014
Received: 22 October 2013
Accepted: 21 May 2014
Published: 30 June 2014
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.
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.
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.
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.
KeywordsDiscrete event simulation Fractures Health care costs Osteoporosis
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 . 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%  (depending on risk profile)) or breast cancer (around 11% ).
Studies have shown that a decrease in bone mineral density (BMD) is a strong predictor for fractures on different sites of the human body [4–6]. 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 . 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” . 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 , health related quality of life  and are faced with an increased mortality . Especially for advanced elderly persons, fracture related functional impairment often necessitates external help from professional home care services, relatives or friends . In some cases institutionalization in a nursing home (NH) is unavoidable . On the societal level fractures have repeatedly been shown to cause a high economic burden [15–21]. As a consequence of the (expected) demographic change in Germany , the economic burden due to osteoporotic fractures is likely to increase strongly in the next decades .
Internationally, the economic burden of osteoporosis and osteoporotic fractures has frequently been investigated in the last decade [16, 18–21, 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 [24–27]. One study used an incidence based approach to estimate lifetime costs caused by hip fractures . 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 . 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.
A DES model [29–31] 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 . The main components in DES are entities, attributes, events, time, resources (optional) and queues (optional) . 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 . 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 . 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) . Optionally, DES can also handle resource constraints using queuing-systems and allow for interaction between individuals .
There is evidence that fractures increase the probability of institutionalization in a NH . 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 ; furthermore it is known that individuals living in a NH have higher mortality rates compared to those who do not at the same age . 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).
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 . 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 . A half-cycle correction was applied to long term care costs .
All other cost categories in the model are tracked as unit costs (rewards) per fracture event, similar to the common cohort model approach . 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 . 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 .
Epidemiological input data
Fracture probabilities – general population
Age-dependent fracture probabilities for all considered fracture types were calculated based on official German population  and hospital discharge data  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) , some of clinical vertebral, humerus, pelvis and wrist fractures are treated exclusively in an outpatient setting [41–44]. 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 , pelvis , wrist ) (see Additional file 1: A.3.a-b).
Total yearly fracture probability of the female general population by age and fracture type
Fracture probabilities for women with and without a previous fracture
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 . 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. . Necessary input data such as fracture specific relative risks by a decrease of one standard deviation in BMD measured at the femoral neck (RRfx) [4–6], the osteoporosis BMD threshold , and age-dependent population BMD values (reference values from NHANES III)  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.  (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  (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 . Moreover, they frequently depend on external care or need to be institutionalized as a result of a fracture . Similar to hip fractures, other types of osteoporotic fractures have been shown to increase the probability of institutionalization . 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  were transformed into NH institutionalization probabilities . 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  (see Additional file 1: A.4.a-b).
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) . 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 . 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
Overview of direct unit costs in € by cost category and fracture type
Hospital treatment (plus outpatient aftercare costs)
Rehabilitation (if required after hospitalization)
Outpatient costs (if no hospitalization required)
Professional home care (age > 65, not in NH)
Informal home care (age > 65, not in NH)
Yearly long term care cost (age > 65, in NH)
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 . 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  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  to obtain DRG costs. Capital costs were calculated by combining the fracture specific length of stay with a daily capital cost rate  (inflated to 2009 ) (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  with a daily cost rate  (see Additional file 1: B.1.b). Long term care costs are based on the official care statistic 2009 . 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  were added to the yearly unit costs (see Additional file 1: B.1.c).
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 , inflated to 2009 . Necessary resource use data in the three outpatient categories were taken from a German cost-effectiveness study . 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  and hourly unit costs . As no information on fracture type specific hours of home care was available for Germany, a study from Austria was used . 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 , fracture type specific amount of hours spend by relatives of a fracture patient were evaluated with the market cost approach (proxy good) . 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 , 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  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).
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 . Productivity loss was evaluated using the human capital  and the friction cost approach . Necessary information, such as employment rate , yearly average gross earnings  corrected by employer share of social contribution , time away from work (due to disability and rehabilitation ), length of friction period  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  (see Additional file 1: B.2.a-b).
Overview of the model assumptions
Assumptions regarding event probabilities
We applied osteoporosis prevalence rates and BMD-values from US-NHANES III reference data
We estimated “total” fracture probabilities by dividing fracture probabilities based on hospital cases with age-independent hospital probabilities
We assumed highest fracture related NH probability when more than one fractures occurs in the same time interval
We modeled fracture related entry in a NH only after a hospital stay
We assumed that only NH entries within 3 months after a fracture may be attributable to the fracture event itself
We applied age-dependent relative fracture risk by one standard deviation decrease in BMD to hip fractures and age-independent relative risks to other fractures
We assumed that osteoporosis risk attributions were calculated exclusively on BMD values measured at the femoral neck
We assumed that osteoporosis prevalence rates do not differ between women living in a NH and women who do not
We applied relative fracture risk and prevalence for previous fractures from an international meta-analysis
We applied relative fracture risk by one standard deviation decrease in BMD from international studies
We applied fracture mortality data from a Canadian study
We assumed the highest fracture excess mortality when more than one fracture occurs in the same time interval
We allowed first entry in a NH firstly for women aged 65 or older
We assumed that individuals in a NH remain there for their remaining lifetime
We assumed that patients with osteoporosis will have osteoporosis for their remaining lifetime
We allowed a maximum possible age of 100 years
Assumptions regarding costing
We assumed that rehabilitation probabilities after a hospital stay do not differ between women living in NH and those who do not.
We applied Austria data for average hours of informal and professional home care by fracture type, also we assumed that the consumed hours are equivalent for hip, other femur and pelvis
We assumed age-dependent fracture unit costs
We assumed that the outpatient costs for humerus and wrist as well as the costs for pelvis, other femur and hip fractures are equivalent
We took outpatient resource use data from a study considering fracture patients with inflammatory bowel disease
We assumed that average informal and professional home care costs are only applicable for individuals not living in NH aged older than 65 years
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%.
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  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 . 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) . Common distributional assumptions were used for costs, probabilities and relative risks  (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) .
Undiscounted and discounted fracture lifetime costs (€) of a 50 year old woman, by fracture type for two risk classes of osteoporosis and excess
By fracture type
Average risk (95%UI)
Not at risk (95%UI)
Average risk (95%UI)
Not at risk (95%UI)
Undiscounted and discounted fracture lifetime costs (€) of a 50 year old woman, by healthcare sector for two risk classes of osteoporosis and excess
By healthcare sector
Average risk (95%UI)
Not at risk (95%UI)
Average risk (95%UI)
Not at risk (95%UI)
Long term care
Physiotherapy, physician, analgesics
Prof. home care
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.
Modeled and expected hip fracture rates and proportions of women with 0 to 6 hip fractures
Fractures per woman
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.
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 , 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  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)  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  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.  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 . 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 . 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 . 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” . 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 . 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” . 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 , 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 ). 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 . 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 . However, this may be due to the low osteoporosis specific treatment prevalence in Germany . 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 . 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 . 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 . 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 . Furthermore, only 14% of fracture patients receive an osteodensitometry, which may indicate an under-detection of osteoporosis . Moreover, the adherence to oral osteoporosis medications is suboptimal . 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  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 . 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 .
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.
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.
- Johnell O, Kanis J: Epidemiology of osteoporotic fractures. Osteoporos Int. 2005, 16 (Suppl 2): S3-S7.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Armstrong K, Eisen A, Weber B: Assessing the risk of breast cancer. N Engl J Med. 2000, 342: 564-571.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- Anonymous: Consensus development conference: prophylaxis and treatment of osteoporosis. Am J Med. 1991, 90: 107-110.View ArticleGoogle Scholar
- Wildner M: Osteoporose [Osteoporosis]. Dtsch Med Wochenschr. 2001, 126: A1170-A1172.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Börgstrom F, Sobocki P, Ström O, Jönsson B: The societal burden of osteoporosis in Sweden. Bone. 2007, 40: 1602-1609.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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 OfficeGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Brecht JG, Schädlich PK: Burden of illness imposed by osteoporosis in Germany. HEPAC Health Economics in Prevention and Care. 2000, 1: 26-32.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Weyler EJ, Gandjour A: Sozioökonomische bedeutung von hüftfrakturen in deutschland [socioeconomic burden of Hip fractures in germany]. Gesundheitswesen. 2007, 69: 601-606.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Caro JJ, Moller J, Getsios D: Discrete event simulation: the preferred technique for health economic evaluations?. Value Health. 2010, 13: 1056-1060.View ArticlePubMedGoogle Scholar
- Sobolev B, Sanchez V, Kuramoto L: Health Care Evaluation Using Computer Simulation: Concepts, Methods, and Applications. 2012, New York, Heidelberg, Dordrecht, London: SpringerView ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Sonnenberg FA, Beck JR: Markov models in medical decision making: a practical guide. Med Decis Making. 1993, 13: 322-338.View ArticlePubMedGoogle Scholar
- Stout NK, Goldie SJ: Keeping the noise down: common random numbers for disease simulation modeling. Health Care Manag Sci. 2008, 11: 399-406.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 OfficeGoogle Scholar
- 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 DiseasesGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Finnern HW, Sykes DP: The hospital cost of vertebral fractures in the EU: estimates using national datasets. Osteoporos Int. 2003, 14: 429-436.View ArticlePubMedGoogle Scholar
- Cooper C, O’neill T, Silman A: The epidemiology of vertebral fractures. Bone. 1993, 14: 89-97.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 UnitGoogle Scholar
- Podgor MJ, Leske MC: Estimating incidence from age-specific prevalence for irreversible diseases with differential mortality. Stat Med. 1986, 5: 573-578.View ArticlePubMedGoogle Scholar
- 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 OfficeGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 OfficeGoogle Scholar
- 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 CareGoogle Scholar
- Rice DP, Cooper BS: The economic value of human life. Am J Public Health Nations Health. 1967, 57: 1954-1966.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- 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
- 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.View ArticlePubMedGoogle Scholar
- 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 OfficeGoogle Scholar
- AOK Bundesverband [AOK Federal Association]: Krankheitsartenstatistik 2008 [Disease statistic 2008]. 2009, Berlin (Germany): AOK Federal AssociationGoogle Scholar
- Deutsche Rentenversicherung [German statutory pension insurance]: Rehabilitation 2009. 2010, Berlin (Germany): German statutory pension insuranceGoogle Scholar
- 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 HealthGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- Statistisches Bundesamt [Federal Statistical Office]: Verdienst und Arbeitskosten 2009 [Earning and labour costs 2009]. 2010, Wiesbaden (Germany): Federal Statistical OfficeGoogle Scholar
- 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.
- 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
- 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 OfficeGoogle Scholar
- 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
- Eurostat: Annual Net Earnings. http://appsso.eurostat.ec.europa.eu/nui/show.do?dataset=earn_nt_net&lang=eng.
- 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.View ArticlePubMedGoogle Scholar
- Briggs A, Claxton K, Sculpher M: Decision Modelling for Health Economic Evaluation. 2006, Oxford: Oxford University Press, 1Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Melton LJ, Chrischilles EA, Cooper C, Lane AW, Riggs BL: How many women have osteoporosis?. J Bone Miner Res. 2005, 20: 886-892.View ArticlePubMedGoogle Scholar
- OECD: OECD Health Data: Economic References. http://dx.doi.org/10.1787/data-00548-en.
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6963/14/284/prepub
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.