 Research
 Open Access
 Published:
Modelbased analysis of experimental data from interconnected, rowconfigured huts elucidates multifaceted effects of a volatile chemical on Aedes aegypti mosquitoes
Parasites & Vectorsvolume 11, Article number: 365 (2018)
Abstract
Background
Insecticides used against Aedes aegypti and other disease vectors can elicit a multitude of dosedependent effects on behavioral and bionomic traits. Estimating the potential epidemiological impact of a product requires thorough understanding of these effects and their interplay at different dosages. Volatile spatial repellent (SR) products come with an additional layer of complexity due to the potential for altered movement of affected mosquitoes and diffusion of volatile particles of the product beyond the treated house. Here, we propose a paired experimental design and statistical inference framework for estimating these nuanced effects of volatile SRs.
Method
We fitted a continuoustime Markov chain model in a Bayesian framework to data on marked mosquitoes released in interconnected experimental huts conducted in Iquitos, Peru. We estimated the effects of two dosages of transfluthrin on Ae. aegypti behaviors associated with humanvector contact: repellency, exiting and knockdown in the treated space and in connected, adjacent huts. We validated the framework using simulated data.
Results
The odds of a female Ae. aegypti being repelled, and thus prevented from entering a treated hut (H_{T}), increased at both dosages (low dosage: odds = 1.64, 95% highest density interval (HDI) = 1.30–2.09; high dosage: odds = 1.35, HDI = 1.04–1.67). The relative risk of exiting from the treated hut was reduced (low: RR = 0.70, HDI = 0.62–1.09; high: RR = 0.70, HDI = 0.40–1.06), with this effect carrying over to untreated spaces two huts away from the treated hut (H_{2}) (low: RR = 0.79, HDI = 0.59–1.01; high: RR = 0.66, HDI = 0.50–0.87). Knockdown rates were increased in both treated and downstream huts, particularly under high dosage (H_{T}: RR = 8.37, HDI = 2.11–17.35; H_{1}: RR = 1.39, HDI = 0.52–2.69; H_{2}: RR = 2.22, HDI = 0.96–3.86).
Conclusions
Our framework is effective at elucidating multiple effects of volatile chemicals used in SR products, as well as their downstream effects. For the examined formulations of transfluthrin, we found notable dosedependent effects on repellency, movement and knockdown that carry over to adjacent, untreated spaces.
Background
Insecticidal strategies against adult mosquitoes have been used extensively in the control of mosquitoborne diseases [1]. However, certain mosquito behaviors, such as outdoor and daytime biting, challenge the efficacy of traditional control tools like insecticide treated nets (ITNs) and indoor residual spraying (IRS) [2]. The evolution of physiological resistance to insecticides [3] and behavioral adaptation of mosquitoes [4, 5] also pose limitations to the effectiveness of such products.
The effect of vector control products often goes beyond their acute lethal effects. For example, ITNs can elicit knockdown with potential for mosquito recovery and can divert mosquitoes away from a protected human to alternate hosts [6,7,8]. Volatile chemicals such as transfluthrin and metofluthrin can be delivered in high dosages and result in high lethality or can be formulated at lower dosages where acute toxicity is attenuated and other sublethal effects, such as incapacitated hostattraction and/or bloodfeeding, are elicited instead, as was described previously for residual pyrethroids [9, 10]. Currently, the term “spatial repellency” is used to describe a range of behaviors that products with volatile chemicals, including spatial repellents (SR), may invoke [11], including repellency (reduced entry), diversion (productinduced movement to an untreated space), knockdown, mortality, irritancy (increase in exiting), attractioninhibition and reduced biting [12,13,14]. These modes of action can have a concerted impact on disease transmission on an individual and community level [15,16,17,18].
Markreleaserecapturetype (MRR) experimental hut studies, in which marked mosquitoes are released in or outside an experimental hut system and monitored using traps and knockdown rates, offer unique opportunities to elucidate dosage and behavioral effects of SRs by measuring lethality, repellency and irritancy of a target vector species [10, 19,20,21]. However, studies such as these have not yet provided the granularity required to disentangle distinct behavioral and bionomic effects. The primary challenge associated with the design and interpretation of these studies is that each mosquito is only observed once: when knocked down or when trapped in entry or exit traps. This leaves movement trajectories in between release and recapture locations unobserved, making it challenging to quantify the relative contributions of multiple competing effects that could account for observed individuallevel outcomes under a multitude of equally plausible scenarios. One recent study [22] showed that even short periods of transient exposure to volatile SRs can have significant, and sometimes delayed, effects on vectors. Such unobserved effects may compromise traditional statistical analyses.
Models used for the analysis of MRR data have a long history in ecology [23,24,25,26,27]. Originally developed to estimate survival probabilities and population sizes [28], they are now increasingly being used to inform spatial processes [29]. These models partition animal movement trajectories into states (e.g. breeding or foraging), with multistate MRR models accounting for the probability of the animal occupying any of the possible states at a given time. Given sufficient information from sampling at multiple points in time and appropriate model constraints, these models can be extended for parameter estimation in the presence of unobserved states [30]. Bayesian methods are increasingly being applied to these types of problems given their treatment of all quantities as random variables [31,32,33]. These methods allow for formal treatment and quantification of parameter uncertainty and they allow researchers to explicitly build on previous studies.
Here, we present a new approach to infer effects of SR products on adult female Ae. aegypti by developing a hierarchical Bayesian model and applying it to a MRR study uniquely designed to measure SR productinduced effects resulting from a single treated home amongst a series of adjacent, untreated homes. The design mimicking an anticipated SR operational scenario in a typical semiurban, dengueendemic environment whereby product is implemented at < 100% household coverage and houses are adjoined, thereby forming shared interior spaces. This analytical approach was necessitated by the fact that the MRR study was conducted in a set of five interconnected huts, meaning that the time spent by individual mosquitoes in each of the five huts could not be directly observed. We first demonstrate the accuracy of this approach using data simulated under the same design as in our field experiments. We then demonstrate the dosedependency of knockdown, repellency and exiting effects of the SR in both treated and untreated huts. We discuss the potential use of this framework to inform the projected impact and implementation of SRs and other vector control tools with volatile chemicals.
Methods
The study was performed at the Instituto Veterinario de Investigaciones in Iquitos, Peru (73.2°W, 7.3°S) during MayJune, 2011.
Product
Technical grade transfluthrin (SigmaAldrich 105 Co. LLC, St. Louis, MO, USA), a volatile pyrethroid insecticide, was applied to cotton strips at 1/16th (0.0025 g/m^{2}) and 1/8th (0.005 g/m^{2}) dilutions of a field application rate (0.04 g/m^{2}) using previously established protocols [34]. Control strips of matched cotton material were treated with acetone alone. Cotton material was applied to the interior walls of a designated experimental hut at 25% surface area coverage using magnets and metal frames [34].
Experimental huts
A unique experimental hut configuration was used for this study. Five independent structures were positioned adjacent to one another in a single row creating adjoining walls (Figs. 1, 2). Side eave gaps of each structure were open, subsequently allowing a continuum of indoor space among all five huts available for mosquito and volatile chemical movement. This design was deliberate to mimic the housing configuration common in the study location in Iquitos, Peru and other dengueendemic areas (i.e. semiurban settings in resourcepoor, tropical areas). The huts measured 4 × 6 m and had 2 m high sidewalls. Each sidewall had an eave space 40 cm high and 6 m long (onefifth of wall surface) (Fig. 2). The eave passage between adjoining huts was fitted with netting on all sides except those facing the interior hut space. In essence, mosquitoes could move between hut interiors but not enter the exterior space between huts. Each hut had two windows (one each on the front and back walls) equipped with exit traps and each hut had two doors (one of each on the front and back walls). All window and door openings were equipped with exit intervention traps. The two outermost huts had additional eave exit intervention traps on sidewalls (Fig. 2). Hut construction materials and structural design were based on previous MRR hut studies [19, 20].
Mosquitoes
Female Ae. aegypti test populations (F_{1}F_{2} generations) of 5–7 daysold were reared from fieldcollected larvae following previously established protocols [35]. Mosquitoes were not bloodfed but were provided with cotton soaked with sucrose solution until 24 h before being released in the experimental huts. Prior to release, five cohorts each with 25 female mosquitoes each were marked with a unique color of fluorescent powder that corresponded to a single, specific experimental hut in which a cohort was released. This marking facilitated monitoring of mosquito movement from release origination to recapture location.
Study design
The study was performed using previously described experimental hut collection protocols [34, 36]. Three trials were performed: i) baseline (no cotton material application); ii) low transfluthrin dosage (0.0025 g/m^{2}); and iii) high transfluthrin dosage (0.005 g/m^{2}). A single trial consisted of five experimental days (i.e. five replicates). Transfluthrintreated cotton was applied to all interior walls of the center hut (H_{T}), while solventonly material (control) was applied to all remaining huts adjacent. Treatments remained fixed throughout a single experimental trial. This treatment allocation was purposive to reflect a semiurban scenario in a dengue endemic setting whereby an SR product is implemented at only one home in a group of connected houses. Movement during the baseline trial was measured prior to transfluthrinintegrated trials to monitor residual impact of treatment across trials. There was a period of at least two weeks between each trial, in which the huts were thoroughly cleaned and left ventilating with traps removed from windows, doors and eaves.
On each experimental day a single, uniquely marked mosquito test cohort was released inside each of the huts (here named according to the side and distance to the treated hut: H_{2L}, H_{1L}, H_{T}, H_{1R}, H_{2R}) at 05:30 hours by five persons. Each individual then remained in the hut, under an untreated bed net, during the 12 h evaluation replicate period to generate hostseeking cues and monitor indoor mosquito knockdown. Knockdown on the floor inside huts was monitored at the top of each hour starting from 06:00 hours until 18:00 hours. Specimens were placed into labeled holding cups, immediately removed from huts and provided access to 10% sugar solution to monitor 24 h mortality. Indoor collectors rotated among huts at the beginning of each replicate day to control for host cue bias.
All window, door and eave exit interception traps were monitored for both alive and knocked down mosquitoes every 30 minutes for a standardized 3minute collection period per trap. Trap monitoring began at 06:00 hours and continued until 18:00 hours by twoperson collection teams (five teams total) positioned outdoors. Collector teams rotated among huts at each sampling period (every 30 min) to control for observer bias and remained at a screened base station positioned 50 m from the experimental huts when not sampling to avoid influencing exit behavior due to outdoor human host cues.
At 18:00 hours, handheld Prokopack aspirators [37] were used by indoor collectors to recapture remaining mosquitoes inside each hut to calculate loss to followup. All recaptured mosquitoes (those from indoor Prokopack aspiration and from exit intervention traps) were held with access to a 10% sugar source to monitor 24 h mortality. Color codes were used to record release origin and location of recapture in a single day.
Model
A continuoustime Markov chain model was developed for the analysis of these data [38]. At any given time, mosquitoes can occupy any one of five huts (transient states: H_{2L}, H_{1L}, H_{T}, H_{1R} or H_{2R}) or have experienced one of 15 events represented by the absorbing states: X_{2L}, X_{1L}, X_{T}, X_{1R} or X_{2R} for the exit traps in each hut, K_{2L}, K_{1L}, K_{T}, K_{1R} or K_{2R} for knockdown in each hut and U_{2L}, U_{1L}, U_{T}, U_{1R} or U_{2R} for mosquitoes that were unaccounted for at the end of the experiment and were thus lost to followup at some unknown time. The infinitesimal generator matrix A contains the rates at which mosquitoes leave one state to move to another, such that a_{ ij } gives the rate at which a mosquito in state i moves to state j. These rates were assumed to be independent of time or previous trajectories; therefore, the time spent in state i before leaving follows an exponential distribution with mean \( {a}_i^{1} \) with \( {a}_i={\sum}_{j=1,j\ne i}^{20}{a}_{ij} \). Note that the rates out of the absorbing states are zero and that, given symmetry in the system, the rates for hut 2L and 2R are equivalent (likewise for 1L and 1R). Subscripts in A indicate the distance from the treatment hut. The 20 × 20 matrix A is defined as
with
and
The rate q_{ i } signifies the movement rate out of a hut. The direction of this movement depends on r_{ i } (proportion of movement directed to outdoors) and, for H_{1}, it further depends on repellency p_{1} (defined as the proportion of indoor movement from H_{ 1 }, adjacent to H_{T,} directed away from H_{T}). The knockdown rate k_{ i } is allowed to vary by hut, whereas the loss to followup rate u is assumed to be the same across huts. Hereafter, we refer to the exit rate q_{i}r_{i} as x_{i} (Fig. 3).
The dynamics of the probabilities P_{ ij }(t) of occupying any of the 20 states are governed by a system of differential equations with rates A and are known as the backward Kolmogorov differential equations [38]
From this, we can derive the rates of change in the probability of occupying a given state
By initializing this system in one hut (e.g. H_{2L} = 1 and all other states are zero at t = 0), solving for this system of differential equations gives the probability that a mosquito released in a given hut occupies a specific state at time t.
The absorbing states (i.e. X_{ i }, K_{ i } and U_{ i }) represent competing endpoints in the sense that an individual who enters one of these states is no longer capable of entering any of the other states at some future time. This is reflective of the nature of the experimental data, wherein captured mosquitoes can only be classified as knocked down or exited, but not both. The Markov chain accounts for these competing endpoints visàvis the property that the states are discrete and mutually exclusive. In addition, a mosquito released in 2L can only be knocked down in 2R conditional on having moved there prior to the knockdown event. The absence of nonzero rates to any of the absorbing states from other huts ensures this conditionality.
Likelihoods
To estimate A, we fitted eqn. (3) to the data using a likelihoodbased approach. The data collected during the experiments consisted of a set of interval and rightcensored timetoevent data. Outcome measures of interest included exiting (i.e. leaving a space), repellency (i.e. reduced entry into a treated space), knockdown, diversion (defined as the movement to a hut other than the release hut) and loss to followup, where exiting, knockdown and loss to followup are competing events (i.e. the occurrence of one event precludes the occurrence of another). The cumulative, conditional probabilities for all events observed in the experiment can be directly obtained from the solutions of eqn. (3), as detailed in eqns. (4), (5) and (6).
Intervalcensored events
Data pertaining to knockdown and exit events are intervalcensored between time points t_{1} and t_{2}, with exit events recorded at 30 minute intervals and knockdown events at hourly intervals. Given model parameter set θ, the probability that a mosquito released in H_{ rel } is observed to be knocked down in hut H at time t_{2} is
where F(t) denotes the probability that a specific event (here knockdown and movement to H) occurred in hut H by time t and S(t) denotes the survival function (i.e. the probability that no knockdown, exit or loss to followup has occurred by time t). Exit and knockdown events contain indirect information on the diversion event, namely that the mosquito has moved from its release location to the hut where the event took place before the event occurred. This condition, as illustrated by F_{div} in the denominator of eqn. (4), is implicitly accounted for within eqn. (3); hence, the absence of conditioning in the second part of eqn. (4).
Loss to followup
Of mosquitoes that are not retrieved at the end of the experiment, we know that they must have been lost to followup at some point between the start and the end of the experiment with probability
Here, the rate of loss to followup u was assumed constant across the huts.
Rightcensored data
Mosquitoes retrieved by the end of the experiment are treated as right censored. Namely, the time before knockdown, exit or loss to followup would have occurred is longer than the duration of the study, but by how much is uncertain. In addition, we know that the mosquito moved from the release hut to the hut where it was retrieved with probability
Likelihood function
The overall likelihood of the parameters given the data is equal to the product of the probabilities of each individual observation conditional on the parameters. These observations include the number of mosquitoes exited or knocked down during specific time intervals during an experiment for different release huts, event huts and experimental day, as well as numbers recaptured or lost to followup at the end of the experiment, resulting in
with i time points, j eventhuts (H_{2L}, H_{2L}, H_{T}, H_{1L}, H_{2R},), k releasehuts and l experiment days (1 to 5). Each k_{exit}, for instance, denotes the number of mosquitoes exiting from hut H observed at time t, by release hut and experiment day. The corresponding probabilities p are derived as detailed in eqns. (4), (5) and (6) and are assumed to be independent of the experiment day.
Model fitting
We used a Bayesian Markov chain Monte Carlo (MCMC) approach for parameter estimation. Using Bayes’ theorem, we define the posterior probability density of the model’s parameters (θ) given the data as
where P(θ) is the prior probability of the parameters. We utilized betadistributed priors with median 0.5 for p_{1} and median 0.25 for r_{ i } (i.e. a mosquito is three times as likely to move to an adjacent hut than to move outside), a gammadistributed prior with mean 0.02 on the movement rates q_{i} (i.e. average time before moving to another hut of 50 minutes) and uniform priors for the remaining parameters (see Table 1 for distribution parameters). Average times before exiting from each hut (1/q_{ i }) were constrained between 5 minutes and 20 hours and the average time until knockdown (1/k_{ i }) between 12 hours and 10 days [10, 19]. We explored the parameter space of θ more broadly using the MetropolisHastings algorithm.
We started from an initial parameter set θ_{1}, which was randomly sampled from uniform distributions with bounds: q: 360^{1}–30^{1}, p_{1}: 0.5–1, r_{i}: 0–0.5, k: 1400^{1}–720^{1} and u: 2000^{1}–1000^{1}. A new parameter was proposed such that λ_{2} = λ_{1}+Λ, where Λ is a random value from a truncated normal proposal distribution g with mean λ_{1} and standard deviation formulated relative to λ_{1} and selected so as to ideally have an acceptance rate between 10% and 50% [39]. Which parameter was updated at a given iteration was determined by taking a random draw from a categorical distribution with 11 categories (i.e. the number of model parameters to be estimated) and equal probabilities for each parameter. The probability for λ_{2} to be accepted depends on the likelihood of both θ_{1} and θ_{2} according to the MetropolisHastings rule as
where θ_{2} differs from θ_{1} only with respect to λ and g denotes truncated normal proposal distributions (between zero and one for each of p and r and from zero to infinity otherwise):
where Λ is normally distributed with mean λ_{1} and standard deviations corresponding to each parameter’s proposal distribution.
In the event that the acceptance probability was larger than a randomly generated uniform value between zero and one, θ_{2} was accepted into the chain. Otherwise, θ_{1} was retained. Multiple iterations of this routine were performed (n = 90,000). This process was repeated five times starting from different initial parameter sets to assess convergence using the GelmanRubin (GR) statistic [40]. The resulting chains of accepted parameters (φ), after discarding a ‘burnin’ period (10,000), were combined to represent our sample from the posterior distribution (π).
Simulation experiments
To validate the accuracy of the modelfitting algorithm, we simulated data with a known datagenerating process corresponding to our likelihood formulation and with known model parameters. Probabilities for released mosquitoes to occupy a specific state over time were derived using eqn. (3). As follows from eqns. (4–6), these probabilities are defined for interval and rightcensored events. Random draws from a multinomial distribution with the simulated probabilities and a given number of released mosquitoes were taken to simulate numbers of mosquitoes occupying each state at the time points at which sampling was simulated to occur. In general, these simulation experiments were designed to mimic features of the empirical experiments.
Ten distinct simulated parameter sets were used to validate the accuracy of our statistical inference framework. These parameter sets were sampled from across the composite parameter space θ using the Sobol algorithm [41, 42], where the same bounds to this sampling space were applied as for the prior distributions (Table 1). Data were simulated for different numbers of released mosquitoes (25: field scenario; 1000: large sample size scenario) for five replicates per parameter set and fitted to eqn. (3) as described before (n = 60,000, of which 10,000 was ‘burnin’).
Results
Validating the inference methodology
We first validated the inference framework against data simulated with the system of ordinary differential equations described in eqn. (3), with an observation process that mimicked the field experiment and with parameters reflecting the range of values in the prior distributions.
Large sample size scenario
In the large sample size scenario (five replicates with 1000 released mosquitoes each), we accurately estimated the values of all parameters used in the simulations. All true parameter values fell within the 95% highest density interval (HDI) of the estimated posterior distributions (Fig. 4). Most posterior medians approximated the true parameter well (Pearson’s r > 0.98), but somewhat less so for knockdown in the treated hut (Pearson’s r = 0.74). Posterior distributions were relatively wider for rate parameters associated with the treated hut (x_{T} and k_{T}). Standard deviations of these parameters were a fraction (i.e. 11% and 12%) of their respective medians, whereas the s.d.:median ratio was below 3.5% for all other parameters. This reduced precision may be a consequence of the fact that rate parameters associated with huts other than the treated hut were informed by twice as much data as were the rate parameters associated with the treated hut, which derives from our assumption of shared parameters for huts a given distance from the treated hut (Fig. 3). GR statistics were below 1.1 for most simulation sets (average 1.04). When simulation sets resulted in parameters with GR statistics above 1.1, these were related to mosquito movement (q_{i}, r_{i} and p_{1}) and were most commonly associated with the untreated huts (Additional file 1: Table S1). This indicates that those parameters may be among the most difficult to estimate.
Field scenario
We also tested the performance of the inference framework on data simulated with the same number of mosquitoes as were released in the experiments (Fig. 5). All true parameter values fell within the 95% HDI, but the posterior medians were less consistent with the simulated values (r > 0.8 for all but x_{T}: 0.68; k_{2}:0.03; k_{T}: 0.29) than under the large sample size scenario (Fig. 5). No systematic underestimation or overestimation was observed based on these simulations, suggesting that the additional discrepancy between simulated and inferred parameter values in the field scenario relative to the large sample size scenario was due to stochasticity associated with the smaller sample size in the field scenario (i.e. n = 25 vs n = 1000). GR statistics were, across all parameters and simulation sets, close to 1 (average GR 1.01) (Additional file 1: Table S1).
Product effects on mosquito behavior
We first fitted the Markov chain model to the experimental hut data allowing all parameters to vary. Strong correlations between r_{ i }, q_{ i } and p_{ i } indicated that these parameters were not identifiable given that a wide range of combinations of values of these parameters explained the data equally well (Additional file 1: Figure S1). To resolve this identifiability issue, we fitted the exit rate x_{ i } as a single composite parameter (q_{ i }r_{ i }). The rate of movement between huts is directly related to the exit rate; namely, it is a proportion (1r_{i}) of the overall movement rate out of a specific hut (q_{ i }). In doing so, we fixed the values of r_{ i } at the medians of the posterior marginal density of the r_{ i } corresponding to each hut that was obtained from the full parameter fit on the baseline data set (r_{T}: 0.040, r_{1}:0.036 , r_{2}: 0.088) (Additional file 1: Figure S1). This reduced the amount of crosscorrelation from Pearson’s r as high as 0.86 in the original (Additional file 1: Figure S1) to as low as 0.72 upon fixing r_{ i } (Additional file 1: Figure S2, Figure S3 and Figure S4). Most importantly, it markedly improved convergence from GR statistics as high as 1.38 (q_{ 2 }, low dosage) to as low as 1.00 for all parameters after fixing r_{ i }, indicating that other parameters became identifiable once this adjustment was made (Additional file 1: Figure S7, Figure S8, Figure S9). Choosing either the 2.5th (r_{T}: 0.031, r_{1}:0.027, r_{2}: 0.065) or 97.5th (r_{T}: 0.054, r_{1}:0.040 , r_{2}: 0.11) percentile of r_{ i } instead did not affect this conclusion (Additional file 1: Figure S5 and Figure S6). Acceptance rates for each chain tended to remain relatively constant following a ‘burnin’ period and varied across chains and parameters within the range of 21–54%.
Exit and movement rates
Under baseline conditions (no chemical), exit rates (x_{ i }) from huts at different distances i from the treatment hut were relatively similar (medians for x_{T}: 2.2 ×10^{3}, x_{1}: 1.6 ×10^{3}, x_{2}: 1.8 ×10^{3}) (Fig. 6ac). In subsequent treatment experiments, exit rates out of the treated hut were reduced relative to the baseline in response to both the low (RR = 0.70, HDI = 0.62–1.09) and the high transfluthrin dosage (RR = 0.70, HDI = 0.40–1.06), with no perceptible difference in the respective effects of the two dosages (Fig. 6c). This effect carried over to the adjacent huts (H_{1}) with exit rates lower than observed in the baseline experiment (low: RR = 0.79, HDI = 0.59–1.01; high: RR = 0.66, HDI = 0.50–0.87) (Fig. 6b). In the huts furthest from the SR application (H_{2}), the low dosage had no effect on exit rates relative to when no SR was applied (RR = 0.94, HDI = 0.72–1.18). In contrast, the high dosage reduced exit rates (RR = 0.71, HDI = 0.54–0.92) in all huts adjacent to the source of transfluthrin, including the furthest adjoining structures (Fig. 6a). Given that the proportion of movement that was directed outdoors (r_{ i }) was held constant in this exercise, these results on exit rates (q_{ i }r_{ i }) are directly proportional to movement rates (q_{ i }).
Repellency
In the baseline experiment, mosquitoes moved away from or towards the treated center hut (H_{T}) with roughly equal probability (p_{1} = 0.54, HDI = 0.48–0.59), although with a possible slight preference for movement away from H_{T} (odds of moving away = 1.16, HDI = 0.92–1.41) (Fig. 6d). In the experiment using lowdosage SR treatment, significant repellency from the treated center hut was observed (odds = 1.64, HDI = 1.30–2.09), with a median probability of moving away from this hut of 0.62 (HDI = 0.57–0.68) (Fig. 6d). In the highdosage treatment, repellency was still clear (odds = 1.35, HDI = 1.04–1.67), but the effect was somewhat smaller (p_{1} = 0.57, HDI = 0.52–0.63) (Fig. 6d).
Knockdown
Knockdown was a very rare event during baseline experiments (2/125 mosquitoes across all five replicates). As a consequence, estimates of knockdown rates in the baseline approached the lower boundary of the prior distribution (medians for H_{ T } = 5.8 × 10^{5}_{,}H_{1} = 4.4 × 10^{5}, H_{2} = 4.0 × 10^{5}) (Fig. 6eg). There was no effect of the low SR dosage on knockdown rates relative to the baseline, both in the treated hut H_{ T } (RR = 1.39, HDI = 0.26–3.84) (Fig. 6g) and in the H_{1} huts directly adjacent (RR = 1.00, HDI = 0.45–1.76) (Fig. 6f). In the H_{2} huts furthest away from the treatment, a somewhat increased knockdown rate was observed in response to the low dosage relative to the baseline (RR = 1.37, HDI = 0.64–2.46) (Fig. 6e). Knockdown rates in the highdosage scenario were elevated in all huts, in particular in the H_{ T } treatment huts (RR = 8.37, HDI = 2.11–17.35) (Fig. 6g) but also in the H_{1} and H_{2} huts (H_{1}: RR = 1.39, HDI = 0.52–2.69; H_{2}: RR = 2.22, HDI = 0.96–3.86) (Fig. 6e, f).
Loss to followup
Rates of loss to followup were similar across the baseline and two SR treatment experiments, although there was a signal for a small increase in these rates with increasing dosage (low: 5%; high: 8%) (Fig. 6h). In comparing posterior samples across dosages, a signal for a positive doseresponse relationship [i.e. u(high) > u(low) > u(baseline)] was confirmed in 61% of samples from the posterior. Loss to followup rates in the baseline scenario were an order of magnitude higher than the estimated knockdown rates. While some of mosquitoes that were lost to followup at the end of the experiment could have been knocked down, the differences in trends between treatment effects on knockdown and loss to followup indicate that this could likely only be true for a small portion of the mosquitoes.
Time spent in a hut
The total amount of time a mosquito spent in each hut results from the composite of treatment effects. By running simulations of the system of ordinary differential equations (eqn. (3)) with the estimated posterior parameter values, we derived a posterior estimate of the proportion of the time a mosquito spent in each hut relative to the total time a mosquito was in the hut system (i.e. before exit, knockdown or loss to followup). This proportion was found to be similar but slightly reduced for the treated hut H_{T} relative to the baseline scenario in either treatment scenario (Fig. 7e) and without any effect in the downstream huts H_{1} and H_{2} (Fig. 7a, c). However, when considering the total duration of the experimental day, the proportion of time spent in the adjacent, downstream huts H_{1} and H_{2} was higher during experiments using both low and high SR dosage than during baseline (Fig. 7b, d). This was a result of reduced exit rates and thus an overall increase in time spent in the hut system as a whole (Fig. 6a, c).
Discussion
Novel Ae. aegypti vector control strategies are currently being evaluated to address challenges related to dengue transmission expansion [2]. Spatial repellent (SR) products, which release volatile chemicals into treated spaces to interrupt hostvector contact, are among these [14]. One challenge for evaluating the efficacy of SRs, and other products that may include nonlethal outcomes, is characterizing the multifaceted, incapacitating effects of a given product on mosquito behavior under field conditions and across a dosage gradient, which is critical for identifying opportunities to exploit new mechanisms of action and guiding product formulation and implementation under operational conditions. To quantify concurrent and downstream (i.e. untreated spaces adjacent to the treated space) effects of a transfluthrinbased SR product on Ae. aegypti behavior, we used a continuoustime Markov chain model informed by experimental data under a Bayesian inference framework.
Examination of posterior estimates of model parameters showed that test mosquitoes were deterred from entering the experimental hut where the product was located and that this effect was stronger for the lower dosage SR, an effect of transfluthrin consistent with other studies [43, 44]. Posterior estimates of model parameters also indicated lower movement rates out of the treatment hut (either to a neighboring hut or out of the hut system) under both treatments, presumably due to confusion effects, similar to previously demonstrated effects of metofluthrin [13]. Under the higher SR dosage, the reduced outofhut movement was noticeable in adjacent untreated huts as far as two huts away from the SR application. This dosedependent effect may be a result of irritancy and disorientation, as reviewed in Bibbs & Kaufman [45]. Similarly strong effects in adjacent connected huts were observed on knockdown rates, which were markedly increased in all huts in the presence of the higher SR dosage, yet little effect on knockdown was observed at the lower dosage. These results are in agreement with dose and distancedependent effects on knockdown observed upon metofluthrin exposure in Ae. aegypti [46].
While metrics of repellency and exit rates from traditional experimental hut studies are often affected by treatment effects on mosquito knockdown, our inference framework explicitly accounts for such confounding effects. We found knockdown and loss to followup rates to be small relative to mosquito movement rates, indicating that for the formulations considered, effects on repellency and exit rates could have meaningful impacts on mosquitohuman contact rates.
We validated our inference method by demonstrating its ability to accurately estimate the model’s parameters given simulated data. This assessment was conditional, however, on the assumption that the model is a realistic representation of reality. Some of the known limiting assumptions of our analysis include (i) effects that depend on distance from the treated hut rather than on each hut individually; (ii) equal loss to followup across huts; and (iii) timeinvariant parameters. Of these, the first may be most problematic when considering that air flow within the hut system could result in asymmetric effects of transfluthrin dispersion to huts of the same distance from the treated hut but on different sides of it [47]. In principle, it would be possible to account for such factors in future studies by measuring air flow and incorporating its effect on the data through an appropriate modification of the model. For example, repellency (p_{ 1 }) could be allowed to vary across huts and treated as a function of readings from a wind gauge. Posterior estimates of the parameters governing the relationship between wind and p_{ 1 } would then allow for inferences about the repellency of the product under varying airflow conditions and beyond those observed in the current experiment. Planning for required sample sizes and sampling schemes for such experiments would benefit from our model and results by using our posterior distributions to inform prior distributions in those future studies [48].
Repellency (i.e. reduced entry into the treated hut from adjacent connected huts) and increased knockdown reduced the overall time that mosquitoes spent transiently in the treated hut, whereas decreased movement rates have the potential to offset this effect. The result of the SR’s impact on the time mosquitoes spent transiently in the treated hut indicates potential for such a product to limit humanvector contact (and thereby reduce the probability of pathogen transmission) in the treated hut. Furthermore, the reduced exit rates that we observed may be a result of confusion and disorientation and may be accompanied by impaired hostseeking [45]. However, because bloodfeeding was prohibited in our study, it is uncertain to what extent hostseeking and bloodfeeding behaviors of these mosquitoes exposed to the SR may have been affected in the current study. Other studies using similar volatile products have shown these effects to also be associated with reduced human landing [13, 49, 50]. The inclusion of bloodfeeding metrics in experiments with volatile pyrethroid products using anophelines under field conditions [51] and against the topical repellent DEET using Ae. aegypti in the laboratory [52] have been valuable in establishing expectations of such synergistic chemical effects.
The effect of SR products on untreated neighboring premises has been a consistent and critical question to the public health value of these products [14, 53]. Three aspects of our results suggest that the risk of diversion (i.e. movement of mosquitoes from a treated space to an adjacent untreated space) may be limited for the formulation used in our experiments. First, SR exposure reduced movement rates between huts. Secondly, there was a marked increase in knockdown in untreated huts at the high SR dosage. At the same time, there was also a marked reduction in exit rates out of untreated huts, which resulted in prolonged time spent in adjacent huts. Evaluating the overall potential for diversion based on these effects will require pairing experimental results such as ours with theory that is capable of translating this range of behavioral effects into estimates of their epidemiological consequences [15, 54, 55].
Under our experimental design, we cannot distinguish between downstream effects caused by volatile particles dispersed into untreated huts or by a residual, postexposure effect of transfluthrin on mosquitoes that are exposed in the treated hut and move elsewhere. Indeed, this highlights that SR postexposure effects remain understudied. A mosquito that remains alive in a treated space, and/or does not exit after being exposed to a SR product, may exhibit attractioninhibition to human hosts inside the space or may be inhibited to bloodfeed due to being chemically incapacitated. Underlying physiological mechanisms causing these responses are yet to be characterized. Correlations between air sampling measurements in experimental huts and mosquito behavior responses have been explored in previous studies using spatial repellents but with limited success due to limits of chemical detection and quantification [56,57,58]. Combining air chemistry inferences of specific active ingredients (i.e. vapor pressure or particle weight) with environmental data (i.e. air current or flow rate) into our new inference framework is therefore warranted and could enable quantification of the extent to which downstream effects result from movement of the volatile chemical or movement of exposed mosquitoes with lingering postexposure effects. The latter possibility has been indicated in other studies to have potential for innovative applications of SRs [59].
Conclusions
The need for development and efficient testing of new vector control products and innovative formulations of existing tools is evident [2, 60]. Advancing the characterization of SRs and other volatile vector control products must include understanding how all productinduced effects, independent or in combination, contribute to potential impact on pathogen transmission. Here, we advance SR product characterization methods by integrating multiple aspects of the complexity of assumed productinduced responses. The complementary experimental design and inferential framework we introduce provides a flexible approach for estimating a product’s effects on mosquito behavior in a quantitative and probabilistic fashion. Using this framework, we disentangled complex, dosedependent effects of transfluthrin on mosquito behavior. While a lowdosage treatment resulted in notable repellency and minor knockdown in Ae. aegypti mosquitoes, higher dosages presented with much higher knockdown rates yet minor repellency effects. Exposure to transfluthrin has the potential to diminish mosquito movement, host attractioninhibition and/or bloodfeeding inhibition, possibly due to confusion or disorientation of the mosquito. The use of an experimental rowhouse hut design allowed for the estimation of transfluthrin effects beyond a treated space, which is of particular interest as adjacent, interconnected houses are common in semiurban environments where SR products are under evaluation as a tool against arbovirus transmission. These findings highlight the need to estimate context and dosagespecific diversion (productinduced movement to an untreated space) and downstream (adjacent spaces) effects under scenarios similar to operational settings for which the product is intended to be used.
Abbreviations
 GR:

