Hostname: page-component-7c8c6479df-94d59 Total loading time: 0 Render date: 2024-03-19T10:19:04.332Z Has data issue: false hasContentIssue false

Cost-utility analysis of antiviral use under pandemic influenza using a novel approach – linking pharmacology, epidemiology and heath economics

Published online by Cambridge University Press:  15 February 2018

D. B. C. Wu
Affiliation:
School of Pharmacy, Monash University Malaysia, Bandar Sunway, Malaysia Asian Centre for Evidence Synthesis in Population, Implementation and Clinical Outcomes (PICO), Health and Well-being Cluster, Global Asia in the 21st Century (GA21) Platform, Monash University Malaysia, Bandar Sunway, Malaysia
N. Chaiyakunapruk*
Affiliation:
School of Pharmacy, Monash University Malaysia, Bandar Sunway, Malaysia Asian Centre for Evidence Synthesis in Population, Implementation and Clinical Outcomes (PICO), Health and Well-being Cluster, Global Asia in the 21st Century (GA21) Platform, Monash University Malaysia, Bandar Sunway, Malaysia Department of Pharmacy Practice, Faculty of Pharmaceutical Sciences, Naresuan University, Center of Pharmaceutical Outcomes Research (CPOR), Phitsanulok, Thailand School of Pharmacy, University of Wisconsin, Madison, USA
C. Pratoomsoot
Affiliation:
Faculty of Public Health, Naresuan University, Phitsanulok, Thailand
K. K. C. Lee
Affiliation:
School of Pharmacy, Monash University Malaysia, Bandar Sunway, Malaysia
H. Y. Chong
Affiliation:
School of Pharmacy, Monash University Malaysia, Bandar Sunway, Malaysia
R. E. Nelson
Affiliation:
University of Utah, Salt Lake City, UT, USA VA Salt Lake City Health Care System, Salt Lake City, UT, USA
P. F. Smith
Affiliation:
d3 Medicine, a Certara Company, Parsippany, New Jersey, USA
C.M. Kirkpatrick
Affiliation:
Centre for Medicine Use and Safety, Faculty of Pharmacy and Pharmaceutical Sciences, Monash University, Melbourne, Australia
M. A. Kamal
Affiliation:
Department of Clinical Pharmacology, Regeneron Pharmaceuticals, Tarrytown, NY, USA
K. Nieforth
Affiliation:
d3 Medicine, a Certara Company, Parsippany, New Jersey, USA
G. Dall
Affiliation:
d3 Medicine, a Certara Company, Parsippany, New Jersey, USA
S. Toovey
Affiliation:
Pegasus Research, Basel, Switzerland
D. C. M. Kong
Affiliation:
Centre for Medicine Use and Safety, Faculty of Pharmacy and Pharmaceutical Sciences, Monash University, Melbourne, Australia
A. Kamauu
Affiliation:
Anonlinx LLC, Salt Lake City, USA
C. R. Rayner
Affiliation:
d3 Medicine, a Certara Company, Parsippany, New Jersey, USA
*
Author for correspondence: N. Chaiyakunapruk, E-mail: nathorn.chaiyakunapruk@monash.edu
Rights & Permissions [Opens in a new window]

Abstract

Simulation models are used widely in pharmacology, epidemiology and health economics (HEs). However, there have been no attempts to incorporate models from these disciplines into a single integrated model. Accordingly, we explored this linkage to evaluate the epidemiological and economic impact of oseltamivir dose optimisation in supporting pandemic influenza planning in the USA. An HE decision analytic model was linked to a pharmacokinetic/pharmacodynamics (PK/PD) – dynamic transmission model simulating the impact of pandemic influenza with low virulence and low transmissibility and, high virulence and high transmissibility. The cost-utility analysis was from the payer and societal perspectives, comparing oseltamivir 75 and 150 mg twice daily (BID) to no treatment over a 1-year time horizon. Model parameters were derived from published studies. Outcomes were measured as cost per quality-adjusted life year (QALY) gained. Sensitivity analyses were performed to examine the integrated model's robustness. Under both pandemic scenarios, compared to no treatment, the use of oseltamivir 75 or 150 mg BID led to a significant reduction of influenza episodes and influenza-related deaths, translating to substantial savings of QALYs. Overall drug costs were offset by the reduction of both direct and indirect costs, making these two interventions cost-saving from both perspectives. The results were sensitive to the proportion of inpatient presentation at the emergency visit and patients’ quality of life. Integrating PK/PD–EPI/HE models is achievable. Whilst further refinement of this novel linkage model to more closely mimic the reality is needed, the current study has generated useful insights to support influenza pandemic planning.

Type
Original Papers
Copyright
Copyright © Cambridge University Press 2018 

Introduction

Simulation models are widely used in pharmacology, epidemiology and health economics (HEs). These models can explicitly link concepts and ideas to data in order to produce outcomes that are useful to healthcare decision makers [Reference Weinstein1]. Importantly, there have been no attempts to incorporate the models from pharmacology, epidemiology and HEs into a single integrated model.

