Received 2 April 2016; accepted 23 May 2016; published 26 May 2016
Malaria is a vector born disease that can be transmitted to people through mosquito bites, prevalent in tropical and subtropical regions of the world. The World Health Organisation estimates that every year nearly one million people die as a result of malaria  . A lot of effort has been put into treatment and prevention strategies, and an ambitious goal of organisations such as the World Health Organisation and the Bill and Melinda Gates foundation is the complete elimination of malaria.
One of the key methods of malaria prevention is the use of Insecticide-Treated Bed-Nets (ITNs) that kill and also repel mosquitoes  .  found that the use of ITNs can reduce malaria cases by 50% and death in children by 20%. The use of ITNs was demonstrated to be very cost-effective, see e.g.  -  . This is important given both the scale of the problem (about 200 million people are at a very high risk  ), and the fact that most of the areas of high malaria prevalence are relatively poor. However, human behavior such as inconsistent usage due to hot weather diminishes the effectiveness of the ITNs.
The impact of ITN usage on the spread of the malaria infection has been recently studied in a series of mathematical models such as  -  . In particular,  developed an ODE model for the effects of ITN on the malaria transmission dynamics and called for a better understanding of the impact of human behavior on the usage.
Game theory is a well suited tool to study human behavior in a quantitative way  . Starting by  , game theory has been long used in biology  , and  used it to study vaccination decisions in situation where individuals face a potentially costly preventive action (such as to use an ITN) that can reduce a risk of con- tracting the disease. Game theory has been applied to vaccination against major public health threats  -  and it is particularly useful because the outcome of the individual decision (to vaccinate or to use an ITN) depends on the decisions of other members of the population because the overall vaccination or ITN usage yields the disease prevalence in the community  .
In particular, game theoretical models can be used to understand why in many cases disease eradication can be very hard  even when the cost of preventive action is very low  . There is often a great incentive to use preventative measures when a disease is at high levels, but when brought below a certain level the incremental benefit of taking the presentative action declines, and so individuals have less incentive to use it. Thus pre- vention can sometimes (but not always) become ineffective before total eradication occurs. In this paper we adapt a non-game theoretical model of ITN usage  to study the impact of individual decisions on the effectiveness of ITN usage.
2. The Model
2.1. Model of Malaria Transmission
We will use the basic model for the transmission dynamics of malaria infection as presented in  . The host population (humans) is divided into two compartments (susceptible and infectious), which are denoted by, respectively. A total host population is given by. The vector population (mosquitoes) is also divided into two compartments (susceptible and infectious) which are denoted, respectively. The total vector population is given by. When a mosquito bites a human, the susceptible human is infected by an infected mosquito with probability and a susceptible mosquito is infected by an infected human with probability. Hosts die with natural mortality rates being and can also die because of malaria with rate. Hosts recover with rate. Vectors can die with rate
where is the natural mortality rate of hosts and is the mortality rate due to ITN use in the population. Here, denotes the average ITN usage in the population. Each individual chooses its own individual ITN usage. We assume that all individuals are born as susceptible (with recruitment rates and respectively) and no infected individuals come from outside. All parameters and their values are summarized in Table 1.
Let be the force of infection for the focal susceptible host.  give the formula
where is the transmission probability per bite from infectious mosquitoes to humans, is the average
number of bites per mosquito per unit of time, is the average number of mosquitoes per host, and thus
is the average number of bites a susceptible host gets per unit time, and finally is the
Table 1. Parameters, their description and value. The cost of malaria infection varies by the country an the used treatment, see for example  -  . Over the life of the ITN, the total cost of infections is more than the cost of the ITN. Moreover, ITNs are often subsidized, so the actual cost of ITN to the individuals is only around $1 per ITN  .
proportion of infectious mosquitoes.  considers
with and which yields and (i.e. ITNs provide complete pro- tection).
Also, the average force of the infection is
In the state when malaria reaches equilibrium,  evaluates as a solution to
In  it is shown that there are potentially two stable steady states of infections and if the malaria was originally common, the stable states yields
Figure 1 shows values of for specific parameter values. Note that eventually becomes 0 (quite abruptly as seen in Figure 1(a)) as increases. This is caused by the fact that for large, the disease can be eradicated according to the model presented in  .
2.2. Game-Theoretical Model
Let be the payoff for a focal individual using () or not using () the ITN in a population where the average use of ITNs is. We shall find the ITN usage equilibrium usage level satisfying
Such a usage level is stable, because if, then and thus individuals are better off not using ITNs, and the average population usage level will thus go down (provided the individuals acts rationally). On the other hand, if, then and thus individuals are better off if they use the ITN and the population usage will thus go up. Here we assume that while individuals cannot accurately assess the overall state of the population (in particular the prevalence of malaria) instantaneously, they can change their behavior at a faster rate than the state of the population changes, in particular when the disease prevalence is either very high or very low. This seems reasonable, since individuals have a complete control over whether they use an ITN and can adapt quickly (though how they act will depend upon the information that they have available, which we consider in the Discussion).
Following  , we will evaluate
Figure 1. Values of for different vector recruitment rates. (a), (b). Values of the remaining parameters are
where is the probability of the focal individual getting infected in a population with average ITN usage and is the relative cost of the ITN and the cost of getting infected by malaria.
To evaluate, we will use and adapt the model developed in  . The probability of a host getting infected is then given by the ratio of the rate of getting infected to the rate of leaving the susceptible state for any reason. Thus,
Given a relatively small value of and a relatively large value of, we get from (15) and (5)
There will of course be values of just above 0 which are not small compared to, but for the parameters in our model (see Table 1) this happens for a very narrow range of values of, and so the above is a good approximation.
It follows from (14) that
and thus, by (16), (17),
Because by (16), for most such that, we get that
and otherwise there is a unique satisfying the equilibrium condition (13). The condition is however not realistic for real populations, and thus in practice we will always get a unique. The plots of are shown in Figure 2.
Moreover, the unique is given exactly as a solution of
where and are given by (7)-(9). As in  , (22) happens for a critical value given by
whenever the term under is nonnegative. The critical ITN usage is then given as a unique solution of a fixed point problem where
where. The plots of for different parameter values are shown at Figure 3.
Figure 2. Values of for different vector recruitment rates. (a), (b). Values of the remaining parameters are
It follows from our analysis that it is in the individual’s best interest to use the ITN as long as there is malaria present in the population. For the parameter values used above, namely, the optimal usage (from the individual’s perspective) coincide with the average usage in the population that eradicates the disease.
However, the model used shows particular sensitivity to the parameter value. This is unfortunate because the effectiveness of the insecticide sprayed on the ITN deteriorates over 3 - 5 years  , the net cannot effectively provide protection for the whole 24 hours, and the mosquitoes may adapt their biting behavior in areas of high ITN coverage  . Consequently, the realistic value of is positive and can be relatively close to. As shown in Figure 4 when (or even for larger values of), the malaria does not get eradicated. Still, the force of infection, , decreases with; and the decrease is significant. The quantity is also always positive (close to 1), thus it is in the individual’s interest to keep using the ITN. Consequently, the ITN usage does not eradicate the disease, but helps to significantly reduce its force.
However, in all of the above, we considered that malaria was transmitted between vectors and host during every bite. If we relax this assumption and consider, then again malaria can be eradicated even for larger values of as in Figure 4(c) and Figure 4(d).
In this paper we have considered a game-theoretical model of ITN usage, allowing individuals in the population to decide whether (and how often) they use ITNs. We have shown that optimal behavior from the individual’s perspective leads to the decrease (or even eradication) of the disease, i.e. an optimal outcome from the popu- lation perspective. In particular if it is not possible to eradicate the disease through ITN usage, individuals use ITNs all the time, which leads to the minimum possible level of malaria. If eradication is possible, individuals select a strategy which is approximately at the level that leads to eradication. Moreover, the optimal behavior does not depend on the exact cost of the ITN use (as long as the cost is smaller than the cost of the disease). This type of outcome is not always the case; see for example  .
We have seen that individual strategies lead to the eradication level effectively independently of the cost in our case, because individuals have a very high likelihood of getting infected by the disease whenever the disease is present in the population, due to the high rate of infection (), which is much higher than the rate of natural death, meaning that people become infected with probability almost equal to 1.
We have assumed that individuals are able to adjust their ITN usage strategy faster than changes in the underlying state of the population, such as the prevalence of the disease, to always obtain the optimal level of
Figure 3. Plots of for different vector recruitment rates. (a), (b). Values of the remaining parameters are
Figure 4. Values of for different vector recruitment rates and minimal contact rates. (a), (b), (c), (d). Values of the remaining parameters are
usage. It is clear that they are physically able to do this, because to put up or remove an ITN takes a relatively short time. However, they clearly need sufficient information about the state of the population to make a rational decision to change. What information would they have on which to base their decision? It is clear from Figure 2 that except for within a very narrow range the disease is either absent or effectively everyone has it; this information is obvious, and so any individual with the wrong strategy would quickly change it. Thus the popu- lation usage would quickly evolve to a level close to the critical level that we have assumed.
We note that because of the high infection rate, a level of ITN usage close to the eradication level but still below it could formally still lead to high levels of disease prevalence. Thus if individuals had perfect infor- mation about the state of the population they could choose their strategies to achieve the precise level of the equilibrium strategy, which would lead to some intermediate level of infection between 0 and 1. It is our assumption in this paper that they cannot do this precisely. In any case natural variations within the real popu- lation parameters, would mean that the exact value of this critical level would be subject to some fluctuation. The consequence is that the chosen strategy is sufficiently close to the eradication level, that there is a good chance of such eradication occurring, on the assumption that individuals play a fixed bed-net usage strategy.
In reality, it may be that the level of malaria will not have enough time to tend to 0 before individuals in the population notice the low level and give up the use of ITNs. Thus if the malaria level is significantly different to the equilibrium level from our analysis, individual behaviour will prevent it from being eliminated. Thus a more complex dynamically coupled model of the evolution of ITN usage and malaria level may be needed to decide the precise behaviour at such intermediate levels of infection.
The research was supported by the Simons Foundation grant #245400 (Jan Rychtář).
 Raghavendra, K., Barik, T.K., Niranjan Reddy, B.P., Sharma, P. and Dash, A.P. (2011) Malaria Vector Control: From Past to Future. Parasitology Research, 108, 757-779.
 Lengeler, C. (2004) Insecticide-Treated Bed Nets and Curtains for Preventing Malaria. Cochrane Database of Systematic Reviews, 2, Article ID: CD000363.
 Goodman, C.A. and Mills, A.J. (1999) The Evidence Base on the Cost-Effectiveness of Malaria Control Measures in Africa. Health Policy and Planning, 14, 301-312.
 Goodman, C.A., Coleman, P.G. and Mills, A.J. (1999) Cost-Effectiveness of Malaria Control in Sub-Saharan Africa. The Lancet, 354, 378-385.
 White, M.T., Conteh, L., Cibulskis, R. and Ghani, A.C. (2011) Costs and Cost-Effectiveness of Malaria Control Interventions—A Systematic Review. Malaria Journal, 10, 1475-2875.
 Miller, J.M., Korenromp, E.L., Nahlen, B.L. and Steketee, R.W. (2007) Estimating the Number of Insecticide-Treated Nets Required by African Households to Reach Continent-Wide Malaria Coverage Targets. JAMA, 297, 2241-2250.
 White, L.J., Maude, R.J., Pongtavornpinyo, W., Saralamba, S., Aguas, R., Van Effelterre, T., Day, N.P.J. and White, N.J. (2009) The Role of Simple Mathematical Models in Malaria Elimination Strategy Design. Malaria Journal, 8, 212.
 Chitnis, N., Schapira, A., Smith, T. and Steketee, R. (2010) Comparing the Effectiveness of Malaria Vector-Control Interventions through a Mathematical Model. The American Journal of Tropical Medicine and Hygiene, 83, 230.
 Agusto, F.B., Del Valle, S.Y., Blayneh, K.W., Ngonghala, C.N., Goncalves, M.J., Li, N.P., Zhao, R.J. and Gong, H.F. (2013) The Impact of Bed-Net Use on Malaria Prevalence. Journal of Theoretical Biology, 320, 58-65.
 Smith, J.M. and Price, G.R. (1973) The Logic of Animal Conflict. Nature, 246, 15-18.
 Bauch, C.T. and Earn, D.J.D. (2004) Vaccination and the Theory of Games. Proceedings of the National Academy of Sciences of the United States of America, 101, 13391-13394.
 Bauch, C.T., Galvani, A.P. and Earn, D.J.D. (2003) Group Interest versus Self-Interest in Smallpox Vaccination Policy. Proceedings of the National Academy of Sciences of the United States of America, 100, 10564-10567.
 Bauch, C.T. (2005) Imitation Dynamics Predict Vaccinating Behaviour. Proceedings of the Royal Society B: Biological Sciences, 272, 1669-1675.
 Galvani, A.P., Reluga, T.C. and Chapman, G.B. (2007) Long-Standing Influenza Vaccination Policy Is in Accord with Individual Self-Interest but Not with the Utilitarian Optimum. Proceedings of the National Academy of Sciences, 104, 5692-5697.
 Shim, E., Kochin, B. and Galvani, A. (2009) Insights from Epidemiological Game Theory into Gender-Specific Vaccination against Rubella. Mathematical Biosciences and Engineering: MBE, 6, 839-854.
 Shim, E., Grefenstette, J.J., Albert, S.M., Cakouros, B.E. and Burke, D.S. (2012) A Game Dynamic Model for Vaccine Skeptics and Vaccine Believers: Measles as an Example. Journal of Theoretical Biology, 295, 194-203.