GelmanRubin
 HDI:

Highest density interval
 IRS:

Indoor residual spraying
 ITN:

Insecticidetreated net
 MCMC:

Markov chain Monte Carlo
 MRR:

Markreleaserecapture
 RR:

Relative risk
 SR:

Spatial repellent
References
 1.
Eldridge BF, Edman J. Medical entomology: a textbook on public health and veterinary problems caused by arthropods. Dordrecht, Netherlands: Springer; 2012.
 2.
Achee NL, Gould F, Perkins TA, Reiner RC Jr, Morrison AC, Ritchie SA, et al. A critical assessment of vector control for dengue prevention. PLoS Negl Trop Dis. 2015;9:e0003655.
 3.
Corbel V. N’Guessan. RDistribution, mechanisms, impact and management of insecticide resistance in malaria vectors: a pragmatic review. In: Manguin S, editor. Anopheles mosquitoes  new insights into malaria vectors. Rijeca: InTech; 2013. p. 633.
 4.
Gatton ML, Chitnis N, Churcher T, Donnelly MJ, Ghani AC, Godfray HCJ, et al. The importance of mosquito behavioural adaptations to malaria control in Africa. Evolution. 2013;67:1218–30.
 5.
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 insecticidetreated nets in rural Tanzania. Malar J. 2011;10:80.
 6.
