Search tips
Search criteria 


Logo of plosbiolPLoS BiologySubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
PLoS Biol. 2007 March; 5(3): e42.
Published online 2007 February 20. doi:  10.1371/journal.pbio.0050042
PMCID: PMC1802755

Revisiting the Basic Reproductive Number for Malaria and Its Implications for Malaria Control

Bryan T Grenfell, Academic Editor


The prospects for the success of malaria control depend, in part, on the basic reproductive number for malaria, R 0. Here, we estimate R 0 in a novel way for 121 African populations, and thereby increase the number of R 0 estimates for malaria by an order of magnitude. The estimates range from around one to more than 3,000. We also consider malaria transmission and control in finite human populations, of size H. We show that classic formulas approximate the expected number of mosquitoes that could trace infection back to one mosquito after one parasite generation, Z 0(H), but they overestimate the expected number of infected humans per infected human, R 0(H). Heterogeneous biting increases R 0 and, as we show, Z 0(H), but we also show that it sometimes reduces R 0(H); those who are bitten most both infect many vectors and absorb infectious bites. The large range of R 0 estimates strongly supports the long-held notion that malaria control presents variable challenges across its transmission spectrum. In populations where R 0 is highest, malaria control will require multiple, integrated methods that target those who are bitten most. Therefore, strategic planning for malaria control should consider R 0, the spatial scale of transmission, human population density, and heterogeneous biting.

Author Summary

Each year malaria results in more than a million deaths. Controlling this disease involves understanding its transmission. For all infectious disease, the basic reproductive number, R 0, describes the most important aspects of transmission. This is the expected number of hosts that can trace their infection directly back to a single host after one disease generation. For vector-borne diseases, such as malaria, R 0 is given by a classic formula.

We made 121 estimates of R 0 for Plasmodium falciparum malaria in African populations. The estimates range from around one to over 3,000, providing much higher estimates than previously thought. We also show that in small human populations, R 0 approximates transmission when counting infections from mosquito to mosquito, but overestimates it from human to human.

Previous studies showed that transmission is amplified if some humans are bitten more than others. We confirm that such heterogeneous biting amplifies transmission counting from mosquito to mosquito, but it can also dampen transmission counting from human to human. Humans who are bitten most both infect a large number of mosquitoes and absorb many infectious bites.

What does this mean for control? When R 0 is in the thousands, eliminating malaria may seem impossible. If transmission from the humans who are bitten the most can be targeted, however, local elimination can still be within reach.


Each year, Plasmodium falciparum causes approximately 515 million clinical malaria cases [1] and over one million deaths [2,3]. Most malaria-related mortality and a large fraction of malaria cases occur in sub-Saharan Africa, where transmission can be very intense [4]. Strategic planning for malaria control should consider the transmission intensity of malaria, which is described by several parasitological and entomological indices (Table 1). The intensity of malaria transmission affects most aspects of malaria epidemiology and control, including the age at first infection, the fraction of a population that is infected (i.e., the parasite rate [PR]), the frequency and type of disease syndromes, the incidence of severe disease, the development and loss of functional immunity (i.e., immunity that reduces the frequency and severity of clinical symptoms), total malaria mortality, and the expected outcome of malaria control [48]. Good estimates of malaria transmission intensity are therefore necessary to compare and interpret malaria interventions conducted in different places and times and to objectively evaluate options for malaria control.

Table 1
Indices of Malaria Transmission

The basic reproductive number, R 0, has played a central role in epidemiological theory for malaria and other infectious diseases because it provides an index of transmission intensity and establishes threshold criteria. R 0 is generally defined as the expected number of hosts who would be infected after one generation of the parasite by a single infectious person who had been introduced into an otherwise naïve population [9,10]. If R 0 is greater than one, the number of people infected by the parasite increases, and if R 0 is less than one, that number declines. Thus, if sustained disease control reduces transmission intensity by a factor that exceeds R 0, the parasite will eventually be eliminated. Alternatively, the fraction of a population that would need to be protected to confer “herd immunity” and interrupt transmission is 1 − 1/R 0.

The classic formula for R 0 is based on a quantitative description of the P. falciparum life cycle [11,12] (Figure 1). It assumes that human populations are effectively infinite and that all humans are bitten at the same rate, but human populations are finite and some people are bitten by vectors more than others [13,14]. In infinite human populations, heterogeneous biting increases R 0 because those humans who are bitten most are also most likely to become infected and then, by infecting a large number of mosquitoes, to amplify transmission [15,16]. Thus, in infinite human populations, the classic formulas underestimate R 0.

Figure 1
The Life Cycle Model and R 0

Classic and neoclassic (i.e., with heterogeneous biting) formulas for R 0 describe idealized populations, where each infectious bite lands on a different host. In reality, some infectious bites land on previously infected hosts because malaria transmission is local. The spatial scale of malaria transmission is affected by vector ecology, especially the distribution of larval habitat and host-seeking behavior, human population density and distribution, and human movement [17,18]. Therefore, we reconsider R 0 in finite human populations with heterogeneous biting, where some bites reinfect humans or mosquitoes. When the number of humans is not effectively infinite, what is the expected number of infected hosts or vectors after one complete generation of the parasite? How are these expectations changed when biting is heterogeneous, and what do these ideas imply about malaria control?

Because R 0 is both an index of how well malaria spreads and a measure of the effort required to eliminate malaria, it would be the ideal index for strategic malaria control planning, but it has not been routinely recorded. Previous estimates of R 0 were made with a variety of methods, and they have a limited spatial coverage. Since each method introduces different sources of potential error and bias, the estimates are not directly comparable [10]. One method estimates each parameter in the classic and neoclassic formulas [19,20]; this is rarely done because it is technically and logistically difficult. A second method is based on the rate of increase in the number of human cases during an epidemic in an uninfected and immunologically naïve population [21,22]. Obviously, this method has limited application in most African populations, where a substantial fraction of people harbor malaria infections. Equilibrium methods, originally suggested by Macdonald and colleagues [23] (see Dietz [10] for a review), rely on mathematical models that describe the relations between R 0 and the population at the steady state. The terms of R 0 are rearranged into a set of indices that can be measured in populations where malaria is endemic, so they are most broadly applicable.

Here, we introduce new equilibrium methods for estimating R 0 that consider heterogeneous biting and factors that introduce a bias, such as sampling issues and immunity. We have used these new ideas to estimate R 0 for 121 African populations. These estimates are based on a common methodology and have a continental spatial coverage, so they provide a more useful index of malaria transmission than previous attempts, and one that is suitable for strategic planning for malaria control.


Estimating R 0