While still performed infrequently, the importance of epidemiological modelling to inform economic analyses of prevention and treatment strategies for infectious diseases is well accepted [Reference Jit and Brisson2]. Cost-effectiveness analyses (CEAs) are typically conducted using static models in which the probability of disease exposure is independent of interventions designed to treat or prevent it. While this is realistic for non-transmissible diseases, it is not suited for communicable diseases given interventions in such circumstance can have a significant impact on disease exposure. For example, vaccines can reduce the susceptible population and antiviral treatments can reduce viral shedding, thus reducing viral transmission. Accordingly, dynamic transmission models are needed when there is interdependence between intervention and disease exposure [Reference Lugner, Mylius and Wallinga3]. These models allow the indirect or herd immunity effect to be predicted. Indeed, static and dynamic models produced dramatically different results [Reference Brisson and Edmunds4, Reference Brisson and Edmunds5]. Importantly, a recent review on economic evaluations of interventions for pandemic influenza planning reported most of the studies have used static models [Reference Lugner and Postma6].

In pharmacology, population pharmacokinetic/pharmacodynamic (PK/PD) models are used to identify optimal dosing regimens [Reference Duffull, Wright and Winter7] and were used recently to demonstrate that elevated doses of oseltamivir may result in reduced viral shedding [Reference Rayner8]. This has important implications for the economic evaluation of influenza pandemic planning as influenza viral transmission is a function of viral shedding [Reference Bulik9]. These indirect effects, which can be larger than the direct effects [Reference Pradas-Velasco, Antonanzas-Villar and Martinez-Zarate10], are often neglected in cost-effectiveness analyses. Any economic analysis failing to consider these indirect effects will have inaccuracies [Reference Pitman11].

Globally, the emergence or re-emergence of influenza outbreaks has been challenging [Reference Girard1214]. Several economic evaluation studies have been conducted to facilitate the management of influenza outbreaks [Reference Kelso15, Reference Milne, Halder and Kelso16]. A recent systematic review reported that most economic studies used models that were incapable of incorporating variability of all parameters and did not consider the viral transmission process [Reference Pradas-Velasco, Antonanzas-Villar and Martinez-Zarate10]. Most importantly, existing models did not consider the indirect effects which are crucial in a pandemic mitigation strategy.

The current study explored the feasibility of linking pharmacology, epidemiology and HE models within the realm of anti-influenza treatment for pandemic influenza. We previously published an overall framework for the linking without a full description of HE models [Reference Kamal17]. To our knowledge, this is the first CEA to directly integrate data from PK/PD and disease transmission models into the economic analysis.

Methods

Linking of pharmacology and epidemiology to HEs models

Previously, a single semi-mechanistic framework assessed quantitatively the impact of oseltamivir pharmacology and treatment approaches on the burden of influenza infection in a hypothetical population of 100 000 individuals across a 1-year flu season [Reference Kamal17]. This model connected three discrete quantitative modules viz. a population pharmacokinetic model for oseltamivir; a PK/PD evaluation of oseltamivir carboxylate (OC, active metabolite) on viral shedding; and a susceptible-exposed-infected-recovered (SEIR) epidemiologic model, adapted to incorporate the impact of antiviral therapy.

Briefly, OC PK variability was incorporated using a published population PK model [Reference Kamal18] for OC incorporating 390 healthy and infected subjects ranging from 1–78 years across a dose range of 20–1000 mg. The final covariate model from this population PK analysis included allometric scaling and a relationship between weight on the OC central volume of distribution and between weight and creatinine clearance on OC clearance. From this final covariate model, OC PK parameter profiles in 5000 70-kg adult patients aged 18–65 with normal renal function receiving the standard 75 mg twice daily (BID) and 150 mg BID oseltamivir for 5 days was simulated. The area under the concentration-time curve (AUCs) for each patient were quantified and the fraction of patients with OC AUCs above a published PK/PD value associated with faster cessation of viral shedding (14 180 ng.h/ml) [Reference Rayner8] was calculated for each dosing regimen; and denoted as F AUChigh.

PD variability was integrated via incorporating viral shedding (T shed) distributions for placebo subjects, and those above and below the PK/PD target as obtained from literature [Reference Rayner8]. These PK/PD distributions and other SEIR input parameters (Table 1), allowed the average number of infected cases per 100 000 population to be determined for a range of scenarios relating to dose, viral transmissibility and treatment uptake (vide infra). The important innovation in the current work is the extension of the aforementioned quantitative framework to a decision analytic model to estimate the cost and outcomes of interventions of interest (Fig. 1). Both virus transmissibility and impact of anti-influenza agents alter the number of infected cases entering into the decision analytic model. As previously described, ß governs transmissibility, and is a composite of both frequencies of individual interactions, and the probability that an interaction will result in a successful influenza infection in a susceptible individual [Reference Kamal17]. The impact of antiviral treatment is via reducing the time an infected subject is infectious; governed by the duration of viral shedding. As an individual disease model is currently unavailable (Fig. 1), the severity of illness associated with influenza virulence is assumed to directly alter the probabilities of complications and associated costs in the decision analytic model.

Fig. 1. Overarching Pharmacology to payer system including ‘modules’. The solid lines indicate that adequate data exists to be able to create semi-mechanistic links to each of adjacent ‘modules’. The dotted lines and light grey describe where significant unknowns remain and are not mature enough to have been incorporated into the current framework. PopPK, population pharmacokinetic; SEIR, susceptible-exposed-infected-recovered; ß, the rate of infectivity; TShed, viral shedding; VK, viral kinetics.