Fanello C, Carneiro I, IlboudoSanogo E, CuzinOuattara N, Badolo A, Curtis C. Comparative evaluation of carbosulfanand permethrinimpregnated curtains for preventing houseentry by the malaria vector Anopheles gambiae in Burkina Faso. Med Vet Entomol. 2003;17:333–8.
 7.
Pleass R, Armstrong J, Curtis C, Jawara M, Lindsay S. Comparison of permethrin treatments for bednets in The Gambia. Bull Entomol Res. 1993;83:133–9.
 8.
Lindsay S, Adiamah J, Armstrong J. The effect of permethrinimpregnated bednets on house entry by mosquitoes (Diptera: Culicidae) in The Gambia. Bull Entomol Res. 1992;82:49–55.
 9.
Achee NL, Sardelis MR, Dusfour I, Chauhan KR, Grieco JP. Characterization of spatial repellent, contact irritant, and toxicant chemical actions of standard vector control compounds. J Am Mosq Control Assoc. 2009;25:156–67.
 10.
Grieco JP, Achee NL, Chareonviriyaphap T, Suwonkerd W, Chauhan K, Sardelis MR, Roberts DR. A new classification system for the actions of IRS chemicals traditionally used for malaria control. PLoS One. 2007;2:e716.
 11.
World Health Organization. Guidelines for efficacy testing of spatial repellents. Geneva: WHO; 2013.
 12.