Our estimates of R 0 are based on two more commonly measured indices called the entomological inoculation rate (EIR) (E in equations), which is the average number of infectious bites received by a person in a year, and the PR (also called the parasite ratio) (X in equations), which is the prevalence of malaria infection in humans. Like other equilibrium methods, our method relies on mathematical models that define the steady state relationships between indices and parameters; these are the EIR, the PR, the vectorial capacity, V, which measures the number of infectious bites that arise from all the mosquitoes that are infected by a single infectious person on a single day [24], the infectivity of humans to mosquitoes, c, and the stability index, S, which measures the number of human bites taken by a vector during its lifetime [25]. The classical parameters and several malaria transmission indices are described in Tables 1 and and2.2. At the equilibrium, the relationship between these indices is given by a simple formula (Methods):

equation image

A simple relationship exists between R 0 and vectorial capacity. R 0 sums vectorial capacity, discounted for imperfect transmission efficiency, over the average infectious period [26,27]. In a population with heterogeneous biting, where the squared coefficient of variation in biting rates is α, R 0 is larger by the factor 1 + α, because the humans who are bitten most amplify transmission [15,16]; we call α the index of biting disparity. The relationship between R 0, vectorial capacity, and the other indices is given by the formula

equation image

These formulas are based on the classic assumptions: mosquito lifespan and the duration of human infections are assumed to be exponentially distributed, and R 0 is computed for a single parasite type (for a longer discussion of the assumptions, see the Methods).

Table 2
The Parameters

Using equation 2, estimates of annual EIR and PR from studies of 121 African populations [3], and parameter estimates from other studies, we generated 121 estimates of R 0 (Figure 2). Parameter estimates for b/r and α were taken from 91 of these studies that included only children less than 15 y old [14]. Published estimates of the stability index range from less than one up to five [9,28]; we use the estimate S ≈ 1, at the low end of published studies. For the infectivity, we use the value c = 0.5, a number that agrees with estimates from direct-feeding experiments [29].

Figure 2
R 0 Estimates for 121 African Populations

The R 0 estimates range from near one to more than 3,000. The median was 115 and the interquartile range was 30−815. These values are consistent with previous estimates, including one estimate of 1,600 [20] in Mngeza, in northwest Tanzania, and another of 2,000–5,000 [19] in Lira township, in central Uganda. Had these studies considered heterogeneous biting, they would have exceeded our highest estimates.

In an area around Madang, Papua New Guinea, where entomological surveys have shown that annual EIR is approximately 150 [30], and where our methods would suggest that R 0 is larger than 500, an estimate based on age seroprevalence was R 0 ≈ 7. The biological basis for the large discrepancy remains unresolved; one possibility is the strain theory of transmission [31].

Immunity and Sampling Bias

Equilibrium methods for estimating R 0 are based on the simple assumptions of mathematical models; the difference between these simple assumptions and variance in real populations can introduce a large bias. When biting rates are heterogeneous, for example, mosquitoes bite infected humans at a different frequency than when humans are sampled in a study. Thus, PR may be a biased measure of the probability a mosquito becomes infected after biting a human. In addition, the intensity of transmission at equilibrium may be lower than it would be in that same population without immunity; immunity would reduce the infectivity of humans to mosquitoes (i.e., transmission-blocking immunity) [29,32], or mosquitoes to humans (i.e., by clearing an infection before the stages that infect red blood cells develop). We have derived new formulas that consider these potential sources of bias, and we use them to modify the previous estimates of R 0.

When biting is heterogeneous and when there is some transmission-blocking immunity, it is necessary to introduce a new term called the sampling bias index, σ, that estimates the bias introduced by assuming that the fraction of mosquitoes that would become infected after biting a human is proportional to PR. σ is the ratio of two proportions. The numerator is the proportion of mosquitoes that become infected after biting a human, in a population at equilibrium; it is determined by EIR, by the index of biting disparity (Table 1), and by the level of transmission-blocking immunity. The denominator is the estimated PR, the proportion of humans that test positive in a study (Methods). Thus, the parameter σ encompasses several complex and poorly quantified processes, including differences in the way that human populations are “sampled” by mosquitoes and scientists, sporadic production of the infectious sexual stages during an infection (PR is an estimate of the prevalence of the noninfectious asexual stages), the reduced infectivity of humans to mosquitoes following the development of transmission-blocking immunity, and the sensitivity of the method used to detect parasites in humans.

When infectivity is estimated in a population where malaria is endemic and where there is some degree of immunity, the average infectivity of mosquitoes and humans, denoted bE and cE, respectively, may vary with EIR. The relevant parameters in the formula for R 0 are taken from populations without immunity, so infectivity estimates would be from naïve populations, b 0 and c 0. The bias introduced by transmission-blocking immunity is included in σ. A correction for infections that are cleared before patency (i.e., before the stages that infect red blood cells are detected) is found by multiplying the formulas for R 0 by the term BE = b 0/bE, which we call the susceptibility bias index.

Thus, we have a new formula for R 0:

equation image

For the same 121 estimates of annual EIR and PR, we generated new estimates of R 0 based on different assumptions about σ and BE (Figure 2). The original estimates effectively assumed that PR is a constant and unbiased index of infectivity (i.e., σ = 1) and that our estimates of susceptibility were not biased (i.e., BE = 1).

Our analysis suggests that σ is a complicated function of EIR (Figure 3; Methods). At low EIR (less than ten per year), mosquitoes sample infected individuals more efficiently than a stratified random sample of the population, so estimates of PR are biased by a factor that equals the product of infectivity and the amplification from heterogeneous biting, i.e., σ ≈ c 0(1 + α). At moderate to high EIR (10–700 per year), transmission-blocking immunity reduces the average infectivity of infectious humans to mosquitoes, and since bites on those who have the most immunity account for a large fraction of bites, PR severely overestimates infectivity at high EIR.

Figure 3
The Index of Sampling Bias, σ

When we assumed that transmission-blocking immunity develops, as illustrated in Figure 3, estimates of R 0 ranged from below one to nearly 11,000, with a median of 86 and an interquartile range of 15–1,000.

R 0 in Finite Human Populations

The extremely high estimates of R 0 raise the question of this index's interpretation in finite human populations; when R 0 exceeds the number of humans, what does R 0 actually describe? To interpret R 0, we simulated transmission in small well-mixed human populations of size H through one complete parasite generation with heterogeneous and homogeneous biting (Figure 4; Methods). Let R 0(H) denote the expected number of humans who could trace an infection back to one human, and Z 0(H) the expected number of mosquitoes who could trace an infection back to one mosquito. (To clarify our notation, R 0 is synonymous with R 0(∞), so when population sizes are effectively infinite, R 0 = R 0(∞) = Z 0(∞).) R 0, R 0(H), and Z 0(H) can differ, depending on the host population size. When these three indices don't differ, the assumption that populations are effectively infinite is reasonably good. When they differ by more than 10%, we call the populations “small.” Small populations are defined by R 0 and H, as well as the index of biting disparity, α, and the stability index, S.

