- Research
- Open Access

# Importance of factors determining the effective lifetime of a mass, long-lasting, insecticidal net distribution: a sensitivity analysis

- Olivier JT Briët
^{1, 2}Email author, - Diggory Hardy
^{1, 2}and - Thomas A Smith
^{1, 2}

**11**:20

https://doi.org/10.1186/1475-2875-11-20

© Briët et al; licensee BioMed Central Ltd. 2012

**Received:**15 November 2011**Accepted:**13 January 2012**Published:**13 January 2012

## Abstract

### Background

Long-lasting insecticidal nets (LLINs) reduce malaria transmission by protecting individuals from infectious bites, and by reducing mosquito survival. In recent years, millions of LLINs have been distributed across sub-Saharan Africa (SSA). Over time, LLINs decay physically and chemically and are destroyed, making repeated interventions necessary to prevent a resurgence of malaria. Because its effects on transmission are important (more so than the effects of individual protection), estimates of the lifetime of mass distribution rounds should be based on the effective length of epidemiological protection.

### Methods

Simulation models, parameterised using available field data, were used to analyse how the distribution's effective lifetime depends on the transmission setting and on LLIN characteristics. Factors considered were the pre-intervention transmission level, initial coverage, net attrition, and both physical and chemical decay. An ensemble of 14 stochastic individual-based model variants for malaria in humans was used, combined with a deterministic model for malaria in mosquitoes.

### Results

The effective lifetime was most sensitive to the pre-intervention transmission level, with a lifetime of almost 10 years at an entomological inoculation rate of two infectious bites per adult per annum (ibpapa), but of little more than 2 years at 256 ibpapa. The LLIN attrition rate and the insecticide decay rate were the next most important parameters. The lifetime was surprisingly insensitive to physical decay parameters, but this could change as physical integrity gains importance with the emergence and spread of pyrethroid resistance.

### Conclusions

The strong dependency of the effective lifetime on the pre-intervention transmission level indicated that the required distribution frequency may vary more with the local entomological situation than with LLIN quality or the characteristics of the distribution system. This highlights the need for malaria monitoring both before and during intervention programmes, particularly since there are likely to be strong variations between years and over short distances. The majority of SSA's population falls into exposure categories where the lifetime is relatively long, but because exposure estimates are highly uncertain, it is necessary to consider subsequent interventions before the end of the expected effective lifetime based on an imprecise transmission measure.

## Keywords

- LLIN
- Simulation
- Life
- Mass
- Distribution
- Lasting
- Insecticidal
- Net

## Background

Over the period 2008-2010, an estimated 290 million long-lasting insecticidal nets (LLINs) were distributed in sub-Saharan Africa [1]. LLINs reduce malaria transmission by protecting individuals from infectious bites, and by reducing the probability that a mosquito survives the extrinsic incubation period. Whereas continuous distribution through antenatal clinics is common, most LLINs are being distributed through mass campaigns, reaching a large proportion of the population at risk of malaria. Over time, after a mass distribution, the proportion of the population sleeping under an LLIN decreases. This is partly due to attrition (the loss of nets available for their intended use, *e.g*. by alternative use), but also due to new births and user fatigue adding to the unprotected population.

The effective protection of LLINs against mosquito bites also wanes as they decay physically (hole formation) and chemically (insecticide loss).

Many LLIN programmes work with the assumptions that there is little variability in the decay among nets and that they last about 3 years, at which time they need replacement. However, variability in net decay appears to be substantial and the average 'lifespan' could be considerably less than 3 years [2]. The number of LLINs remaining in households does not take the physical and chemical state of the nets into account, and a proportion of those nets may have lost considerable functionality [3]. The World Health Organization (WHO) recommends tracking the physical integrity of nets (number, size and location of holes) and the insecticidal activity, measured by knock down and killing in standard WHO cone and tunnel tests [2]. Unfortunately, little is known about how these quantities, alone or in interaction, affect personal protection, and how they could be used to define when a net is worn out and at the end of its 'useful life'.

Even if individual nets do not adequately prevent mosquitoes from inoculating the user, they may still reduce mosquito survival and thus affect transmission at the population level. This community effect is likely to be more important than personal protection in preventing inoculations [4, 5]. The timing of repeat LLIN distributions may also depend on the characteristics of the human population, in particular the transmission level and immune status. Whereas knowing the 'useful life' of individual LLINs (for which a cut-off minimum functionality would need to be defined, below which an LLIN would be declared 'dead') might facilitate planning in continuous distribution programmes. For a round of mass distributed LLINs, the 'effective lifetime', based on the duration of the malaria preventive effect at population level capturing all the effects described above, might be more useful for planning the timing of subsequent rounds.

This paper describes a simulation experiment to predict the duration of epidemiological protection offered by a mass LLIN distribution targeting the general population and identifies the factors that are important in determining it.

## Methods

The OpenMalaria modelling platform [6] is an open source C++ programme and takes scenario specification inputs in eXtensible Markup Language (XML). In this platform, stochastic individual-based models for malaria in humans are combined with a deterministic model for malaria in mosquitoes, which have been fitted to multiple field data sets [7]. For each five-day time step, data on a human population is updated via components representing new infections, parasite densities, acquired immunity, uncomplicated and severe malaria episodes, direct and indirect mortality, infectiousness to mosquitoes, and case management. Each simulated malaria infection has a distinct parasite density that varies by time step, while the malaria transmission level varies seasonally. The models can accommodate multiple mosquito species with varying periodical emergence rates, and non-human hosts [8, 9]. An ensemble of 14 model variants [10] is currently available, capturing a range of possibilities for the dynamics of malaria in humans.

For this experiment, the existing insecticide treated net (ITN) intervention model component [8, 9] was developed to include capability to model physical and chemical decay of LLINs. The effect of LLINs, depending on their physical and chemical state, on deterrence, pre-prandial and post-prandial killing of malaria vector mosquitoes (see Appendix) was parameterised, using published experimental hut data [11–14].

*i.e*. the year with the minimum number of malaria episodes), as compared to a scenario without any intervention. A sensitivity analysis was done on how the effective lifetime of a mass LLIN distribution depended on the pre-intervention entomological inoculation rate (EIR) (which was varied by scaling vector emergence), initial coverage, attrition rate, LLIN effects, and physical and chemical decay rates and other LLIN-related parameters. Each of 14 parameters or parameter groups (listed in Figure 1 and discussed in the Appendix) were varied over three values (low, central and high), while keeping all other parameters and parameter groups at their central values (see Additional file 1 and the link [15] for an intervention scenario with all parameters at central value).

Also, selected parameter and parameter group combinations that seemed likely to be interdependent or to have multiplicative effects (Figure 1), were varied. For each combination, parameter values were chosen that acted on the outcome in the same direction. This was done both with those parameter values in the combination associated with lower outcomes, and with those values associated with higher outcomes. The parameters outside the combinations were kept at central values.

The relationship between the effective lifetime and pre-intervention EIR was studied in more detail by varying it over a wider range, together with LLIN attrition half-life. The effect of initial coverage was also looked at in more detail. Because the cut-off at half the impact on malaria episodes in the definition of the effective lifetime is arbitrary, sensitivity to this was studied. Finally, whether or not the protective epidemiological impact could be sustained in hypothetical situations where coverage is sustained throughout and where nets do not decay, depending on the pre-intervention EIR and coverage, was explored.