Ogoma SB, Lorenz LM, Ngonyani H, Sangusangu R, Kitumbukile M, Kilalangongono M, et al. An experimental hut study to quantify the effect of DDT and airborne pyrethroids on entomological parameters of malaria transmission. Malar J. 2014;13:131.
 13.
Ritchie SA, Devine GJ. Confusion, knockdown and kill of Aedes aegypti using metofluthrin in domestic settings: a powerful tool to prevent dengue transmission. Parasit Vectors. 2013;6:262.
 14.
Achee NL, Bangs MJ, Farlow R, Killeen GF, Lindsay S, Logan JG, et al. Spatial repellents: from discovery and development to evidencebased validation. Malar J. 2012;11:164.
 15.
Killeen GF, Smith TA. Exploring the contributions of bed nets, cattle, insecticides and excitorepellency to malaria control: a deterministic model of mosquito hostseeking behaviour and mortality. Trans R Soc Trop Med Hyg. 2007;101:867–80.
 16.
Kiware SS, Chitnis N, Moore SJ, Devine GJ, Majambere S, Merrill S, Killeen GF. Simplified models of vector control impact upon malaria transmission by zoophagic mosquitoes. PLoS One. 2012;7:e37661.
 17.
Briët OJ, Smith TA, Chitnis N. Measurement of overall insecticidal effects in experimental hut trials. Parasit Vectors. 2012;5:256.
 18.