Figure 4
R 0, R 0(H), and Z 0(H) in Finite Populations

When the size of the human population was small and malaria transmission was very intense, R 0(H) was limited by the number of humans; obviously, R 0(H) ≤ H. If every human received exactly the same number of bites, some of them would have remained uninfected, by chance. With stochastic biting, there would be some variance in the number of bites received by each individual, even if the expected biting rates were uniformly distributed. Since humans received multiple bites, this tended to increase the proportion of bites that were absorbed by already infected humans, thereby reducing R 0(H).

When human population sizes were effectively infinite, each infectious bite landed on a different human. In finite populations, heterogeneous biting amplifies transmission, as measured by Z 0(H), just as it does for infinite populations, because those who are bitten most infect a large number of mosquitoes [15,16]. Surprisingly, heterogeneous biting reduced R 0(H) below the expected number for homogeneous biting, especially when R 0 was large and H was low (Figure 5). Heterogeneous biting reduced R 0(H), i.e., the 20% of individuals who were bitten most also absorbed 80% of the infectious bites. Thus, a larger fraction of infectious bites were “reinfection” events; the transmission amplification associated with heterogeneous biting was nullified by a “superabsorbing” effect when those same individuals received most of the infectious bites.

Figure 5
R 0(H) and Z 0(H) in Finite Populations with Heterogeneous and Uniform Biting at Three Biting Intensities

The range of human population sizes that would be considered “small” differed for Z 0(H) and R 0(H) (Figure 5). Z 0(H) rises to R 0 much more rapidly than does R 0(H), when considered as a function of human population size, H. Some mosquitoes become infected and return to bite the same human again; reinfection of mosquitoes affects both R 0(H) and Z 0(H). The fraction of mosquitoes that are reinfected depends mainly on the stability index, S, the index of biting disparity, α, and human population size. For reasonable estimates of S and α, Z 0(H) approaches R 0 when the neighborhood includes less than 100 humans (Methods).

“Small” for R 0(H) depends on the ratio of R 0 to H and the index of biting disparity. Some reinfection of mosquitoes does reduce R 0(H), but this is a relatively unimportant effect for H > 25 (Methods). The much larger effect is reinfection of humans. Obviously, when R 0 and H are of comparable size, repeat infection of humans substantially reduces R 0(H), but when the human population is several times larger than R 0, R 0(H) ≈ R 0, because very few people receive multiple bites. As a rule of thumb, R 0(H) approaches R 0 when H > 2R 0. When 20% of people get 80% of the bites, the two measures are not close to one another until human population sizes are much larger: R 0(H) ≈ R 0 when H > 100R 0.

The asymmetry between R 0(H) and Z 0(H) as a function of R 0 and H arises because of the large difference in the number of humans infected by each mosquito and the number of mosquitoes infected by each human. Mosquitoes have short lives, typically 1–2 wk. The expected number of humans infected per mosquito—cS by our assumptions—is typically much less than three. The infectious period in humans, by contrast, stretches out over several months. The number of mosquitoes that bite a human during that time can range upwards to several thousand, limited mainly by the ratio of mosquitoes to humans. The number of mosquitoes infected by a single human can be so large that it exceeds the number of humans available to be bitten. When a large number of bites are distributed back on a limited number of humans, a substantial fraction result in reinfection.

Control in Finite Populations

The large range of R 0 estimates suggests that malaria control presents a variable challenge across Africa. At low transmission intensities, local elimination of malaria might be a practical goal. At the highest transmission intensities, classic theory suggests that transmission would need to be reduced by a factor of thousands, or that greater than 99% of hosts would need to be protected from infection. The amplification asymmetry that defines the relationship between R 0, H, R 0(H), and Z 0(H) suggests that malaria control measures set different targets depending on the control method deployed. Here, we consider the implications of the extreme variation in R 0 for control in finite populations with heterogeneous biting, where a few humans might account for a very large fraction of all infectious bites. In such populations, control measures that target those who are bitten most will tend to disproportionately reduce transmission. To explore these ideas, we simulated malaria control.

Because of differences in the way that control measures scale with human population size and alter transmission, we considered three categories of malaria control: host-based, vector-based, and mixed. Host-based methods, including antimalarial drugs or vaccines, reduce or completely neutralize transmission from hosts. Vector-based methods target vector populations in a general way: they lower the intensity of malaria transmission by reducing total vector density or adult lifespan. Mixed methods include insecticide-treated nets (ITNs) and indoor residual spraying (IRS). Like vector-based methods, they achieve their greatest effects by killing vectors, but like host-based methods, they are deployed around hosts to whom vectors are attracted.

Host-based methods include chemotherapy, chemoprophylaxis, and vaccines. Chemotherapy to clear infections would shorten the infectious period and reduce transmission. Obviously, case management does reduce the number of infectious individuals, but much larger reductions could be achieved through active detection of asymptomatically infected individuals followed by chemotherapy to clear infection. Since a person can become reinfected immediately after clearing an infection, more durable reductions would be achieved through chemoprophylaxis that completely neutralizes a host's ability to transmit. Similar effects would be also achieved through a vaccine that prevented infection, but no commercial vaccine for malaria is currently available or registered for public health use.

For perfect targeting, we simulated neutralizing that fraction of the individuals who were bitten most (Methods). With perfect targeting, herd immunity was achieved by neutralizing a relatively small fraction of hosts (Figure 6); neutralizing transmission from those who are bitten most makes the most of superabsorbing. The threshold population coverage required to confer herd immunity increased approximately linearly with the logarithm of R 0, rising from around 20% of the human population when R 0 was 50, to 50% when R 0 was 2,000, much lower than the 98% and 99.95% coverage predicted to be necessary assuming homogeneous biting in infinite populations. When hosts were randomly neutralized, much higher coverage was required to achieve herd immunity; classic theory provided a useful guide, although the coverage required to achieve herd immunity was slightly lower in small populations (Figure 6).

Figure 6
Achieving Herd Immunity with Random and Targeted Intervention

Vector-based methods include mass spraying and other methods that target adult mosquitoes or larvae with pesticides or that reduce larval habitat. Our analysis suggests that R 0 provides a fairly good measure of the factor by which transmission would have to be reduced by vector control to eliminate malaria. In very small populations, there is some repeat infection of mosquitoes regardless of R 0, but as vector control reduces R 0, repeat infection of humans becomes much less common.