## Results

*R0000*). At the beginning of the sixth year of the simulation, LLINs were assigned to 70% of the population. Over time, the coverage declined, the insecticide in the remaining LLINs declined, and the hole index in the remaining LLINs increased (Figure 2a). Because of variation in the rates of insecticide loss and hole formation, the curves for the mean hole index and mean insecticide content became erratic as fewer LLINs remained. In the absence of intervention, the number of episodes per person per five-day time step (Figure 2b) reflects the annual periodicity of the vector emergence. The effect of an LLIN distribution is best illustrated by the de-seasonalised trend. Immediately after LLIN assignment, the number of episodes declined, reaching a minimum after two years. The number then rose with an S-shaped curve, reaching a slightly higher level than pre-intervention, and declined gradually to the pre-intervention level. The red arrow in Figure 2b illustrates the approximate point where the impact of the LLIN distribution round is half of its maximum, and the time from distribution until that point.

The results of the sensitivity analysis are plotted in a 'skeleton' diagram (Figure 1), which is an adaptation and expansion of a tornado diagram. A skeleton diagram not only displays how strongly the outcome varies with each parameter over its tested range, with all other parameters at their central values; it also indicates, for selected parameter combinations only, how the sensitivity to a parameter is altered if other parameters on which the parameter in question might be contingent are at their extremes.

Effective lifetime appears to be particularly sensitive to three parameters: pre-intervention EIR (*annualEIR*), attrition of nets (*attritionOfNets half-life*) and half-life of insecticide decay (*insecticideDecay L*).

The skeleton diagram shows that, for example, the bar lengths for *coverage*, *deterrency*, *preprandialKillingEffect*, and *postprandrialKillingEffect* are shorter for the centre-to-low (red) bars connected to the lower extreme for *attritionOfNets half-life* (red cross), than for those connected to the central level of *attritionOfNets half-life* (blue cross).

Because the attrition rate determines in part how many nets are left (together with the initial receipt proportion, defined by the parameter *coverage*), it strongly interacts with parameters that describe the LLIN effects on the vector population. At the lower extremes of both *attritionOfNets half-life* and *preprandialKillingEffect*, the centre-to-low (red) bar for the *insecticideScalingFactor* is only slightly shorter than the centre-to-low (red) bar with both *attritionOfNets half-life* and *preprandialKillingEffect* (and all other parameters) at central level. Further, one level deeper, at the lower extreme of the *insecticideScalingFactor*, the length of the centre-to-low (red) bar for *insecticideDecay L* is shorter than that with all other parameters at central level, yet still important. From Figure 1, it is not possible to establish how much the variation in each of the parameters for *attritionOfNets half-life*, *preprandialKillingEffect*, and *insecticideScalingFactor* contributes to the shortening of the *insecticideDecay L* bar. However, it is likely that *attritionOfNets half-life* has an important role, as insecticide decay only acts on remaining nets and is thus strongly contingent on the attrition rate. At a long attrition half-life, the sensitivity to the mean insecticide decay rate and its 'downstream' parameter (the parameter for variation in insecticide decay is contingent on the mean insecticide decay rate) will be much stronger than at a shorter attrition half-life because most nets will disappear in the latter scenario before the decay rate matters much. This is similar for the mean hole formation rate and its 'downstream' parameters (the parameters for variation in the hole formation rate and the rip factor are contingent on the mean hole formation rate). The *holeRate mean* is more important at the higher extreme of *attritionOfNets half-life*, *preprandialKillingEffect*, and *holeScalingFactor* than with all parameters at central level. Nevertheless, the outcome appears to be insensitive to parameters determining physical decay (hole formation).

*R0063*,

*R0065*and

*R0068*, which vary heterogeneity in human exposure [10], yield results that are different from the 11 other variants, with longer lasting LLIN effects. This is illustrated by Figure 4, which shows how the effective lifetime outcome depends on pre-intervention EIR and model variant, at an attrition half-life of four years. If effective lifetime is defined based on the incidence of infections (Figure 4c & 4d) instead of episodes (Figure 4a & 4b), then lifetime is consistently higher for these model variants than for the base model. For episodes, however, the effective lifetime with model

*R0068*converged with that of the base model

*R0000*at high pre-intervention EIR levels, and the lifetime of models

*R0063*and

*R0065*decreased more sharply with increasing pre-intervention EIR than with the base model

*R0000*.

## Discussion

### Effective lifetime

To plan the timing of malaria control interventions after a mass LLIN distribution, it would be useful to understand how long the protective effect of the LLIN distribution will last. Following intervention, the impact on clinical malaria incidence reaches a maximum within a short period and is sustained for some time, after which the number of episodes rises steeply. Clearly, re-intervention (LLIN distributions or other) should take place before case incidence exceeds pre-intervention levels (due to the loss of acquired immunity). The 'effective lifetime', defined as the period from intervention until the reduction in incidence falls to a set cut-off proportion of the maximum impact, can be used to decide when to re-intervene. Because of the sharp rise in incidence it matters little whether the cut-off is set at 40, 50 or 60% of the maximum, and there is little interaction between the cut-off value and the effect of the pre-intervention EIR, indicating that the sensitivity analysis would probably give similar results if a cut-off other than 50% had been used. Even though optimal criteria for determining when to re-intervene would include economic factors, the use of the 'effective lifetime' thus provides a reasonably robust alternative.

### Pre-intervention EIR level

LLINs have approximately^{a} the same proportionate effect on vectorial capacity at different pre-intervention EIRs. However, the apparent rate at which malaria transmission resurges from a level lower than its steady state is strongly positively related to vectorial capacity. High resurgence rates appear to be associated with shallower minima in the annual number of episodes and in the annual number of (asymptomatic) infections, but this observation results mainly from temporal smoothing.

Depending on the pre-intervention EIR, three model variants showed very different effective lifetimes from the base model *R0000*. These are models *R0063*, *R0065*, and *R0068*, in which the age-specific susceptibility is independent of exposure (contrasting with the sigmoidal relationships with exposure in the other variants), and that include extra-Poisson variation in the probability of being bitten. Model *R0063* assigns most of this to inter-host variation; model *R0065* is intermediate, and model *R0068* assigns the variation predominantly to 'within host variation', assuming that some individuals are always more likely to get bitten than others [10]. An explanation for why these model variants have different results is outside the scope of this paper. The models R0063 and R0065, both with inter-host variation in exposure, showed a stronger relationship between effective lifetime and pre-intervention EIR than did the base model. Heterogeneity in exposure to infectious bites between individuals is highly likely, as there are differences in mosquito access to houses [16] (depending on house quality and geographical situation) and in individual attractiveness to mosquitoes [17, 18]. The effect of pre-intervention EIR is thus possibly even stronger than shown in the sensitivity analysis, based on averages of 14 model variants, 12 of which ignore inter-host variation in exposure.