Chitnis N, Schapira A, Smith T, Steketee R. Comparing the effectiveness of malaria vectorcontrol interventions through a mathematical model. Am J Trop Med Hyg. 2010;83:230–40.
 19.
Chareonviriyaphap T, Grieco JP, Suwonkerd W, Prabaripai A, Polsomboon S, Thainchum K, et al. An improved experimental hut design for the study of Aedes aegypti (Diptera: Culicidae) movement patterns in Thailand. J Vector Ecol. 2010;35:428–31.
 20.
Achee NL, Grieco JP, Andre RG, Rejmankova E, Roberts DR. A markreleaserecapture study using a novel portable hut design to define the flight behavior of Anopheles darlingi in Belize, Central America. J Am Mosq Control Assoc. 2005;21:366–79.
 21.
Randriamaherijaona S, Briët OJ, Boyer S, Bouraima A, N’Guessan R, Rogier C, Corbel V. Do holes in longlasting insecticidal nets compromise their efficacy against pyrethroid resistant Anopheles gambiae and Culex quinquefasciatus? Results from a releaserecapture study in experimental huts. Malar. J. 2015;14:332.
 22.
Buhagiar TS, Devine GJ, Ritchie SA. Effects of sublethal exposure to metofluthrin on the fitness of Aedes aegypti in a domestic setting in Cairns, Queensland. Parasit Vectors. 2017;10:274.
 23.