Mixed methods merit a separate consideration from host- or vector-based methods because their success is often measured in terms of the proportion of hosts covered, and the effects on vector populations are more complicated than for purely vector-based control. Mixed methods reduce transmission from some hosts, but some vectors successfully feed despite ITNs or IRS, some mosquitoes are killed, and some vectors are repelled and attempt to bite again. We simulated targeted and random control with mixed methods (Methods; Figure 7). Like host-based methods, ITNs and IRS were very effective at reducing transmission when they were targeted, but the benefits also saturated after those who accounted for most of the bites were protected.

Figure 7
Changes in Transmission in Finite Populations with Heterogeneous Biting under Control by ITNs or IRS

Despite the promise of enormous reductions in R 0 through reductions in the lifespan of vectors [27], the total reductions in transmission from ITNs or IRS were limited. For the parameters we considered here, ITNs or IRS did not confer herd immunity, even with 100% coverage, for values of R 0 well below our median estimate. The maximum reductions in transmission depended on the fraction of mosquitoes that were killed or deflected by ITNs or IRS, and there is substantial uncertainty about these parameters under field conditions. The maximum reductions were also affected by the stability index, a measure of an individual vector's transmission efficiency. The same level of transmission can be generated by a very large number of inefficient vectors, or a lesser number of efficient ones. ITNs and IRS were most effective at reducing transmission from very efficient vectors, i.e., vector populations with a high stability index. The lower the stability index, the lower the potential proportional reductions in transmission.

Since it was possible to achieve most of the reductions in transmission by targeting those who are bitten most, it might be possible to reduce costs by targeting. A side effect of ITNs or IRS was that the deflected bites were redistributed, so biting increased on those members of the population who were not protected (Figure 7). ITNs and IRS lower the risk of infection to unprotected individuals in the surrounding population by depleting vector populations, deflecting bites onto nonhuman hosts, or shortening vector lifespan. Despite the lower risk of infection overall, increased biting on unprotected hosts could increase their risk. Our analysis was focused on changes in R 0(H) and Z 0(H), so it did not explicitly consider the risk of infection, as measured by either EIR or PR. To evaluate these questions, a different sort of analysis would be required.


Estimates of the basic reproductive number (R 0), the factor by which malaria transmission must be reduced through vector control in order to eliminate malaria, ranged from near one to more than 3,000 in a sample of 121 African populations. Revised estimates that considered other factors, such as sampling biases and immunity, that are potentially important but difficult to estimate suggest that the true range of R 0 is even larger.

To put these R 0 estimates in a broader context, the highest estimates of R 0 are up to a thousand times higher than estimates of R 0 for acute, directly transmitted infectious diseases [9]. However, R 0 measures the number of new cases through one complete generation of the parasite, not the rate of increase in the number of cases per day. The time for malaria to complete one generation is more than 200 d [14]. During that time, diseases with an R 0 of around two and a generation time of about 10 d, such as flu, for example, would have doubled 20 times in an effectively infinite host population and generated a million cases. Malaria generations overlap, so the number of expected cases after one disease generation would be higher than R 0, but these extremely high R 0 values do not necessarily represent a faster daily rate of increase for malaria compared to acute diseases with a much smaller R 0. The goals that R 0 values set for malaria control are high, but the longer generation times imply that there is more time for control.

Strategic Planning

R 0 is an important metric for strategic planning for malaria control because it helps to set priorities and define realistic expectations about the outcome of control. Despite the importance of R 0, it has not been commonly estimated; the new estimates presented here increase the total number of published R 0 estimates for malaria by an order of magnitude. The extremely large range of these R 0 estimates suggests that a globally defined “one-size-fits-all” malaria control strategy would be inefficient. Where R 0 is low, local elimination of malaria may be practical, even optimal. Where R 0 is in the thousands, malaria may resist elimination even after heavy investments in multiple control measures [33]. In such populations, focused research to identify important aspects of local transmission would help to target control and achieve larger reductions.

Mathematical modeling and R 0 provide a quantitative framework for strategic planning, one that can be modified to suit the local micro-epidemiology [34]. Important factors for control include the density and distribution of humans, the distribution of larval habitat, the vector species and their biting habits, and the seasonal patterns of transmission. Our analysis here suggests that the size of the local human population is also an important factor to consider, and that different methods may be effective (or cost-effective), depending on the distribution of humans and vectors.

Thus, an important factor in evaluating the success of malaria control is the spatial scale of malaria transmission, which is determined by several factors. Mosquito flight distances may be shorter when human blood meals are close to oviposition sites, so the spatial scale of transmission is codetermined by human population density, the distribution of humans and vector habitat, vector ecology, and vector behavior [17,18]. The spatial scale is also affected by the movement of humans. The formulas that link commonly measured entomological and parasitological indices to transmission intensity, and that correct these estimates for vector ecology and human population density, provide obvious opportunities for extensive mapping of malaria endemicity to help guide and rationalize control. These opportunities are explored in detail elsewhere [35].

Targeting Intervention

The large reductions in transmission from targeting control are only possible if those who are bitten most can be identified, as has been done for some vector-borne diseases [36]. The feasibility of targeting depends strongly on the underlying causes of heterogeneous biting. Potential causes include mosquito aggregation around places where adult mosquitoes emerge [17] or vectors oviposit [18]; also, some components of breath and sweat [37] and dirty linen [38,39] make some humans inherently more attractive to mosquitoes [39,40]. Other causes of differential biting include the use of bed nets, protective clothing, and repellants [41], housing quality and design [42], pregnancy [43], alcohol consumption [44], body size [45], and defensive behavior [46]. With research, some of these may be exploited to identify and target those who are bitten most, and thereby improve malaria control.

One practical idea is to target those with clinical malaria and presumptively treat their families and nearest neighbors with efficacious antimalarial drugs with antigametocidal properties (i.e., that clear the infectious stages) [47] to clear infection and reduce the local reservoir. In low transmission areas, where a large fraction of new malaria infections result in clinical malaria, such targeting has demonstrably reduced transmission [48,49]. In high transmission areas, where a lower fraction of new cases result in clinical malaria, clinical malaria in young children may provide some indication of where drug treatment would be most effectively targeted. In such areas, the required reductions in transmission intensity are unlikely to be achieved by any single control measure. Where R 0 exceeds a thousand, the additional widespread use of ITNs and supplementary targeted IRS may be required to achieve desired reductions in morbidity and mortality [33].