For the sake of simplicity, results in the sensitivity analysis were averaged over 14 model variants. This can be problematic if there is strong interaction between the effect and the model over the range studied, as was the case with the pre-intervention EIR. Sensitivity to coverage, however, was similar for all model variants. Although the 14 model variants reflect a plausible range of models, they are not necessarily evenly spread out over this range. In order to calculate a mean effect, models should ideally be weighted based on their overall fit, and on correlations between both structure and parameter values with other variants in the ensemble. Model variants with a poor fit, and/or similar to other variants included in the analysis, should receive a low weight. To perform such a weighting is a challenging task and there is a need for more methodological development in this area to address the problem appropriately.

In all the model variants, pre-intervention EIR was varied by scaling vector emergence rates. The pre-intervention EIR could possibly also be varied by changing other variables influencing vectorial capacity, such as pre-intervention survival rate, and such a pre-intervention EIR may have a somewhat different relationship with the LLIN effective lifetime.

Vector emergence was modelled as a fixed repeating seasonal pattern, independent of the adult vector population size. If the models were to include feedback from the reduced vector population due to interventions, leading to fewer emerging mosquitoes, then longer effective lifetimes would be expected, especially at locations with lower emergence rates (where local extinction of vectors might occur). Therefore, these models give a conservative estimate of effective lifetime and its dependency on emergence rates.

### LLIN decay

Effective lifetime is highly sensitive to the attrition rate. With a short attrition half-life, nets disappear before net decay can have much impact. Hole formation and insecticidal content decay are therefore of more importance with slower attrition rates. In turn, variability in hole rate and insecticide decay rate, hole size (rip factor) and the effects of decayed nets on mosquito biting and survival, are only important when hole formation and insecticidal content decay are themselves important. Thus, the shorter the attrition half-life, the more important is its accuracy relative to that of the other net decay parameters.

Insecticide decay rate is conditional on slow attrition and one of the most important factors determining the effective lifetime of an LLIN distribution. Lifetime varies strongly between the lower half-life estimate (0.5 years, characteristic for first generation LLINs) and the central value (1.5 years, characteristic for second generation LLINs). Insecticide half-life is strongly product-dependent and the values for specific products are relatively well-defined, so, despite its importance, insecticide decay is not a major contributor to uncertainty in effective lifetime, provided that mosquitoes remain sensitive to the active compound. It remains to be studied how sensitive the effective lifetime of LLINs will be to this parameter in the presence of pyrethroid resistance [19].

Effective lifetime of an LLIN distribution was surprisingly insensitive to parameters specifying the hole formation process in the nets. It is possible that the guessed values for the *holeScalingFactor*, not based on any data, were unreasonable, and that the effect of LLINs on mosquitoes wanes much faster than presumed as the hole index increases. Even though insensitive to hole parameters, it would still be useful to have evidence-based estimates as, with the emergence and spread of pyrethroid resistance [20–22], tear resistance is expected to gain in relative importance.

### Coverage

Coverage targets can be varied relatively easily in mass distribution campaigns, thus the effects of coverage were examined in more detail. But, effective lifetime was found to be insensitive to the initial coverage, expressed as the percentage of people that had access to an LLIN at the time of distribution. Since simulated people were not grouped into households, where they could share commonly owned nets, each covered person was simulated independently with a distinct net, used every night. This precluded explicit modelling of the distinction between household ownership of nets and personal net use, net transfer patterns within families, or local protective effects shared between net users and non-users [23]. Thus the simulated coverage is equivalent to the percentage of people that used an LLIN. However, the insensitivity of effective lifetime to coverage implies that it would also be insensitive to measures of usage or familial correlations in ownership or usage.

The simulation experiment with sustained net coverage sheds some light onto what might happen if mass LLIN distributions are repeated at regular intervals, or supplemented by LLIN distribution through continuous delivery channels, such as the Expanded Program on Immunization (EPI). Although the assumption of no net attrition or decay is unrealistic, this experiment allows to distinguish the effects of attrition and decay from the transient dynamics induced by reducing exposure.

In settings with a medium to high pre-intervention EIR, after initial reduction the episode incidence will increase over time even in the absence of decay in number and state of LLINs, and reach a new equilibrium that is higher than its minimum. This is due to a reduction in the equilibrium level of acquired immunity caused by the reduction in exposure. The rate of decline in prevented episodes with sustained coverage is relatively slower than with a mass distribution with a four-year attrition half-life and central values for decay of both net coverage and physical and chemical states of the nets. For instance, at pre-intervention EIR of 128 ibpapa, with 70% coverage, for sustained coverage and for a mass distribution, the effective lifetimes were 15.8 and 3.0 years, respectively. This suggests that a decline in natural immunity has little impact on the effective lifetime of a single mass LLIN distribution.

Over a long period, lower acquired immunity levels will reduce the number of clinical episodes prevented by LLINs and will be accompanied by a shift in the age distribution of the clinical episodes. Some of the simulations suggest that prolonged low coverage levels in areas with a high pre-intervention EIR could lead to an increase in incidence that exceeds pre-intervention levels. However, many of the episodes in these simulations occur in older children or adults, and may well be milder than had they occurred at younger ages, resulting in overall lower case fatality. Nevertheless, these results suggest that a sustained high level of coverage of 80% could continue to suppress episodes at medium pre-intervention EIRs of 128 ibpapa and below, despite (complete) loss of acquired immunity. This supports high coverage targets for long-term sustained intervention planning. However, a separate simulation study involving multiple distribution mechanisms, such as antenatal care and EPI continuous distribution and repeated mass distributions, would be required before answering the question "what target coverage is most cost effective?".

## Conclusions

The strong dependency of the effective lifetime of an LLIN mass distribution on pre-intervention transmission indicates that the required distribution frequency may vary more with the local entomological situation than with LLIN quality or the characteristics of the distribution system. This highlights the need for monitoring malaria before and during intervention programmes, particularly since there are likely to be strong variations between years and over short distances. The majority of sub-Saharan Africa's population probably falls into exposure categories where the effective lifetime is relatively long [24], but because exposure estimates are highly uncertain, it is necessary to consider subsequent control measures sooner than at the end of the expected effective lifetime based on an imprecise measure of transmission.

## Endnotes

^{a}Only different death rates might influence LLIN ownership, and thus LLIN effect on the vectorial capacity.

## Appendix

### Experiment parameterisation of ITN effects in OpenMalaria and experiment parameter values

In this Appendix, for selected parameters (italicised) and parameter groups (italicised) important for this study, detailed information is given on the choice of the parameter values.

The discussion of the parameters and parameter groups is organised according to the hierarchical organisation of a scenario script. The experiment's 'central' scenario specification in the machine readable language XML is given as Additional file 1, and also on the OpenMalaria site [15]. This XML scenario is richly annotated to explain what function the parameters have. Additional information on the function of the parameters is documented in the wiki section of the project webpage [6]

#### demography

The 'Ifakara' *demography*[25] was used. The population is stationary and approximately stable: individuals move up in age group with time, and because this structure is monotonically decreasing with age, surplus individuals are out-migrated (also above *maximumAgeYrs*).

#### popSize