Chapman DG, Junge CO Jr. The estimation of the size of a stratified animal population. Ann Math Stat. 1956;27:375–89.
 24.
Jolly GM. Explicit estimates from capturerecapture data with both death and immigrationstochastic model. Biometrika. 1965;52:225–47.
 25.
Arnason AN. Parameter estimates from markrecapture experiments on two populations subject to migration and death. Res Popul Ecol. 1972;13:97–113.
 26.
Arnason AN. The estimation of population size, migration rates and survival in a stratified population. Res Popul Ecol. 1973;15:1–8.
 27.
Seber GA. A note on the multiplerecapture census. Biometrika. 1965;52:249–59.
 28.
Lebreton J, Nichols JD, Barker RJ, Pradel R, Spendelow JA. Modeling individual animal histories with multistate capturerecapture models. Adv Ecol Res. 2009;41:87–173.
 29.
Lebreton J, Cefe RP. Multistate recapture models: modelling incomplete individual histories. J Appl Stat. 2002;29:353–69.
 30.
Kendall WL, Nichols JD. Estimating statetransition probabilities for unobservable states using capturerecapture/resighting data. Ecology. 2002;83:3276–84.
 31.
Ergon T, Gardner B. Separating mortality and emigration: modelling space use, dispersal and survival with robustdesign spatial capturerecapture data. Methods Ecol Evol. 2014;5:1327–36.
 32.