In small human populations, transmission may be effectively controlled by identifying those individuals who are most important for transmission and neutralizing their potential to transmit malaria. For example, consider an island that has only a few people, but many vectors. If one additional person came ashore infected with malaria, an epidemic would tend to ensue, if R 0(H) > 1. It may not be possible to control the epidemic with ITNs (i.e., because Z 0(H) [dbl greater-than sign] 1), but malaria could be rapidly eliminated by clearing the infection from these individuals and preventing new infections with chemoprophylaxis. In large human populations, malaria could be controlled by targeting the same fraction of humans, but this might represent a very large number of people, so the costs may differ dramatically relative to control measures in small populations.

Our analysis suggests that R 0 provides a reasonably good estimate of the reductions in transmission intensity that would be required to eliminate malaria through vector-based control. Obviously, the decision to invest in vector-based control depends on many considerations. Like heterogeneous biting, the heterogeneous distribution of adult emergence rates from larval habitats would affect the benefits of larval control. If most of the adult mosquitoes could be eliminated by removing a few larval breeding sites, targeting larval habitats might produce a large gain for little effort. In the extreme case, if all the mosquitoes emerged from a well, the easiest solution might be to cover the well. Since the benefits are related to the number of humans who would benefit, vector-based control will be more cost-effective when there are many humans. In large, urban populations, it might be more cost-effective to target vector populations for control, because of the simple fact that there is much less area to treat and many more people who benefit [50].

Reiterating Basic Principles

The effectiveness of various malaria control methods depends on the context of local transmission, but several general principles derived from the classic modeling efforts are germane. First, since the infectious period for malaria can be extremely long and a substantial fraction of the P. falciparum reservoir resides in asymptomatic cases, the infectious period can be shortened and the reservoir of parasites reduced by the use of antimalarial drugs. Thus, effective antimalarial drugs can be important tools for malaria control as well as for treating clinical malaria, although this does raise concerns about the spread of resistance. Second, although the intensity of malaria transmission is exquisitely sensitive to the mortality rate of adult mosquitoes [27], potential reductions in transmission intensity via manipulations of this parameter are limited by the fact that ITNs and IRS are not completely efficient; the maximum benefits depend on the fraction of mosquitoes that are killed or repelled and on aspects of the vector populations, especially the stability index. Because most of the reductions in transmission come from protecting a few humans, it is far more important to improve the killing effects of ITNs or IRS around those who are bitten most than to improve coverage on those who are bitten least; however, complete coverage and improved killing effects may be necessary to reach control goals. Finally, when host population sizes are small or transmission is very localized, targeted neutralization may be an extremely effective way to protect other people in the community from getting malaria. In some places, vector control may be an effective and cost-effective way to reduce the burden of malaria [2], and it has had some historical success [51], but it may not be cost-effective everywhere.

In some of the African populations described here, where malaria transmission is very intense, no single control measure will be sufficient. Nevertheless, if the suite of interventions appropriate for the transmission regime could be implemented at the appropriate targeted scale in many malaria-endemic nations, the malaria-related millennium development goals could be achieved well before an effective vaccine is available. Clarifying the optimal mix of interventions and how these can be mapped and optimally targeted at scale thus remains an important direction for our collective future research.

Materials and Methods

Estimate 1: The life-cycle model.

Ross developed and Macdonald modified a mathematical model for the transmission of a vector-borne disease that is a simplified quantitative description of the parasite life cycle [11,12]. The parameter names, following Macdonald's notation, are given in Table 2. The life-cycle model tracks the fraction of infected humans, X, and the fraction of infectious mosquitoes, Y, over time:

equation image

In this system of equations, the parasite persists if R 0 > 1, where

equation image

If R 0 > 1, the equilibria are given by the expressions

equation image

Since the average mosquito lifespan is short (i.e., 1/g ≈ 10–20 d), but the malaria infections in humans last months (i.e., b/r ≈ 170 d [14]), the proportion of infectious mosquitoes adjusts rapidly to the proportion of infectious humans, i.e., the sporozoite rate tracks PR when mosquito populations are constant (but see the discussions by Aron and May [52] and by Smith et al. [17]).

Thus, EIR is given by the formula

equation image

where V denotes vectorial capacity, following the original definition (see Table 2) [24]. Solving for V, we get

equation image

By our notation R 0 = bcV/r, so we can compute R 0 by solving for vectorial capacity:

equation image

Dietz [15] and Dye and Hasibeder [16] have demonstrated that R 0 is higher because of heterogeneous biting:

equation image

where α is the squared coefficient of variation of the human biting rate.

In these equations, mortality during sporogony is counted, but the delay for sporogony is not [17]. These equations give expressions for R 0 and equilibria, Xand Y, that are consistent with the simple assumptions of the classic model. These equations differ slightly from those given by Anderson and May, who write An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex001.jpg [9], but the equilibrium An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex002.jpg would not be consistent with the standard assumptions when mortality during sporogony is incorporated by setting c′ = ce gn [27]. Closely related delay equations are given by Aron and May [52]. An alternative approach incorporating a realistic incubation period was modeled by Smith et al. [17]. All these models assume constant per capita mortality for mosquitoes, and so they ignore important factors such as temperature-dependent mortality and senescence.

Macdonald et al.'s equilibrium method estimates R 0 from the force of infection [23]; usually, these estimates of h are based on the change in PR with age in cross-sectional surveys:

equation image


equation image


The estimates of b/r and α come from a nonlinear regression analysis using a model with superinfection (i.e., multiple infections) [14]. Here, the connection between that model and the life-cycle model is explained.

A generalized form of the life-cycle model tracks the fraction of the human population with some number of parasite “broods” [5355], denoted i. New broods are introduced by new infections at the happenings rate, which might depend on the number of broods present, hi, and these broods are cleared naturally, also depending on the number of broods present, ρi. The change in the fraction of uninfected humans is described by an equation:

equation image

The change in the fraction of humans that are infected with i broods is given by

equation image

This is an extremely general formulation of a model for infection, although the idea of a “brood” remains poorly defined. For different assumptions about hi and ρi, and for explicit assumptions about transmission of different broods by mosquitoes, it is possible to generate a very large number of models for infection in humans; some of these have been worked out by Dietz [56].

With a single brood, the dynamics reduce to the classical formulation. If there are a very large (effectively infinite) number of broods, then the force of infection is constant, hi = bE. For an infinite number of broods that clear independently, i.e., ρi = ir, the distribution of brood number at equilibrium is Poisson with mean bE/r [55], and the fraction infected is given by

equation image

These estimates of R 0 are based on Smith et al.'s estimate of b/r, which is based on the infinite brood and independent clearance model [14]. In turn, the formulas for R 0 consider the invasion of a population by a single brood.