Table 1. Parameters used in the SEIR model

κ, the delay rate between exposure to influenza and symptom development; F AUChigh, the mean (SD) fraction of the simulated population receiving oseltamivir with an AUC > 14 180 ng.h/ml, obtained from the pharmacology module; γ0, the duration of viral shedding under no treatment; γlow, the mean (SD) duration of viral shedding if OC AUC is <14 180 ng.h/ml; γhigh, the mean (SD) duration of viral shedding if OC AUC is >14 180 ng.h/ml, ß, the rate of infectivity; AUC, area under the concentration–time curve; BID, twice daily.

Health economics model structure and data inputs

The decision analytic model (Fig. 2) was developed to evaluate the cost and outcomes of oseltamivir standard (75 mg BID) or high (150 mg BID) dose or no treatment under pandemic influenza scenarios. The number of infected patients entering into the model was simulated from the SEIR model. A cost-utility analysis was undertaken based on the US population of healthy adults, aged 18–64 years old from the payer and societal perspective. Societal perspective included direct and indirect costs. Total costs, quality-adjusted life years (QALYs), and incremental cost-effectiveness ratios (ICERs) were determined over a 1-year time horizon.

Fig. 2. Decision analytic model. Influenza patients entered the decision analytic model from epidemiology model. They received treatment in outpatient or inpatient setting. OPD, outpatient; ED, emergency department; AVR, antiviral; GW, general ward; ICU, intensive care unit; ARDS, acute respiratory distress syndrome.

The infected individual entered the HE model either as an outpatient or inpatient (Fig. 2). Those admitted as inpatient would be admitted to a general ward or an intensive care unit, may experience either pneumonia, sepsis and acute respiratory distress syndrome (ARDS) [Reference Jain19Reference Skarbinski21]. We assumed that patients could only experience one influenza-related complication. The infected patient either recovered from the infection or died. We assumed that oseltamivir was prescribed within 48 h of influenza symptoms and all patients were 100% adherent to treatment received.

A conservative approach involving oseltamivir acting only by reducing the duration of viral shedding and not directly on transmission rate was adopted. As such, only drug treatment effects and not prophylaxis effects were captured in the model. We assumed that oseltamivir reduces the time of symptom alleviation by 21 h and conservatively has no effects on complication, hospitalisation and mortality rate [Reference Burch22Reference Michiels25]. Data about influenza disease progression, medical resource utilisation, cost of treatment and health state utilities were obtained from the literature (Medline searches) and Healthcare Cost and Utilization Project Nationwide Inpatient Sample database [26]. Where possible, data related to the 2009 pandemic H1N1 influenza was used. Meta-analyses were performed when applicable. All costs were converted to 2013 US dollars using the Consumer Price Index [27]. Base-case estimates and ranges for probabilities, costs, utilities and length of stay for hospitalisation are shown in Table 2; detailed descriptions available in Appendix 1.

Table 2. Input parameters, values and data sources used in the health economics model

CI, confidence interval; SE, standard error; ED, emergency department; GP, general practitioner; GW, general ward; ICU, intensive care unit; ARDS, acute respiratory distress syndrome; USD, United States Dollar; AWP, average wholesale price.

a Oseltamivir is available as 75 mg/capsule and in a pack of 10 capsules in the USA. The AWP for one pack of oseltamivir is USD 132.77.

b Weighted probability or cost estimate was calculated using probability or cost reported for two age groups, i.e. 18–49 and 50–64 years, assuming the age-distribution of the influenza adults is similar to that of the U.S. population in 2012 [50].

c Weighted cost or length of stay estimate was calculated using cost or length of stay reported for two age groups, i.e. 18–44 and 45–64 years.

d Reported as standard deviation.

e Calculated by multiplying 5 miles round trip and the standard business mileage rate 2013 for the use of a car.

f Calculated using the age- and sex-stratified median weekly earnings in the U.S. Bureau of Labor Statistics for the age group 18–64 years, taking account of the unemployment rate of 7.3%.

g Calculated using the length of stay reported and adding one additional day for convalescence [Reference Meltzer, Cox and Fukuda40].

h Khazeni et al. [Reference Khazeni41] adjusted for possible on-going disability following recovery from severe (all patients who received inpatient ICU care) influenza, which assumed an estimated quality of life of 0.9 for the remainder of those patients’ lifetimes. In our HE model, we assumed the value of 0.90 for the rest of the year for those patients.

Simulation scenarios and treatment comparison

We investigated the HE impact of two extreme hypothetical pandemic scenarios, namely, influenza virus with low virulence and low transmissibility and virus with high virulence and high transmissibility [Reference Russell51, Reference Kawaoka52]. The ß values for transmissibility were from Kamal et al. [Reference Kamal17], (see Table 1) whereas virulence was based on health care utilisation; low virulence was based on pH1N1 experience [Reference Girard12]; and the high virulence scenario involved doubling the probability of hospitalisation [Reference Lee53] for the low virulence scenario within the decision analytic model.

We evaluated three different management strategies: (1) no treatment, (2) the standard approved oseltamivir dose (75 mg BID) and (3) high dose (150 mg BID). We also evaluated the treatment strategies according to varying treatment uptake rates (25, 50 and 80%), defined as the percentage of infected subjects treated, across an infected population. Detailed descriptions of all values are available in Appendix 1. Analyses were performed using Berkeley Madonna™ version 8.3.18, R version 2.15.3 and Microsoft Excel 2010 (Redmond, Washington).