A population size (*popSize*) of 10,000 was used. This is a balance between computational effort (which increases with larger population size) and the level of stochasticity (which decreases with larger population size). At a size of 10,000, the effect of rounding integers in the population demography (which is noticeable below a size of 5,000) is minimal.

#### monitoring

For this study, the following output variables are relevant: *nets owned*: the total number of nets (irrespective of physical and chemical state) present in the population; *nUncomp*: the number of uncomplicated malaria episodes; *nSevere*: the number of complicated, severe malaria episodes; *nNewInfections*: the number of new infections.

#### interventions > ITN > usage value

Usage represents the proportion of time during the night that a net is used by a simulated individual. It is not the proportion of people that use a net conditional on ownership. Because, in the current model, mosquito species bite homogenously throughout the night, the *usage value* can also be interpreted as the probability that host searching occurs during the time that people who own a net are using the net. In literature, this is called the 'π_{i} value' [26]. Govella and colleagues [27] define it as "the proportion of normal exposure of unprotected humans lacking nets that occurs at times and places when net users would be protected by sleeping under them". The current model version (schema 29) allows only one *usage value* to be set, thus it is not possible to vary the π_{i} value for different species within the same scenario through the *usage value*, nor can the *usage value* be varied over net users.

In this experiment, all mosquitoes were assumed to display the same host searching behaviour, with a fixed probability to search in places during times when people are protected by a net (if they own one). This, as opposed to situations where the mosquito population might be, to a degree, divided into sub-populations which either always or never search only during times and in places when people are protected by a net. Such separate behaviour could be caused by genetics, and by learning; repeating the behaviour of whatever happened in the first feeding cycle.

Comparisons of indoor *versus* outdoor human landing catches throughout the night, combined with studies of human behaviour and the source of blood meals, give insights into the proportion of host searching mosquitoes that would encounter an ITN protected host (given ownership). However, the degree to which sub-populations exist that display different behaviours is unknown.

This parameter was varied for the sensitivity analysis. A *usage value* of 0.75 was used as the central value, assuming high endophagy (the propensity to bite indoors) and biting peaks after average bed-time of the population. The extreme low parameter value was taken as 0.5 and the extreme high value was taken as 1.0. These values are based on the π_{i} values reported by Govella and colleagues [27] and Russell and colleagues [28].

#### interventions > ITN > holeRate mean

*holeRate mean*, was set at 1.8 holes per net per year. This value was based on re-analysis of the data on distribution of the total number of holes in Olyset nets after seven years of use [29], provided by Christian Lengeler. An outline normalised histogram of this distribution is shown in Figure 11. This parameter was varied for the sensitivity analysis. The extreme low parameter value was taken as 0.9 (half of the central value) and the extreme high value was taken as 3.6 (double of the central value). The effects of these parameters are also shown in Figure 11.

#### interventions > ITN > holeRate sigma

The value of the hole formation rate is varied among nets by multiplying with a distribution factor which is log normally distributed with mean one and the standard deviation of the log transformed variable *sigma*). The distribution factor is generated by taking one sample per net from a Gaussian distribution with mean zero and standard deviation one. For each parameter (*holeRate*, *ripRate*, *insecticideDecay* rate), the same sample is multiplied by the respective *sigma* and a constant (*mu*) added such that, once exponentiated, the mean of the variable over nets is one. For *insecticideDecay* rate, this constant can be chosen freely. The transformed sample is then exponentiated to obtain the respective distribution factor. This procedure implies that the distribution of *holeRate*, *ripRate* and insecticide decay rate are supposed to be covariant: nets that are heavily used decay fast both chemically and physically, whereas nets that are gently used decay slowly both chemically and physically. There is some evidence that these are indeed associated [30].

*sigma*parameter of the distribution factor for hole formation rates was set to 0.8. This value was also based on re-analysis of the raw data on distribution of the total number of holes in Olyset nets after seven years of use [29], provided by Christian Lengeler. The

*holeRate sigma*parameter was varied for the sensitivity analysis (Figure 12). The extreme low parameter value was taken as 0.60 and the extreme high value was taken as 1.00.

#### interventions > ITN > ripRate mean

The *ripRate mean* value was set equal to the value of the *holeRate mean*. (The ripping process was assumed to be similar to the hole formation process). This parameter was thus varied for the sensitivity analysis, but not independently.

#### interventions > ITN > ripRate sigma

The *riprate sigma* value was set equal to the value of the *holeRate sigma*. (The ripping process was assumed to be similar to the hole formation process). This parameter was thus varied for the sensitivity analysis, but not independently.

#### interventions > ITN > ripFactor value

*ripFactor value*expresses how important rips are in increasing the (proportionate) hole index. A net's hole index is the hole count plus the

*ripFactor value*multiplied with the cumulative number of rips. With the central values for

*holeRate mean*,

*ripRate mean*,

*holeRate sigma*and

*ripRate sigma*, a

*ripFactor value*of 0.30 allowed to approximate the upward curve in the mean hole index shown by Kilian and colleagues [31] (See Figure 13). Based on this, the central level of the

*ripFactor value*was set to 0.30. Extreme low and extreme high values were chosen of 0.15 and 0.60, respectively.

#### interventions > ITN > initialInsecticide mu

The mean insecticide content of new nets (*initialInsecticide mu*) is set to 68.4 corresponding to the declared deltamethrin content of 68.4 mg.m^{-2} for long-lasting (incorporated into filaments) insecticidal nets according to WHO interim specification 333/LN/3 [32].

#### interventions > ITN > initialInsecticide sigma

The insecticide concentration of new nets is Gaussian distributed. The standard deviation (*sigma*) was set to 14, based on the interquartile range observed by Kilian and colleagues [31], for Permanet 2nd generation.

#### interventions > ITN > insecticideDecay L and function

*insecticideDecay function*chosen was "exponential", ${\varphi}_{t}=exp\left(-\frac{tln\left(2\right)}{L}\right)$ with

*φ*

_{ t }the proportion of the initial insecticide concentration remaining at time

*t*(in years). The

*insecticideDecay L*parameter then directly translates into the insecticide half-life in years. However, if the decay rate

*λ*= ln(2)/

*L*is heterogeneous, the mean half-life is longer (Figure 14). The

*insecticideDecay L*parameter was varied for the sensitivity analysis. The central level of the

*insecticideDecay L*for the decay rate of the insecticide in the nets was taken as 1.5, which, if combined with a central distribution factor

*insecticideDecay sigma*of 0.8, yields a mean half-life of about two years. This roughly corresponds to the decay of second generation LLINs [30, 31], The extreme low parameter value was taken as 0.5, roughly corresponding to first generation LLINs. The extreme high value was taken as 2.5.

#### interventions > ITN > insecticideDecay sigma (and mu)

*insecticideDecay mu*and

*insecticideDecay sigma*are for the distribution factor (same samples as for the

*holeRate*distribution factor). Figure 15 shows how the variation in the insecticide increases over time due to the heterogeneity in the insecticide decay rate

*λ*. Such behaviour is also apparent from data presented by Killian and colleagues [31]. The

*insecticideDecay sigma*was varied for the sensitivity analysis. The central level of

*insecticideDecay sigma*was chosen at 0.8 and

