key: cord-0004236-k3z30aaw authors: Arefin, Md. Rajib; Masaki, Tanaka; Kabir, K. M. Ariful; Tanimoto, Jun title: Interplay between cost and effectiveness in influenza vaccine uptake: a vaccination game approach date: 2019-12-18 journal: Proc Math Phys Eng Sci DOI: 10.1098/rspa.2019.0608 sha: 39860401de24cbae315f33d3e9566552ac45e352 doc_id: 4236 cord_uid: k3z30aaw Pre-emptive vaccination is regarded as one of the most protective measures to control influenza outbreak. There are mainly two types of influenza viruses—influenza A and B with several subtypes—that are commonly found to circulate among humans. The traditional trivalent (TIV) flu vaccine targets two strains of influenza A and one strain of influenza B. The quadrivalent (QIV) vaccine targets one extra B virus strain that ensures better protection against influenza; however, the use of QIV vaccine can be costly, hence impose an extra financial burden to society. This scenario might create a dilemma in choosing vaccine types at the individual level. This article endeavours to explain such a dilemma through the framework of a vaccination game, where individuals can opt for one of the three options: choose either of QIV or TIV vaccine or none. Our approach presumes a mean-field framework of a vaccination game in an infinite and well-mixed population, entangling the disease spreading process of influenza with the coevolution of two types of vaccination decision-making processes taking place before an epidemic season. We conduct a series of numerical simulations as an attempt to illustrate different scenarios. The framework has been validated by the so-called multi-agent simulation (MAS) approach. Pre-emptive vaccination is regarded as one of the most protective measures to control influenza outbreak. There are mainly two types of influenza virusesinfluenza A and B with several subtypes-that are commonly found to circulate among humans. The traditional trivalent (TIV) flu vaccine targets two strains of influenza A and one strain of influenza B. The quadrivalent (QIV) vaccine targets one extra B virus strain that ensures better protection against influenza; however, the use of QIV vaccine can be costly, hence impose an extra financial burden to society. This scenario might create a dilemma in choosing vaccine types at the individual level. This article endeavours to explain such a dilemma through the framework of a vaccination game, where individuals can opt for one of the three options: choose either of QIV or TIV vaccine or none. Our approach presumes a mean-field framework of a vaccination game in an infinite and well-mixed population, entangling the disease spreading process of influenza with the coevolution of two types of vaccination decision-making processes taking place before an epidemic season. We conduct a series of numerical simulations as an attempt to illustrate different scenarios. The framework has been validated by the so-called multi-agent simulation (MAS) approach. inspiration from [40] , which mainly focused on a networked population. A similar work related to [40] can be found in [41] . The number of vaccinees remains the same during an epidemic season as we are considering seasonal influenza. Moreover, as influenza vaccines are not perfectly effective [42] , we take into account vaccine efficacy in our modelling. The assumption of imperfect vaccine in the vaccination game can tackle more realistic situations, and lead to interesting results, which otherwise can be intractable (for instance, see [43, 44] ). In our current setting, at the end of each season, individuals assess their perceived payoffs based on the previous season's experience in terms of vaccine efficacy, cost, etc., and decide whether to vaccinate or not or if vaccinating then which vaccine to take for the next season. The strategy adaptation after each season is based on imitating others' strategy with better payoffs. The group who chooses QIV vaccine are expected to get better protection against the disease than that of the group taking TIV vaccine, although this may depend on the circulation of influenza B virus. However, the cost for TIV vaccine should be no more than the cost for QIV vaccine. Taking into account all these factors, we combine the disease spreading dynamics of influenza with the vaccination decision-making process through an evolutionary framework. We perform our simulations for repeated seasons until we reach a social equilibrium (steady state). By varying vaccine effectiveness and cost level, our numerical simulations intend to show an overall picture of influenza epidemics. At equilibrium, we calculate to what extent the average social payoff attained from the vaccination game falls short from the desired or social optimum payoff, which we name as, 'social efficiency deficit'. Moreover, our approach has been validated with the multi-agent simulation (MAS) approach. We also explore the robustness of our framework in explaining the situation with influenza A or B virus-dominant circulation. The rest of this paper is organized as follows: §2 provides the detailed illustration of the model covering the epidemic spreading, payoff structure, strategy update and the evolutionary dynamics of vaccination coverage. Section 3 delivers a comprehensive discussion resulting from a series of numerical simulations. Section 4 concludes the findings of this work and provides some future indications. We entangle the simultaneous spreading of two influenza viruses (A and B) and the evolution of vaccination (QIV or TIV or none) decision by constructing a repetitive sequence of a two-stage process in an infinite and well-mixed population. The disease spreading process is governed by susceptible-infected-recovered (SIR)-like dynamics [45] coupled with vaccination that allow us to estimate several fractions of individuals at the end of each epidemic period. These different groups then evaluate their payoffs and revisit their vaccination strategies by copying other individuals' strategy who achieved higher payoffs in the previous season. Finally, we update the fraction of vaccinees of both types (QIV and TIV) using evolutionary equations given in §2e. Let us note that the whole process is repeated for several steps (generations) until we reach a steady state. Figure 1 illustrates the whole dynamical set-up. The first stage of our dynamical system is the disease spreading process that is based on a SIR-like model coupled with vaccination. Initially, we divide the whole population into three categories: susceptible (unvaccinated) denoted by S, QIV vaccinees denoted by V Q and TIV vaccinees denoted by V T . It is not usually recommended to take both vaccines (TIV and QIV) in a single season [46] ; therefore, we disregard any group taking both vaccines in the same season. As influenza vaccines are not 100% perfect, we consider imperfect vaccinations in our disease modelling, where both vaccines are presumed to provide the same level of efficiency against A virus but different efficacies against B virus. The inclusion of these parameters allows us to consider some individuals from vaccinated groups who fail to get immunity from vaccination (QIV or TIV) and still face the risk of infection like susceptible. The unvaccinated or vaccinated people who fail to get immunity may become infected with either influenza virus (A or B), and then these infected people transfer to the recovery state after suffering from illness. Therefore, at time t, an individual can be at any of the six compartments: susceptible, S(t); QIV vaccinees, V Q (t); TIV vaccinees, V T (t); infected with influenza A, I A (t); infected with influenza B, I B (t); recovered state, R(t). We presume e Q (0 ≤ e Q ≤ 1) as the vaccine efficacy of QIV against both influenza A and B viruses. As TIV (QIV) vaccine contains two strains of influenza A virus and one (two) strain(s) of influenza B virus, we assume that both QIV and TIV vaccines bestow the same level of effectiveness (e Q (0 ≤ e Q ≤ 1)) against influenza A but TIV vaccine provides relatively lower effectiveness, e T (0 ≤ e T ≤ e Q ) against influenza B than that of QIV vaccine. If V Q (0) is the fraction of vaccinees at the beginning of an epidemic season, then a proportion e Q of V Q (0) i.e. e Q V Q (0) gets perfect immunity from QIV vaccine and, the remaining fraction, V Q (t) − e Q V Q (0) becomes vulnerable to the disease, where V Q (t) is the fraction of vaccinators in the V Q compartment at time t. A similar logic has been used for TIV vaccinees. The transmission rates and recovery rates for influenza A and B viruses are assumed β A , γ A and β B , γ B , respectively. It is worth mentioning that the state space variables represent several fractions of the total population; therefore, we are not considering a population of any particular size. Although coinfection and superinfection with two influenza viruses are possible in some cases [47, 48] , the incidence rate is not so significant. Therefore, our model disregards the incidence of coinfection and superinfection with two influenza viruses (A and B) and presumes long-term cross-immunity between them. Although we assume long-term crossimmunity between two viruses, short-term immunity is possible in reality. The incidence of short-term cross-immunity between two viruses may allow some individuals recovered from one virus to be infected with the other one after a short time interval (e.g. [49] ). The schematic of our epidemic model is given in the disease spreading part of figure 1. The model equations are given by and x + y ≤ 1) are the fractions of QIV and TIV vaccinees at the beginning of a season, i.e. V Q (0) = x, V T (0) = y, then where ε A ( → 0 + ) and ε B ( → 0 + ) are initial infections with influenza A and influenza B, respectively. After each epidemic season, we estimate several groups by numerically calculating fluxes from one state to another. We mainly estimate nine fractions of individuals that are, QIV vaccinees and healthy-HV Q , QIV vaccinees but infected with influenza A(B) virus-V Q I A (V Q I B ), TIV vaccinees and healthy-HV T , TIV vaccinees but infected with influenza A(B) virus-V T I A (V T I B ), successful free riders-SFR (who are unvaccinated but did not suffer from any influenza virus), failed free riders-FFR A or FFR B (unvaccinated individuals infected by either influenza virus). Figure 1 . The layout of the whole dynamical set-up. The vaccine efficacy of TIV vaccine against influenza B virus is assumed e T , which is lower than that of QIV vaccine. However, both vaccines are assumed to have same efficacy (e Q ) against influenza A virus. Once the disease spreading process ends, we estimate several fractions such as vaccinated and healthy, vaccinated but infected, infected with influenza A(B), etc., and evaluate their payoffs. These fractions then update their strategies for the next season. This process is repeated until we reach a steady state. Arrows depict the sequence of the evolutionary process. (Online version in colour.) who remain at S, V Q and V T state at equilibrium (t → ∞), that is S(∞), V Q (∞) and V T (∞), are termed successful free riders, healthy QIV vaccinees and healthy TIV vaccinees, respectively. Suppose ϕ P→Q is the total proportion transferred from state P to state Q during an epidemic season. Thus, QIV vaccinees who get infected with influenza A(B) virus can be evaluated from the flux, In a similar way, we can estimate this for TIV vaccinees. The failed free riders can be estimated by fluxes, ϕ S→I A and ϕ S→I B . Once an epidemic season ends, several groups estimate their payoffs prior to the onset of the next epidemic season. Suppose the cost of vaccination and infection are, respectively, C V and C I . The vaccination cost includes the vaccine price with possible side effects due to vaccination. Since the infection cost is higher than the vaccination cost, we define the relative cost of vaccination as C r = C V /C I , (0 ≤ C r ≤ 1). Without loss of generality, we can choose C I = 1. Then the relative cost for QIV (TIV) vaccine can be assumed as C Q (C T ). Since the cost for TIV is assumed to be no more than QIV, we can presume 0 ≤ C Q ≤ 1 and 0 ≤ C T ≤ C Q . Therefore, QIV(TIV) vaccinees who remain healthy during the epidemic period can have a payoff, − C Q (−C T ); however, QIV(TIV) vaccinees who become infected with either influenza virus (A or B) have a payoff, −C Q − 1(−C T − 1). In addition, payoffs for successful free riders and failed free riders are assigned as 0 and −1, respectively. All fractions along with their payoffs are summarized in table 1. Using table 1, we estimate the average social payoff π ; average payoff for QIV vaccinees ( π V Q ), TIV vaccinees ( π V T ) and non-vaccinators ( π NV ) as follows: and Individuals assess their payoffs after an epidemic season and decide whether to imitate others' strategy or stay with their previous season's strategy. In a seminal work regarding the agent-based vaccination game model, Fu et al. [40] first used the pairwise Fermi function [50, 51] to calculate the likelihood of copying others' strategy by comparing payoffs, where a focal agent i having the strategy S i with the payoff π i adopts the strategy S j of a randomly chosen agent j having the payoff π j with the probability where the parameter, k > 0 measures the strength of selection; smaller k signifies the strong selection, that is players are more responsive to the payoff difference. A typical choice of k is 0.1, which has been used in many previous studies such as [16, [18] [19] [20] 29] . Fukuda et al. [18] later slightly modified the formula (2.6), where rather than comparing with a single agent, a focal player i compares its payoff π i with an average payoff π j of all neighbours having the same strategy S j of a randomly selected neighbour j. The modified formula becomes Although formulae (2.6) and (2.7) have been mainly used for agent-based modelling, these can also be used for the mean-field approach (such as [16, 20, 28, 52] ). Inspired by that, our approach uses formula (2.7) to update vaccination strategies for the next season. In our case, every fraction (each of nine fractions in table 1) compares its payoff with the average payoff of a mutually exclusive group (QIV, TIV, non-vaccinated in table 1); for instance, the fraction, HV Q can, respectively, compare its payoff with the average payoff of TIV vaccinees group (V T ) and non-vaccinators group (NV) by the formula Similarly, we estimate all necessary transition probabilities. As mentioned before, individuals update their vaccination strategies prior to onset of the next epidemic season. The evolution of vaccination strategies (QIV or TIV or unvaccinated) is estimated using the following evolutionary equations that have been derived by extending the master equation of the mean-field framework in a well-mixed population [52, 53] Here additions and subtractions, respectively, indicate inflows to and outflows from a state. We solve equations (2.8) and (2.9) numerically, allowing the system to reach the steady state as, t → ∞. This section deliberately explains results of the whole dynamical process. We vary several parameters such as transmission rates, vaccine effectiveness, vaccination cost, etc., to illustrate different scenarios. The symptoms of usual influenza sickness resolve after a period of 3-7 days, so the average durability in the infected stage is 5 days [49] . Therefore, we presume γ A = γ B = 1/5 = 0.2 throughout our numerical experiments. Let us first briefly discuss the relative dynamics of infections due to influenza A and B viruses for different transmission rates and degree of initial infections in a single season without considering game aspect (figure 2). The vaccination coverage for this case is chosen as 2a,b) ). However, the virus with higher transmission rate is found to dominate irrespective to the initial infection. Hence, the sensitivity of initial infections seems insignificant comparing to the transmission rates, i.e. transmission rates show more impact on virus dominance. We estimate the coevolution of both types of vaccination at the end of each season. Owing to the complexity of equations (2.8) and (2.9), it is difficult to theoretically derive all possible equilibria; however, it is still possible to derive that numerically by tuning different parameters. Considering equal transmission rate for both viruses, we estimate six possible combinations of evolutionary outcomes (figure 3), which are (x, y) ≡ (1, 0);(x, y) ≡ (0, 0);(x, y) ≡ (x*, y*);(x, y) ≡ (x*, 0);(x, y) ≡ (0, 1);(x, y) ≡ (0, y*); where 0 < x*, y* < 1. Figure 3 suggests that there is none of the bi-stability arising from equations (2.8) and (2.9). One possible reason is that the present model as a to be dominant whenever both vaccination costs are comparable. Also, we can perceive the coexistence of both vaccinees in the borderline region between blue and red colours (yellowish regime enclosed by dotted lines in (a-iii) and (a-iv)). Results obtained from the mean-field framework have been justified by a sequence of numerical simulations based on the MAS approach [18, 40] presuming a complete graph as an underlying network (since we assume well-mixed population) with the population size, N = 10 000. The transmission rate parameters are chosen as β A = β B = 5.19957 × 10 −5 , which is estimated as the minimum transmission rate that surpasses the preset threshold final epidemic size of 0.9 without any vaccination [16] . The corresponding initial infections are I A (0) = 2 agents and I B (0) = 4 agents. Note that we estimate ensemble average for each 100 realizations. Generally, results obtained from the mean-field framework and MAS approach are showing overall the same tendency, albeit having subtle discrepancies in terms of colour scaling that comes from the fact that the MAS approach undertakes a finite population. Now let us focus on analysing similar phase diagrams by varying vaccine efficacies. To this aim, we generate several heatmaps as above as a function of (e Q , e T ), 0 ≤ e T ≤ e Q ≤ 1, for the case of similar and different cost levels ( figure 5) . Obviously, similar cost for both vaccines would encourage people to choose QIV vaccine instead of TIV because the former one bestows better protection against both viruses. This situation illustrates why the sensitivity comes only along the direction of e Q (upper panels in figure 5) . A lower degree of vaccine efficacy would not entice individuals to take vaccine (red region in figure 5a-iii); however, if e Q passes a threshold level, the fraction of QIV vaccinees mounts to the highest level (blue region in figure 5a-iii)). A further increase of e Q seems to suppress both infections significantly, which allows some people to avoid vaccination by free riding on herd immunity. This is why, after passing a transient regime having maximum vaccination coverage, the fraction of vaccinees decreases monotonically even with the increase of e Q (yellowish region in figure 5a-iii) . Nevertheless, if the cost for QIV is higher than TIV (lower panels in figure 5 ), then individuals' vaccination choice differs depending upon the degrees of e Q and e T . More specifically, under the current setting, QIV vaccine is favoured over TIV vaccine if e Q is above 50% (approximately) and e T is below 40% (approximately); on the other hand, if e T is above 40% (approximately), then TIV vaccine is preferred over QIV vaccine regardless of e Q (of course with the condition, e T ≤ e Q ). Notably, in this case, the fraction of QIV vaccinees never reaches the maximum as the corresponding cost is higher; contrarily, TIV vaccinees are seen to reach the highest level for a mid-range of e T , although it starts decreasing with the further increase of e T that arises with the prevalence of free riders. Remarkably, the average payoff of society in case of different costs seems lower than that of the equal cost (see heatmaps for ASP in figure 5 ) as the former case imposes a higher financial burden to society. The case β A > β B (β A = 0.6, β B = 0.4) leads to infection dominance of influenza A over B virus ( figure 6a,b) . This situation can be controlled by TIV vaccine alone as it targets two strains of A virus. Individuals would then mostly prefer TIV vaccine over QIV vaccine as the price of TIV is lower than the price of QIV ( figure 6c,d) . As a consequence, we perceive the sensitivity of choosing TIV vaccine only along the direction of C T . The predominance of influenza B virus over A virus creates a dilemma of choosing vaccine types. Individuals prefer QIV vaccine over TIV whenever both costs are comparable or the cost difference between QIV and TIV is not so high (region enclosed by black dotted lines in figure 6g) ; however, TIV vaccine seems mostly favourable in case of higher vaccination cost for QIV or higher cost difference between QIV and TIV vaccine (region enclosed by red dotted lines in figure 6h ). Figure 6g ,h also depicts the coexistence of both types of vaccinees in the yellowish transient regions. Furthermore, comparing the light blue region (enclosed with black dotted lines) with the greenish region (enclosed by red dotted lines) in figure 6f , it is observed that the QIV vaccine-dominant region imposes more effect on disease suppression compared to TIV vaccine. (iv) β A versus β B phase plane We intend to observe an overall impact of the transmission rates on disease propagation and vaccine uptake. To this aim, we draw heatmaps for total infection (FES) due to both viruses, the basic reproduction number (R 0 ) and vaccination coverage for both vaccines (QIV, TIV), as a function of (β A , β B ) ( figure 7) . The vaccine-dependent basic reproduction number R 0 [54, 55] . With fixed (C * Q ,C * T ) and y*, we vary x (0 ≤ x ≤ 1) and estimate the maximum average payoff, which we term as social optimum (SO) payoff for x at (C * Q ,C * T ) and y*, i.e. ASP SO x (C * Q ,C * T |y * ) (b). After that we derive SED for x at (C * Q , C * T ) by subtracting ASP at equilibrium from the payoff at social optimum. We follow the same procedure to derive SED for y using (a) and (c). (Online version in colour.) where R A 0 and R B 0 are the vaccine-dependent basic reproduction number for influenza A and B viruses, respectively. We derive R A 0 and R B 0 as and Other Consequently, we observe no vaccination coverage inside that box (figure 7c,d). We perceive a red dotted line in FES-heatmap (figure 7a) below(above) which influenza A(B) predominates. As QIV vaccine provides better protection against B virus than TIV vaccine, we perceive the predominance of QIV vaccinees in case of higher β B ; contrarily, TIV vaccine is favoured over QIV for higher β A because its price is lower than QIV. Now we intend to explore how much the payoff attained at equilibrium falls short from the desired payoff or social optimum (SO) payoff. We name this payoff shortfall as the social efficiency deficit (SED) and define as SED, δ = ASP EQM − ASP SO , (3.1) where ASP EQM and ASP SO denote the average social payoff at equilibrium and social optimum, respectively. We use this parameter in explaining the degree of dilemma associated with each vaccine choice. The whiteout regions in (c,d) depict that there is no gap between SO payoff (for x or y) and ASP at equilibrium. More areas in SED-heatmap for x appears to have darker regime than that of the case for y, which indicates the higher gap between SO payoff for x and ASP at equilibrium. (Online version in colour.) for both vaccines. Suppose we wish to derive SED as a function of (C Q , C T ). Note that it is also possible to express SED as a function of (e Q , e T ). Step (I) At first, we estimate ASP at equilibrium (EQM) as a function of (C Q , C T ), for instance consider the heatmap for ASP given in figure 4a-v. Step (II) Each point on this heatmap corresponds to a pair, (C * Q , C * T ). Also, there is a corresponding pair (x*, y*), i.e. vaccination coverage for each vaccine, for each pair (C * Q , C * T ). Now choosing y*(x*) as fixed and varying x(y) from 0 up to 1, we estimate the maximum average payoff (using formula (2.2)), what we call a social optimum payoff for x(y) at the point (C * Q , C * T ) (figure 8a-c). Step (III) Next, we derive SED for x(y) at the point (C * Q , C * T ) by simply subtracting ASP attained in step I from the social optimum payoff attained in step II and denote it by δ x Figure 8 schematically explains the whole procedure. Now let us explain the heatmaps for SED that are generated following the above procedure. As an archetype, we choose the case illustrated in the upper panel of figure 4 , where we assumed β A = β B = 0.5 and e Q = 0.6, e T = 0.4. In that case, we observed that most of the people are inclined to choose TIV vaccine whenever the cost difference between C Q and C T gets higher but individuals prefer QIV if the cost difference becomes lower. It is conceivable that the overall desired payoff or social optimum payoff for QIV vaccinees (x) would be relatively higher than the TIV vaccinees (y), which is appeared in figure 9a ,b, but in this case, the vaccination coverage for QIV seems much lower than TIV (figure 4a-iii-a-iv), which indicates that the dilemma of choosing QIV vaccine is higher than the other one. Obviously, this comes from the fact of higher cost associated with QIV vaccine. This factor results in less contribution of QIV vaccine to ASP ( figure 4a-v) . As a result, we perceive relatively higher payoff-gap between social optimum payoff for x and ASP (at equilibrium); hence, more areas of the heatmap for SED x are seemed dark compared to its counterpart (SED y ). The whiteout region in heatmaps (figure 9c,d) depicts the scenario with no SED, that is there is no gap between social optimum payoff and ASP. There are two white regions in the SED yheatmap (figure 9d), which either come from 100% TIV vaccination coverage (region enclosed by blue dotted lines) resulting from lower C T and higher C Q (blue region in figure 4a-iii), or from the scenario when everyone becomes a free rider (region enclosed by red dotted lines in figure 9b,d). The latter situation can arise for higher vaccination cost or for very low vaccine effectiveness. In this case, avoiding any type of vaccination is the best strategy (desired strategy), which agrees with the D(unvaccinated)-dominant equilibrium (figure 4a-iii,a-iv). As a result, the payoff at this equilibrium also agrees with the social optimum payoff; hence no SED or no dilemma arises in this case. In this paper, we consider a simultaneous spreading of two influenza viruses (A and B) coupled with two types of vaccinations. The QIV vaccine contains one extra strain of influenza B virus and hence provides better protection compared to the TIV vaccine. However, the cost of the QIV vaccine is higher than TIV, which may create a dilemma for individuals' vaccination choice. Here we design a mean-field vaccination game model in an infinite and well-mixed population to explain people's vaccination dilemma, entangling an epidemic spreading process of influenza with the coevolution of two types of vaccinations. Our evolutionary framework assumes three strategies-choose QIV vaccine, TIV vaccine or none-that evolve at the end of each season by imitating strategies based on the previous season's experience. We perform a series of numerical simulations by varying vaccination cost, vaccine effectiveness, transmission rates, etc., to depict different scenarios. Our results show that individuals are more inclined to take QIV vaccine whenever both vaccination costs are comparable; however, they prefer TIV vaccine if the cost difference gets higher. Notably, the framework has been validated by the so-called MAS approach. Generally, people's vaccination choice is influenced by the costs as well as vaccine efficacies; however, our results show that individuals' vaccination options (QIV or TIV) are more affected by the cost difference of both vaccines as long as the vaccine with the lower price is able to bestow a considerable level of efficiency, especially against B virus. Another interesting result concerns the variation of transmission rates for both viruses and its resulting effect on disease propagation and vaccination choice. The higher transmission rate of influenza A virus leads to a maximum coverage for TIV vaccine since it possesses the same degree of effectiveness against influenza A virus as that of QIV vaccine but with a lower cost; contrarily, the higher transmission rate of influenza B virus results in almost complete dominance of QIV vaccinees that is quite conceivable. Our analysis also incorporates how much the overall payoff or the average social payoff (ASP) attained by the process of voluntary vaccination (at equilibrium) falls short from the desired payoff or social optimum payoff, what we have termed as social efficiency deficit (SED)-the payoff gap between social optimum and equilibrium. This parameter also helps us elucidate the degree of dilemma associated with each vaccination choice. It appears that for the case of equal transmission rates, people face a bigger dilemma to choose QIV vaccine rather than TIV that arises from the fact of the higher vaccination cost for QIV even though it confers relatively better protection. In sum, this work endeavours to explore individuals' dilemma on choosing influenza vaccines through the lens of a vaccination game, where the disease spreading process is modelled by a simple SIR-like model coupled with two types of vaccinations. Our approach does not consider any strain-specific infection, rather it relies on virus-specific infection (influenza A and B) with long-term cross-immunity. It would be interesting to explore strain-specific infection, where an individual infected with a particular strain can have short-term or long-term cross-immunity against the other one (e.g. [49] ). The inclusion of age structure (for instance [49, 56] ) in the current framework could be another aspect to explore. The influenza epidemic varies from one country to another country, thus it would be also interesting to explore the outlined framework for a specific country with proper estimation of relevant parameters [56, 57] . In addition to a voluntary vaccination campaign, subsidies [14, 58] provided by the government (especially for QIV vaccine) can enhance vaccination coverage to prevent disease-outbreak more efficiently. Therefore, the current framework with a subsidy provision could be another avenue to explore. In the present study, we have presumed a well-mixed situation. The next step for the current framework may be the inclusion of an underlying network. The focal point of our present model is plural vaccinations having different capabilities and costs for plural viruses of the same disease. Such a specific but quite realistic situation would be significantly influenced by a spatial structure connecting individuals. The well-mixed situation can be said, in a sense, to be an ideal assumption because the social information helping individuals to decide for vaccination is quickly shared and can also be said to be an extreme assumption because disease spreading occurs quickly vis-à-vis the case with the spatial structure. The time delay between information and disease spreading on the underlying network might be more significant for the model with plural vaccinations and plural viruses (or plural strains) than the conventional setting with a single vaccination and a single strain. Data accessibility. All data of this research work have been generated from a series of numerical simulations. The source code (c++) of our numerical simulation has been uploaded as electronic supplementary material in the 'File Upload' section. A systematic review of the social and economic burden of influenza in low-and middle-income countries The cost-effectiveness of trivalent and quadrivalent influenza vaccination in communities in South Africa, Vietnam and Australia World Health Organization. Factsheet Influenza (seasonal) Rationale for vaccination with trivalent or quadrivalent live attenuated influenza vaccines: protective vaccine efficacy in the ferret model Update: influenza Activity in the United states during the 2017-18 season and composition of the 2018-19 influenza vaccine Phylogenetic and evolutionary history of influenza B Viruses Reappearance and global spread of variants of influenza B/Victoria/2/87 lineage viruses in the Prevention and control of seasonal influenza with vaccines: recommendations of the advisory committee on immunization practices-United States Impact of Influenza B lineage-level mismatch between trivalent seasonal influenza vaccines and circulating viruses Simulation studies to assess the long-term effects of Japan's change from trivalent to quadrivalent influenza vaccination Influenza cost and cost-effectiveness studies globally -a review Cost-effectiveness of quadrivalent versus trivalent influenza vaccine in the United States. Value Heal Cost-utility of quadrivalent versus trivalent influenza vaccine in Germany, using an individual-based dynamic transmission model To vaccinate or not to vaccinate: a comprehensive study of vaccination-subsidizing policies with multi-agent simulations and mean-field modeling Analysis of SIR epidemic model with information spreading of awareness Which is more effective for suppressing an infectious disease: imperfect vaccination or defense against contagion? Analysis of epidemic outbreaks in two-layer networks with different structures for information spreading and disease diffusion Risk assessment for infectious disease and its impact on voluntary vaccination behavior in social networks Impact of imperfect vaccination and defense against contagion on vaccination behavior in complex networks A game theoretic approach to discuss the positive secondary effect of vaccination scheme in an infinite and well-mixed population Vaccination and the theory of games Statistical physics of vaccination Dynamical behaviors for vaccination can suppress infectious disease -a game theoretical approach Influence of breaking the symmetry between disease transmission and information propagation networks on stepwise decisions concerning vaccination Effect of information spreading to suppress the disease contagion on the epidemic vaccination game Realistic decision-making processes in a vaccination game Vaccination strategies in a two-layer SIR/V-UA epidemic model with costly information and buzz effect Evolutionary vaccination game approach in metapopulation migration model with information spreading on different graphs Three-strategy and four-strategy model of vaccination game introducing an intermediate protecting measure Effects of stubborn decision-makers on vaccination and disease propagation in social networks Evolutionary dynamics: exploring the equations of life Evolutionary game theory with sociophysics: analysis of traffic flow and epidemics Fundamentals of evolutionary game theory and its applications Evolutionary games and population dynamics Vaccination and epidemics in networked populations-an introduction Saving human lives: what complexity science and information systems can contribute Dilemma game structure observed in traffic flow at a 2-to-1 lane junction Social dilemma structures hidden behind traffic flow with lane changes Effect of intermediate defense measures in voluntary vaccination games Imitation dynamics of vaccination behaviour on social networks Evolutionary vaccination dilemma in complex networks The future of influenza vaccines: a historical and clinical perspective. Vaccines 6 Interplay between cost and benefits triggers nontrivial vaccination uptake Imperfect vaccine and hysteresis A contribution to the mathematical theory of epidemics Immunization Action Coalition Dual infections by influenza A/H3N2 and B viruses and by influenza A/H3N2 and A/H1N1 viruses during winter Occurrence of a case of influenza A(H1N1)pdm09 and B co-infection during the epidemic season 2012-2013 Intermediate levels of vaccination coverage may minimize seasonal influenza outbreaks The statistical mechanics of strategic interaction Evolutionary Prisoner's dilemma game on a square lattice Heterogeneous update mechanisms in evolutionary games: mixing innovative and imitative dynamics Evolutionary games on graphs Reproduction numbers of infectious disease models Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission Prioritization of vaccine strategy using an age-dependent mathematical model for 2009 A/H1N1 influenza in the Republic of Mathematical models of transmission dynamics and vaccine strategies in Hong Kong during the 2017-2018 winter influenza season Game dynamic model of optimal budget allocation under individual vaccination choice Competing interests. We declare we have no competing interests. Funding. This study was partially supported by Grant-in-Aid for Scientific Research from JSPS, Japan, KAKENHI (grant no. 18K18924), SCAT (Support Center for Advanced Telecommunications Technology) Research Foundation and Pfizer Health Research Foundation awarded to Prof. Tanimoto. Also, the computation was mainly carried out using the computer facilities at Research Institute for Information Technology, Kyushu University. We would like to express our gratitude to them. Authors' contributions. M.R.A. developed the model, performed numerical simulations, analysed results and drafted the manuscript. T.M. carried out the multi-agent simulation (MAS) approach to validate the model and critically revised the manuscript. K.M.A.K. assisted in model formulation, result analyses and critically revised the manuscript. J.T. helped designing the study, coordinated the study, also helped draft the manuscript. All authors gave final approval for publication and agree to be held accountable for the work performed therein.