Sensitivity analyses

One-way sensitivity analyses were performed to investigate the effects of altering parameters within the plausible ranges (Table 2). Additional sensitivity analysis was undertaken to determine the robustness of estimates when a different approach for productivity loss estimation was adopted. In the base-case analyses, we adopted the approach in Meltzer et al. [Reference Meltzer, Cox and Fukuda40] Alternatively, we calculated productivity loss using the relationship between ‘time to return to normal activity’ and ‘illness duration’ (Appendix 1). We also varied the probability of developing ARDS in highly virulent and highly transmissible pandemic influenza case. Probabilistic sensitivity analysis (PSA) was undertaken to address uncertainty in the assumptions underlying the model by allowing each input parameter values to vary simultaneously over their respective feasible ranges within the model; 5000 second-order Monte Carlo simulations were performed in which the parameter values were drawn from pre-defined distributions. Results from the PSA are presented as scatterplots in cost-effectiveness plane.

Results

Low virulence and low transmissibility

Based on a hypothetical cohort of 100 000 individuals, under low virulence and low transmissibility, standard-dose oseltamivir (75 mg BID) was cost-saving from both perspectives compared to no treatment (Table 3). It could save USD30,246,490 (USD81,272,885), USD33,352,767 (USD93,993,932) and USD30,406,894 (USD92,614,326) for oseltamivir uptakes of 25%, 50% and 80% by averting 706, 412 and 783 deaths, respectively, from payer (societal) perspective. This translates to 361, 430 and 441 QALYs gained for uptakes of 25%, 50% and 80%, respectively.

Table 3. Base-case analyses: high-dose vs. no treatment, and standard dose vs. no treatment

a All costs are expressed in 2013 USD.

High-dose oseltamivir (150 mg BID) when compared to no treatment was also cost-saving from payer (societal) perspective (Table 3), leading to cost-saving of USD29,839,058 (USD85,192,106), USD27,742,305 (USD89,785,170) and USD20,483,671 (USD83,687,273) for uptakes of 25%, 50% and 80% by preventing 376, 423 and 430 deaths, respectively, leading to the associated QALY gains of 392, 441 and 449, respectively.

Compared to low-dose, high-dose oseltamivir would avert −30, −11, −7 deaths at 25%, 50% and 80% uptakes, respectively, leading to the associated QALY gains of 31, 11, 7, respectively. It led to overall increased cost of USD407,432, USD5,610,462, USD9,923,223 from payer perspective for uptakes of 25%, 50%, 80% uptakes, respectively, resulting in USD13,117, USD546,753, USD1,353,188 per QALY gained. From societal perspective, using high-dose oseltamivir could save USD3,919,221 at 25% uptake (cost-saving) but incurred extra cost of USD4,110,702 and USD8,927,052 at 50% and 80% uptakes, leading to USD377,524 and USD1,217,344 QALY gained, respectively.

High virulence and high transmissibility

Under the scenario of high virulence and high transmissibility, compared with no treatment, oseltamivir 75 mg BID was dominant (Table 3). It could avert 168, 617 and 1098 deaths, translating to 174, 598 and 1063 QALYs gained for the uptakes of 25%, 50% and 80%, respectively (Table 3). The standard-dose could save USD11,885,337 (USD27,413,072), USD49,309,677 (USD101,904,191) and USD88,899,525 (USD179,850,525) for 25%, 50% and 80% uptakes from payer (societal) perspective.

When high-dose oseltamivir was compared to no treatment (Table 3), it was shown to save USD24,304,917 (USD53,631,551), USD63,252,396 (USD133,577,886) and USD95,501,955 (USD200,976,888) for uptakes of 25%, 50% and 80%, by preventing 341, 844 and 1288 deaths, respectively, leading to the associated QALY gains of 330, 818 and 1247, respectively.

Compared to low-dose, our model indicated that high-dose oseltamivir would save USD12,419,581 (USD26,218,479), USD13,942,719 (USD31,673,695) and USD6,602,430 (USD21,126,363) from the payer (societal) perspective for uptakes of 25%, 50% and 80% (Table 3). In addition, 174, 227 and 190 deaths could be prevented, leading to QALY gains of 174, 227 and 189, respectively.

Sensitivity analyses

The tornado diagrams (Fig. 3) summarise the results of 1-way sensitivity analyses. Under low virulence and low transmissibility, at 80% oseltamivir uptake, the ICER was most sensitive to the probability of hospitalisation for ED cases and baseline average quality of life. Under the high virulence and high transmissibility scenario, the baseline average quality of life and median weekly earnings in the USA were the two most influential factors. The results were also very similar to the other uptakes and thus not shown.

Fig. 3. Tornado diagram (150 mg vs. no treatment with 80% uptake of oseltamivir): One-way sensitivity analysis under (1a) low virulence and low transmissibility and (1b) high virulence and high transmissibility. (a) Low virulence and low transmissibility. #+/−: The higher value of the parameter leads to lower ICER. $−/+: The higher value of the parameter leads to higher ICER. (b) High virulence and high transmissibility. #+/−: The higher value of the parameter leads to lower ICER.