Villela DA, Codeço CT, Figueiredo F, Garcia GA, MacieldeFreitas R, Struchiner CJ. A Bayesian hierarchical model for estimation of abundance and spatial density of Aedes aegypti. PLoS One. 2015;10:e0123794.
 33.
Ovaskainen O. Habitatspecific movement parameters estimated using markrecapture data and a diffusion model. Ecology. 2004;85:242–57.
 34.
Manda H, Shah P, Polsomboon S, Chareonviriyaphap T, CastroLlanos F, Morrison A, et al. Contact irritant responses of Aedes aegypti using sublethal concentration and focal application of pyrethroid chemicals. PLoS Negl Trop Dis. 2013;7:e2074.
 35.
McLeanCooper N, Achee N, Foggie T, Grieco J, Williams J. Space optimizing methods for laboratory rearing of Aedes aegypti. J Am Mosq Control Assoc. 2008;24:460–2.
 36.
Grieco JP, Achee NL, Andre RG, Roberts DR. A comparison study of house entering and exiting behavior of Anopheles vestitipennis (Diptera: Culicidae) using experimental huts sprayed with DDT or deltamethrin in the southern district of Toledo, Belize, C.A. J Vector Ecol. 2000;25:62–73.
 37.
VazquezProkopec GM, Galvin WA, Kelly R, Kitron U. A new, costeffective, batterypowered aspirator for adult mosquito collections. J Med Entomol. 2009;46:1256–9.
 38.
Taylor HM, Karlin S. An introduction to stochastic modeling. San Diego: Academic Press; 2014.
 39.
Gilks WR. Markov Chain Monte Carlo. In: Encyclopedia of Biostatistics. Wiley Online Library; 2005. https://doi.org/10.1002/0470011815.b2a14021.
 40.
Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–72.
 41.
Sobol IM. On the distribution of points in a cube and the approximate evaluation of integrals. Zh Vychisl Mat Mat Fiz. 1967;7:784–802.
 42.
King A, Ionides E, Bretó C, Ellner S, Kendall B, Wearing H, et al. pomp: Statistical inference for partially observed Markov processes (R packageversion~1.15); 2010. https://kingaa.github.io/pomp/.
 43.
McPhatter LP, Mischler PD, Webb MZ, Chauhan K, Lindroth EJ, Richardson AG, Debboun M. Laboratory and semifield evaluations of two (transfluthrin) spatial repellent devices against Aedes aegypti (L.) (Diptera: Culicidae). US Army Med Dep J. 2017;1–17:13–22.
 44.
Nentwig G, Frohberger S, Sonneck R. Evaluation of clove oil, icaridin, and transfluthrin for spatial repellent effects in three tests systems against the Aedes aegypti (Diptera: Culicidae). J Med Entomol. 2016;54:150–8.
 45.
Bibbs CS, Kaufman PE. Volatile pyrethroids as a potential mosquito abatement tool: a review of pyrethroidcontaining spatial repellents. J Integr Pest Manag. 2017;8:21.
 46.
Darbro JM, Muzari MO, Giblin A, Adamczyk RM, Ritchie SA, Devine GJ. Reducing biting rates of Aedes aegypti with metofluthrin: investigations in time and space. Parasit Vectors. 2017;10:69.
 47.
Hoffmann EJ, Miller JR. Reduction of mosquito (Diptera: Culicidae) attacks on a human subject by combination of wind and vaporphase DEET repellent. J Med Entomol. 2002;39:935–8.
 48.
Restif O, Hayman DT, Pulliam JR, Plowright RK, George DB, Luis AD, et al. Modelguided fieldwork: practical guidelines for multidisciplinary research on wildlife ecological and epidemiological dynamics. Ecol Lett. 2012;15:1083–94.
 49.
Ogoma SB, Mmando AS, Swai JK, Horstmann S, Malone D, Killeen GF. A low technology emanator treated with the volatile pyrethroid transfluthrin confers long term protection against outdoor biting vectors of lymphatic filariasis, arboviruses and malaria. PLoS Negl Trop Dis. 2017;11:e0005455.
 50.
Govella NJ, Ogoma SB, Paliga J, Chaki PP, Killeen G. Impregnating hessian strips with the volatile pyrethroid transfluthrin prevents outdoor exposure to vectors of malaria and lymphatic filariasis in urban Dar es Salaam, Tanzania. Parasit Vectors. 2015;8:322.
 51.
Ogoma SB, Ngonyani H, Simfukwe ET, Mseka A, Moore J, Killeen GF. Spatial repellency of transfluthrintreated hessian strips against laboratoryreared Anopheles arabiensis mosquitoes in a semifield tunnel cage. Parasit Vectors. 2012;5:54.
 52.
Sugiharto VA, Grieco JP, Murphy JR, Olsen CH, ColaciccoMayhugh MG, Stewart V, et al. Effects of preexposure to DEET on the downstream bloodfeeding behaviors of Aedes aegypti (Diptera: Culicidae) mosquitoes. J Med Entomol. 2016;53:1100–4.
 53.
Maia MF, Onyango SP, Thele M, Simfukwe ET, Turner EL, Moore SJ. Do topical repellents divert mosquitoes within a community? Health equity implications of topical repellents as a mosquito bite prevention tool. PLoS One. 2013;8:e84875.
 54.
Brady OJ, Godfray HC, Tatem AJ, Gething PW, Cohen JM, McKenzie FE, et al. Vectorial capacity and vector control: reconsidering sensitivity to parameters for malaria elimination. Trans R Soc Trop Med Hyg. 2016;110:107–17.
 55.
Briët OJ, Penny MA, Hardy D, Awolola TS, Van Bortel W, Corbel V, et al. Effects of pyrethroid resistance on the cost effectiveness of a mass distribution of longlasting insecticidal nets: a modelling study. Malar J. 2013;12:77.
 56.