*insecticideDecay mu*was chosen such that the mean was equal to one (

*insecticideDecay mu*= -0.32 for the central value). The extreme low parameter value was taken as 0.60 and the extreme high value was taken as 1.00. Figure 16 shows the distribution of the percentage insecticide remaining at 38 months. Figure 16a looks similar to Figure 6 presented by Smith and colleagues [33] for first generation LLINs. Figure 17 shows the density of insecticide concentration in nets over time, with all parameters at central values.

#### interventions > ITN > attritionOfNets L and function and k

*φ*

_{ t }the proportion of the initial net coverage remaining at time

*t*(in years). A

*k*value of 18 was used. The smooth-compact function with this

*k*value was applied by Nakul Chitnis to data on net ownership provided by Albert Kilian (Chitnis and Kilian, personal communications). The

*L*parameter was varied for the sensitivity analysis and chosen such that 50% of nets initially distributed had disappeared after 3 (low extreme), 4 (central level) or 5 (high extreme) years. This was at L values of 15.579, 20.773 and 25.966, respectively. It should be noted that from the simulated population, which is kept at a stationary size, people are out-migrated (with their nets) due to population growth. Therefore, the attrition rate of nets per person in the simulated population is slightly higher than the attrition of nets (Figure 18); if the half-life of the attrition of nets would be infinity, with a population growth of 3.47%, the half-life of nets per person in the simulated population would be about 20 years. Population growth may thus explain part of the observed difference in attrition rates between prospective studies (cohort based) and population wide surveys.

#### interventions > ITN > anophelesParams

*PA*, die in the process with probability

*PA μ*, or encounter a host of type

*i*(

*PA*

_{ i }).

After encountering a host of type *i*, they determine whether to attack or not (in which case they continue host tracking). These terms are not separately modelled, and included in *PA*_{
i
} and *PA*, respectively. In this model, once a mosquito is determined to attack, it will either successfully feed, or die in the process. Unfed alive (UA) mosquitoes found in experimental huts are thus regarded as those that encountered a host and entered the hut but decided not to attack. This is a simplification of reality, where a mosquito may survive after unsuccessfully trying an attack. Deterrency acts on the determining phase. Deterrency is defined as one minus the relative number of affected mosquitoes (RA_{1 vs 2}) of a host of type 1, as compared to another host of type 2. The number of affected mosquitoes is calculated as the sum of fed alive (FA), fed dead (FD) and unfed dead (UD) mosquitoes. A host type that is protected by an ITN will likely have a ratio below one, relative to a similar host type without ITN protection.

A mosquito determined to attack host *i* will either succeed in inserting its proboscis with probability *PB*_{
i
} , or die during the process without inserting its proboscis, with probability *PB μ*_{
i
} . For transmission from mosquito to human, it is important that the proboscis is inserted, and not if blood feeding was successful. However, in OpenMalaria, for simplicity, only blood-fed mosquitoes are assumed to have potentially inoculated hosts with sporozoites. ITNs will have an effect on the pre-prandial killing probability *PB μ*_{
i
} , which can be approximated by the proportion of UD mosquitoes out of the total of determined (unfed dead, and fed dead or alive) mosquitoes.

After proboscis insertion, feeding takes place and the mosquito tries to escape from the host's vicinity, which is successful with probability *PC*_{
i
} , or unsuccessful and the mosquito dies in the process. ITNs will have an effect on the probability of successfully escaping the host after a blood meal, called the post-prandial killing effect. This can be approximated by the proportion of FD mosquitoes out of the total of fed (dead or alive) mosquitoes.

Having escaped the host's vicinity, the mosquito will search for an appropriate resting spot. An ITN might interfere with this (through an excito-repellent effect), but this was not modelled specifically.

After a mosquito finds a resting spot it rests, and survives with probability *PD*_{
i
} . Whereas indoor residual spraying acts during this phase, ITNs were assumed to have no influence here. Experimental hut study procedures will typically collect mosquitoes in the morning (before the resting is complete), and observe alive mosquitoes for an extended period to account for deaths during the resting phase. This is because contact with insecticide picked up during earlier phases may have a delayed effect on mortality.

Rested mosquitoes will search for an oviposition site, oviposit and start another gonotrophic cycle by host tracking.

ITNs give both personal protection by reducing the number of (infectious) bites, and reduce transmission by reducing the mosquito survival per gonotrophic cycle through increased mortality during attack and escape phases, and more time spent host tracking (with associated mortality) due to deterrency.

*e.g*.: 8 holes of 10 × 20 cm [12]; 6 holes of 10 × 5 cm [11]; 6 holes of 10 × 10 cm [14].

Host types used

Abbreviation | Host type |
---|---|

NN | Person without a net |

IT | Person with an intact treated net |

IU | Person with an intact untreated net |

HT | Person with a holed treated net |

HU | Person with a holed untreated net |

HsatT | Person with a treated net, holed to saturation |

HsatU | Person with an untreated net, holed to saturation |

θT | Person with a treated net, with undefined hole index |

θU | Person with an untreated net, with undefined hole index |

#### interventions > ITN > anophelesParams > deterrency

Relative number of affected mosquitoes by host type

Gokool and colleagues [14] and Curtis and colleagues [11] (in their second experiment) found more fed and/or dead *Anopheles gambiae* s.l. in huts with untreated, holed mosquito nets than without nets. Clearly this is an artefact caused by imperfect hut traps, with more mosquitoes escaping from huts without nets, than from huts with nets, where mosquitoes got trapped under the net. It is extremely unlikely that a holed untreated net makes a person more attractive to mosquitoes than when without a net.

The geometric mean data suggest that 0.5 might be a reasonable central value for the RA of IU and IT host types *versus* NN, and 0.67 a good central value for an HT host type. As explained above, if the HU host type is considered not to increase the number of affected mosquitoes, a central value of 1 seems reasonable. These values were used to compute the 'medium' *deterrency* parameter group setting values.

The highest RA values were 0.78, 0.79 and 0.77 for IU, HT, and IT host types, respectively. These values were used to compute the 'low' *deterrency* parameter group setting values.

The lowest RA values were 0.41, 0.55 and 0.35 for IU, HT, and IT host types, respectively. These values were used to compute the 'high' *deterrency* parameter group setting values.

*p*) and hole index (

*h*), where

*p*is scaled from zero (no insecticide) to one (maximum active surface concentration for a new ITN) and

*h*(which could be expressed as the number of holes, a composite hole index based on number and size of holes, or holed surface in cm

^{2}or as a percentage) is scaled from zero (intact) to one (badly torn net), these could be written as follows:

*holeScalingFactor*and

*insecticideScalingFactor*are scaling factors for the number of holes, and the insecticide concentration, respectively. The

*holeScalingFactor*and

*insecticideScalingFactor*describe how fast the effect of the hole index, and insecticide concentration, respectively, plateaus. For the 'medium'

*deterrency*level this is:

*deterrency*for all three modelled levels (low, medium, and high).

Although the *holeScalingFactor* and the *insecticideScalingFactor* can be specified for each ITN effect (*deterrency*, *preprandialKillingEffect* and *postprandialKillingEffect*) separately and for each mosquito (sub) population or species separately, the same *holeScalingFactor* and *insecticideScalingFactor* were used for ITN effects and species in these simulations.