Estimate 2: Immunity and heterogeneous biting.

The probability that a mosquito becomes infected, per bite, in the life-cycle model is denoted cX. In reality, transmission-blocking immunity and heterogeneous biting skew the probability that a mosquito becomes infected, per bite. Let X~ denote the probability that a mosquito becomes infected after biting a human (i.e., in the life-cycle model X~ = cX); then, infection in mosquitoes follows the equation

equation image

Following similar arguments as before, we get that vectorial capacity is given by the formula

equation image

Because of transmission-blocking immunity, infectivity of humans declines as a function of EIR, denoted cE. Similarly, immunity at the liver stage affects the average infectivity of mosquitoes, denoted bE. Since R 0 is defined for naïve populations, the formulas are based on infectivity in naïve hosts, c 0 and b 0. Following similar arguments as before:

equation image

Since our estimate of b/r may actually be an estimate of bE/r, we need to correct the estimate by the ratio BE = b 0/bE.

The bias introduced by transmission-blocking immunity depends implicitly on heterogeneous biting. With heterogeneous biting, mosquitoes bite individuals with index s at the rate sE; s is called a biting weight. Let X(s) denote the fraction of individuals with biting weight s that are infected, and let Γ(s,α) be the fraction of the population that has index s [14]. Finally, let c(sE) denote the average infectivity of humans who have a personal expected biting rate, sE. It follows that the probability a mosquito becomes infected after biting a human is

equation image

We let c(sE) = c 0 e −γsE , so that, because of the development of transmission-blocking immunity, infectivity declines in those who are bitten most. Using the Γ distribution and the equations for superinfection, as in [14], equation 19 can be solved:

equation image

Similarly, prevalence, Y, is given by [14]:

equation image

We assume that a well-designed study would estimate Y, while a mosquito sees X~.

The sampling bias index is An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex003.jpg . Using this formula, we can estimate R 0 as a function of EIR and PR:

equation image

We note that when EIR is low, σ ≈ c 0(1 + α), so this formula simplifies to the following:

equation image

Human to human in finite populations.

Here, we explore the interpretation of R 0 in finite populations, of size H. This approach is motivated by the extremely high estimates of R 0 (or R 0(∞)), which in some cases may even exceed the local human population size. Here, R 0(H) is defined as the expected number of different individual humans that can trace an infection back to a single human after one complete generation of the parasite, and Z 0(H) is the number of mosquitoes that can trace an infection back to a single mosquito.

Mathematical approaches to R 0 have evolved since Macdonald [12], and so have the definitions. We maintain the connection to Macdonald's original definition, in part, for historical continuity. Nowadays, R 0 is computed using next-generation approaches [10,57]. By those definitions, the quantity that we compute is called An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex004.jpg . Next-generation approaches are linearized approximations, and R 0 is an eigenvalue associated with asymptotic growth rates. Our reevaluation of R 0 is motivated by a different case—when R 0 and H are of comparable size—so asymptotic growth rates are not our primary interest. Since R 0(H) and Z 0(H) differ, it is possible that R 0(H) < 1, but Z 0(H) > 1. In finite populations where R 0(H) and Z 0(H) are both near one, malaria would be likely to random walk to extinction, in any case.

To compute R 0(H) or Z 0(H) in heterogeneous populations, let i = 1...H index humans in a population, and let masi denote their individual biting rates, where the distribution of biting weights, {s}, is constrained to have a mean of 1; An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex005.jpg . The proportion of bites that land on the ith individual is therefore si/H.

First, we compute the number of infected humans, per human. While infectious, the ith human receives masi/r bites. Each mosquito biting the individual becomes infected with probability c, but some fraction of these mosquitoes return to bite the ith human again, so we need to discount multiple infection of mosquitoes. The fraction of bites on the ith human is si/H, so a short time after the ith human has become infected, the fraction of mosquitoes that were already infected by that human is

equation image

The proportion of those bites that infect a different mosquito is

equation image

Note that more than 90% of bites are new infections when H > 9cS, so reinfection of mosquitoes is a relatively small effect when H > 25. Thereafter, the mosquito survives to become infectious with probability ce−gn, and then is expected to give a/g infectious bites.

Thus, the total number of infectious bites that arise from the ith human is

equation image

The jth human in that population is expected to be bitten at the rate Zisj/H, and each bite causes an infection with probability b. Thus, the probability that the jth individual remains uninfected is

equation image

If the ith person is the index case, the expected number of infected humans is

equation image

There are two reasonable expectations to be computed. First is the unweighted expectation:

equation image

The second is the weighted expectation:

equation image

We prefer this second, weighted expectation because it reflects heterogeneous biting, because those who are bitten most are most likely to be the index case, and because in the infinite human population limit, it converges to the formula for R 0 derived by Dietz [15] and Dye and Hasibeder [16]. Note that α is the squared CV of {s} and that An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex006.jpg .

Mosquito to mosquito in finite populations.

From a single infectious mosquito, the expected number of bites that produce an infection is bS. The probability that the ith person becomes infected is

equation image

Thereafter, that person gets masi/r bites before clearing an infection.

The number of infected humans is

equation image

A fraction c of all bites infect uninfected mosquitoes. As before, some fraction of mosquitoes are already infected. We consider only those infected mosquitoes that can trace their infection back to the index mosquito, so following the previous argument, the fraction of mosquito infections that are not reinfections is An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex007.jpg

Thereafter, e−gn infected mosquitoes survive to become infectious.

Therefore, the total number of infectious mosquitoes per infectious mosquito is given by the formula

equation image

The fraction of newly infected mosquitoes increases rapidly as a function of H. In a very large population, T is less than bS; more than 90% of bites are new infections when H > 9cST > 9bcS 2.

Simulated control in finite populations: Human-based methods.

When transmission from humans is neutralized by a perfect vaccine or by chemoprophylaxis, infected humans continue to absorb bites, but don't infect any mosquitoes. We construct a vector of length H where Vj = 0 if an individual is protected, and Vj = 1 otherwise. With targeted protection,

equation image

Note that Vj appears in the exponent to account for bites on neutralized individuals. To compute R 0(H) with neutralization, we compute the weighted expectation:

equation image

Here, Vi removes protected individuals from the computation—if a person is protected, then it is not possible for him to be the index case, by assumption.

Simulated control in finite populations: Pure mosquito-based methods.

After controlling vector populations, estimates of R 0(H) and Z 0(H) would be computed as before, but with different estimates of m or g. It is also possible that vector control would change the distribution of biting weights, but this is not a question that we have addressed here.

Simulated control in finite populations: Mixed methods.