Martin NJ, Smith PA, Achee NL, DeLong GT. Determining airborne concentrations of spatial repellent chemicals in mosquito behavior assay systems. PLoS One. 2013;8:e71884.
 57.
Martin NJ, Smith PA, Brown CW, Achee NL, DeLong GT. Dichlorodiphenyltrichloroethane determination in air by thermal desorption gas chromatographymass spectrometry. Pest Manag Sci. 2012;68:1360–7.
 58.
Achee N, Masuoka P, Smith P, Martin N, Chareonviryiphap T, Polsomboon S, et al. Identifying the effective concentration for spatial repellency of the dengue vector Aedes aegypti. Parasit Vectors. 2012;5:300.
 59.
Choi DB, Grieco JP, Apperson CS, Schal C, Ponnusamy L, Wesson DM, Achee NL. Effect of spatial repellent exposure on dengue vector attraction to oviposition sites. PLoS Negl Trop Dis. 2016;10:e0004850.
 60.
Bowman LR, Donegan S, McCall PJ. Is dengue vector control deficient in effectiveness or evidence?: systematic review and metaanalysis. PLoS Negl Trop Dis. 2016;10:e0004551.
Acknowledgements
We express our gratitude to Carlos Ique, Director of the Institute Veterinario de Investigaciones (IVITA), Iquitos, Peru for the use of land for experimental hut studies. Thanks to Roxanne Burrus, Kirk Mundal and Victor Lopez (NAMRU6) for their logistic support during the experiments and to Neil Lobo for valuable insights. Special thanks to Edwin Requena, Hugo Jaba and to the team of mosquito collectors for their dedicated effort and to reviewers for valuable comments.
Funding
This work was supported and funded by the Bill & Melinda Gates Foundation (Grant #48513): "A pushpull strategy for Aedes aegypti control." The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. QAtB was supported by a graduate student fellowship from the Eck Institute for Global Health at the University of Notre Dame.
Availability of data and materials
The data supporting the conclusions of this article are included within the article and its additional file. The datasets generated during the current study are available in the Open Science Framework repository, https://osf.io/xtmy7/ as part of the OSF project https://osf.io/5hcpf/. All code is available on Github.com https://github.com/quirine/ExperimentalHuts.
Disclaimer
The views expressed in this work are those of the authors and do not reflect the official policy or position of the Department of the Navy, Department of Defense or U.S. Government.
Copyright statement
FCL and ACM are employees of the U.S. Government. This work was prepared as part of their official duties (PJTNMRCD 0.24). Title 17 U.S.C. § 105 provides that ‘Copyright protection under this title is not available for any work of the United States Government’. Title 17 U.S.C. § 101 defines a U.S. Government work as a work prepared by a military service member or employee of the U.S. Government as part of that person’s official duties.
Author information
Affiliations
Contributions
FCL, HM, ACM, JPG and NLA designed the hut experiments. FCL and HM performed the experiments. QAtB and TAP developed the modeling and simulation framework. QAtB performed the analyses. QAtB, JPG, NLA and TAP interpreted the results. QAtB, NLA and TAP wrote the manuscript. All authors read and approved the final manuscript.
Corresponding authors
Correspondence to Quirine A. ten Bosch or T. Alex Perkins.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1:
Table S1. Average GelmanRubin statistics across simulated data sets (median and the upper bound of the 95% confidence interval). Figure S1. Correlations between parameter posteriors of model fit on baseline scenario with all parameters estimated. Marginal posteriors are depicted on the diagonals. The numbers on the right of the diagonal depict the Spearman rank correlation coefficients for each side by side comparison. Figure S2. Correlations between parameter posteriors of model fit on baseline scenario with r_{ i } fixed (r_{T}: 0.040, r_{1}:0.036, r_{2}: 0.088). Marginal posteriors are depicted on the diagonals. The numbers on the right of the diagonal depict the Spearman rank correlation coefficients for each side by side comparison. Under this parameterization, the movement rate q_{ i } is exactly equal to the product x_{ i }/ r_{ i }. Figure S3. Correlations between parameter posteriors of model fit on low dosage scenario with r_{ i } fixed (r_{T}: 0.040, r_{1}:0.036, r_{2}: 0.088). Marginal posteriors are depicted on the diagonals. The numbers on the right of the diagonal depict the Spearman rank correlation coefficients for each side by side comparison. Under this parameterization, the movement rate q_{ i } is exactly equal to the product x_{ i }/ r_{ i }. Figure S4. Correlations between parameter posteriors of model fit on high dosage scenario with r_{ i } fixed (r_{T}: 0.040, r_{1}:0.036, r_{2}: 0.088). Marginal posteriors are depicted on the diagonals. The numbers on the right of the diagonal depict the Spearman rank correlation coefficients for each side by side comparison. Under this parameterization, the movement rate q_{ i } is exactly equal to the product x_{ i }/ r_{ i }. Figure S5. Posterior distributions of model parameters fitted to experimental data while fixing the values of r_{ i } at the 2.5th percentile of the posterior from the full parameter fit to the baseline data (r_{T}: 0.031, r_{1}:0.027, r_{2}: 0.065). Posteriors are shown for the baseline (gray), low dosage (orange) and high dosage (pink) for the SRhut (subscript 0) and huts 2 or 1 removed (subscript 2 and 1, respectively). ac rates at which mosquitoes exit the huts, d proportion of movement from H_{1} (hut directly adjacent to the treatment hut) away from the SRproduct. eg knockdown rates and h loss to followup rates. Under this parameterization, the movement rate q_{ i } is exactly equal to the product x_{ i }/ r_{ i }.The algorithm was run for 25,000 iterations with a ‘burnin’ period of 10,000. Figure S6. Posterior distributions of model parameters fitted to experimental data while fixing the values of r_{ i } at the 97.5th percentile of the posterior from the full parameter fit to the baseline data (r_{T}: 0.054, r_{1}:0.040, r_{2}: 0.11). Posteriors are shown for the baseline (gray), low dosage (orange) and high dosage (pink) for the SRhut (subscript 0) and huts 2 or 1 removed (subscript 2 and 1, respectively). ac rates at which mosquitoes exit the huts, d proportion of movement from H_{1} (hut directly adjacent to the treatment hut) away from the SRproduct. eg knockdown rates and h loss to followup rates. Under this parameterization, the movement rate q_{ i } is exactly equal to the product x_{ i }/ r_{ i }. The algorithm was run for 25,000 iterations with a ‘burnin’ period of 10,000. Figure S7. GelmanRubin convergence diagnostics by iteration for the baseline scenario. Figure S8. GelmanRubin convergence diagnostics by iteration for the low dosage scenario. Figure S9. GelmanRubin convergence diagnostics by iteration for the high dosage scenario. Figure S10. Trace plots for the baseline scenario. Figure S11. Trace plots for the low dosage scenario. Figure S12. Trace plots for the high dosage scenario. (DOCX 50818 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Aedes aegypti
 Bayesian parameter estimation
 Continuoustime Markovchain models
 Dengue
 Spatial repellent
 Transfluthrin
 Vector control
 Zika
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.