Monte Carlo simulations (5000) for different levels of viral virulence, transmissibility and oseltamivir uptake from the societal perspective are presented in Fig. 4. In both scenarios, a very large proportion of the simulated results fell in the south-east quadrant, suggesting that 75 mg BID vs. no treatment at 25%, 50% and 80% uptakes were highly likely to be cost-saving. The results of 150 mg BID from a societal perspective were also very similar (data not provided).

Fig. 4. Scatter plots (incremental cost vs. incremental QALY) of 75 mg vs. no treatment under societal perspective. (a) Low virulence and low transmissibility. (b) High virulence and high transmissibility.

Discussion

This is the first study to link pharmacology, epidemiology and HE models into a single integrated model. We examined the costs and outcomes of different levels of oseltamivir dosing for pandemic influenza planning using the information from the PK/PD model to ultimately inform an HE model by way of an SEIR model. While we applied this novel linked model within the context of pandemic influenza, its usefulness is not limited to this arena. Our approach can also be useful for policy-oriented decision making surrounding the use of other antimicrobials agents.

Indeed, our results show that using standard- or high-dose oseltamivir will afford a significant reduction in deaths due to three major influenza-related complications over a 1-year time horizon, thus leading to substantial savings in both life years and QALYs under various degrees of virus virulence and transmissibility. The extra cost for oseltamivir was offset by reductions in the number of disease episodes and influenza-related deaths, resulting in substantial savings.

Several economic evaluations on influenza pandemic preparedness strategies have been published, including studies evaluating vaccine, pharmaceutical and non-pharmaceutical (e.g. social distancing) interventions with the majority focusing on anti-influenza drugs [Reference Pradas-Velasco, Antonanzas-Villar and Martinez-Zarate10]. Pharmaceutical interventions combined with non-pharmaceutical interventions have been found to be relatively cost-effective in comparison to vaccines and/or antivirals alone. However, pharmaceutical interventions vary from cost-saving to being not cost-effective [Reference Pradas-Velasco, Antonanzas-Villar and Martinez-Zarate10]. The wide range of cost-effectiveness results is due to several factors such as different dosages and durations for prophylaxis as well as the comparator. Our results showing cost-savings of oseltamivir under all scenarios and dosages were possibly due to the absence of consideration of concomitant treatments and the social contact structure among various subpopulations.

Despite the implementation of the World Health Organization's (WHO's) pandemic preparedness plan and response guide during the 2009 pandemic H1N1, experience and evidence suggest that the global community is ill-prepared to respond to a severe influenza pandemic or to any similarly global, life-threatening public-health emergency [54]. The possibility of the avian H5N1 becoming pandemic has raised concern and has been anticipated by the public health authorities [Reference Fineberg55]. This strain is a significant threat given its high fatality rate in humans, the large reservoir of poultry in South-East Asia and that H5N1 influenza may be genetically modified to become more transmissible [Reference Russell51, Reference Herfst56, Reference Imai57]. Importantly, the global community should be prepared for the next pandemic by having a full understanding of the consequences of influenza pandemics and their related interventions. In the current study, by using the 2009 pandemic data in the linked models, we explored the cost-effectiveness of oseltamivir as a pandemic mitigation strategy in two hypothetical pandemic scenarios with different severity, inspired by the 2009 pandemic H1N1 and H5N1. Our study did not consider the reported case fatality ratio (CFR) of H5N1 by the WHO due to inconclusive findings on the real CFR estimate by different sources [Reference Li58, Reference Cowling59]. Furthermore, the CFR may decrease as more information becomes available or as the epidemic progresses. Several recent studies from the Cochrane collaboration [Reference Jefferson60] and Muthuri et al. [Reference Muthuri61] have called into question the efficacy of oseltamivir. Therefore, we conservatively estimated the drug effect of oseltamivir based solely on viral shedding in our pharmacology component of the integrated model. Subsequently, in the HEs model, assumptions were made that oseltamivir shortened the duration of influenza symptoms by 21–22 h, but had no effect on the progression to secondary complications, hospitalisation or mortality [Reference Burch22Reference Hsu24]. These assumptions may underestimate the true effects of oseltamivir during an influenza pandemic.

Our work has several limitations. First, we did not take into consideration other interventions such as the use of masks, school closure and influenza vaccine. The concurrent use of other interventions may reduce the impact of oseltamivir. Second, we assumed that all infected cases seek medical care, which might not be true in reality. Third, our model was confined to only a healthy population between 18 and 65 years of age. However, the model could easily be extended to include other patient populations. Finally, it is noteworthy that patients’ viral shedding duration for 75 or 150 mg BID was based upon detailed shedding data from a human inoculation study, where adults receive oseltamivir at 28-h post-inoculation. In reality, the impact on viral shedding will decrease, as patients initiate treatment later in their disease course. Therefore, even though many assumptions above may underestimate the true effects of oseltamivir during an influenza pandemic, delayed initiation of treatment would result in smaller clinical and HE benefits.