*holeScalingFactor*and the

*insecticideScalingFactor*. For the

*insecticideScalingFactor*, the dose-response curve of insecticide on killing in WHO cone tests [31, 34] was used. Figure 22 shows such a relationship for deltamethrin (adapted from [34], Figure 4A, page 40), with the deltamethrin concentration expressed in mg.m

^{-2}(by multiplying the g.kg

^{-1}scale with 85/2.8). Based on the shape of the relationship between the proportion of dead mosquitoes after contact with insecticide in the WHO cone tests, depending on deltamethrin concentration, the value of 0.1 was chosen as the central value, 0.2 was chosen as the high extreme value, and 0.05 was chosen as the low extreme value. This is largely in agreement with experimental hut data presented by Lindsay and colleagues [35] for killing and deterrency. However, the solvent proved to be highly deterrent in that study, making it difficult to use the results to predict the deterrency of aged nets

*versus*new nets.

For holes, very little information was available. Carnevale and colleagues [36] published data over a range of physical damage to nets (0.5, 1 and 2% of the net surface) and number of bites, but the relationships varied over the experiments. For this work, the same range of values for the *holeScalingFactor* as used for the *insecticideScalingFactor* was chosen: 0.05, 0.1 and 0.2.

*deterrency*level, with the three modeled parameter values for the

*holeScalingFactor*and the

*insecticideScalingFactor*varied simultaneously.

#### interventions > ITN > anophelesParams > preprandialKillingEffect

In order to calculate a standardised effect, corrected for different mortality in the NN host type arm, probably due to varying environmental and experimental conditions, the following process was adopted:

Standardised pre-prandial killing

The mean proportion killed for NN host types was lowest (0.09), and the mean proportion for HU host types was somewhat higher (0.23). IU host types had a much larger proportion killed (0.66), followed by HT host types (0.73) and IT host types (0.84).

The fact that the average proportion killed for HU host types (0.23) was higher than that of NN host types (0.09) indicates that holed untreated nets do have a small effect on pre-prandial killing, despite the (severe) damage done to the nets. If the nets were damaged even more, the effect would presumably be smaller and eventually, at a saturation point, such a net (host type HsatU), would no longer impede flight of mosquitoes and the pre-prandial mortality would be the same as for NN. It is of interest to estimate the effect of the insecticide at this saturation point, thus for treated nets holed to saturation (HsatT). Presumably, the effect of an HsatT net would be similar to sleeping next to, not under an HT net. If the effect decays with the same shape function (linear or exponential) for both untreated nets and treated nets, this can be calculated.

Let *h* denote the hole index in a net in arbitrary units. The host type *i* can be described by the following letter combinations: NN, HU, IU, IT, HT, HsatT, HsatU, θT or θU (Table 1), depending on the net type the host uses. The letter θ indicates here an undefined hole state. For example, if a linear decay in effect is presumed, if *h*_{HsatU} = 1 and *h*_{IU} = 0, then *PB μ*_{θU} = *PB μ*_{IU} - (*PB μ*_{IU} - *PB μ*_{
NN
} )*h*_{θU} = 0.66 - (0.66 - 0.09)*h*_{θU}. Thus, if *PB μ*_{HsatU} = 0.23, then ${h}_{\mathsf{\text{HU}}}=\frac{-\left(0.23-0.66\right)}{0.66-0.09}=0.7544$. Then, if *h*_{HsaT} = 1, *h*_{IT} = 0, and *h*_{HT} = *h*_{HU} = 0.7544, $PB{\mu}_{\mathsf{\text{HsatU}}}=PB{\mu}_{\mathsf{\text{IT}}}-\frac{PB{\mu}_{\mathsf{\text{IT}}}-PB{\mu}_{\mathsf{\text{HT}}}}{0.7544}=0.694$. Thus the values are 0.66, 0.69, 0.84, and 0.09 for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'medium' pre-prandial killing parameter group values.

Similarly, the highest values are 0.83, 0.89, 0.93, and 0.09 for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'high' pre-prandial killing parameter group values.

The lowest values are 0.44, 0.46, 0.72, and 0.09 for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'low' pre-prandial killing parameter group values.

*PB μ*

_{ i }the probability of dying before feeding as a result of being committed to biting a host of type

*i*, and

*p*the active insecticide concentration, with

*p*= 1 the maximum active surface concentration for a new ITN, and

*p*= 0 for an untreated net, if linear relationships between

*PB μ*

_{ i }and

*p*and

*h*are assumed, these can be written as follows:

For the medium *preprandialKillingEffect* level, this is thus: $PB{\mu}_{i}=0.09+0.57\left(1-h\right)+0.60p-0.42\left(1-h\right)p$.

Note that for NN, *PB μ*_{
NN
} = 0.09, and that the preprandial killing due to a net is the difference between *PB μ*_{
1
} (where the subscript 1 indicates a host type with a net) and *PB μ*_{
NN
} .

*preprandialKillingEffect*and insecticide and holes are not linear, but decrease asymptotically as the number of holes increase, and increase asymptotically as the insecticide concentration increases. Thus, for the

*preprandialKillingEffect*, the following relationship was used:

*preprandialKillingEffect*level this is:

*holeScalingFactor*and the

*insecticideScalingFactor*, here both varied simultaneously.

#### interventions > ITN > anophelesParams > postprandialKillingEffect

Post-prandial killing

Standardised post-prandial killing

The mean post-prandial killing for the HU host type (0.08) was similar to that of the IU host type (0.09), Similarly, the mean post-prandial killing for the HT host type (0.64) was similar to that of the IT host type (0.65). For the mean post-prandial killing effect, the physical net state (intact or badly holed) thus has a minimal effect on post-prandial killing.

The HU and IU host type values were slightly lower than the value for NN, which is counter-intuitive, because nets are not expected to protect against mosquito mortality.

The values 0.10, 0.65, 0.0.65, and 0.10 are proposed for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'medium' *postprandialKillingEffect* parameter group setting values.

Similarly, the values 0.36, 0.70, 0.70, and 0.10 for IU, HsatT, IT, and HsatU host types, respectively, were used to compute the 'high' *postprandialKillingEffect* parameter group setting values.

The values 0.10, 0.50, 0.50, and 0.10 for IU, HsatT, IT, and HsatU host types, respectively, were used to compute the 'low' *postprandialKillingEffect* parameter group setting values.

*preprandialKillingEffect*, the following relationship was used:

with *PC μ*_{
i
} the probability of dying (for a mosquito) after biting as a result of the biting process on host type *i*. For example for the medium *postprandialKillingEffect* level this is:

*PC μ*

_{ i }= 0.10 + 0 + 0.55(1-exp(-

*p × insecticideScalingFactor*))+0. This is illustrated in Figure 27.

*holeScalingFactor*and the

*insecticideScalingFactor*, here both varied simultaneously.

#### interventions > ITN > timed coverage