When humans are protected from infection by ITNs or by IRS, some fraction of the mosquitoes that attempt to bite a protected human are killed, and some fraction are diverted onto other hosts. To model both effects, we assume that the biting weights describe the probability of finding a host during each attempt, that a fraction of biting attempts on protected humans kill the mosquito each visit (denoted δ), that a fraction of mosquitoes successfully feed (ψ), and that those mosquitoes that neither die nor successfully feed fly off to begin a new search. Of these, a fraction Q finds a human, again. Let N denote the set of people who are protected, then the fraction of visits that find a protected human is An external file that holds a picture, illustration, etc.
Object name is pbio.0050042.ex008.jpg .

We ignore the delay required to find another host, and assume that the vectors instantaneously reassort themselves onto hosts until they have either died or successfully fed. The fraction of mosquitoes that die is δP at the first attempt, plus δP times all those who failed to feed the first time and again find a protected human, and so on:

equation image

Thus, [var phi] is the fraction of human feeding attempts by vectors that result in mosquito death. With ITN use, the mosquito death rate increases to g′ = g + [var phi]a. By a similar argument, the feeding rate on the ith protected host is

equation image

And the proportion of bites on the jth unprotected hosts increases to

equation image

In a finite population, we compute R 0(H) and Z 0(H) as before, with new parameters describing human feeding, mosquito mortality, and biting weights (which may not sum to one). Obviously, the success of ITNs depends on the baseline parameters, Q, δ, and ψ. Here, we simulate control for Q = 0.9, δ = 0.3, and ψ = 0.2.