Despite these limitations, our study has many, important strengths. First, because we used a dynamic transmission model, we were able to incorporate important indirect effects of oseltamivir treatment. Second, we undertook substantial effort through comprehensive literature review and meta-analysis to synthesise the best available evidence in generating input parameters. We, therefore, feel that our model inputs are of highest quality. Third, not only is our integration of pharmacology into an epidemiological model of influenza transmission novel, it also allows for a more realistic drug effect variability than previous efforts which have relied simply on blunt assumptions as to the efficacy of antivirals in reducing the infectiveness of the virus [Reference Ferguson62]. Because of the increased realism of the drug effect, the transmission dynamics in our model more closely mirror those that would occur in the real-world. Therefore, the results from our economic analysis will be more reflective of real-life, thus, improving the usability of this economic model for decision makers. Finally, including pharmacology, epidemiology and HEs into one integrated model may facilitate earlier meaningful dialogue between the key stakeholders, such as sponsors, regulators and payers. This concept provides a crucial element for so-called adaptive licensing approaches of drug development such as staggered approval, managed entry and progressive authorisation [Reference Eichler63].

Conclusions

This is the first study demonstrating the feasibility of linking of PK/PD–EPI/HE models. Our linked model, which explored the use of oseltamivir in pandemic influenza as an example, allows for direct integration of antiviral effectiveness on viral shedding to the disease transmission model. This, in turn, leads to better understanding of the effects of the different doses of oseltamivir on mortality, morbidity and economic outcomes that make up the outputs of our economic model. Simulations that better reflect reality lead to more informed policy making, which, in turn, can save healthcare costs as well as patient lives. However, efforts are needed to further refine the model to better represent the reality. It enables us to better inform multiple stakeholders even at the early phase of drug development.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268818000158.

Acknowledgement

This work was supported by the funding from F. Hoffmann-La Roche Ltd.

References