*coverage*, which describes the proportion of people that receive a net during mass distribution of nets, was set at 0.7 (70%). This parameter was varied for the sensitivity analysis. The low extreme level was set at 0.6, and the high extreme level at 0.8. Eighty percent is a reasonable figure for mass distributions targeting the entire population (de Savigny, personal communication). The nets were not randomly distributed over the population, but spread over the population according to age-use curves (Figure 30). Data from seven countries [37] was logit-transformed, and the average over all countries was taken for each age group. A constant was added such that, when the sum of the age group specific value and constant was back transformed and multiplied with the size of the age group, the total

*coverage*was at the required level over the whole population (all age groups). The mass delivery was done in the last five-day step of the fifth year into the simulation run. The 0

*coverage*level was used as comparator.

#### interventions > importedInfections

From time step zero onwards, 10 infections per 1,000 population per year were imported by stochastically infecting individuals in the population. This was done to ensure that malaria would not be eliminated from the simulated population, which might overestimate the protective effect of an intervention.

Thus, even if an intervention provides full protection to the entire population, 1% of the population will be infected once per year. These infections do not necessarily develop into disease episodes. These could be seen as infections obtained while travelling to a malarious area.

#### healthSystem

The "Tanzania ACT" health system was used, described elsewhere [38].

#### entomology > annualEIR

In the OpenMalaria schema version used (Schema 29), during the warm-up phase of the simulation run, the mosquito emergence rate is scaled such that the average annual malaria transmission to an adult (expressed as the number infectious bites per adult per annum) is approximately equal to this parameter value. Post warm-up, in the absence of interventions, this annual EIR is approximately constant. However, interventions such as ITN distributions will affect the transmission. Hence, the *annualEIR* is referred to is the pre-intervention EIR. For the sensitivity analysis, the *annualEIR* was varied between 8, 16, and 32 ibpapa. For the in-depth analysis, also values of 2, 4, 64, 128 and 256 ibpapa were used, covering a low to high range of stable malaria transmission.

#### entomology > mode and name

