Fitting the Incidence Data from the City of Campinas, Brazil, Based on Dengue Transmission Modellings Considering Time-Dependent Entomological Parameters

Four time-dependent dengue transmission models are considered in order to fit the incidence data from the City of Campinas, Brazil, recorded from October 1st 1995 to September 30th 2012. The entomological parameters are allowed to depend on temperature and precipitation, while the carrying capacity and the hatching of eggs depend only on precipitation. The whole period of incidence of dengue is split into four periods, due to the fact that the model is formulated considering the circulation of only one serotype. Dengue transmission parameters from human to mosquito and mosquito to human are fitted for each one of the periods. The time varying partial and overall effective reproduction numbers are obtained to explain the incidence of dengue provided by the models.


Introduction
Dengue virus, a flavivirus transmitted by arthropod of the genus Aedes, is prevalent in different parts of the world. As a result of being pathogenic for humans and capable of transmission in heavily populated areas, dengue virus (arbovirus) can cause widespread and serious epidemics, which constitute one of the major public health problems in many tropical and subtropical regions of the world where Aedes aegypti and other appropriate mosquito vectors are present [1]. Additionally, the seasonality in these regions influences on the size of mosquito population, and also on the extrinsic incubation period due to the fact that mosquitoes belong to coldblooded (poikilothermic) animals.
With respect to the seasonality of abiotic conditions, entomological parameters such as transition and mortality rates in aquatic phase, mortality rate of adult mosquito and oviposition rate vary with temperature [2][3] [4]. On the other hand, there are scarce knowledge about [13] for a model including eggs compartment), is determined by the oviposition rate φ(M) = ϕ m M, where ϕ m is the per-capita oviposition rate and M, the number of female mosquitoes at time t. Defining L as the number of larvae (female) at time t, the effective larvae production rate is given by qf[1−L/(DC)]ϕ m M, where we assume that a fraction q of eggs will hatch to the larva stage and a fraction f of these eggs become female mosquitoes, and DC is the total (carrying) capacity of the breeding sites. The constant parameter D represents the magnitude of the breeding sites, while C carries on the time varying abiotic conditions and demography of humans. The number of larvae decreases according to the change of larvae to pupae and death, described, respectively, by the changing σ l and the mortality μ l rates. The number of pupae at time t, designated by P, increases with the change of larvae to pupae (σ l ) and decreases according to the transformation of pupae to adult mosquitoes and death, described, respectively, by the emerging σ p and the mortality μ p rates. Finally, the number of female mosquitoes increases according to the emerging of pupae (σ p ) and decreases according to the mortality rate μ f . Table 1 presents summary of parameters related to mosquito population and representative values. The value of D was chosen to be of the order of the magnitude of human population in the City of Campinas.
The entomological parameters ϕ m , σ l , σ p , μ l , μ p and μ f clearly depend on temperature [3]. But, besides the temperature, the rainfall influences on the entomological parameters μ l and μ p (μ f is practically insensitive [14]). The rainfall affects also on the carrying capacity C and the hatching fraction q. Let us describe the dependency of the model parameters on temperature and rainfall.
From a series of daily recorded maximum (T max ) and minimum (T min ) temperatures [12], we interpolated any temperature T at time t as TðtÞ ¼ T minðjÞ þ ½T maxðjÞ À T minðjÞ sin ðpt 0 Þ; ð1Þ where the index j, which is the integer part of time t, refers to the j-th calendar day, and t 0 is the fractional part of time t, that is, 0 t 0 < 1. Notably, t 0 = 0 and t 0 = 1 correspond to midnight, whereas t 0 = 0.5 corresponds to midday.
To incorporate the rainfall in the model, we define W(j) as the amount of rainfall in j-th calendar day. The amount of rain precipitation influences on the mortality rates μ l and μ p in the aquatic phase and, also, on the carrying capacity (C) and the capacity of eggs being hatched (q). During heavy rain periods, larvae and pupae can be flushed, thereby cleansing breeding sites of aquatic forms [15]. We refer to this situation as physically induced mortality, and we do not consider this type of mortality among adult mosquitoes [16].
From a series of daily recorded precipitation [12], we define the simplest physically induced mortality rates m c l and m c p as ( m c l ¼ m l f1 þ g l ½WðjÞ À V c y½WðjÞ À V c g m c p ¼ m p f1 þ g p ½WðjÞ À V c y½WðjÞ À V c g; where μ l and μ p are temperature-dependent mortality rates, and g l and g p measure how the precipitation affects the additional physical mortality rates. Additionally, V c is the critical rain volume sufficiently high to originate the overflow of breeding sites; θ(x) is the Heaviside function, that is, θ(x) = 1, if x ! 0 and is otherwise θ(x) = 0. The fraction of eggs that are hatching (q) and the carrying capacity (C) are assumed to depend only on precipitation. We assume that the rain that has fallen in the past few days also influences both parameters.
For the fraction of eggs that hatch into larvae q, we follow similar argumentation given for the carrying capacity (see below). Because the relative humidity is almost 100% when there is abundant rain, we assume that where q 0 is the capacity of eggs hatching with rainfall, q 1 is the the critical amount of rain to promote hatching (the lower q 0 , the higher the influence of humidity on hatching). Notably, q = 1 when W ! 1. This type of Hill function was chosen to describe the saturating behavior [17]. The term W M represents the rainfall of past days defined by Wðj À oÞ fw 1 ½T max ðj À oÞ þ T min ðj À oÞg o ; where W(j − o) is the rainfall in j-o-th days before, w 1 is the residual effect of past rainfall. Observe that in the proposed expression for W M the effect of past rain depends on the temperature since evaporation increases at higher temperatures; the parameter w 1 regulates the past contributions in the sense that the higher w 1 , the lower the contribution of past rain due to evaporation. If w 1 = 1/2, then the evaporation is a function of the daily mean temperature. Heavy and continuous rainfall may be a regulating factor for A. aegypti proliferation, as immature forms might have been washed away and thus significantly reduced the adult production. However, heavy rainfall prepares the habitat for A. aegypti proliferation in the post rainy season [7]. During the dry seasons, the number of breeding sites is dramatically reduced. We assume that the amount of rainfall always increases the breeding sites, even in the case of heavy rainfall (cleansing of breeding sites). Hence, for carrying capacity, we describe the influence of rainfall on its number as where C 2 is the rain independent variation in the breeding sites but can depend on time due to demographic changes in human population (not considered here), C 0 is the capacity of rain in producing new breeding sites due to rain, and C 1 is the the critical amount of rain in the formation of breeding sites (the lower C 1 , the higher the influence of rain on increasing breeding sites). Saturating behavior is again observed: C = C 2 + C 0 when W ! 1.
We now want to consider control efforts, such as insecticide application (chemical control) and the remotion of breeding sites (mechanical control). But, breeding sites can be removed in two ways: (1) remotion of breeding sites targeting houses surrounding case notification, and (2) preventive remotion to avoid dengue transmission. The total size of inhabitations in the City of Campinas [18] is approximately n T = 350,000. We assume that the number of houses surveyed for breeding sites, and for further remotion, is proportional to decreasing in the carrying capacity, represented by À a a n a n T ; where n b and n a are the numbers of houses surveyed surrounding dengue cases, with α b and α a corresponding to the efficacy of mechanical control. Because the preventive visit is a random event (subscript a), whereas the dengue-targeting house visitation is case-search event (subscript b), we assumed that α b > α a . Spraying insecticide aims to kill adult mosquitoes, which is represented by where n i is the number of houses where insecticide was sprayed surrounding dengue cases, and μ f is the temperature-dependent mortality rate. Susceptible mosquitoes are distributed in all region, while the infectious mosquitoes are located in the neighborhood of infectious individuals. Hence, we assume that susceptible and infectious mosquitoes are under different additional mortalities. Defining α s as the insecticide induced mortality among susceptible mosquitoes, and α d , the induced mortality among infectious mosquitoes, then α d > α s . Summarizing, the expressions for the carrying capacity, the hatching fraction and the entomological parameters are Two extreme weather conditions are not considered in this modelling: monsoon season and long period of absence of rainfall. During monsoon season significant reduction of adult and washing away of immature stages are observed, which significantly reduces the breeding sites.
In tropical and subtropical regions like the City of Campinas, heavy rains occur during few days in comparison with monsoon season. In our assumption, breeding sites are always increased with rain (which occurs during few hours or days), but immature stages are under physically induced mortality. During severe and prolonged drought period, people tend to storage water at home, sometimes in not well covered containers, and small bodies of water tend to dry and form multiple breeding sites. We are not considering this situation neither, which indeed occurred in the City of Campinas in the summer of 2015. Table 2 presents a summary of parameters and representative values related to precipitation, and Table 3 presents a summary of parameters and representative values related to vector control. With respect to the vector control parameters in Table 3, the number of houses visited by public health agents was small. For this reason, they visit neighborhood of reported cases of dengue. Hence the effects on susceptible mosquitoes (distributed in a very large area) is assumed to be very small in comparison with that on infectious mosquitoes, which are located in the neighborhood of dengue cases. Hence, we attributed to the efficacy parameters plausible values based on the previous reasonings.

Hypotheses on the dengue transmission
With respect to dengue transmission, the human population is divided into four compartments according to the natural history of the disease: S, E, I and R, which are the numbers at time t of, respectively, susceptible, exposed, infectious and recovered humans, with S + E + I + R = N, where N is the size of the human population. The fractions of human subpopulations are defined by the size of compartment divided by N, or s = S/N, e = E/N, i = I/N and r = R/N. The female mosquito population is divided into three compartments: M 1 , M 2 and M 3 , which are the numbers at time t of, respectively, susceptible, exposed and infectious mosquitoes. The size of mosquito population is given by M = M 1 + M 2 + M 3 . By dividing M 1 , M 2 and M 3 by the magnitude D, for instance m 1 = M 1 /D, we consider densities in order to have the same order of magnitude of fractions of sub-populations of humans. However, the fractions are defined by dividing by M; for instance, the fraction of susceptible mosquitoes is M 1 /M = m 1 /m. Hence, we stress the fact that m 1 is the density of susceptible mosquitoes, while m 1 /m is the fraction, with m = m 1 + m 2 + m 3 . Dengue transmission is sustained by the flows among human and mosquito compartments according to the dengue epidemics cycle. Susceptible humans are infected during the blood meal by infectious mosquitoes, with the force of infection (or per-capita incidence rate) being designated by B h . The exposed persons are, then, transferred to an infectious class by rate γ h , where 1/γ h is the intrinsic incubation period. These infectious persons progress to recovered (immune) class at rate σ h . Neither loss of immunity nor induced mortality due to the disease are considered (a unique serotype infection). With respect to the vector, the susceptible mosquitoes are infected at a force of infection B m . These exposed mosquitoes are transferred to infectious class at a rate γ m , where 1/γ m is the extrinsic incubation period, and remain infective until death.
The intrinsic incubation rate γ h does not depend on ambient temperature, while the extrinsic incubation rate γ m does. Following Lindsay and Birley [19], the extrinsic incubation rate depends on temperature according to where T is the ambient temperature (in degrees Celsius°C), T s is the thermic sum (in°C ×day), and T m is the threshold of temperature below which dengue virus cannot multiply, hence T > T m . The forces of infection (or per-capita incidence rates) B h and B m depend on the frequency of bites on humans by mosquitoes (assumed to be dependent on the per-capita oviposition rate ϕ m ) and on the number (or fraction) of infective populations (designated by U h for humans and U m for mosquitoes). Hence, the forces of infection B h and B m can be written as B h = β h ϕ m U m and B m = β m ϕ m U h , where β h and β m are the transmission coefficients (dimensionless with respect to time). Then, the transmission of dengue infections among susceptible individuals are given by where t is in days. The demographic parameters ϕ h and μ h are, respectively, the natality and mortality rates. Different dependencies of the force of infection with infective populations U h and U m are studied. First, the frequency dependent (FD) model is considered. The rate at which one mosquito bites one human can be assumed being directly proportional to the average biting rate of mosquitoes (intrinsic behavior of mosquitoes), and inversely proportional to the number of humans (one specific person has the chance of being bitten by one mosquito is reduced in a large population). In this situation, U m = M 3 /N and U h = I/N, and the forces of infection B h and B m are given by ðFDÞ : ( ( and the coefficients β h and β m are dimensionless with respect to populations, and they are total transmission coefficients. Finally, a special incidence rate (SIR), or special mass action law, model is considered. In this modelling, the propagation of dengue is proportional to both fractions of susceptible and infectious populations. In this case, U m = M 3 /(MN),U h = I/(MN), and the per-capita incidence rates B h and B m are ðSIRÞ :  Table 4 presents a summary of parameters and representative values related to dengue transmission and human population. The value of γ m was obtained using Eq (3).

Dynamics of dengue transmission
We assume that the size of human population at time t, designated by N, obeys Malthus law. Hence, its size varies according to the difference between natality and mortality rates designated by ϕ h and μ h , respectively, or Taking into account this equation, the equations for the size of sub-populations are related to the corresponding fractions. For instance, for the fraction of susceptibles s, we have and using Eq (4), it is linked to the number of susceptibles S through where equation for S was previously defined. The same procedure must be done for equations related to fractions e and i. But, for mosquito population, the numbers L, P, M 1 , M 2 and M 3 are divided by the magnitude D, for instance, the density l = L/D. Fractions and densities are dimensionless variables with respect to populations. Based on the foregoing descriptions of the model parameters and variables, dengue transmission is described by two systems of differential equations. First system considers sub-populations where forces of infection do not appear, that is, where the decoupled fraction of immune humans is given by r = 1 − s − e − i. The second system describes the susceptible and exposed populations according to the hypothesis related to the disease transmission: (1) for frequency dependent model (FD), > > > > > > < > > > > > > : (2) for pseudo mass action law model (PMAL), (3) for true mass action law model (TMAL), and (4) for special incidence rate model (SIR), In all models, β h and β m are assumed to be constant. The entomological parameters ϕ m , σ l , σ p , μ l , μ p and μ f depend on temperature T, which is a function of time t according to Eq (1). The parameters m c l , m c p , m s f and m d f are the entomological parameters μ l , μ p and μ f , respectively, encompassing rainfall and controlling efforts, whose expressions is given in Eq (2). The parameters C and q are dependent on precipitation. Therefore, all these parameters are function of time, and the dynamical systems are non-autonomous.
Summarizing, the first model (Model FD) considers the force of infection determined by biting probability, leading to a system formed by Eqs (5) and (6); the second model (Model PMAL) considers the pseudo mass action law, which gives the system formed by Eqs (5) and (7); the third model (Model TMAL) considers true mass action law, given the Eqs (5) and (8); the fourth model (Model SIR) considers the special incidence rate, which leads to Eqs (5) and (9). In all models, the dynamics of the human population is described by Eq (4).

Effective reproduction numbers
In non-autonomous modelling, there are not well established epidemiological parameters to evaluate the risk of epidemics, such as the basic reproduction number. However, to a nonautonomous system, we can associate an autonomous by assuming constant parameters given by the mean values of the respective time-dependent parameters. The behaviors of the solutions of this associated autonomous system are then approximations of the dynamical behaviors of the solutions of the original system, at least for small time intervals; their long time behavior are also interesting to know since they inform what would be the dynamical tendencies if the parameters do not change in time. Therefore, by analyzing the special case where all model parameters are considered constant in time, we can borrow from autonomous modeling concepts, such as the basic reproduction number, and have at least some related information for the non-autonomous case.
In the particular case of constant model parameters and constant human and mosquito populations, the steady state of the autonomous model can be obtained. In S1 Text, from the steady states analysis of the autonomous dynamical system, the basic reproduction number R 0 , equation (S.10), is defined as the product of partial reproduction numbers R m 0 and R h 0 , which, depending on the model, are given by equations (S.6), (S.7), (S.8) and (S.9). The values shown in Tables 1 to 4 can be taken as the representative values for the model parameters, which could be the mean values over one year.
From R 0 and its partial reproduction numbers R m 0 and R h 0 , let us defined the time varying effective reproduction number R ef in terms of the partial effective reproduction numbers in humans (R h ef ) and mosquitoes (R m ef ) as (In the definitions for R h ef and R m ef below, consider that the model parameters and size of populations N and m are constant in time, but susceptible populations s and m 1 are varying in time.) Suppose that at t = 0 one infectious case (does not matter if human or mosquito) is introduced in a completely susceptible populations of humans and mosquitoes (then, controls are absent). For t 0, before and just at the time of the beginning of epidemics, R ef = R 0 because s = 1 and m 1 = m (then, Dm Ã 1 ¼ M Ã , corresponding to the trivial equilibrium). Notwithstanding, when t ! 1, the epidemics reaches a steady state, which occurs due to R ef = 1, and equation (S.5) can be obtained, assuming that controls are applied.
In an autonomous modelling, the dynamical trajectory can be determined when one infectious human is introduced in a community free of dengue. The equilibrium point before the introduction of infectious case is given by P 0 , with equilibrium values of mosquito population being given by equation (S.1). Depending on the value of the basic reproduction number, the dengue disease fades out after a small epidemics (R 0 < 1), or attains an endemic level (R 0 > 1). Now, in non-autonomous models developed here, the majority of model parameters (exceptions are f, β m , β h , and human related parameters ϕ h , μ h , σ h and γ h ) and the size of human and mosquito populations vary with time. Following the idea of the partial reproduction numbers obtained from the autonomous modelling, let us define the time varying effective reproduction number R ef in terms of the partial effective reproduction numbers in humans (R h ef ) and mosquitoes (R m ef ) for our non-autonomous models. For model FD, the partial effective reproduction numbers at any time t are defined by, using M = Dm, ðFDÞ : For model PMAL, the partial effective reproduction numbers at time t are ðPMALÞ : For model TMAL, the partial effective reproduction numbers at time t are ðTMALÞ : For model SIR, the partial effective reproduction numbers at time t are ðSIRÞ : The partial effective reproduction numbers are evaluated at each time t.
The non-autonomous systems of equations depend on the temperature and rainfall, which assume different values in each time t. For this reason, the time dependent model parameters are presented below, which will be used in the numerical simulations of the four models named FD, TMAL, PMAL and SIR.

Data collected from the City of Campinas
The daily maximum and minimum temperatures and precipitation [12] from October 1 st , 1995 to September 30 th , 2012 in the City of Campinas, São Paulo State, Brazil, are used to numerically simulate the previously described dynamics systems. From Fig 1(a), the maximum temperatures oscillate between 16 and 37°C, while the minimum, between 4 and 24°C. A mean upper amplitude of the precipitation situates around 50 mm, and in each year some peaks of rainfall were observed, with the highest peak situating near 145 mm (see Fig 1(b)).
From a series of daily recorded maximum (T max ) and minimum (T min ) temperatures [12], shown in Fig 1(a), we interpolate any temperature T at time t according to Eq (1). . The worst epidemics occurred in the third period. The possible circulating serotypes can be determined in these four periods based on serological screen in a small fraction of total cases. In the first period of time, serotype 1 predominated (96.4%), and serotype 2 was 3.6%. In the second period of time, serotype 1 predominated (72.2%), and serotype 3 was 27.8%. In the third period of time, only serotype 3 circulated (100%). In the last period, serotype 1 predominated (76.7%), followed by serotype 2 (16.7%) and by serotypes 3 and 4 (3.3% each). No serotype were found in years 2000, 2004 and 2008, which was another reason to grouping in the previous four periods. Fig 2(b) shows the number of houses visited to spray insecticides, and Fig 2(c), remotion of breeding sites. In the second period, the decreased cases of dengue could be explained by  intense controlling efforts, while the third period presented highest epidemics due to the controlling efforts have been applied after the beginning of epidemics. Another factor to explain this elevated cases of dengue is the massive infection by serotype 3, which was already circulating in the second period.
Population in the City of Campinas varied from 1991 (847,595) to 2010 (1,080,999) [10]. This varying population was taken into account estimating the values of ϕ h and μ h using Eq (4). The estimated values were given in Table 4. The total size of inhabitations in City of Campinas is approximately n T = 350,000 [18].

Time dependent entomological parameters
Based on the temperature and rainfall shown in Fig 1, time dependent parameters related to mosquito population are presented. The temperature dependent entomological parameters ϕ m , μ l , μ p , μ f , σ l and σ p of A. aegypti are obtained using fittings given in Table 5 [2][3] [4]. An n-th degree polynomial P n (T) = b 0 + b 1 T + Á Á Á + b n T n was used to fit experimental data. The carrying capacity C and the fraction of eggs hatching q shown in Fig 3 were obtained using Eq (2) and the volume of rainfall shown in Fig 1(b). The values of the parameters dependent on the precipitation were given in Table 2. The transition rates σ l and σ p were obtained using polynomials fitted given in Table 5, and the temperature at each time t was calculated using Eq (1), with maximum and minimum daily temperatures being shown in Fig 1(a). The temperature dependent mortality rates of aquatic forms μ l and μ p were obtained using polynomials fitted given in Table 5, and the temperature at each time t was calculated using Eq (1), with maximum and minimum daily temperatures being shown in Fig 1(a). The Table 5. Estimated entomological parameters of female mosquito and aquatic phase.
These parameters were estimated using an n-th degree polynomial, with unit of coefficients b i being days −1 × (°C) −i . Standard deviations are not shown.
doi:10.1371/journal.pone.0152186.t005   Mathematical Modelling and Fitting of Dengue Incidence contribution of the rainfall was calculated by Eq (2) and the volume of rainfall is shown in Fig  1(b). The values of the parameters dependent on the precipitation were given in Table 2. Fig 6 shows the mortality rate of adult mosquitoes μ f (a) and the oviposition rate ϕ m (b) during the period of time (October) 1995-2012 (September). These parameters are allowed to depend only on temperature.
The temperature dependent mortality rate of adult mosquitoes μ f and the oviposition rate ϕ m were obtained using polynomials fitted given in Table 5, and the temperature at each time t was calculated using Eq (1), with maximum and minimum daily temperatures being shown in Fig 1(a). Figs 3 to 6 show clearly the model parameters depending on temperature and also on precipitation. In unfavorable seasons, there is a dramatic decrease in these parameters.

Estimation method
In the foregoing subsection, the model parameters related to mosquito (time dependent) and human (time independent) populations were obtained. Time-dependent model parameters were shown graphically. The simulations done in Yang et al. [14] only gave us certain intervals where probably the values of β m and β h are located. In the next section, the transmission coefficients β m and β h are fitted for the dengue incidence recorded in the City of Campinas. To obtain the estimated values of β m and β h , it is necessary to appeal to a suitable numerical algorithm for parameter identification.
In order to estimate β m and β h , it is natural to start by considering a weighted discrete L 2error; thus, a possible functional Er to be minimized could be: 4. i obs (t j ) is the fraction of the expected infected individuals, which is computed as where cases(t j ), for t 1 , Á Á Á, t N , are the recorded cases of dengue in the City of Campinas, and N(t j ) is the size of the human population. The additional parameter prop encompasses all cases of dengue which were not detected (asymptomatic individuals) [21]; the value of this parameter must also be found during the fitting process.
5. p(λ, t j ) is the weighting factor that depends on a pre-established parameter λ ! 1; such weighting factor could be, for instance: In this case, we observe that by choosing large value of λ, one can favor the days with more infected people in the process of fitting the data; λ = 1 attributes equal weights to the given data; The transmission coefficients β m and β h , and the term prop are then estimated using confirmed dengue cases. In the procedure of the estimation method, we need to provide the intervals where these parameters are located as Interval where b m is probably located : ½a; b Interval where b h is probably located : ½c; d Interval where prop is probably located : ½A; B: In S1 Text, the estimation method is described in detail.
With respect to observed data, the recorded cases of dengue in the City of Campinas are cumulative number of dengue cases I obs (t j ) in each month. Hence, i obs (t j ) = I obs (t j )/30, where I obs (t j ) is given in Fig 2(a), and t j is the mean calendar time of each month.
The Due to the characteristics of our data, in the actual fitting procedure, we used a weighting factor distinct of the one described by Eq (16). In fact, it is observed a pattern in the incidence cases in each year (see Fig 2(a)): in October initiates the increasing phase, which is followed by decreasing up to September. Hence, let us define the annual epidemiological cycle comprised between October 1 st and September 30 th of the next year. In each annual epidemiological cycle, the peak of incidence can vary broadly. Then, instead of a maximum incidence in the whole period of estimation, maximum in each annual cycle max j |i obs (t j )| g , for g = 1, Á Á Á, u, where u is the number of annual cycles in this period, is determined. In this case, a possible weighting factor could be p g ðl; t j Þ ¼ 1 þ ðl À 1Þ i obs ðt j Þ max j ji obs ðt j Þj g ; for g ¼ 1; Á Á Á ; u; and j = 1, Á Á Á, j g , where j g is the number of recorded dengue cases in the annual cycle g, with P u g¼1 j g ¼ k l , where k l is the number of recorded dengue cases in the l-th period. Additional modifications to the weighting factor p g (λ, t j ) are done as follows. In an annual epidemiological cycle, all observed points contribute with the same weight, then λ = 1 and p g (λ, t j ) = 1. However, it is possible to enhance the importance of the maximum dengue case in each annual cycles for the fitting process. Let us definet g , for g = 1, Á Á Á, u, as the calendar time when occurred the maximum recorded dengue cases in the g-th annual cycle. But this maximum may occur at different time than that yielded by the model, which is denoted by t g , for g = 1, Á Á Á, u. In the error function, besides the equal contribution of all samples, we include the difference between the peaks of theoretical and observed number of cases (jiðt g Þ À i obs ðt g Þj) and, also, the time when both peaks occurred (j t g À t g j). In the minimization process, we impose that the difference between the peaks i(t g ) and i obs ðt g Þ must be very close, then weight λ 1 > 1 is associated to it. On the other hand, the times when the peaks occur must be also very close, then we introduce the penalty factor λ 2 . In this case, the error function Eq (15) becomes where the second term is the penalty when the times between observed and theoretical peaks are distant. During the estimation procedure, the dynamical system of equations are simulated by providing initial conditions. The goal of an arbitrary initial condition is to eliminate the transient effects on the dynamics of mosquito population. For this reason the initial conditions at t = t 0 for all models are the same, given by where t 0 corresponds to January 1 st , 1991, and N 0 = 847,595. Dengue infection is assumed to be introduced in October 1 st of years 1995, 2000, 2004 and 2008, respectively named periods 1, 2, 3 and 4. The introduction of infectious individuals is described by, for instance for FD model, ðFDÞ : where δ(t − t d ) is the Dirac delta function, and i d = I d /N d , with N d being the size of human population at time t d . In other words, a number I d is transferred from susceptible population to infectious population at exactly t = t d . This instantaneous transference is captured by the following boundary conditions at t = t d given by where l d = l(t d ), p d = p(t d ) and m 1d = m 1 (t d ). The boundary conditions given by Eq (20), however, are taken at different times. The boundary conditions correspond to four distinct and uncorrelated epidemics (this strong assumption can be weakened by introducing more than one serotypes in a future modeling). For each period, we estimate the parameters β m , β h and prop considering the four models. We show, using the estimated parameters β m and β h , the following figures: (1) fraction of cases of dengue multiplied by the term prop, that is, prop × cases(t j ), and model based fraction i; (2) the proportion of infectious mosquitoes m 3 /m; (3) the partial effective reproduction numbers R h ef and R m ef ; and (4) the effective reproduction number R ef . The effective reproduction number R ef is given by Eq (10), with the partial effective reproduction numbers R h ef and R m ef being given by Eqs (11), (12), (13) and (14), respectively for models FD, TMAL, PMAL and SIR.
The mean value for the effective reproduction number R ef over the period of time [t 1 , t 2 ], denoted by R ef , is calculated by where t 1 and t 2 are the calendar year when one of the periods of estimation begins and ends, respectively. Similarly, the mean values for the partial effective reproduction numbers R h ef and R m ef , denoted respectively by R h ef and R m ef , can be calculate.

Data and data sources
We present a brief description of data and data source.

Results
We present the fittings of the incidence of dengue observed in the City of Campinas applying the four models. The tables that summarize the fittings provide the estimated β m , β h and prop, with the corresponding value of error Er, using Eq (18) for the period denoted by [t 1 , t 2 ], and the weighting factors λ 1 and λ 2 ; the value of the fraction of recovered individuals obtained at the end of the period of estimation r f = r(t 2 ); and the mean values for the partial and overall effective reproduction numbers R ef , R h ef and R m ef , using Eq (21). The four models are fitted against dengue incidence data recorded during four periods in the City of

Frequency dependent model-FD
The frequency dependent model, given by Eqs (4), (5) and (6), is considered to fit data from the City of Campinas. The transmission terms are such that only m 3 (density) is multiplied by D/N. In Table 6 we provide the estimated β m , β h and prop, and other parameters calculated using these estimated values. The value for λ 2 was chosen to the penalty factor be proportional to the sum of the difference between observed and calculated fraction of infectious individuals.
The figures corresponding to the four periods are presented separately. For the four periods, we show, using the estimated parameters β m and β h , the following figures: (1) fraction of cases of dengue multiplied by the term prop, that is, prop × cases(t j ) (labelled by O, red), and model based fraction i (labelled by E, black); (2) the fraction of infectious mosquitoes m 3 /m; (3) the partial effective reproduction numbers R h ef (labelled by h, black) and R m ef (labelled by m, red); and (4) the effective reproduction number R ef . Fig 7 shows the fittings for Period 1 (1995-2000). The third cycle (year 1998) is well adjusted (the fitted curve is slightly delayed), and for other cycles, the fitted epidemics are more severe than real epidemics. The fraction of recovered individuals reaches 22.5% at the end of this period. The maximum peak of the proportion of infectious mosquitoes occurs in year 1998, with value reaching 0.8%. The average partial effective reproduction numbers are R m ef ¼ 0:33 and R h ef ¼ 4:24, and for the overall number, R ef ¼ 2:15.      All epidemics cycles are well adjusted, except the last cycle, when the estimated epidemics is more severe than the real one. The fraction of recovered individuals reaches 11.9% at the end of this period. The maximum peak of the proportion of infectious mosquitoes occurs in year 2011, with value reaching 0.24%. The average partial effective reproduction numbers are R m ef ¼ 0:14 and R h ef ¼ 13:2, and for the overall number, R ef ¼ 2:71. The FD model fitted better the third and fourth periods, while in the first and second periods, only one and two, respectively, of the epidemics peaks are well adjusted. In the second period, the order of magnitude between peaks of R h ef and R m ef situates at 10 2 , while for other three periods, the difference was 10. Notice that the mean partial numbers follow these order of magnitude. The fraction of recovered individuals reached around 10%, except for period 1, when reached 22%. Comparing fittings from periods 1 (the lowest R ef ) and 3 (the highest R ef ), we observe that: (1) a prominent epidemics in the third cycle, (2) higher peak of infectious mosquitoes and prop, but lower R h ef , in the period 1 than period 3, and (3) period 1 presented the highest error E r and initial infectious individuals I d .
The FD model showed a plateau in curves for R h ef , but triangular shape for R m ef curves. Triangle shaped curve is also presented by the overall R ef . These findings show strong dependency of infective term among mosquito population with abiotic variations, but quite insensitive among humans. However, the strong seasonal dependency of mosquito population prevails in the overall transmission of dengue.

Pseudo mass action law model-PMAL
The pseudo mass action law model, given by Eqs (4), (5) and (7), is considered to fit data from the City of Campinas. The transmission terms are such that i and m 3 (density) are multiplied by N and D, respectively. In Table 7 we provide the estimated β m , β h and prop, and other parameters calculated using these estimated values.
The figures corresponding to the four periods are presented separately. For the four periods, we show, using the estimated parameters β m and β h , the following figures: (1) fraction of cases of dengue multiplied by the term prop, that is, prop × cases(t j ) (labelled by O, red), and model based fraction i (labelled by E, black); (2) the fraction of infectious mosquitoes m 3 /m; (3) the partial effective reproduction numbers R h ef (labelled by h, black) and R m ef (labelled by m, red); and (4) the effective reproduction number R ef .  Table 7. Model PMAL: Summary of fitted parameters and epidemiological values for four periods, with r f = r(t 2 ) and I d = i(t 1 )N(t 1 ).   1 (1995-2000). The fitting of incidence data is quite similar than that obtained for FD model. The fraction of recovered individuals reaches 23.8% at the end of this period. The maximum peak of the proportion of infectious mosquitoes occurs in year 1999, with value reaching 0.3% (in FD model occurred at 1998, with peak around 2.5-fold higher). The average partial effective reproduction numbers are R m ef ¼ 0:12 and R h ef ¼ 11:7, and for the overall number, R ef ¼ 2:23.      The PMAL model fitted better the third and fourth periods, while in the first and second periods, only one and two, respectively, of the epidemics peaks are well adjusted. In all periods, the order of magnitude between peaks of R h ef and R m ef situates at 10 2 . The fraction of recovered individuals reached around 25% (periods 1 and 2) and 15% (periods 3 and 4). Comparing fittings form periods 1 (the lowest R ef ) and 3 (the highest R ef ), the results are similar than that observed in FD model.
Both PMAL and FD models presented quite the same pattern of fittings for all periods. Also, both PMAL and FD models showed same characteristics with respect to effective reproduction numbers.

True mass action law model-TMAL
The true mass action law model, given by Eqs (4), (5) and (8), is considered to fit data from the City of Campinas. The transmission terms i and m 3 /m (fraction) appear alone. In Table 8 we provide the estimated β m , β h and prop, and other parameters calculated using these estimated values.
The figures corresponding to the four periods are presented separately. For the four periods, we show, using the estimated parameters β m and β h , the following figures: (1) fraction of cases of dengue multiplied by the term prop, that is, prop × cases(t j ) (labelled by O, red), and model based fraction i (labelled by E, black); (2) the fraction of infectious mosquitoes m 3 /m; (3) the partial effective reproduction numbers R h ef (labelled by h, black) and R m ef (labelled by m, red); and (4) the effective reproduction number R ef . Fig 15 shows the fittings for Period 1 (1995-2000). The fitting is poor, with the most severe epidemics overlapping the next epidemics. All estimated epidemics cycles are severe. The fraction of recovered individuals reaches 51.3% at the end of this period. The maximum peak of   Table 8. Model TMAL: Summary of fitted parameters and epidemiological values for four periods, with r f = r(t 2 ) and I d = i(t 1 )N(t 1 ).    , and for the overall number, R ef ¼ 2:12. The TMAL model fitted reasonably the third period. The epidemics overlap each other, eliminating the inter-epidemics observed in previous FD and PMAL models. Other feature is the increasing phase of epidemics occurs during winter (middle of the year) in some epidemics cycles. In the first and fourth periods, the quotient of the magnitude of the peaks of R h ef and R m ef is both of order 10 2 , while for second and third periods, the order of magnitudes were 10 3 and 10, respectively. Notice that the mean of the partial numbers ( R h ef and R m ef ) follow these order of magnitudes. The fraction of recovered individuals reached between 50% and 80%. Comparing fittings form periods 1 (the lowest R ef ) and 3 (the highest R ef ), we observe that: (1) a prominent epidemics in the third cycle, (2) lower peak of infectious mosquitoes and lower R h ef , but higher   prop, in the period 1 than period 3, and (3) period 1 presented the highest error E r and initial infectious individuals I d .
The TMAL model showed a plateau in both curves for R h ef and R m ef curves. This shape is also presented by the overall R ef . These findings show quite insensitive dengue transmission among humans and mosquitoes with respect to abiotic variations. However, there are measurable decreasing in the partial effective number in humans R h ef and the overall effective number R ef after epidemics. This decreasing is much more pronounced than in the previous two models.

Special incidence rate model-SIR
The special incidence rate model, given by Eqs (4), (5) and (9), is considered to fit data from the City of Campinas. The transmission terms are such that i and m 3 /m (fraction) are divided by Dm and N, respectively. In Table 9, we provide the estimated β m , β h and prop, and other parameters calculated using these estimated values.
The figures corresponding to the four periods are presented separately. For the four periods, we show, using the estimated parameters β m and β h , the following figures: (1) fraction of cases of dengue multiplied by the term prop, that is, prop × cases(t j ) (labelled by O, red), and model based fraction i (labelled by E, black); (2) the fraction of infectious mosquitoes m 3 /m; (3) the partial effective reproduction numbers R h ef (labelled by h, black) and R m ef (labelled by m, red); and (4) the effective reproduction number R ef . Fig 19 shows the fittings for Period 1 (1995Period 1 ( -2000. The fitting is poor, with the most severe epidemics overleaping the next epidemics. All estimated epidemics cycles are severe. The fraction of recovered individuals reaches 39.6% at the end of this period. The maximum peak of the proportion of infectious mosquitoes occurs between years 1997 and 1998, with value reaching 0.16%. The average partial effective reproduction numbers are R m ef ¼ 0:04 and R h ef ¼ 52:6, and for the overall number, R ef ¼ 2:73. Fig 20 shows the fittings for Period 2 (2000Period 2 ( -2004. The fitted second and third cycles (years 2002 and 2003) show much more severe epidemics than the real ones. The first and fourth peaks are absent. Again, the fitting is poor. The fraction of recovered individuals reaches 78.5% at the end of this period. The maximum peak of the proportion of infectious mosquitoes occurs Table 9. Model SIR: Summary of fitted parameters and epidemiological values for four periods, with r f = r(t 2 ) and I d = i(t 1 )N(t 1 ).  The SIR model fitted very poorly in all periods. In the first and the second periods, the quotients of the magnitudes of the peaks of R h ef and R m ef are of order 10 4 and 10 3 , respectively, while for other two periods, the order was 10 2 . The fraction of recovered individuals reached around 40% in first and fourth periods, with the highest 78% in the second period, and the lowest 30% in the third period. In this model, the beginning of annual epidemics occurs at the middle of the year.
Both TMAL and SIR models presented similar pattern of fittings for all periods. TMAL model presented plateau for both R h ef and R m ef , but SIR model reverted the pattern found in FD and PMAL models: triangle shape for R h ef and plateau for R m ef curves. This change practically eliminated the overlapping of epidemics obtained from TMAL model, however, in all epidemics cycles, the increasing phase occurred in the winter season.  then all models are equivalent.
Here, the non-autonomous models have parameters as well as the human and mosquito populations being time dependent. Two types of models were considered depending on the way the mosquito and human populations interaction is described. Type 1 modelling (FD and PMAL) considered the infection among susceptible humans being proportional to the density of infectious mosquitoes, that is, m 3 . In the Type 2 modelling (TMAL and SIR), the infection term was proportional to the fraction of infectious mosquitoes, that is, m 3 /m.
Let us discuss Type 1 and 2 modellings. In Type 1 modelling, variables D and N appear differently in the dynamical systems: in FD model, the quotient D/N appears multiplying the density of the infectious mosquitoes m 3 , while in PMAL model, m 3 is multiplied by only D, but the fraction of infectious individuals i comes multiplied by N. This difference in the equations, however, resulted in a quite similar fittings.
The total transmission coefficients estimated using FD model were 10 6 order of magnitude in comparison with the coefficients fitted from PMAL model. In the latter model, the transmission coefficients are per-capita, but the order of D and N are 10 6 , and the total coefficients are per-capita coefficients multiplied by D and N.
FD model showed a better fittings than PMAL (Er was lower, except period 4, which is quite similar). Estimation with FD model resulted in lower risk of dengue, presenting lower final fraction of recovered individuals and mean reproduction number, in the latter case only period 4 resulted in higher R ef . In FD and PMAL models, epidemics cycle approximately initiates at the beginning of the calendar year and ends in September, encompassing an epidemic period of around 9 months. There is clearly inter-epidemics periods of 3 months, with practically absent dengue cases.
In both models, when a severe epidemics (period 3) or epidemics with relatively same intensity (period 4) are present, then dengue cases are well adjusted. Periods 1 (other peaks are not negligible in comparison with period 3) and 2 (epidemics with decreasing intensity in comparison with period 4) did not fit so well. In general, both models presented good fittings.
Both models showed curves for R h ef and R m ef with similar pattern: curves with plateau for R h ef and triangle shaped curves for R m ef . This is an expected finding, since the temperature and rain dependent parameters have little influence in human population (s), while mosquito population is strongly influenced (m 1 ). Notice that the plateau in R h ef is slightly decreased after epidemics, due to decrease in susceptible humans. The overall effective reproduction number R ef follows the trend of R m ef . The value of R ef is the potential risk of dengue epidemics, with the peak of the risk situating in March-April. The susceptible populations of humans and mosquitoes plus the entomological parameters determine the intensity of the maximum hazard.
In Type 2 modelling, variables Dm(= M) and N appear only in SIR model: m 3 /m is divided by N, but the fraction of infectious individuals i comes divided by Dm. In TMAL model, m 3 /m and i do not come accompanied by M and N. This difference in the equations resulted in different fittings.
In TMAL and SIR models, the incidence data were poorly fitted. When epidemics start, the annual cycles are overlapped in TMAL model, eliminating the inter-epidemics periods [22], but SIR model persistently initiated epidemics in unfavorable season. Type 2 modelling yielded fittings with more intense dengue epidemics than Type 1 modelling.
In TMAL model, the shape of the epidemics is generally similar than models of Type 1 (annual cycle begins at the beginning of the year and range of epidemics around 8 months), but there are non-zero transmission in winter seasons. The annual epidemics periods are highly overlapped, with the absence of inter-epidemic periods. However, in SIR model, epidemics cycle approximately begins in September/October of the calendar year and ends in April/May of the next year, encompassing an epidemic period of around 9 months, but practically eliminated the overlapping phenomenon.
The transmission of dengue epidemics in unfavorable seasons are due to the use of fraction of infectious mosquitoes. In unfavorable seasons, the number of infectious mosquitoes M 3 (or its density m 3 = M 3 /D) goes to zero. However, the fraction m 3 /m does not tend to zero, which results in non-zero risk of dengue transmission. This occurs in TMAL model. But, in SIR model, besides the fraction m 3 /m, there is another fraction m 1 /m (in the transmission term among humans). This term allowed humans to be strongly dependent on the abiotic variation with respect to dengue infection. This second fraction of mosquito population results in a more possibility of dengue transmission in unfavorable seasons. Hence, in SIR model, the dengue transmission during unfavorable seasons (inter-epidemics) is enhanced in comparison with TMAL model, but the overlapping is practically eliminated.
Notice that the triangle shape in Type 1 modelling for R m ef was changed to a plateau in Type 2 modelling. The division m 1 /m, see Eq (13), somehow diminished the strong variation of mosquito population with temperature and rain, leading its quotient to a constant value. The plateau in Type 1 modelling for R h ef was changed by a sharp triangle shape only in SIR model. The division s/m, see Eq (14), in somehow introduced a strong variation with temperature and rain in R h ef , and by the fact that R h ef appears in the denominator, its effect is much enhanced, and the triangle shape is very sharp.
Both Type 1 and Type 2 modellings yielded similar mean values for the effective reproduction number, except periods 3 and 4 of SIR model. Type 1 modelling resulted in a triangle shape with relatively large basis for the risk of overall dengue (R ef ). For this reason the peak of R ef situated in the range 10-20. However, in Type 2 modelling, we obtained two shapes: a plateau in TMAL model, while a very sharp triangle shape is present in SIR model. For this reason, the peak of R ef situated around 8 for TMAL model, and in the range 60-180 for SIR model. For this reason, SIR model presented the highest peaks of infection among humans and mosquitoes. In general, Type 2 modelling resulted in more severe dengue epidemics to explain the incidence data from the City of Campinas.
FD and PMAL (Type 1 modelling) reproduced abiotic fluctuation mainly in mosquito population (R m ef ), while human population (R h ef ) was less affected. Additionally, the risk of dengue (R ef ) increases at the end of the year, and attains its maximum at month February. However, in Type 2 modelling, the risk of dengue is practically constant during entire year for TMAL model, but for SIR model, the highest risks of dengue situated in the middle or at the end of the year, with very narrow basis for R ef curve.
Let us now discuss the fitting of the parameters β m , β h and prop.
The quotient between the estimated values for β m and β h have magnitudes of order 6 according to the models considered. For PMAL model, where β m and β h are per-capita transmission coefficients, estimations of both parameters have order of 10 −9 /10 −8 . For FD and TMAL models, where β m and β h are total transmission coefficients, estimations of both parameters have order of 10 −3 /10 −2 . For SIR model, where β m and β h are bi-total transmission coefficients, estimations of both parameters have order of 10 3 /10 5 . The quotient of order of magnitude of 10 6 results due to the fact that the human population is around 10 6 , as well as the value D for mosquito population. This difference is in accordance with Eq (22).
The estimated parameter prop and an arbitrary initial number of infectious individuals I d play important role in the estimation process. The functional given by Eq (18) presents several minima, with values very close. Besides the best fittings for β m and β h , we allowed to prop to vary in order to obtain the fraction of recovered individuals in the range 20-50% [23], and the peak of the fraction of infectious mosquitoes in the range 0.5-5% [24]. We believe that these ranges are biologically acceptable, than for instance, 95% or 2% of recovered individuals at the end of the epidemics, or a peak of 25% of infectious mosquitoes. In the schistosomes transmission modelling, biological acceptance was adopted to chose better fitting [25] [26]. Finally, I d was allowed to vary in order to adjust the initial phase of the epidemics (larger I d , earlier occurs the ascending phase of epidemics). This degree of freedom was introduced by the assumption that the entire population is under a new epidemics, hence the number of recovered individuals due to previous circulating serotype of dengue is set zero.

Discussion
The incidence data recorded in the City of Campinas were split in four periods: Period 1 (1995Period 1 ( -2000, Period 2 (2000Period 2 ( -2004, Period 3 (2004Period 3 ( -2008 and Period 4 (2008Period 4 ( -2012, where each period begins at month October and ends at September. Each one of these four periods was fitted considering models with unique serotype of dengue. Type 1 modelling (FD and PMAL) fitted well all periods of incidence of dengue in the City of Campinas, differently from the fittings yielded by Type 2 modelling (TMAL and SIR).
Type 1 modelling showed a clear inter-epidemic period between consecutive annual epidemics cycles. The epidemics began at the end of the previous year, and ended at the middle of the next year. However, Type 2 modelling presented overlapping of epidemic cycles, and the initial phase of an epidemics could occur in unfavorable season.
By observing the partial effective reproduction numbers (R h ef and R m ef ), Type 1 modelling showed a strong seasonality in the risk of mosquitoes to be infected by dengue, but the risk of infection among humans was practically insensitive to seasonal variations. However, Type 2 modelling presented mosquitoes being insensitive to seasonal variations, but humans can be strongly dependent (SIR model) or not (TMAL model) on the seasonal variations with respect to dengue infection. Comparing the overall effective reproduction number (R ef ), in Type 1 modelling there is a triangle shaped curve, however, a plateau (TMAL model) and very sharped triangle shape (SIR model) in Type 2 modelling.
We conclude that models FD and PMAL fitted better incidence data from the city of Campinas. The results yielded by these models are in better agreement with biological reasonings (periodicity of epidemics in favorable seasons, and well established inter-epidemics period during unfavorable seasons). This conclusion is valid for regions where there are strong variations in abiotic factors. When small abiotic variations occur, all models could be used indistinctly.
The prediction of next dengue epidemics is a hard task. We address this question, choosing the model FD estimated with respect to Period 1, considering the following factors: (1) temperature and precipitation, (2) initial number of infectious humans introduced in a community, and (3) the time of this introduction.
To evaluate the influence of temperature and precipitation on the incidence of dengue, we fix the estimated parameters β m and β h for Period 1 of Model FD and the initial number I d of infectious humans (see first column of Table 6). Using these fixed values, we obtain the fraction of infectious individuals i for Periods 2, 3 and 4. The time of introduction of infectious individuals in these periods is October 1 st of the initial year of the corresponding period. Fig 23 shows the strong influence of temperature and precipitation. Notably, period 2 presents lower risk of dengue incidence than period 1 (approximately one fourth), while periods 3 and 4 are practically safe from dengue incidence. In these last periods, temperature and precipitation affected quite similarly, being the last period slightly more hazardous. Another source of uncertainty is the number of infectious individuals that trigger an epidemics. To evaluate the influence of this number, we consider only Model FD with estimations given in Period 1, by varying I d . Fig 24 shows the fraction of infectious individuals i for: I d = 1 (a), I d = 25 (b), I d = 150 (c), and I d = 1000 (d). As I d increases, the risk of dengue incidence increases, anticipating the highest peak of epidemics. For small I d (Fig 24(a) and 24(b)), the most severe epidemics occurs after 4 years of the introduction of dengue cases, and as I d increases, the highest peak is anticipated: In Fig 24(c) and also Fig 7(a), the highest peak is anticipated in one year, while in Fig 24(d), in two years.
Finally, let us evaluate how the time of the introduction of infectious humans affects the incidence. Again, we consider only Model FD with estimations given in Period 1, now varying t 1 . Fig 25 shows the fraction of infectious individuals i for t 1 corresponding to year 1995, in the day: July 1 st (a), September 1 st (b), December 1 st (c), and December 31 st (d). We observe a pattern similar than the one in the previous case. Notably, if dengue cases are introduced in the beginning of December (Fig 25(c)), the first epidemics is as severe as the highest peak occurred in Fig 25(a) and 25(b), cases of early introduction of dengue. By delaying the introduction of  dengue cases, we observe that initial peaks are decreasing, but latter epidemics are more severe ( Compare Fig 25(c) and 25(d)).
In future works, the previous models will be improved by introducing the concomitant action of two or more serotypes of dengue instead of just one serotype as in the present work, in order to fit the entire period of dengue incidence (1995-2012) recorded from the City of Campinas. A feature of the estimation procedure used in the present work was the introduction of a parameter prop in order to encompass the unnotified cases of dengue. However, another possibility, also to be considered in a future work, is to modify the model by introducing a compartment of mild cases of dengue, which could be considered asymptomatic [27].
In this paper, the temperature dependency of the entomological parameters were well defined from experiments conducted in laboratory. However, the parameters depending on precipitation were chosen based on simulations [14]. Sensitivity analysis of these parameters must be done to evaluate their effects on the dengue transmission (especially in extreme conditions).
Supporting Information S1 Text. Appendix containing mathematical analysis of autonomous model and description of the estimation method. (PDF)