1. Weinstein, MC, et al. (2003) Principles of good practice for decision analytic modeling in health-care evaluation: report of the ISPOR task force on good research practices – modeling studies. Value in Health 6(1), 917.Google Scholar
2. Jit, M and Brisson, M (2011) Modelling the epidemiology of infectious diseases for decision analysis: a primer. PharmacoEconomics 29(5), 371386.Google Scholar
3. Lugner, AK, Mylius, SD and Wallinga, J (2010) Dynamic versus static models in cost-effectiveness analyses of anti-viral drug therapy to mitigate an influenza pandemic. Health Economics 19(5), 518531.Google Scholar
4. Brisson, M and Edmunds, WJ (2003) Economic evaluation of vaccination programs: the impact of herd-immunity. Medical Decision Making 23(1), 7682.Google Scholar
5. Brisson, M and Edmunds, WJ (2006) Impact of model, methodological, and parameter uncertainty in the economic analysis of vaccination programs. Medical Decision Making 26(5), 434446.CrossRefGoogle ScholarPubMed
6. Lugner, AK and Postma, MJ (2009) Mitigation of pandemic influenza: review of cost-effectiveness studies. Expert Review of Pharmacoeconomics and Outcomes Research 9(6), 547558.CrossRefGoogle ScholarPubMed
7. Duffull, SB, Wright, DF and Winter, HR (2011) Interpreting population pharmacokinetic-pharmacodynamic analyses – a clinical viewpoint. British Journal of Clinical Pharmacology 71(6), 807814.Google Scholar
8. Rayner, CR, et al. (2013) Pharmacokinetic-pharmacodynamic determinants of oseltamivir efficacy using data from phase 2 inoculation studies. Antimicrobial Agents and Chemotherapy 57(8), 34783487.Google Scholar
9. Bulik, C, et al. Pharmacokinetic-pharmacodynamic (PK-PD) evaluation of the impact of oseltamivir on influenza viral endpoints. 54th Interscience Conference on Antimicrobial Agents and Chemotherapy (ICAAC). Washington DC, US.Google Scholar
10. Pradas-Velasco, R, Antonanzas-Villar, F and Martinez-Zarate, MP (2008) Dynamic modelling of infectious diseases: an application to the economic evaluation of influenza vaccination. PharmacoEconomics 26(1), 4556.Google Scholar
11. Pitman, R, et al. (2012) Dynamic transmission modeling: a report of the ISPOR-SMDM modeling good research practices task force working group-5. Medical Decision Making 32(5), 712721.Google Scholar
12. Girard, MP, et al. (2010) The 2009 A (H1N1) influenza virus pandemic: a review. Vaccine 28(31), 48954902.Google Scholar
13. World Health Organisation (WHO) (2013) Cumulative Number of Confirmed Human Cases for Avian Influenza A (H5N1) Reported to WHO, 2003–2013. Geneva: WHO.Google Scholar
14. World Health Organization. Global alert and response (GAR). Avian influenza. Geneva: WHO.Google Scholar
15. Kelso, JK, et al. (2013) Economic analysis of pandemic influenza mitigation strategies for five pandemic severity categories. BMC Public Health 13, 211.Google Scholar
16. Milne, GJ, Halder, N and Kelso, JK (2013) The cost effectiveness of pandemic influenza interventions: a pandemic severity based analysis. PLoS ONE 8(4), e61504.Google Scholar
17. Kamal, MA, et al. (2017) Interdisciplinary pharmacometrics linking oseltamivir pharmacology, influenza epidemiology and health economics to inform antiviral use in pandemics. British Journal of Clinical Pharmacology 83(7), 15801594.CrossRefGoogle Scholar
18. Kamal, MA, et al. (2013) Population pharmacokinetics of oseltamivir: pediatrics through geriatrics. Antimicrobial Agent and Chemotherapy 57(8), 34703477.Google Scholar
19. Jain, S, et al. (2012) Influenza-associated pneumonia among hospitalized patients with 2009 pandemic influenza A (H1N1) virus – United States, 2009. Clinical Infectious Diseases 54(9), 12211229.Google Scholar
20. Jain, S, et al. (2009) Hospitalized patients with 2009 H1N1 influenza in the United States, April–June 2009. New England Journal of Medicine 361(20), 19351944.Google Scholar
21. Skarbinski, J, et al. (2011) Hospitalized patients with 2009 pandemic influenza A (H1N1) virus infection in the United States – September–October 2009. Clinical Infectious Diseases 52(suppl. 1), S50S59.Google Scholar
22. Burch, J, et al. (2009) Antiviral drugs for the treatment of influenza: a systematic review and economic evaluation. Health Technology Assessment 13(58), 1265, iii–iv.Google Scholar
23. Jefferson, T, et al. (2012) Neuraminidase inhibitors for preventing and treating influenza in healthy adults and children. Cochrane Database of Systematic Reviews 1, CD008965.Google ScholarPubMed
24. Hsu, J, et al. (2012) Antivirals for treatment of influenza: a systematic review and meta-analysis of observational studies. Annals Internal Medicine 156(7), 512524.Google Scholar
25. Michiels, B, et al. (2013) The value of neuraminidase inhibitors for the prevention and treatment of seasonal influenza: a systematic review of systematic reviews. PLoS ONE 8(4), e60348.CrossRefGoogle ScholarPubMed
26. HCUP Nationwide Inpatient Sample (NIS) (2009) Healthcare Cost and Utilization Project (HCUP). Rockville, MD: Agency for Healthcare Research and Quality.Google Scholar
27. U.S. Bureau of Labor Statistics (2013) CPI Inflation Calculator. Washington DC, US: U.S. Bureau of Labor Statistics.Google Scholar
28. Molinari, NA, et al. (2007) The annual impact of seasonal influenza in the US: measuring disease burden and costs. Vaccine 25(27), 50865096.Google Scholar
29. Gonzalez, MK, et al. (2013) The Evolving Role of Emergency Departments in the United States. Santa Monica, CA: RAND Corporation.Google Scholar
30. Angus, DC, et al. (2001) Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care. Critical Care Medicine 29(7), 13031310.Google Scholar
31. Presanis, AM, et al. (2009) The severity of pandemic H1N1 influenza in the United States, from April to July 2009: a Bayesian analysis. PLoS Medicine 6(12), e1000207.Google Scholar
32. Physician Desk Reference (2013) PDR Red Book: Pharmacy's Fundamental Reference. Montvale: Thompson Healthcare Inc.Google Scholar
33. Caldwell, N, et al. (2013) “How much will I get charged for this?” patient charges for top ten diagnoses in the emergency department. PLoS ONE 8(2), e55491.Google Scholar
34. Wiesen, J, et al. (2012) Relative cost and outcomes in the intensive care unit of acute lung injury (ALI) due to pandemic influenza compared with other etiologies: a single-center study. Annals Intensive Care 2(1), 41.Google Scholar
35. Yu, H, et al. (2012) Clinical and economic burden of community-acquired pneumonia in the medicare fee-for-service population. Journal of the American Geriatrics Society 60(11), 21372143.Google Scholar
36. MacLaren, R, et al. (2008) Clinical and economic outcomes of involving pharmacists in the direct care of critically ill patients with infections. Critical Care Medicine 36(12), 31843189.Google Scholar
37. Internal Revenue Service (2013) Standard Mileage Rates for 2013. Washington DC, US: IRS.Google Scholar
38. Center of Disease Prevention (2009) Distance to Nearest Hospital Files, NAMCS and NHAMCS. Atlanta: CDC.Google Scholar
39. U.S. Bureau of Labor Statistics (2013) Median Weekly Earnings by Age, Sex, Race and Hispanic or Latino Ethnicity, First Quarter 2013. Washington DC, US: U.S. Bureau of Labor Statistics.Google Scholar
40. Meltzer, MI, Cox, NJ and Fukuda, K (1999) The economic impact of pandemic influenza in the United States: priorities for intervention. Emerging Infectious Diseases 5(5), 659671.Google Scholar
41. Khazeni, N, et al. (2009) Effectiveness and cost-effectiveness of expanded antiviral prophylaxis and adjuvanted vaccination strategies for an influenza A (H5N1) pandemic. Annals Internal Medicine 151(12), 840853.Google Scholar
42. Fryback, DG, et al. (1993) The beaver dam health outcomes study: initial catalog of health-state quality factors. Medical Decision Making 13(2), 89102.CrossRefGoogle ScholarPubMed
43. Turner, DA, et al. (2006) The cost-effectiveness of influenza vaccination of healthy adults 50–64 years of age. Vaccine 24(7), 10351043.Google Scholar
44. Song, Y, et al. (2012) The potential economic value of a Staphylococcus aureus vaccine among hemodialysis patients. Vaccine 30(24), 36753682.Google Scholar
45. Angus, DC, et al. (2001) Quality-adjusted survival in the first year after the acute respiratory distress syndrome. American Journal of Respiratory and Critical Care Medicine 163(6), 13891394.Google Scholar
46. Talmor, D, et al. (2008) The costs and cost-effectiveness of an integrated sepsis treatment protocol. Critical Care Medicine 36(4), 11681174.CrossRefGoogle ScholarPubMed
47. Davies, A, et al. (2005) Cost effectiveness of drotrecogin alfa (activated) for the treatment of severe sepsis in the United Kingdom. Anaesthesia 60(2), 155162.Google Scholar
48. Drabinski, AWG and Formica, C (2001) Observational evaluation of health state utilities among a cohort of sepsis patients. Value in Health 4(2), 128129.CrossRefGoogle Scholar
49. Macario, A, Chow, JL and Dexter, F (2006) A Markov computer simulation model of the economics of neuromuscular blockade in patients with acute respiratory distress syndrome. BMC Medical Informatics and Decision Making 6, 15.Google Scholar
50. U.S. Bureau of Labor Statistics (2012) Household Data, Annual Averages. Washington DC, US: U.S. Bureau of Labor Statistics.Google Scholar
51. Russell, CA, et al. (2012) The potential for respiratory droplet-transmissible A/H5N1 influenza virus to evolve in a mammalian host. Science 336(6088), 15411547.Google Scholar
52. Kawaoka, Y (2012) H5n1: flu transmission work is urgent. Nature 482(7384), 155.Google Scholar
53. Lee, BY, et al. (2010) To test or to treat? An analysis of influenza testing and antiviral treatment strategies using economic computer modeling. PLoS ONE 5(6), e11284.Google Scholar
54. World Health Organisation (WHO) (2011) Report of the Review Committee on the Functioning of the International Health Regulations (2005) in Relation to Pandemic (H1N1) 2009. Geneva: WHO.Google Scholar
55. Fineberg, HV (2014) Pandemic preparedness and response – lessons from the H1N1 influenza of 2009. New England Journal of Medicine 370(14), 13351342.Google Scholar
56. Herfst, S, et al. (2012) Airborne transmission of influenza A/H5N1 virus between ferrets. Science 336(6088), 15341541.Google Scholar
57. Imai, M, et al. (2012) Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature 486(7403), 420428.CrossRefGoogle ScholarPubMed
58. Li, FC, et al. (2008) Finding the real case-fatality rate of H5N1 avian influenza. Journal of Epidemiology and Community Health 62(6), 555559.Google Scholar
59. Cowling, BJ, et al. (2013) Comparative epidemiology of human infections with avian influenza A H7N9 and H5N1 viruses in China: a population-based study of laboratory-confirmed cases. Lancet 382(9887), 129137.CrossRefGoogle Scholar
60. Jefferson, T, et al. (2014) Oseltamivir for influenza in adults and children: systematic review of clinical study reports and summary of regulatory comments. British Medical Journal 348, g2545.Google Scholar
61. Muthuri, SG, et al. (2014) Effectiveness of neuraminidase inhibitors in reducing mortality in patients admitted to hospital with influenza A H1N1pdm09 virus infection: a meta-analysis of individual participant data. Lancet Respiratory Medicine 2(5), 395404.Google Scholar
62. Ferguson, NM, et al. (2006) Strategies for mitigating an influenza pandemic. Nature 442(7101), 448452.Google Scholar
63. Eichler, HG, et al. (2012) Adaptive licensing: taking the next step in the evolution of drug approval. Clinical Pharmacology and Therapeutics 91(3), 426437.Google Scholar
Figure 0