The "Namawala" seasonality and relative species abundance was used [39], with a dynamic transmission mode. Briefly, there are three species (*Anopheles gambiae* s.s., *Anopheles arabiensis* and *Anopheles funestus*. *An. gambiae* s.s. and *An. arabiensis* had the same seasonality. The seasonality of the natural logarithm of the observed density was fitted with a third order Fourier series, thus with a mean term, two terms for the annual frequency and two terms for the semi-annual frequency. The other parameter values in this section are described elsewhere [8].

#### model

Here, the specification for the "base" model [7, 25, 40–42] is given. All scenarios included in the experiment were run with 14 different malaria model variants [10].

## Declarations

### Acknowledgements

We acknowledge our colleagues in the Swiss TPH Biostatistics & Computational Sciences Group, particularly Michael Tarantino, for help with the models and simulations. We also thank Christian Lengeler, Don de Savigny, Jo Lines, Sarah Hoibak and Albert Kilian for helpful discussions, and Amena Briët for language editing. Simulations were run with the help of volunteers providing computational power. This work was supported through the Malaria Modelling Project #OPP1032350 funded by the Bill and Melinda Gates Foundation (BMGF) and through the NetWorks project funded by USAID in partnership with John Hopkins Bloomberg School of Public Health.

## Authors’ Affiliations

## References

- WHO Global Malaria Programme: World Malaria Report 2010. 2010, Geneva: World Health OrganizationGoogle Scholar
- Control of Neglected Tropical Diseases WHO Pesticide Evaluation Scheme, Global Malaria Programme Vector Control Unit: Guidelines for monitoring the durability of long-lasting insecticidal mosquito nets under operational conditions. 2011, Geneva: World Health Organization, WHO/HTM/NTD/WHOPES/2011.5Google Scholar
- Kilian A: How long does a long-lasting insecticidal net last in the field?. Public Health Journal. 2010, 21: 43-47.Google Scholar
- Hawley WA, Phillips-Howard P, ter Kuile F, Terlouw DJ, Vulule JM, Ombok M, Nahlen B, Gimnig JE, Kariuki SK, Kolczak MS, Hightower AW: Community-wide effects of permethrin-treated bed nets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003, 68 (Suppl. 4): 121-127.PubMedGoogle Scholar
- Killeen GF, Smith TA: Exploring the contributions of bed nets, cattle, insecticides and excitorepellency to malaria control: a deterministic model of mosquito host-seeking behaviour and mortality. Trans R Soc Trop Med Hyg. 2007, 101: 867-880. 10.1016/j.trstmh.2007.04.022.PubMed CentralView ArticlePubMedGoogle Scholar
- OpenMalaria: a simulator of malaria epidemiology and control. [http://code.google.com/p/openmalaria/]
- Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M: Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of
*Plasmodium falciparum*malaria: overview. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 1-10.PubMedGoogle Scholar - Chitnis N, Smith T, Steketee R: A mathematical model for the dynamics of malaria in mosquitoes feeding on a heterogeneous host population. J Biol Dyn. 2008, 2: 259-285. 10.1080/17513750701769857.View ArticlePubMedGoogle Scholar
- Chitnis N, Hardy D, Smith T: A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012, [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/pubmed/22218880]Google Scholar
- Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, Brooks A, Penny M, Tanner M: Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012, 9: e1001157-10.1371/journal.pmed.1001157.PubMed CentralView ArticlePubMedGoogle Scholar
- Curtis CF, Myamba J, Wilkes TJ: Various pyrethroids on bednets and curtains. Mem Inst Oswaldo Cruz. 1992, 87 (Suppl 3): 363-370.View ArticlePubMedGoogle Scholar
- Lines JD, Curtis CF, Myamba J, Njau R: Tests of repellent or insecticide impregnated curtains, bednets and anklets against malaria vectors in Tanzania. 1985, World Health Organization, WHO/VBC/85.920Google Scholar
- Lines JD, Myamba J, Curtis CF: Experimental hut trials of permethrin-impregnated mosquito nets and eave curtains against malaria vectors in Tanzania. Med Vet Entomol. 1987, 1: 37-51. 10.1111/j.1365-2915.1987.tb00321.x.View ArticlePubMedGoogle Scholar
- Gokool S, Smith DF, Curtis CF: The use of PCR to help quantify the protection provided by impregnated bednets. Parasitol Today. 1992, 8: 347-350. 10.1016/0169-4758(92)90072-A.View ArticlePubMedGoogle Scholar
- Central scenario. [http://openmalaria.googlecode.com/svn/application_deployment/examples/centre_scenario.xml]
- Smith T, Charlwood JD, Takken W, Tanner M, Spiegelhalter DJ: Mapping the densities of malaria vectors within a single village. Acta Trop. 1995, 59: 1-18. 10.1016/0001-706X(94)00082-C.View ArticlePubMedGoogle Scholar
- Qiu YT, Smallegange RC, Van Loon JJ, Ter Braak CJ, Takken W: Interindividual variation in the attractiveness of human odours to the malaria mosquito
*Anopheles gambiae*s.s. Med Vet Entomol. 2006, 20: 280-287. 10.1111/j.1365-2915.2006.00627.x.View ArticlePubMedGoogle Scholar - Lindsay SW, Adiamah JH, Miller JE, Pleass RJ, Armstrong JR: Variation in attractiveness of human subjects to malaria mosquitoes (Diptera: Culicidae) in The Gambia. J Med Entomol. 1993, 30: 368-373.View ArticlePubMedGoogle Scholar
- Briët O, N'Guessan R, Hardy D, Rowland M: Epidemiological effectiveness of insecticide treated nets in the presence of pyrethroid resistance [abstract]. Trop Med Int Health. 2011, 16 (Suppl. 1): 140-Google Scholar
- Kawada H, Dida GO, Ohashi K, Komagata O, Kasai S, Tomita T, Sonye G, Maekawa Y, Mwatele C, Njenga SM, Mwandawiro C, Minakawa N, Takagi M: Multimodal pyrethroid resistance in malaria vectors,
*Anopheles gambiae*s.s.,*Anopheles arabiensis*, and*Anopheles funestus*s.s. in western Kenya. PLoS ONE. 2011, 6: e22574-10.1371/journal.pone.0022574.PubMed CentralView ArticlePubMedGoogle Scholar - Chandre F, Darrier F, Manga L, Akogbeto M, Faye O, Mouchet J, Guillet P: Status of pyrethroid resistance in
*Anopheles gambiae*sensu lato. Bull World Health Organ. 1999, 77: 230-234.PubMed CentralPubMedGoogle Scholar - Doannio JM, sou-Yovo J, Diarrassouba S, Chauvancy G, Darriet F, Chandre F, Henry MC, Nzeyimana I, Guillet P, Carnevale P: Efficacité des moustiquaires pré-imprégnées de perméthrine Olyset Net
^{®}en zone de résistance des vecteurs aux pyréthrinoïdes. I--Evaluation entomologique. Med Trop (Mars). 1999, 59: 349-354.Google Scholar - Lines JD, Myamba J, Curtis CF: Experimental hut trials of permethrin-impregnated mosquito nets and eave curtains against malaria vectors in Tanzania. Med Vet Entomol. 1987, 1: 37-51. 10.1111/j.1365-2915.1987.tb00321.x.View ArticlePubMedGoogle Scholar
- Maire N, Shillcutt S, Walker DG, Tediosi F, Smith T: Cost effectiveness of the introduction of a pre-erythrocytic malaria vaccine into the Expanded Program on Immunization in sub-Saharan Africa: analysis of uncertainties using a stochastic individual-based simulation model of
*Plasmodium falciparum*malaria. Value Health. 2012, 14 (8): 1028-1038.View ArticleGoogle Scholar - Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L: A model for natural immunity to asexual blood stages of
*Plasmodium falciparum*malaria in endemic areas. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 19-31.PubMedGoogle Scholar - Killeen GF, Kihonda J, Lyimo E, Oketch FR, Kotas ME, Mathenge E, Schellenberg JA, Lengeler C, Smith TA, Drakeley CJ: Quantifying behavioural interactions between humans and mosquitoes: evaluating the protective efficacy of insecticidal nets against malaria transmission in rural Tanzania. BMC Infect Dis. 2006, 6: 161-10.1186/1471-2334-6-161.PubMed CentralView ArticlePubMedGoogle Scholar
- Govella NJ, Okumu FO, Killeen GF: Insecticide-treated nets can reduce malaria transmission by mosquitoes which feed outdoors. Am J Trop Med Hyg. 2010, 82: 415-419. 10.4269/ajtmh.2010.09-0579.PubMed CentralView ArticlePubMedGoogle Scholar
- Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF: Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania. Malar J. 2011, 10: 80-10.1186/1475-2875-10-80.PubMed CentralView ArticlePubMedGoogle Scholar
- Tami A, Mubyazi G, Talbert A, Mshinda H, Duchon S, Lengeler C: Evaluation of Olyset insecticide-treated nets distributed seven years previously in Tanzania. Malar J. 2004, 3: 19-10.1186/1475-2875-3-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Kilian A, Byamukama W, Pigeon O, Gimnig J, Atieli F, Koekemoer LL, Protopopoff N: Evidence for a useful life of more than three years for a polyester-based long-lasting insecticidal mosquito net in Western Uganda. Malar J. 2011, 10: 299-10.1186/1475-2875-10-299.PubMed CentralView ArticlePubMedGoogle Scholar
- Kilian A, Byamukama W, Pigeon O, Atieli F, Duchon S, Phan C: Long-term field performance of a polyester-based long-lasting insecticidal mosquito net in rural Uganda. Malar J. 2008, 7: 49-10.1186/1475-2875-7-49.PubMed CentralView ArticlePubMedGoogle Scholar
- Deltamethrin long-lasting (incorporated into filaments) insecticidal net. WHO specifications and evaluations for public health pesticides. 2010, World Health OrganizationGoogle Scholar
- Smith SC, Joshi UB, Grabowsky M, Selanikio J, Nobiya T, Aapore T: Evaluation of bednets after 38 months of household use in northwest Ghana. Am J Trop Med Hyg. 2007, 77 (6 Suppl): 243-248.PubMedGoogle Scholar
- Report of the 12th WHOPES working group meeting 8-11 December 2008. 2009, World Health Organization, WHO/HTM/NTD/WHOPES/2009.1Google Scholar
- Lindsay SW, Adiamah JH, Miller JE, Armstrong-Schellenberg JR: Pyrethroid-treated bednet effects on mosquitoes of the
*Anopheles gambiae*complex in The Gambia. Med Vet Entomol. 1991, 5: 477-483. 10.1111/j.1365-2915.1991.tb00576.x.View ArticlePubMedGoogle Scholar - Carnevale P, Bitsindou P, Diomande L, Robert V: Insecticide impregnation can restore the efficiency of torn bed nets and reduce man-vector contact in malaria endemic areas. Trans R Soc Trop Med Hyg. 1992, 86: 362-364. 10.1016/0035-9203(92)90219-3.View ArticlePubMedGoogle Scholar
- Otten M, Lines JD: Where did the LLINs go?. An analysis of data from 7 countries with the most recent surveys (2008-2009). 2011, Slide 22.Google Scholar
- Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, Tanner M: An approach to model the costs and effects of case management of
*Plasmodium falciparum*malaria in sub-Saharan Africa. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 90-103.PubMedGoogle Scholar - Smith T, Charlwood JD, Kihonda J, Mwankusye S, Billingsley P, Meuwissen J, Lyimo E, Takken W, Teuscher T, Tanner M: Absence of seasonal variation in malaria parasitaemia in an area of intense seasonal transmission. Acta Trop. 1993, 54: 55-72. 10.1016/0001-706X(93)90068-M.View ArticlePubMedGoogle Scholar
- Ross A, Maire N, Molineaux L, Smith T: An epidemiologic model of severe morbidity and mortality caused by
*Plasmodium falciparum*. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 63-73.PubMedGoogle Scholar - Maire N, Tediosi F, Ross A, Smith T: Predictions of the epidemiologic impact of introducing a pre-erythrocytic vaccine into the expanded program on immunization in sub-Saharan Africa. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 111-118.PubMedGoogle Scholar
- Penny MA, Maire N, Studer A, Schapira A, Smith TA: What should vaccine developers ask? Simulation of the effectiveness of malaria vaccines. PLoS ONE. 2008, 3: e3193-10.1371/journal.pone.0003193.PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.