- Research article
- Open Access
- Open Peer Review
Health economic evaluation of vaccination strategies for the prevention of herpes zoster and postherpetic neuralgia in Germany
© Ultsch et al.; licensee BioMed Central Ltd. 2013
- Received: 15 April 2013
- Accepted: 18 September 2013
- Published: 26 September 2013
Herpes zoster (HZ) is a self-limiting painful skin rash affecting mostly individuals from 50 years of age. The main complication is postherpetic neuralgia (PHN), a long-lasting pain after rash has resolved. A HZ-vaccine has recently been licensed in Europe for individuals older than 50 years. To support an informed decision-making for a potential vaccination recommendation, we conducted a health economic evaluation to identify the most cost-effective vaccination strategy.
We developed a static Markov-cohort model, which compared a vaccine-scenario with no vaccination. The cohort entering the model was 50 years of age, vaccinated at age 60, and stayed over life-time in the model. Transition probabilities were based on HZ/PHN-epidemiology and demographic data from Germany, as well as vaccine efficacy (VE) data from clinical trials. Costs for vaccination and HZ/PHN-treatment (in Euros; 2010), as well as outcomes were discounted equally with 3% p.a. We accounted results from both, payer and societal perspective. We calculated benefit-cost-ratio (BCR), number-needed-to-vaccinate (NNV), and incremental cost-effectiveness ratios (ICERs) for costs per HZ-case avoided, per PHN-case avoided, and per quality-adjusted life-year (QALY) gained. Different target age-groups were compared to identify the most cost-effective vaccination strategy. Base-case-analysis as well as structural, descriptive-, and probabilistic-sensitivity-analyses (DSA, PSA) were performed.
When vaccinating 20% of a cohort of 1 million 50 year old individuals at the age of 60 years, approximately 20,000 HZ-cases will be avoided over life-time. The NNV to avoid one HZ (PHN)-case was 10 (144). However, with a BCR of 0.34 this vaccination-strategy did not save costs. The base-case-analysis yielded an ICER of 1,419 (20,809) Euros per avoided HZ (PHN)-case and 28,146 Euros per QALY gained. Vaccination at the age of 60 was identified in most (sensitivity) analyses to be the most cost-effective vaccination strategy. In DSA, vaccine price and VE were shown to be the most critical input-data.
According to our evaluation, HZ-vaccination is expected to avoid HZ/PHN-cases and gain QALYs to higher costs. However, the vaccine price had the highest impact on the ICERs. Among different scenarios, targeting individuals aged 60 years seems to represent the most cost-effective vaccination-strategy.
- Herpes zoster
- Postherpetic neuralgia
Before implementation of routine childhood varicella vaccination in 2004, the lifetime risk of acquiring a varicella-zoster-virus (VZV) infection in Germany was almost 100% [1, 2]. Children under 10 years of age were predominantly affected by the infection, which clinically manifests as chickenpox (varicella). After recovery from varicella, the virus remains latent life-long in individual’s dorsal root ganglia . Later in life, the VZV can reactivate and manifest as shingles (herpes zoster, HZ) as a result of decreasing VZV-specific T-cell-immunity [4, 5]. Natural waning immunity and other causes like psychological stress or immunosuppression can contribute to a VZV-reactivation, too .
HZ is a painful and self-limiting skin-rash that lasts approximately four weeks [7–9]. The main complication of HZ is postherpetic neuralgia (PHN) a long-lasting pain in the formerly HZ-affected skin region after rash has resolved [5, 10–12]. Both, HZ and PHN cause limitations on quality of life (QoL) [13–16]. Approximately 20-30% of individuals experience at least one HZ-episode in their life . The risk and burden of HZ increase with age, most cases being 50 years and older [4, 18–21]. In Germany, the risk in this age-group is estimated at 9 per 1,000 persons per year (PY), which is comparable to other industrialized countries [21–23]. Estimates for the proportion of HZ-patients developing PHN range between 6 and 7% [21, 23]. The options for therapy and prevention of HZ and PHN are limited [24–28].
In 2006, the first vaccine for the prevention of HZ and PHN for people 50 years of age and older was licensed in Europe . This live-attenuated vaccine demonstrated its efficacy in two double-blinded, placebo-controlled randomized trials (RCT), which included healthy individuals aged >50 years [29, 30]. Another HZ-vaccine candidate is currently tested in clinical phase III trials and might be available in the near future . As of June 2013, the licensed vaccine was not yet available on the market in most European countries. However, the market launch in Germany might be expected soon. To support an informed decision-making of a potential vaccine-recommendation by relevant authorities, we performed a health economic evaluation of routine HZ-vaccination among the elderly and compared various vaccination strategies.
State (From state → To state)
Death → Death
Any state → Death
Status quo [no vaccination]
Healthy → Herpes zoster
Healthy → Healthy
= 1 - α i - β i
Herpes Zoster → Death
Herpes Zoster → PHN
Herpes Zoster → Healthy after disease
= 1 - α i - δ i - ϵi
PHN → PHN
PHN → Healthy after disease
= 1 - α i - ηi
Healthy after disease
Healthy after disease → Healthy after disease
Healthy after disease → Herpes zoster
= 1 - α i - κ i
Vaccine-Scenario [adapted transition probability due to vaccine efficacy (VE)]
Healthy → Herpes zoster
Healthy → Healthy
Herpes zoster → PHN
Healthy after disease
Healthy after Disease → Healthy after disease
Healthy after Disease → Herpes zoster
% of HZ-cases
95% CI (Cases/1,000 PY)
95% CI (% of HZ-cases)
Treatment costs per case (€)
Payer 95% CI
Societal 95% CI
Payer 95% CI
Societal 95% CI
Background mortality [probability (%) p.a.]***
HZ related mortality [deaths/100,000 PY]***
Other input parameter
DSA [lower/higher value]
Vaccine Price/dose (€)
Vaccine administration costs (€)**
Discount rate (% p.a.)
Waning immunity rate (% p.a.)
Coverage (%) assumption
HZ susceptibility**** (%)
EQ-5D utilities for healthy individuals (age in years)
Several studies have assessed the risk of HZ-recurrence, i.e. the occurrence of a subsequent HZ-episode [4, 20, 44–47]. There are discrepant findings on the frequency of this phenomenon in the literature. However, studies suggesting that the risk of HZ-recurrence is low have a rather short follow-up time and/or a low number of HZ-patients included [20, 44, 46]. In contrast, studies with a follow-up period of several years conclude consistently that the risk of a subsequent HZ-episode is a critical factor when considering HZ-related disease burden [4, 45, 47]. A large study addressing specifically the risk of HZ-recurrence, which involved 1,669 HZ-cases with an average follow-up of 7.3 years, suggests that the incidence of HZ-recurrence is similar to the incidence of a first HZ-episode . We have accounted for this risk by using same input-data for HZ-incidence and HZ-recurrence (see Sections ‘Base-case analysis’ and ‘Structural sensitivity analysis’ for further exemplifications). Incidence of HZ-associated death was available from a retrospective data-analysis in Germany . The probability of death due to all causes by age (‘background mortality’) in the general population living in Germany was derived from public statistics of the federal statistical office in Germany  (Table 2).
For cost-utility analyses, quality of life (QoL) data were considered. Since HZ-related QoL and pain intensity data among HZ- or PHN-patients is lacking for Germany, we used data from a study by Drolet and colleagues, which was conducted in a Canadian population of individuals 50 years of age and older without stratification by HZ severity . This fitted well to our input-data (incidence and costs), which were derived from all patients without specification of pain intensity. We informed our model with the corresponding ‘EQ-5D scores’ (utilities). We fitted the data by calculating the age-specific average utilities based on values of day 0 and day 30 as HZ-related utilities, and the age-specific average based on values of day 90 and day 180 as PHN-related utilities, as defined by Drolet et al. . For individuals located in the states ‘Healthy’ and ‘healthy after disease’ we chose 1 as a value for baseline utility in the base-case. To account for vaccine-related adverse reactions on the QoL, we assumed a utility of 0.99 for two days per vaccinee, since the vaccine seems to be generally well tolerated [29, 49].
Cost data was derived from the same study where the incidence data was utilized. Data was stratified by HZ and PHN, by age-groups, and by perspectives (Table 2). As all cost values in our model were at Euro 2010 price levels, we considered the German pharmacy retail price of HZ-vaccine for the same year reduced by the obligatory pharmacy and manufacturer discounts. Thus, the HZ-vaccine price per dose in our model was set to 140.48 Euros . We used the mean administration costs per pneumococcal vaccine dose in Germany in 2010 (7 Euros) as a proxy for HZ-vaccine administration in the base-case assumption and considered the regional minimum and maximum values for sensitivity analyses (SA, Table 2).
All epidemiological, demographic, and VE data were, if necessary, converted into age-group specific probabilities or relative risks with respect to the model’s cycle-length (see Table 1 and Figure 1) .
In base-case analysis, the model’s start age was set to 50 years and the age of vaccination to 60 years. For incidence and recurrence rates as well as for treatment costs per case we used the respective point estimators from studies as presented in ‘Input data’ section and Table 2. For the tunnel state ‘Healthy after Disease’, we defined that a recurrence cannot occur within three years (12 cycles) after preceding the first HZ-episode . Furthermore we adjusted the tunnel state ‘PHN’ to keep individuals three cycles, since the average duration of PHN is 9 month (confer ‘Model design’ section). The utility values, waning rate, vaccine price, and vaccine honorary were set as described in ‘Input data’ section. For ICER calculation, costs and outcomes were discounted equally with 3% per annum as suggested in Germany  (Table 2). Additionally we conducted a ‘Break-even’ analysis. For this we used the base-case conditions to identify the vaccine price at which the ICERs represent cost-savings.
Univariate deterministic sensitivity analyses
With univariate deterministic sensitivity analyses (DSA) we investigated the impact of different age-specific incidences and treatment costs, VE, waning immunity rates, discount rates, and vaccine prices as well as administration costs and utility values (Table 2). For measuring the impact of an absent HZ-recurrence we set the respective value to zero. To measure the impact of a high HZ-recurrence, we used the upper limit 95% confidence interval (CI) of HZ-incidence values and removed the ‘tunnel’-state condition. Hence, in the latter scenario a HZ-episode was theoretically possible right after the previous one. Since in a static model a change in coverage would not affect the ICERs, we neglected this variation in the SA. In contrast to the base-case we changed the age of vaccination from 60 years to 50, 55, 65, 70, 75, and 80 years and stratified the analyses by outcome and perspective.
Structural sensitivity analysis
In a ‘Best Case Scenario’ we assumed a situation in favor of the new HZ-vaccine and changed input-data accordingly. We used the higher values of 95% CI of the HZ and PHN disease-incidences and costs per case. Furthermore, we set the recurrence rate as high as described before and used the lower utilities for HZ and PHN, respectively. Higher 95% CI values were taken for VE whereas the waning rate was set to zero. Vaccine administration costs were lowered to 5.20 Euros and the vaccine price by 50%.
The ‘Worst case Scenario’ represented possible least preferable results. In terms of epidemiology we used the lower 95% CI estimates of the HZ- and PHN-incidence. For HZ and PHN-treatment costs we took the lower 95% CI values and neglected recurrence by setting its probability to zero. We considered higher utilities, an annual waning rate of 20% and the lower 95% CI values of the VE. We increased the vaccine price by 50% and set the administration costs to 10.81 Euros per vaccination.
The variation of ‘Duration of stable VE’ aims to reduce the potential uncertainty associated with the model’s base-case that considers 10 years of stable VE following vaccination. Since this is an assumption based on trial data (see section ‘Vaccine characteristics’) we conducted analyses where we changed the period of stable VE following vaccination from 0, 5, 15, to 20 years.
A combined analysis of the impact of different annual waning immunity rates (1%, 5%, 8.3%, and 20%), duration of stable VE (0, 5, 10, and 15 years), and the age at vaccination (50, 60, and 70 years) on the costs per QALY gained was performed to identify the optimal age at vaccination.
Furthermore, we implemented different ‘booster scenarios’ into the model. Therefore, HZ-vaccinated individuals were administrated a second dose 20 years after first HZ-vaccination. This booster was given to 50% of these formerly HZ-vaccinated individuals located in the states ‘Healthy’ or ‘Healthy after disease’ within the vaccine-scenario (Figure 1). Besides additional costs due to vaccination and a two day reduction of individual’s QoL because of vaccine related adverse reactions the booster shot increased the vaccine-induced immunity. We modeled two oppositional booster scenarios: In favor for the HZ-vaccine and its booster we assumed in the ‘Best Case Booster Scenario’ an increase of the vaccine-induced immunity to the same level as after the first vaccination (Table 2), and a constant VE for 10 years before the waning with a rate of 1% p.a. began (Figure 2). In the ‘Worst Case Booster Scenario’ the booster increased the vaccine induced immunity only to the respective level of the actual age as presented in RCTs considering VE (Table 2) and this immunity diminished immediately by the annual waning rate of 20%. In order to account for the influence of age at vaccination on the ICERs, we modeled both booster scenarios not only for the base-case (60 years of age at vaccination) but also for 50, 70, and 80 years at vaccination.
To measure the impact of varying ‘duration of PHN’ on ICER (cost per QALY gained), we changed PHN duration from 9 months (base-case) to 6, 15, and 36 month, respectively.
In order to assess the impact of a ‘change in baseline utility’ values for healthy individuals, we used data from a study conducted in Germany . In this study 2,022 healthy individuals between 16 and 93 years of age were interviewed using the EQ-5D questionnaire. Furthermore three different accounting models (sum model, linear model, and a multiplicative model) were used to calculated values resulting from EQ-5D items. For this SA we used the age-specific utilities from the linear model . Since our model uses throughout gender unspecific input-data, we calculated age-specific average values from both gender (Table 2).
Multivariate probabilistic sensitivity analyses
To identify the robustness of the base-case and the probability of extreme favorable or unfavorable results and to get more insights into uncertainties, we performed multivariate probabilistic sensitivity analyses (PSA) . We conducted a Monte Carlo simulation with 10,000 runs for each outcome (HZ- and PHN-cases avoided, QALY gained). The respective input-data were drawn simultaneously (multivariate) from probability distributions within the ranges considered in DSA (Table 2). For guidance purposes four thresholds (20,000, 30,000, 50,000, and 100,000 Euros per QALY gained) representing different fictitious willingness to pay (WTP) levels were implemented within the scatter plot of the PSA considering costs per QALY gained.
External comparison and software
For a proper external comparison of our model results we considered findings from other international health economic evaluations. For a better comparison, we inflated the cost data from other studies to price levels of 2010 (reference year) within each study’s country using consumer prices indices and converted these respective amounts by the gross domestic product purchasing power parity (GDP PPP) into Euro amounts with Germany as reference country based on data from the Organization for Economic Co-operation and Development (OECD) [63, 64]. Hence, amounts were presented in PPP-Euros. Furthermore, we considered those analyses in respective studies that were similar to our base-case. Thus, we did not necessarily refer to base-case results from the other studies in the ‘Discussion’ section.
Software used was: TreeAge® Pro 2012 (TreeAge Software, Williamstown, MA), Microsoft® Excel 2010 (Microsoft, Redmond, WA), and STATA® 12.1 (StataCorp, College Station, TX).
Results from base-case analyses
a: Overall effects [cohort of 1 million individuals not discounted]
Number of HZ-cases
Costs of treatment (€)
Costs of vaccination (€)
Total costs (€)
b: Number needed to vaccinate (NNV)
Avoiding one HZ-case
Avoiding one PHN-case
Gaining one QALY
c: Costs and effects in both scenarios of the model
Results from structural sensitivity analyses (discounted)
a: Extreme scenarios and variations in the duration of stable VE
ICER by outcome
b: Duration of stable VE [years]
ICER by outcome
c: Booster scenarios [by age at vaccination in years]
ICER by outcome
ICER by outcome
d: Variations in the length of PHN [month]
ICER by outcome
Univariate deterministic sensitivity analyses
Considering the age at vaccination, cost-effectiveness is decreasing (indicated by higher ICERs) between the age of 55 and 75 years. Due to low VE, vaccination of 80 year old individuals seems least cost-effective, especially when considering costs per HZ-case avoided (Figure 3a). In terms of PHN-cases avoided the most cost-effective age at vaccination is 65 years. However, the base-case results are only slightly higher (hatched bar in Figure 3b). The lowest ICER with respect to QALYs gained is represented by the base-case (hatched bar in Figure 3c). Due to indirect costs occurring in the working population, the difference between PP and SP is higher between 50 to 60 year old individuals and decreases with increasing age.
Structural sensitivity analysis
The best case booster scenario is more cost-effective than the worst-case scenario independently from age at vaccination. However, the age-specific differences between both scenarios become smaller with increasing age: While ICERs at age 50 at vaccination differ between 30 to over 40% between both scenarios, the ICERs at age 80 at vaccination are almost equal (<0.1%, Table 4c). In Table 4d ICERs are shown by different duration of PHN. An increase from 6 to 9 month reduces ICERs by over 10%; an increase from 9 to 15 month causes a reduction by almost 20%. A PHN-duration of 3 years instead of 9 month decreases the ICER by almost 85% (Table 4d). Using age-specific baseline utility values causes an increase of about 30% compared to base-case figures. The ICER results in 36,629 Euros per QALY gained.
Multivariate probabilistic sensitivity analyses
We developed a static Markov-model and conducted various SA to assess the cost-effectiveness of routine HZ-vaccination in Germany and to identify the most cost-effective strategy when targeting specific age-groups for vaccination. HZ-vaccination of 20% individuals of a cohort of 1 million 50 year old persons at age 60 would avoid 20,000 HZ-cases and reduce the overall HZ-treatment costs by 10 million Euros. However, with a BCR of 0.34 for the base-case, our results show that vaccination against HZ and PHN is not a cost-saving measure. Univariate analyses revealed that the price had the greatest impact on ICERs in the DSA. In the break-even analysis with all other variable set at base-case scenario, a price per vaccine dose of <26.50 Euros caused cost-savings.
Varying several variables in favor for the vaccine (including a 50% reduction in the price per vaccine dose to 70 Euros) the intervention becomes almost cost-saving. However, under least favorable conditions (low estimates for HZ/PHN-related incidence and treatment costs, low HZ-recurrence rate, high utilities, low VE, a high waning immunity rate, high vaccine administration costs, and a high vaccine dose price) the model accounts very high ICERs. In general, health economic evaluations of new vaccines are often impacted by two important input factors: incidence of the target disease and VE. However, since the incidence data are rather reliable and range within narrow CIs, the wide ranges of VE (especially for the protection against PHN), which are based on trial data, mainly force the differences of these two extreme scenarios. Furthermore, the variation of the vaccine price (-50% vs. +50%) and the variation of disease recurrence seem to influence the ICERs, too. However, these two extreme scenarios of varying several input parameter either in favor for the vaccine or under least favorable conditions, can be seen as unlikely outer boundaries of the whole SA’s scope.
Due to decreasing VE by age, an increasing age at vaccination usually leads to higher ICERs indicating a lower cost-effectiveness. However, before the age of 60 the ICERs are not necessarily more favorable, due to lower HZ and PHN-related incidence and assumed waning immunity. Hence, targeting persons aged 60 years is likely the most cost-effective vaccination strategy. The variation in the annual waning rate of vaccine-induced immunity had rather limited impact on the ICERs. Altering the period of stable VE following vaccination from 0 to 20 years produces a relative wide array of ICERs. The combined analysis of varying waning immunity rates and durations of stable VE has a rather high impact on ICERs by age at vaccination. However, this analysis confirms that age 60 seems to be the optimal age at vaccination if the annual waning rate is ≥5%. The booster scenarios showed an increasing similarity between best and worst booster scenario with increasing age at vaccination. This confirms the assumption that independently from how a fictitious booster scenario is designed age of vaccination is one of the most relevant factors. The most cost-effective age of vaccination changes from 50 to 60, when moving from best to worst booster scenario. A variation of PHN-duration (6–15 month) had a moderate impact in ICERs. Whereas an extreme extension of PHN-duration up to 3 years downsizes the costs per QALY gained enormously. However, since recent studies confirm an average PHN-duration of several month, our base-case assumption of 9 month seems to be justified [29, 65, 66].
In terms of internal validity, we compared the epidemiology reported in the literature for Germany (used as model input parameter) and the model results. The model slightly overestimated HZ and accounted less than one percent more HZ-cases per age-group than represented by input-data. Regarding PHN-cases, the model calculated 4% more PHN-cases in age-groups 50 to 60 years. In older age-groups the overestimation was less than 1%. When implementing vaccination into the model, the accounted cases were reduced according to the VE implemented into the model. Thus, the model’s internal validity can be considered as good.
To date there is no other health economic evaluation of the HZ-vaccine from Germany. However, we identified 14 studies from other European countries [17, 34, 40, 54, 67–71] and North-America [56, 72–75]. We found a wide range in vaccine prices from 43.85 to 438.5 PPP-Euros in one US study  and 81 to 147.32 PPP-Euros in the other studies. From SP for vaccination of individuals of 60 years results ranged from cost-savings  to 130,097 PPP-Euros per QALY gained . From PP ICERs ranged from 6,809 , 1,200–46,968  to 40,050 PPP-Euros per QALY gained . Our base-case results range well within these international findings. However country-specific issues like vaccine-price, disease epidemiology, price levels, and treatment pathways as well as model-specific issues like model structure, specific assumptions, and dealing with uncertainty hinder a full comparability. Hence, more uniform methods in studies are needed to make results more comparable .
One limitation of our health economic model is the absence of utility-data considering the impact on health related QoL caused by HZ and PHN specific for Germany. Instead, we used data derived from a Canadian study without country-specific adaptation, which might not necessarily represent the real impact of HZ/PHN on QoL in Germany . Other studies have reported a higher limitation on QoL caused by PHN [13, 77]. Thus, the utilized values in our model might underestimate the impact on QoL due to PHN, whereas the HZ utilities used in our model might overestimate the impact of HZ on QoL [29, 35, 78]. However, we reduced this uncertainty by varying respective utility values within SAs. Furthermore, we set the baseline utility value for healthy individuals in base-case to 1. Thus, the impact of HZ and PHN and therefore the effect of the vaccine might be overestimated, since age-specific utilities among elderly tend to be lower than 1. Even though for health economic evaluation the utilization of accurate utility data for QoL assessments is critical, in Germany age-specific utility data for healthy but also for the most indications are scarce since cost-utility analyses do not have the relevance in decision-making in Germany as in other countries. We identified in a literature search in total five studies evaluating the QoL among healthy individuals in Germany [61, 79–82]. Since all studies have certain limitations (e.g. only visual analog scale values was reported, no age-specific values reported, small study sample) we decided against using values of one of these studies for baseline utility values within the base-case analysis. However, by considering age-specific baseline utility values in SA we provided insights into this factor. The cost data used in the model were derived from a database from one large regional SHI in Germany and might not necessarily be representative for the whole population living in Germany. However, since countrywide treatment guidelines for HZ and PHN exist and since prices are mandatory for all SHIs in Germany, we believe that the utilization of these input-data is justified and representative. Furthermore our incidence data included both immune-competent and incompetent individuals. Since the live-attenuated vaccine is licensed for immune-competent individuals only, the incidence figures used in the model might be slightly overestimated. Our treatment cost input-data did not include privately covered expenditures for health-care services and over the counter drugs. Hence, costs from SP might be underestimated. However, since HZ and PHN-treatments are usually covered by the SHI, we suspect that the level of underestimation is low. As evidence was lacking concerning the duration and waning of vaccine-induced immunity, we had to make a few assumptions but included them in the SA. Finally, a real cost-effectiveness threshold concerning the WTP does not exist for the German health-care system. Therefore especially the ICERs and PSA results on costs per QALY gained have to be interpreted with caution. However, when comparing different age-groups to be targeted by routine HZ-vaccination and when comparing different scenarios (e.g. with and without booster vaccination), the lack of a cost-effectiveness threshold for Germany does not constitute a relevant limitation.
Our model provides valuable analyses and insights when considering the implementation of an efficient strategy for the prevention of HZ and PHN, and it will guide decision-making when developing a vaccination recommendation for Germany. First, this model reflects the efficacy of the vaccine quite well, since the HZ- and PHN-definitions as well as age-strata used in this model were similar to those used in the clinical-trials [21, 29]. HZ-cases with a clinical diagnosis but also with a ‘suspected’ diagnosis were included in both data-sets. Second, the definition of PHN in the clinical trial was pain persisting more than 90 days after rash onset. This matches exactly with the definition in our input-data, in which HZ-cases became PHN-cases when they were diagnosed or received a PHN-specific medication at least three month after HZ-diagnosis. Second, a further strength of our analysis is the intensive parameterization during modeling. The base-case and the SAs demonstrate that on the one side a careful selection of reliable input-data is important, but on the other side a wide range of SAs has to be conducted to reduce uncertainty within model results. Especially HZ- and PHN-related incidence and VE have to be incorporated with caution, since these factors influence the analyses considerably. However, for our model incidence data were utilized from a large study recently conducted in Germany. These data correspond well with results from a nationwide incidence study in Germany and another retrospective data analysis in Germany, but also with study results from other countries [17, 20, 21, 23, 36, 37, 83–85]. Thus, this incidence input-data can be regarded as rather robust. SAs with variations of these data within the respective CIs have little impact on our results. However, vaccine characteristics are based on a less rich data fundament. While data on the VE in individuals from the age of 50 years is available, data on the duration of vaccine-induced protection and waning rates of vaccine-induced immunity is lacking and assumptions had to be made. Therefore, VE data were analyzed in SAs to assess the associated uncertainty. Based on one clinical-trial we assumed the period of stable VE following vaccination to be ten years. Since this assumption carries some uncertainty, we conducted a structural SA in which we varied the period of stable VE in order to assess the associated impact on results. Furthermore, evidence on the exact relative annual waning rate was utilized from literature. To analyze the impact of waning on ICERs we neglected the existence of a waning rate in one scenario and then subsequently increased it within SAs. A combined sensitivity analysis of varying waning immunity rates and durations of stable VE illuminated their combined impact on the optimal age at vaccination. Going one step beyond, we analyzed also the potential impact of differently designed fictitious booster scenarios on the results. This enables to assess booster scenarios way before evidence on the potential need for booster is established . Finally, the SA in respect to HZ-recurrence shows a considerable impact of this parameter on the results. However neglecting HZ-recurrence has a lower impact than considering a high HZ-recurrence. This analysis shows that there is an urgent need to establish more evidence on HZ-recurrence on the long term view. Since HZ-recurrence but also HZ-booster vaccination are important parameters when conducting a health economic evaluation of HZ-vaccination in a given country, we believe that our results are also of high interest to other countries that consider the introduction of routine HZ-vaccination in their health-care systems. For our study we were able to use numerous country-specific input-data of high quality.
Routine HZ-vaccination in the elderly is estimated to avoid a substantial number of HZ- and PHN-cases in Germany. However, under base-case assumptions the costs of vaccination are unlikely to be compensated by lower treatment costs. The vaccine price was identified to be crucial when considering overall cost and ICERs. In view of these findings, a responsible pricing seems needed. Targeting individuals aged 60 years is likely to be the most cost-effective vaccination strategy for the prevention of HZ and PHN in Germany. If another HZ-vaccine product becomes available on the market in the future, a comparative health-economic evaluation can be conducted by utilizing our model. However, results from head-to-head VE comparisons would be desirable for such an analysis.
The authors gratefully acknowledge the supporting consultancy of Oliver Damm [University of Bielefeld, Germany].
- Klein S, Schöneberg I, Krause G: Vom Zwang zur Pockenschutzimpfung zum Nationalen Impfplan. Bundesgesundheitsbl. 2012, 55 (11–12): 1512-1523.View ArticleGoogle Scholar
- Wutzler P, Färber I, Wagenpfeil S, Bisanz H, Tischer A: Seroprevalence of varicella-zoster virus in the German population. Vaccine. 2001, 20 (1–2): 121-124.View ArticlePubMedGoogle Scholar
- Breuer J, Whitley R: Varicella zoster virus: natural history and current therapies of varicella and herpes zoster. Herpes. 2007, 14 (Suppl 2): 25-29.PubMedGoogle Scholar
- Hope-Simpson RE: The nature of Herpes Zoster: a long-term study and a new hypothesis. Proc R Soc Med. 1965, 58: 9-20.PubMedPubMed CentralGoogle Scholar
- Wittek M, Doerr HW, Allwinn R: Varicella and herpes zoster. Part 1: virology, epidemiology, clinical picture, laboratory diagnostics. Med Klin (Munich). 2010, 105 (5): 334-338.View ArticleGoogle Scholar
- Schmader K, Gnann JW, Watson CP: The epidemiological, clinical, and pathological rationale for the herpes zoster vaccine. J Infect Dis. 2008, 197 (Suppl 2): S207-S215.View ArticlePubMedGoogle Scholar
- Drolet M, Brisson M, Levin MJ, Schmader KE, Oxman MN, Johnson RW, Camden S, Mansi JA: A prospective study of the herpes zoster severity of illness. Clin J Pain. 2010, 26 (8): 656-666.PubMedGoogle Scholar
- Opstelten W, Eekhof J, Neven AK, Verheij T: Treatment of herpes zoster. Can Fam Physician. 2008, 54 (3): 373-377.PubMedPubMed CentralGoogle Scholar
- Forbes H, Thomas S, Langan S: The Epidemiology and Prevention of Herpes Zoster. Curr Dermat Rep. 2012, 1 (1): 39-47.View ArticleGoogle Scholar
- Hope-Simpson RE: Postherpetic neuralgia. J R Coll Gen Pract. 1975, 25 (157): 571-575.PubMedPubMed CentralGoogle Scholar
- Jericho BG: Postherpetic neuralgia: a review. Internet Journal of Pain, Symptom Control and Palliative Care. 2010, 8 (1).Google Scholar
- Dworkin RH, Portenoy RK: Pain and its persistence in herpes zoster. Pain. 1996, 67 (2–3): 241-251.View ArticlePubMedGoogle Scholar
- Oster G, Harding G, Dukes E, Edelsberg J, Cleary PD: Pain, medication use, and health-related quality of life in older persons with postherpetic neuralgia: results from a population-based survey. J Pain. 2005, 6 (6): 356-363.View ArticlePubMedGoogle Scholar
- Drolet M, Brisson M, Schmader KE, Levin MJ, Johnson R, Oxman MN, Patrick D, Blanchette C, Mansi JA: The impact of herpes zoster and postherpetic neuralgia on health-related quality of life: a prospective study. CMAJ. 2010, 182 (16): 1731-1736.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmader KE, Sloane R, Pieper C, Coplan PM, Nikas A, Saddier P, Chan IS, Choo P, Levin MJ, Johnson G, et al: The impact of acute herpes zoster pain and discomfort on functional status and quality of life in older adults. Clin J Pain. 2007, 23 (6): 490-496.View ArticlePubMedGoogle Scholar
- Katz J, Cooper EM, Walther RR, Sweeney EW, Dworkin RH: Acute pain in herpes zoster and its impact on health-related quality of life. Clin Infect Dis. 2004, 39 (3): 342-348.View ArticlePubMedGoogle Scholar
- Edmunds WJ, Brisson M, Rose JD: The epidemiology of herpes zoster and potential cost-effectiveness of vaccination in England and Wales. Vaccine. 2001, 19 (23–24): 3076-3090.View ArticlePubMedGoogle Scholar
- Miller E, Marshall R, Vurdien J: Epidemiology, outcome and control of varicella-zoster infection. Rev Med Microbiol. 1993, 4 (4): 222-230.View ArticleGoogle Scholar
- Brisson M, Edmunds WJ, Law B, Gay NJ, Walld R, Brownell M, Roos L, De Serres G: Epidemiology of varicella zoster virus infection in Canada and the United Kingdom. Epidemiol Infect. 2001, 127 (2): 305-314.View ArticlePubMedPubMed CentralGoogle Scholar
- Donahue JG, Choo PW, Manson JE, Platt R: The incidence of herpes zoster. Arch Intern Med. 1995, 155 (15): 1605-1609.View ArticlePubMedGoogle Scholar
- Ultsch B, Köster I, Reinhold T, Siedler A, Krause G, Icks A, Schubert I, Wichmann O: Epidemiology and cost of herpes zoster and postherpetic neuralgia in Germany. Eur J Health Econ. 2012, 1-12.Google Scholar
- Ultsch B, Siedler A, Rieck T, Reinhold T, Krause G, Wichmann O: Herpes zoster in Germany: quantifying the burden of disease. BMC Infect Dis. 2011, 11 (1): 173.View ArticlePubMedPubMed CentralGoogle Scholar
- Schiffner-Rohe J, Jow S, Lilie HM, Köster I, Schubert I: Herpes zoster in Deutschland - Eine retrospektive Analyse von GKV-Daten. MMW Fortschr Med Originalien. 2009, 151: 193-197.Google Scholar
- Volpi A, Gross G, Hercogova J, Johnson RW: Current management of herpes zoster: the European view. Am J Clin Dermatol. 2005, 6 (5): 317-325.View ArticlePubMedGoogle Scholar
- Johnson RW, Wasner G, Saddier P, Baron R: Herpes zoster and postherpetic neuralgia: optimizing management in the elderly patient. Drugs Aging. 2008, 25 (12): 991-1006.View ArticlePubMedGoogle Scholar
- Alper BS, Lewis PR: Does treatment of acute herpes zoster prevent or shorten postherpetic neuralgia?. J Fam Pract. 2000, 49 (3): 255-264.PubMedGoogle Scholar
- Wittek M, Doerr HW, Allwinn R: Varicella and herpes zoster. Part 2: therapy and prevention. Med Klin (Munich). 2010, 105 (6): 399-403.View ArticleGoogle Scholar
- Gross G, Schofer H, Wassilew S, Friese K, Timm A, Guthoff R, Pau HW, Malin JP, Wutzler P, Doerr HW: Herpes zoster guideline of the German Dermatology Society (DDG). J Clin Virol. 2003, 26 (3): 277-289. discussion 291–273View ArticlePubMedGoogle Scholar
- Oxman MN, Levin MJ, Johnson GR, Schmader KE, Straus SE, Gelb LD, Arbeit RD, Simberkoff MS, Gershon AA, Davis LE, et al: A vaccine to prevent herpes zoster and postherpetic neuralgia in older adults. N Engl J Med. 2005, 352 (22): 2271-2284.View ArticlePubMedGoogle Scholar
- Schmader KE, Levin MJ, Gnann JW, McNeil SA, Vesikari T, Betts RF, Keay S, Stek JE, Bundick ND, Su S-C: Efficacy, safety, and tolerability of herpes zoster vaccine in persons aged 50–59 years. Clin Infect Dis. 2012, 54 (7): 922-928.View ArticlePubMedPubMed CentralGoogle Scholar
- Study to Evaluate GSK Biologicals' Herpes Zoster Vaccine GSK1437173A in Adults Aged >= 50 Years [ClinicalTrials.gov Identifier NCT01165177] last update February 7, 2013. 2013, http://clinicaltrials.gov/ct2/show/NCT01165177?term=herpes+zoster&spons=GlaxosmithKline&age=12&phase=2&rank=1.
- Drummond M, Chevat C, Lothgren M: Do we fully understand the economic value of vaccines?. Vaccine. 2007, 25 (32): 5945-5957.View ArticlePubMedGoogle Scholar
- Gesetzliche Krankenversicherung Mitglieder,mitversicherte Angehörige und Krankenstand Monatswerte Januar-Mai 2013. 2013, http://bundesgesundheitsministerium.de/fileadmin/dateien/Downloads/Statistiken/GKV/Mitglieder_Versicherte/KM1_Januar_Mai_2013.pdf.
- Annemans L, Bresse X, Gobbo C, Papageorgiou M: Health economic evaluation of a vaccine for the prevention of herpes zoster (shingles) and post-herpetic neuralgia in adults in Belgium. J Med Econ. 2010, 13 (3): 537-551.View ArticlePubMedGoogle Scholar
- Brisson M: Estimating the number needed to vaccinate to prevent herpes zoster-related disease, health care resource use and mortality. Can J Public Health. 2008, 99 (5): 383-386.PubMedGoogle Scholar
- Gauthier A, Breuer J, Carrington D, Martin M, Remy V: Epidemiology and cost of herpes zoster and post-herpetic neuralgia in the United Kingdom. Epidemiol Infect. 2009, 137 (1): 38-47.View ArticlePubMedGoogle Scholar
- Gonzalez Chiappe S, Sarazin M, Turbelin C, Lasserre A, Pelat C, Bonmarin I, Chosidow O, Blanchon T, Hanslik T: Herpes zoster: burden of disease in France. Vaccine. 2010, 28 (50): 7933-7938.View ArticlePubMedGoogle Scholar
- Max M, Schafer S, Culnane M, Smoller B, Dubner R, Gracely R: Amitriptyline, but not lorazepam, relieves postherpetic neuralgia. Neurology. 1988, 38 (9): 1427.View ArticlePubMedGoogle Scholar
- Opstelten W, Mauritz JW, de Wit NJ, van Wijck AJ, Stalman WA, van Essen GA: Herpes zoster and postherpetic neuralgia: incidence and risk indicators using a general practice research database. Fam Pract. 2002, 19 (5): 471-475.View ArticlePubMedGoogle Scholar
- van Hoek AJ, Gay N, Melegaro A, Opstelten W, Edmunds WJ: Estimating the cost-effectiveness of vaccination against herpes zoster in England and Wales. Vaccine. 2009, 27 (9): 1454-1467.View ArticlePubMedGoogle Scholar
- Yawn BP: Post-shingles neuralgia by any definition is painful, but is it PHN?. Mayo Clin Proc. 2011, 86 (12): 1141-1142.View ArticlePubMedPubMed CentralGoogle Scholar
- General Methods for Assessment of the Relation of Benefits to Costs. https://www.iqwig.de/download/General_Methods_for_the_Assessment_of_the_Relation_of_Benefits_to_Costs.pdf.
- Cost estimation - working paper, Version 1.0. https://www.iqwig.de/download/Working_Paper_Cost_Estimation.pdf.
- Gnann JW, Whitley RJ: Clinical practice. Herpes zoster. N Engl J Med. 2002, 347 (5): 340-346.View ArticlePubMedGoogle Scholar
- Ragozzino MW, Melton LJ, Kurland LT, Chu CP, Perry HO: Population-based study of herpes zoster and its sequelae. Medicine (Baltimore). 1982, 61 (5): 310-316.View ArticleGoogle Scholar
- Tseng HF, Chi M, Smith N, Marcy SM, Sy LS, Jacobsen SJ: Herpes zoster Vaccine and the incidence of Recurrent Herpes Zoster in an immunocompetent elderly population. J Infect Dis. 2012, 55 (10): 1320-1328.Google Scholar
- Yawn BP, Wollan PC, Kurland MJ, St Sauver JL, Saddier P: Herpes zoster recurrences more frequent than previously reported. Mayo Clin Proc. 2011, 86 (2): 88-93.View ArticlePubMedPubMed CentralGoogle Scholar
- Bevölkerung und Erwerbstätigkeit - Sterbetafel Deutschland. https://www.destatis.de/DE/ZahlenFakten/GesellschaftStaat/Bevoelkerung/Sterbefaelle/Tabellen/SterbetafelDeutschland.xls?__blob=publicationFile.
- Gagliardi AM, Gomes Silva BN, Torloni MR, Soares BG: Vaccines for preventing herpes zoster in older adults. Cochrane Database Syst Rev. 2012, 10: CD008858.PubMedGoogle Scholar
- Oxman MN, Levin MJ: Vaccination against Herpes Zoster and Postherpetic Neuralgia. J Infect Dis. 2008, 197 (Suppl 2): S228-S236.View ArticlePubMedPubMed CentralGoogle Scholar
- FDA: FDA Clinical Briefing Document for Merck & Co., Inc. Zoster vaccine live (Oka/Merck) Zostavax™. Vaccines and Related Biological Products Advisory Committee. Edited by: US-Food-and-Drug-Administration. 2005Google Scholar
- Schmader K, Oxman M, Levin M, Betts R, Morrison V, Gelb L, Johnson G, Zhang J, Li X, Annunziato P: Persistence of Zoster Vaccine Efficacy. ICAAC/IDSA: 2008; Durham. 2008Google Scholar
- Schmader KE, Oxman MN, Levin MJ, Johnson G, Zhang JH, Betts R, Morrison VA, Gelb L, Guatelli JC, Harbecke R, et al: Persistence of the Efficacy of Zoster Vaccine in the Shingles prevention study and the short-term persistence substudy. Clin Infect Dis. 2012Google Scholar
- Bilcke J, Marais C, Ogunjimi B, Willem L, Hens N, Beutels P: Cost-effectiveness of vaccination against herpes zoster in adults aged over 60 years in Belgium. Vaccine. 2012, 30 (3): 675-684.View ArticlePubMedGoogle Scholar
- Bilcke J, Ogunjimi B, Hulstaert F, Van Damme P, Hens N, Beutels P: Estimating the age-specific duration of herpes zoster vaccine protection: a matter of model choice?. Vaccine. 2012, 30 (17): 2795-2800.View ArticlePubMedGoogle Scholar
- Pellissier JM, Brisson M, Levin MJ: Evaluation of the cost-effectiveness in the United States of a vaccine to prevent herpes zoster and postherpetic neuralgia in older adults. Vaccine. 2007, 25 (49): 8326-8337.View ArticlePubMedGoogle Scholar
- MMWR: Noninfluenza Vaccination Coverage Among Adults — United States. 2011, http://www.cdc.gov/mmwr/pdf/wk/mm62e0129.pdf.Google Scholar
- Epidemiologisches Bulletin. https://www.rki.de/DE/Content/Infekt/EpidBull/Archiv/2002/Ausgabenlinks/16_02.pdf?__blob=publicationFile.
- Lauer-Taxe online - Preisinfo: ZOSTAVAX. 2010, http://www2.lauer-fischer.de/home/.
- Briggs A, Claxton K, Sculpher M: Decision modelling for health economic evaluation. 2008, Oxford: Oxford University PressGoogle Scholar
- Hinz A, Klaiberg A, Brahler E, Konig HH: [The quality of life questionnaire EQ-5D: modelling and norm values for the general population]. Psychother Psychosom Med Psychol. 2006, 56 (2): 42-48.View ArticlePubMedGoogle Scholar
- Claxton K, Sculpher M, McCabe C, Briggs A, Akehurst R, Buxton M, Brazier J, O'Hagan T: Probabilistic sensitivity analysis for NICE technology assessment: not an optional extra. Health Econ. 2005, 14 (4): 339-347.View ArticlePubMedGoogle Scholar
- Welte R, Feenstra T, Jager H, Leidl R: A decision chart for assessing and improving the transferability of economic evaluation results between countries. Pharmacoeconomics. 2004, 22 (13): 857-876.View ArticlePubMedGoogle Scholar
- Jacobzone S: Pharmaceutical Policies in OECD Countries. Reconciling Social and Industrial Goals. Edited by: 40 OLMaSPOPN. 2000, Paris: OECDGoogle Scholar
- Mick G, Hans G: Postherpetic neuralgia in Europe: The scale of the problem and outlook for the future. J Clin Gerontol Geriatr. 2013Google Scholar
- Reda H, Greene K, Rice FL, Rowbotham MC, Petersen KL: Natural History of Herpes Zoster: Late follow-up 3.9 yrs (n= 43) and 7.7 yrs (n= 10). PAIN. 2013, 10: 2227-2233.View ArticleGoogle Scholar
- van Lier A, van Hoek AJ, Opstelten W, Boot H, de Melker H: Assessing the potential effects and cost-effectiveness of programmatic herpes zoster vaccination of elderly in the Netherlands. BMC Health Serv Res. 2010, 10 (1): 237.View ArticlePubMedPubMed CentralGoogle Scholar
- De Boer PT, Pouwels KB, Cox JM, Hak E, Wilschut JC, Postma MJ: Cost-effectiveness of vaccination of the elderly against herpes zoster in The Netherlands. Vaccine. 2013, 31 (9): 1276-1283.View ArticlePubMedGoogle Scholar
- Szucs TD, Kressig RW, Papageorgiou M, Kempf W, Michel JP, Fendl A, Bresse X: Economic evaluation of a vaccine for the prevention of herpes zoster and post-herpetic neuralgia in older adults in Switzerland. Hum Vaccin. 2011, 7 (7): 749-756.View ArticlePubMedPubMed CentralGoogle Scholar
- Bresse X, Annemans L, Preaud E, Bloch K, Duru G, Gauthier A: Vaccination against herpes zoster and postherpetic neuralgia in France: a cost-effectiveness analysis. Expert Rev Pharmacoecon Outcomes Res. 2013, 13 (3): 393-406.View ArticlePubMedGoogle Scholar
- Moore L, Remy V, Martin M, Beillat M, McGuire A: A health economic model for evaluating a vaccine for the prevention of herpes zoster and post-herpetic neuralgia in the UK. Cost Effectiveness and Resource Allocation. 2010, 8 (1): 7.View ArticlePubMedPubMed CentralGoogle Scholar
- Hornberger J, Robertus K: Cost-effectiveness of a vaccine to prevent herpes zoster and postherpetic neuralgia in older adults. Ann Intern Med. 2006, 145 (5): 317-325.View ArticlePubMedGoogle Scholar
- Rothberg MB, Virapongse A, Smith KJ: Cost-effectiveness of a vaccine to prevent herpes zoster and postherpetic neuralgia in older adults. Clin Infect Dis. 2007, 44 (10): 1280-1288.View ArticlePubMedGoogle Scholar
- Brisson M, Pellissier JM, Camden S, Quach C, De Wals P: The potential cost-effectiveness of vaccination against herpes zoster and post-herpetic neuralgia. Hum Vaccin. 2008, 4 (3): 238-245.View ArticlePubMedGoogle Scholar
- Najafzadeh M, Marra CA, Galanis E, Patrick DM: Cost effectiveness of herpes zoster vaccine in Canada. Pharmacoeconomics. 2009, 27 (12): 991-1004.View ArticlePubMedGoogle Scholar
- Szucs TD, Pfeil AM: A systematic review of the cost effectiveness of herpes zoster vaccination. Pharmacoeconomics. 2013, 31 (2): 125-136.View ArticlePubMedGoogle Scholar
- Coplan PM, Schmader K, Nikas A, Chan IS, Choo P, Levin MJ, Johnson G, Bauer M, Williams HM, Kaplan KM, et al: Development of a measure of the burden of pain due to herpes zoster and postherpetic neuralgia for prevention trials: adaptation of the brief pain inventory. J Pain. 2004, 5 (6): 344-356.View ArticlePubMedGoogle Scholar
- Bala MV, Wood LL, Zarkin GA, Norton EC, Gafni A, O'Brien B: Valuing outcomes in health care: a comparison of willingness to pay and quality-adjusted life-years. Journal of clinical epidemiology. 1998, 51 (8): 667-676.View ArticlePubMedGoogle Scholar
- Graf Von Der Schulenburg JM, Claes C, Greiner W, et al: Die deutsche Version des EuroQol-Fragebogens. Zeitschrift für Gesundheitswissenschaften. 1998, 6: 3-20.Google Scholar
- Hoebel J, Kuntz B, Muters S, Lampert T: Subjective Social Status and Health-Related Quality of Life among Adults in Germany. Gesundheitswesen. 2013Google Scholar
- Konig HH, Bernert S, Angermeyer MC: Health status of the German population: results of a representative survey using the EuroQol questionnaire. Gesundheitswesen. 2005, 67 (3): 173-182.View ArticlePubMedGoogle Scholar
- Mielck A, Vogelmann M, Schweikert B, Leidl R: Health status of adults in Germany: results from a representative survey using the EuroQol 5D (EQ-5D). Gesundheitswesen. 2010, 72 (8–9): 476-486.View ArticlePubMedGoogle Scholar
- Pierik J, Gumbs P, Fortanier S, Van Steenwijk P, Postma M: Epidemiological characteristics and societal burden of varicella zoster virus in the Netherlands. BMC Infect Dis. 2012, 12 (1): 110.View ArticlePubMedPubMed CentralGoogle Scholar
- Gialloreti L, Merito M, Pezzotti P, Naldi L, Gatti A, Beillat M, Serradell L, di Marzo R, Volpi A: Epidemiology and economic burden of herpes zoster and post-herpetic neuralgia in Italy: a retrospective, population-based study. BMC Infect Dis. 2010, 10 (1): 230.View ArticlePubMedPubMed CentralGoogle Scholar
- Cebrian-Cuenca A, Diez-Domingo J, San-Martin-Rodriguez M, Puig-Barbera J, Navarro-Perez J: Epidemiology and cost of herpes zoster and postherpetic neuralgia among patients treated in primary care centres in the Valencian community of Spain. BMC Infect Dis. 2011, 11 (1): 302.View ArticlePubMedPubMed CentralGoogle Scholar
- Vesikari T, Hardt R, Rümke HC, Icardi G, Montero J, Thomas S, Sadorge C, Fiquet A: Immunogenicity and safety of a live attenuated shingles (herpes zoster) vaccine (Zostavax®) in individuals aged ≥ 70 years: a randomized study of a single dose vs. two different two-dose schedules. Human Vaccines & Immunotherapeutics. 2013, 9 (4): 854-860.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6963/13/359/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 cited.