Fig. 1. Overarching Pharmacology to payer system including ‘modules’. The solid lines indicate that adequate data exists to be able to create semi-mechanistic links to each of adjacent ‘modules’. The dotted lines and light grey describe where significant unknowns remain and are not mature enough to have been incorporated into the current framework. PopPK, population pharmacokinetic; SEIR, susceptible-exposed-infected-recovered; ß, the rate of infectivity; TShed, viral shedding; VK, viral kinetics.

Figure 1

Table 1. Parameters used in the SEIR model

Figure 2

Fig. 2. Decision analytic model. Influenza patients entered the decision analytic model from epidemiology model. They received treatment in outpatient or inpatient setting. OPD, outpatient; ED, emergency department; AVR, antiviral; GW, general ward; ICU, intensive care unit; ARDS, acute respiratory distress syndrome.

Figure 3

Table 2. Input parameters, values and data sources used in the health economics model

Figure 4

Table 3. Base-case analyses: high-dose vs. no treatment, and standard dose vs. no treatment

Figure 5

Fig. 3. Tornado diagram (150 mg vs. no treatment with 80% uptake of oseltamivir): One-way sensitivity analysis under (1a) low virulence and low transmissibility and (1b) high virulence and high transmissibility. (a) Low virulence and low transmissibility. #+/−: The higher value of the parameter leads to lower ICER. $−/+: The higher value of the parameter leads to higher ICER. (b) High virulence and high transmissibility. #+/−: The higher value of the parameter leads to lower ICER.

Figure 6

Fig. 4. Scatter plots (incremental cost vs. incremental QALY) of 75 mg vs. no treatment under societal perspective. (a) Low virulence and low transmissibility. (b) High virulence and high transmissibility.

Supplementary material: File

Wu et al. supplementary material

Appendices 1-2

Download Wu et al. supplementary material(File)
File 24.4 KB