We would like to thank Kevin Marsh, Pete Bull, Simon Levin, Sunetra Gupta, Jonathan Dushoff, and Leslie Real for their help and comments on earlier drafts. DLS and FEM conducted this research at the Fogarty International Center, National Institutes of Health. SIH is funded by a Senior Fellowship from the Wellcome Trust (#079091). RWS is a Wellcome Trust Principal Research Fellow (#079080) and acknowledges the support of the Kenyan Medical Research Institute.


entomological inoculation rate
indoor residual spraying
insecticide-treated net
parasite rate


Competing interests. The authors have declared that no competing interests exist.

Author contributions. DLS conceived and designed the experiments. DLS and SIH analyzed the data. DLS, FEM, RWS, and SIH wrote the paper.

Funding. This work forms part of the Malaria Atlas Project (, principally funded by the Wellcome Trust, United Kingdom.


  • Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI. The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005;434:214–217. [PMC free article] [PubMed]
  • Snow RW, Omumbo JA. Malaria. In: Jamison DT, Feachem RG, Makgoba MW, Bos ER, Baingana FK, et al., editors. Disease and mortality in sub-Saharan Africa, 2nd edition. Washington (D. C.): World Bank; 2006. pp. 195–231.
  • Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW. Urbanization, malaria transmission, and disease burden in Africa. Nat Rev Microbiol. 2005;3:81–90. [PMC free article] [PubMed]
  • Snow RW, Omumbo JA, Lowe B, Molyneux CS, Obiero JO, et al. Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa. Lancet. 1997;349:1650–1654. [PubMed]
  • Marsh K, Snow RW. Malaria transmission and morbidity. Parassitologia. 1999;41:241–246. [PubMed]
  • Snow RW, Marsh K. The consequences of reducing transmission of Plasmodium falciparum in Africa. Adv Parasitol. 2002;52:235–264. [PubMed]
  • Struik SS, Riley EM. Does malaria suffer from lack of memory? Immunol Rev. 2004;201:268–290. [PubMed]
  • Reyburn H, Mbatia R, Drakeley C, Bruce J, Carneiro I, et al. Association of transmission intensity and age with clinical manifestations and case fatality of severe Plasmodium falciparum malaria. JAMA. 2005;293:1461–1470. [PubMed]
  • Anderson RM, May RM. Infectious diseases of humans. Oxford: Oxford University Press; 1991. 757
  • Dietz K. The estimation of the basic reproduction number for infectious diseases. Stat Methods Med Res. 1993;2:23–41. [PubMed]
  • Ross R. The prevention of malaria. London: John Murray; 1910. 669
  • Macdonald G. The epidemiology and control of malaria. Oxford: Oxford University Press; 1957. 201
  • Woolhouse ME, Dye C, Etard JF, Smith T, Charlwood JD, et al. Heterogeneities in the transmission of infectious agents: Implications for the design of control programs. Proc Natl Acad Sci U S A. 1997;94:338–342. [PubMed]
  • Smith DL, Dushoff J, Snow RW, Hay SI. The entomological inoculation rate and its relation to the prevalence of Plasmodium falciparum infection in African children. Nature. 2005;438:492–495. [PMC free article] [PubMed]
  • Dietz K. Models for vector-borne parasitic diseases. Lect Notes Biomath. 1980;39:264–277.
  • Dye C, Hasibeder G. Population dynamics of mosquito-borne disease: Effects of flies which bite some people more frequently than others. Trans R Soc Trop Med Hyg. 1986;80:69–77. [PubMed]
  • Smith DL, Dushoff J, McKenzie FE. The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Biol. 2004;2:e368. doi: 10.1371/journal.pbio.0020368. [PMC free article] [PubMed]
  • Le Menach A, McKenzie FE, Flahault A, Smith DL. The unexpected importance of mosquito oviposition behaviour for malaria: Non-productive larval habitats can be sources for malaria transmission. Malar J. 2005;4:23. [PMC free article] [PubMed]
  • Davidson G, Draper CC. Field studies on some of the basic factors concerned in the transmission of malaria. Trans R Soc Trop Med Hyg. 1953;47:522–535. [PubMed]
  • Davidson G. Further studies of the basic factors concerned in the transmission of malaria. Trans R Soc Trop Med Hyg. 1955;49:339–350. [PubMed]
  • Macdonald G. Theory of the eradication of malaria. Bull World Health Organ. 1956;15:369–387. [PubMed]
  • Freeman J, Laserson KF, Petralanda I, Spielman A. Effect of chemotherapy on malaria transmission among Yanomami Amerindians: Simulated consequences of placebo treatment. Am J Trop Med Hyg. 1999;60:774–780. [PubMed]
  • Macdonald G, Cuellar CB, Foll CV. The dynamics of malaria. Bull World Health Organ. 1968;38:743–755. [PubMed]
  • Garrett-Jones C. Prognosis for interruption of malaria transmission through assessment of the mosquito's vectorial capacity. Nature. 1964;204:1173–1175. [PubMed]
  • Davidson G. Estimation of the survival of Anopheline mosquitoes in nature. Nature. 1954;174:792–793. [PubMed]
  • Dye C. Vectorial capacity: Must we measure all its components. Parasitol Today. 1986;2:203–209. [PubMed]
  • Smith DL, McKenzie FE. Statics and dynamics of malaria infection in anopheles mosquitoes. Malar J. 2004;3:13. [PMC free article] [PubMed]
  • Killeen GF, McKenzie FE, Foy BD, Schieffelin C, Billingsley PF, et al. A simplified model for predicting malaria entomologic inoculation rates based on entomologic and parasitologic parameters relevant to control. Am J Trop Med Hyg. 2000;62:535–544. [PMC free article] [PubMed]
  • Githeko AK, Brandling-Bennett AD, Beier M, Atieli F, Owaga M, et al. The reservoir of Plasmodium falciparum malaria in a holoendemic area of western Kenya. Trans R Soc Trop Med Hyg. 1992;86:355–358. [PubMed]
  • Burkot TR, Graves PM, Paru R, Wirtz RA, Heywood PF. Human malaria transmission studies in the Anopheles punctulatus complex in Papua New Guinea: Sporozoite rates, inoculation rates, and sporozoite densities. Am J Trop Med Hyg. 1988;39:135–144. [PubMed]
  • Gupta S, Trenholme K, Anderson RM, Day KP. Antigenic diversity and the transmission dynamics of Plasmodium falciparum . Science. 1994;263:961–963. [PubMed]
  • Muirhead-Thomson RC. Factors determining the true reservoir of infection of Plasmodium falciparum and Wuchereria bancrofti in a West African village. Trans R Soc Trop Med Hyg. 1954;48:208–209. [PubMed]
  • Molineaux L, Gramiccia G. The Garki project: Research on the epidemiology and control of malaria in the Sudan savanna of West Africa. Geneva: World Health Organization; 1980. 311
  • Greenwood B. The microepidemiology of malaria and its importance to malaria control. Trans R Soc Trop Med Hyg. 1989;83:S25–S29. [PubMed]
  • Hay SI, Snow RW. The Malaria Atlas Project (MAP): Developing global maps of malaria risk. PLoS Med. 2006;3:e347. doi: 10.1371/journal.pmed.0030473. [PMC free article] [PubMed]
  • Perkins SE, Cattadori IM, Tagliapietra V, Rizzoli AP, Hudson PJ. Empirical evidence for key hosts in persistence of a tick-borne disease. Int J Parasitol. 2003;33:909–917. [PubMed]
  • Mukabana WR, Takken W, Coe R, Knols BG. Host-specific cues cause differential attractiveness of Kenyan men to the African malaria vector Anopheles gambiae . Malar J. 2002;1:17. [PMC free article] [PubMed]
  • Knols BG. On human odour, malaria mosquitoes, and limburger cheese. Lancet. 1996;348:1322. [PubMed]
  • Murphy MW, Dunton RF, Perich MJ, Rowley WA. Attraction of Anopheles (Diptera: culicidae) to volatile chemicals in Western Kenya. J Med Entomol. 2001;38:242–244. [PubMed]
  • Takken W, Knols BGJ. Odor-mediated behavior of Afrotropical malaria mosquitoes. Annu Rev Entomol. 1999;44:131–157. [PubMed]
  • Srinivas G, Edwin Amalraj R, Dhanraj B. The use of personal protection measures against malaria in an urban population. Public Health. 2005;119:415–417. [PubMed]
  • Lindsay SW, Snow RW. The trouble with eaves; house entry by vectors of malaria. Trans R Soc Trop Med Hyg. 1988;82:645–646. [PubMed]
  • Ansell J, Hamilton KA, Pinder M, Walraven GE, Lindsay SW. Short-range attractiveness of pregnant women to Anopheles gambiae mosquitoes. Trans R Soc Trop Med Hyg. 2002;96:113–116. [PubMed]
  • Shirai O, Tsuda T, Kitagawa S, Naitoh K, Seki T, et al. Alcohol ingestion stimulates mosquito attraction. J Am Mosq Control Assoc. 2002;18:91–96. [PubMed]
  • Port GR, Boreham PFL, Bryan JH. The relationship of host size to feeding by mosquitoes of the Anopheles gambiae giles complex (Diptera: Culicidae) Bull Entomol Res. 1980;70:133–144.
  • Kelly DW. Why are some people bitten more than others? Trends Parasitol. 2001;17:578–581. [PubMed]
  • Sutherland CJ, Ord R, Dunyo S, Jawara M, Drakeley CJ, et al. Reduction of malaria transmission to Anopheles mosquitoes with a six-dose regimen of co-artemether. PLoS Med. 2005;2:e92. doi: 10.1371/journal.pmed.0020092. [PMC free article] [PubMed]
  • Nosten F, van Vugt M, Price R, Luxemburger C, Thway KL, et al. Effects of artesunate-mefloquine combination on incidence of Plasmodium falciparum malaria and mefloquine resistance in western Thailand: A prospective study. Lancet. 2000;356:297–302. [PubMed]
  • Barnes KI, Durrheim DN, Little F, Jackson A, Mehta U, et al. Effect of artemether-lumefantrine policy and improved vector control on malaria burden in KwaZulu–Natal, South Africa. PLoS Med. 2005;2:e330. doi: 10.1371/journal.pmed.0020330. [PubMed]
  • Conteh L, Sharp BL, Streat E, Barreto A, Konar S. The cost and cost-effectiveness of malaria vector control by residual insecticide house-spraying in southern Mozambique: A rural and urban analysis. Trop Med Int Health. 2004;9:125–132. [PubMed]
  • Killeen GF, Knols BG, Gu W. Taking malaria transmission out of the bottle: Implications of mosquito dispersal for vector-control interventions. Lancet Infect Dis. 2003;3:297–303. [PubMed]
  • Aron JL, May RM. The population dynamics of malaria. In: Anderson RM, editor. Population dynamics and infectious disease. London: Chapman and Hall; 1982. pp. 139–179.
  • Walton GA. On the control of malaria in Freetown, Sierra Leone. I. Plasmodium falciparum and Anopheles gambiae in relation to malaria occurring in infants. Ann Trop Med Parasitol. 1947;41:380–407. [PubMed]
  • Dietz K, Molineaux L, Thomas A. A malaria model tested in the African savannah. Bull World Health Organ. 1974;50:347–357. [PubMed]
  • Bailey NTJ. The biomathematics of malaria. Oxford: Oxford University Press; 1982. 210
  • Dietz K. Density dependence in parasite transmission dynamics. Parasitol Today. 1988;4:91–97. [PubMed]
  • Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990;28:365–382. [PubMed]

Articles from PLoS Biology are provided here courtesy of Public Library of Science