Discrete Dynamics in Nature and Society
Volume 2018 (2018), Article ID 6456107, 20 pages
https://doi.org/10.1155/2018/6456107

Research Article

Dynamics and Stability Analysis of a Brucellosis Model with Two Discrete Delays

Paride O. Lolika1,2 and Steady Mushayabasa¹

¹University of Zimbabwe, Department of Mathematics, P.O. Box MP 167, Harare, Zimbabwe
²University of Juba, Department of Mathematics, P.O. Box 82 Juba, Central Equatoria, South Sudan

Correspondence should be addressed to Steady Mushayabasa; steadymushaya@gmail.com

Received 23 April 2018; Revised 29 June 2018; Accepted 17 July 2018; Published 7 August 2018

Academic Editor: Abdul Qadeer Khan

Copyright © 2018 Paride O. Lolika and Steady Mushayabasa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We present a mathematical model for brucellosis transmission that incorporates two discrete delays and culling of infected animals displaying signs of brucellosis infection. The first delay represents the incubation period while the second account for the time needed to detect and cull infectious animals. Feasibility and stability of the model steady states have been determined analytically and numerically. Further, the occurrence of Hopf bifurcation has been established. Overall the findings from the study, both analytical and numerical, suggest that the two delays can destabilize the system and periodic solutions can arise through Hopf bifurcation.

1. Introduction

Brucellosis is one of the neglected zoonotic diseases that remains a major public health problem world over, especially in Middle Eastern countries, southern Europe and North Africa, countries in South and Central Asia, sub-Saharan Africa, Mexico, the Caribbean, and countries in South and Central America [1], with an annual occurrence of more than 500 000 cases [2].

In animals, brucellosis is usually transmitted through direct contact between a susceptible and an infectious animal or indirectly, i.e., when a susceptible animal ingest contaminated forages or the excrement containing large quantities of bacteria, generally discharged by infected animals [3]. In humans, however, the majority of the infections result from direct or indirect exposure to infected animals or consumption of raw animal products such as unpasteurized milk or cheese [4]. Since human-to-human transmission of the disease is extremely rare [5], the ultimate management of human brucellosis can be achieved through effective control of brucellosis in livestock. Some researchers postulate that eradication of brucellosis in animals can be attained by combining vaccination with test-and-slaughter programs [1].

Mathematical models have proved to be essential guiding tools for epidemiologists, biologists, and policy makers. Models can provide solutions to phenomena which are difficult to measure practically. Recently, a number of mathematical models have been proposed to study the spread and control of brucellosis (see, for example, [3, 6–15] and references therein). A limitation of these previous studies, however, is the noninclusion of the time taken before an infectious animal is detected and culled, despite the fact that in many countries where the disease is endemic lack of financial and human resources often results in delay of detection and culling of infectious animals. The size of this delay may play an important role in minimizing the spread of the disease in the community.

It is therefore essential to gain a better and more comprehensive understanding of the effects of time delay on brucellosis transmission and control. Prior studies have shown that epidemic models with time delay often exhibit periodic solutions and as a consequence understanding the nature of these periodic outbreaks plays a crucial role in designing policies that can successfully control the disease. In fact, a recent analysis of brucellosis dataset in countries where the disease is endemic has shown that the disease incidences exhibit a strong periodic behavior with mortality and morbidity of the disease concentrated in a few months each year [16, 17]. Understanding the impact of such seasonal variations is crucial in managing the spread of the disease in the community.

Our main goal in this study is to explore the dynamics and stability analysis of a brucellosis model with two discrete delays. Hence we formulated a mathematical model that incorporates two discrete delays. The first delay represents the incubation period while the second delay accounts for the time taken to detect and cull infectious animals. In addition, we subdivide the total animal population into classes of susceptible, vaccinated, infectious undetected, and infectious detected animals. In certain situations immediate slaughter of detected animals may not be feasible and more often these animals are isolated from the rest. However, due to lack of financial and human resources, in addition to lack of knowledge and attitude of farmers, isolation of detected animals has not been a successful practice in most developing nations where animal infections are rampant. Thus in our modelling process we assume that a proportion of detected animals that are not immediately culled are also responsible for disease transmission. Utilizing both analytical and numerical results we have demonstrated that the two delays can destabilize the system and lead to Hopf bifurcation.

The article is organised as follows. The model description is given in Section 2. Analytical and numerical results are given in Sections 3 and 4, respectively. We end with Section 5 of conclusions.

2. Mathematical Model

We subdivide the total animal population into compartments of susceptible , vaccinated , undetected infectious animals , and infectious detected and unculled . Although brucellosis can be transmitted indirectly (environmental transmission), prior studies [6, 8] suggest that indirect transmission plays a relatively small role on the spread of brucellosis, and as such we have ignored this aspect in our study. Brucellosis dynamics in this study are governed by the following autonomous system:

where is the recruitment rate through birth, is the natural death rate, is the disease direct transmission rate, is fraction of detected animals that have been culled, is the vaccination rate, is the vaccination waning rate, represents the latency delay, is the delay in detection, is the modification factor, is the rate at which animals are detected and quarantined, is the culling rate of detected animals, and is disease-induced death.

3. Analytical Results

3.1. Initial Conditions

The appropriate space for system (1) is the Banach space of continuous functions mapping the interval into equipped with subnorm where . From the standard results of functional differential equations [18] it follows that, given any initial conditions , there exists a unique solution of system (1), which satisfies , the initial conditions are given by

where , with and .

3.2. The Basic Reproduction Number

By direct calculation, we find that system (1) has a disease-free equilibrium , given by , with

By utilizing the next-generation matrix method [19], one can deduce that the basic reproduction number of model (1) is

with , and

3.3. Stability of the Disease-Free Equilibrium

In this section, we study the local and global stability of the disease-free equilibrium.

Theorem 1. The disease-free equilibrium of model (1) is locally asymptotically stable when and unstable when .

Proof. To study the local stability of the disease-free equilibrium , we linearized system (1) about this point and obtained the characteristic equation, given by the following determinant:

where is the eigenvalue.
From (5) the characteristic equation is

Clearly, and are eigenvalues and the other two can be determined from the following equation:

Let

Through direct calculation one can easily verify that is an analytic function and it follows that

Now we proceed to investigate the distribution of the solutions of (7) in the following cases:

(a)If , then . Since the derivative for , and , (7) has no zero root and positive real roots for all positive and . Now we assume that the solution of (7) does not have any purely imaginary roots , for some , . Then by computation, must be positive real root of

 If , (10) has no positive roots. Hence (7) does not have any purely imaginary roots. We can easily see that the roots of all have negative real parts when . By the implicit function theorem and the continuity of , we know that all roots of (7) have negative real parts for positive and which implies that is stable.(b)If then . Since the derivative for , and , (7) has a simple zero and no positive root for all positive and . By the same argument in case , we can obtain that all roots of (7) have negative real parts for positive and except a zero root. Thus is a degenerate equilibrium of codimension and is stable except in one direction.(c)If , then . Since we have and for and , (7) has a unique positive real root for all positive and and is unstable.

Theorem 2. The disease-free equilibrium of model (1) is globally asymptotically stable when and unstable when .

Proof. We denote by the translation of the solution of the system (1); that is, , where and . We consider the following Lyapunov functional:

Taking the derivative of along the solutions of (1) gives

Therefore, holds if . Furthermore, if . Therefore, the largest invariant set of is a singleton such that , , . It follows from LaSalle’s invariance principle [20] that the disease-free equilibrium of system (1) denoted by is globally asymptotically stable whenever . This completes the proof of Theorem 2.

3.4. Disease Persistence

System (1) is said to be uniformly persistent if there exists a constant such that any solution of (1) satisfies

Now we give a result on the uniform persistence of system (1). To proceed we introduce the following notation and terminology. Denote by the family of solution operators corresponding to (1). -limit set of consists of such that there exists a sequence as with as .

Theorem 3. System (1) is uniformly persistent, if it satisfies .

Proof. Let , which is relatively closed in .
In what follows we demonstrate that is positively invariant for . From the third and fourth equations of (1) we have

Since , we have , and it follows from (14) that

Thus is positively invariant for .
We set

We claim that

Assuming , ; it suffices to show that , . If it is not true, then there exists such that either (a) , ; or (b) , . For case (a), from the fourth equation of (1), we have

Hence there is such that , . On the other hand, from there exists such that , . Thus we have , , which contradicts the assumption that . Similarly, we can obtain a contradiction for case (b). This proves the claim (17).
Let , where is the global attractor of restricted to . We show that . In fact, from and the first and second equation of (1), we have and . Thus, is the isolated invariant set in .
Next we show that . If this is not true, then there exists a solution such that

For any sufficiently small constant , there exists a positive constant such that , , . For the constant given above, it follows from the third and fourth equations of (1) that

If , then, by a standard comparison argument and the nonnegativity, the solution of the following monotone system

with initial condition converges to as well. Thus , where is defined by

Differentiating with respect to time yields

Since , we have for a sufficiently small . Therefore goes to either infinity or a positive number as , which leads to a contradiction with . Thus we have . Define by . It is clear that and . Thus by [21] theorem 3 we have for some constant . Let where is the constant such that , . We showed that , , , . This completes the proof of Theorem 3.

3.5. Existence of the Endemic Equilibrium

Theorem 4. If , model (1) admits a unique endemic equilibrium.

Proof. The endemic equilibrium of model (1) is determined by the following:

From the last equation in (24) we have

The first two equations in (24) give

For , substituting (25) into the third equation in (24) gives

Substituting (26) into (27) yields

Direct calculations show

where

then function is monotonic decreasing for , and then we can define the following function:

Therefore, by monotonicity of a function there exists a unique endemic equilibrium

3.6. Stability of the Endemic Equilibrium

In this section, we will investigate the local and global stability of the endemic equilibrium point.

Theorem 5. The endemic equilibrium of the system (1) is locally asymptotically stable if and conditions (36) are satisfied.

Proof. The characteristics equation of system (1) on the infected equilibrium is given by the following determinant:

with

After some algebraic manipulations one can show that the characteristic equation has the following form:

with

By the Routh-Hurwitz criterion, all roots of the characteristics equation (34) have negative real parts and the endemic equilibrium of system (1) is locally asymptotically stable if , if and only if , ,

Now, we wish to explore if there is a possibility of having complex roots with positive real part for , and , . We now proceed to explore the above cases as follows.
(a) If , , then the characteristics equation (34) becomes

with

Now we need to show that all roots of (37) have negative real parts for all . To do so, we show that (37) does not have any purely imaginary roots for all . We assume that with is a root of (37). Then must satisfy the following system:

Now, we square both sides of each equation above and add the resulting equations; thus must be a positive root of

where

Let ; then (40) becomes

One can observe that if , , then (42) has no positive roots. Therefore (37) does not have any purely imaginary roots for all so that all roots of the characteristic equation (37) have negative real parts and the endemic equilibrium of (1) is stable for all .
(b) If , then the characteristics equation (34) becomes

Using the same discussion as in the above case then (43) can be written as

with

It follows that all roots of (43) have negative real parts for when , and this implies that endemic equilibrium is locally asymptotically stable for . This completes the proof.

We now explore the global stability of the endemic equilibrium.

Theorem 6. If , then is globally asymptotically stable.

Proof. Let us consider the Lyapunov function:

Here,

The derivatives of are given by

Substituting the appropriate differentials from (1) we have

At endemic equilibrium, we have

Using the above constants, we have

The derivatives of are given by

The derivatives of are given by

Combining the derivatives of , for , we have

Note that

for all and , because the arithmetic mean is greater than or equal to the geometric mean. Further, note that is always nonpositive for any function , and if and only if . Hence, it follows that and, consequently, . Moreover, the largest invariant set of is a singleton, where , , , and . Using LaSalle’s invariance principle [20], we conclude that the endemic equilibrium point is globally asymptotically stable if

3.7. Hopf Bifurcation Analysis

In this section we determine criteria for Hopf bifurcation to occur using the time delay and as the bifurcation parameters to find the interval in which the infected equilibria is stable and unstable out of the same margins. Now to consider Hopf bifurcation we consider the cases , and and . Our analysis is as follows:

(a) When and we need to show that differentiating (37) with respect to we get

This gives

Thus,

with

Lemma 7. Suppose that , , are the roots of equation and is the largest positive root; then

In our case considering defined in (42), and assuming and as the largest positive root, we have

The above analysis can be summarized into the following theorem.

Theorem 8. Suppose that . If either or and is satisfied, and is the largest positive simple root of (42) then the infected equilibrium of model (1) is locally asymptotically stable when and unstable when where

when , a Hopf bifurcation occurs; that is, a family of periodic solutions bifurcates from as passes through the critical value .

(b) When and we also need to show that differentiating (43) with respect to we get

This gives

Thus,

Lemma 9. Suppose that , , are the roots of equation and is the largest positive root; then

In our case considering defined in (45), and assuming and as the largest positive root, we have

The above analysis can be summarized into the following theorem.

Theorem 10. Suppose that . If either or and is satisfied, and is the largest positive simple root of (45) then the infected equilibrium of model (1) is locally asymptotically stable when and unstable when where

when , a Hopf bifurcation occurs; that is, a family of periodic solutions bifurcates from as passes through the critical value .

From the analysis above, we can deduce that Hopf bifurcations may arise if conditions in Theorems 8 and 10 are satisfied. Thus, the introduction of time delay in system (1) can destabilize the system.

4. Numerical Results

In order to explore the behavior of system (1) and illustrate the stability of the equilibria solutions, we numerically solve system (1) using MATLAB and parameter values adopted from Table 1.

Table 1: Model parameters and variables and their baseline values.

Symbol
Definition
Value
Unit
Source

Elimination rate due to brucellosis
0.15

[8]

Fraction of infectious animals culled upon detection
0.5
-
[22]

Direct transmission rate

[6]

Vaccination waning rate
0.4

[8]

Natural elimination rate
0.25

[6]

Modification factor
0.18
-
[8]

Recruitment rate
76434
animals
[6]

Vaccination rate
0.316
year
[8]

Detection rate
Varied
year
-

Culling rate
0.15
year
[8]

In Figure 1 we illustrate the effects of varying the delay ( on the dynamics of system (1). Figures 1(a) and 1(b) demonstrate that the system approaches the stable disease-free or endemic equilibrium for and , respectively. One should note that, according to Theorems 2 and 6, the stability of the model steady states does not depend on the value of the time delays, but rather on the basic reproduction number , only. In addition, we observe that the range of values for the two time delays does not lead to periodic solutions but an increase in both delays translate to an increase in the infectious population, both detected and undetected.

(a)

(b)

Figure 1: Stability of the infected and free-infected equilibrium of model system (1) showing plots of and with varying delay . The direction of the arrow depicts an increase in delay with a step size of starting from 2.0 to 10. The blue patterns in both (a) and (b) highlight brucellosis dynamics when while the red pattern is for Initial population levels were assumed as follows: animals, animals, animals, and animals.

Figure 2 depicts the numbers of infectious undetected and infectious detected animals over time with varying delays. The results clearly show that the incubation related delay has more influence on shaping the dynamics of brucellosis compared to the culling related delay More precisely, the incubation period delay significantly increases the infectious population (both detected and undetected) for and thereafter its impact will be the same as that of detection

(a)

(b)

Figure 2: Numerical solutions of model system (1) illustrating the effects of different time delay on brucellosis infection level in the community. Initial population levels were assumed as follows: animals, animals, animals, and animals.

In Figure 3 we illustrate the stability of the disease-free equilibrium with and (note that ). As we can observe, for certain parameter values and initial population levels, system (1) exhibits some periodic oscillation. Precisely, we note that the infected population ( and ) oscillates with a reduced amplitude from the start till when is approximately 400; thereafter the oscillations dies off the solutions converges to the disease-free equilibrium. These simulation results demonstrate the occurrence of periodic solutions through Hopf bifurcation for delay values and . In contrast, we can observe that there are no periodic oscillations for the uninfected populations and

(a)

(b)

(c)

(d)

Figure 3: Numerical solutions demonstrating the stability of equilibrium of model system (1) with . We set , , , , , animals , and the remainder retained the baseline values in Table 1. Further, we set the initial conditions as follows: animals, animals, animals, and animals.

In Figure 4, we demonstrate the dynamic for model system (1) with respect to the stability of infection-free equilibrium for different pair of delay values and from the simulation results we can conclude that both delays do not have a huge influence on the stability of disease-free equilibrium.

(a)

(b)

(c)

(d)

Figure 4: Dynamics of model system (1) for different values of , which illustrate the stability of infection-free equilibrium at . We set , , , animals , and the remainder retained the baseline values in Table 1. Further, we set the initial conditions as follows: animals, animals, animals, and animals.

In Figure 5 we observe that, for certain parameter values and initial population levels, system (1) may admit periodic oscillations when . As we can observe, when both the solutions of the infected and uninfected populations exhibit periodic oscillation for a certain period, before stability at endemic point is attained.

(a)

(b)

(c)

(d)

Figure 5: Stability of equilibrium of model system (1) with . The time delay was fixed to be and was fixed to be . We set the model parameters and variables as follows: , , , animals, animals, animals, and animals while the other parameter values are as in Table 1.

In Figure 6, we illustrate the dynamics for model system (1) with respect to the stability of endemic equilibrium for several pair of delay values . The results confirm that the incubation related delay has more influence on shaping the dynamic of brucellosis compared to the culling related delay .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: Numerical results of model system (1) for different values of , which demonstrate the stability of infected equilibrium at . We set the model parameters and variables as follows: , , , animals, animals, animals, and animals while the other parameter values are as in Table 1.

To explore influence of model parameters on the reproduction number , we perform a local sensitivity analysis of the basic reproduction number following the approach in [23]. The local sensitivity analysis will be useful on identifying parameters with greatest influence to change To this end, denoting by the generic parameter of system (1), we evaluate the normalised sensitivity index:

which indicates how sensitive is to a change of parameter Model parameters with positive index increase the value of whenever they are increased while those with a negative index decrease the value of whenever they are increased. We consider the parameter values in Table 1, and we set in order to evaluate the normalized sensitivity index and the results are depicted in Figure 7. Here, we observe that parameters , , , have a positive correlation with , such that increasing these parameters will increase . However, it is the increase of and that has the greatest influence to change . Precisely, increasing either or by 50% will increase by 50%. We also note that increasing parameters, , , , , , and , will lower the reproduction number.

Figure 7: Sensitivity index for with respect to model parameters that define it.

5. Conclusion

Zoonotic brucellosis remains a major public health problem in many developing nations. This is mainly attributed to several challenges associated with effective disease control in these nations. The challenges for effective control of brucellosis in developing nations range from inadequate veterinary personnel and vaccines as well as the failure by farmers to adhere to some of the aforementioned brucellosis control and eradication program activities. Furthermore, these challenges often lead to delay in detection and culling of infectious animals. In this study, we developed and analysed a mathematical model for brucellosis infection that incorporates two discrete delays. The first delay accounts for the latent period and the second delay represents the time taken to detect infectious animals. We computed the basic reproduction number and demonstrated that it is an important threshold quantity for stability of equilibria. By constructing suitable Lyapunov functionals, it has been shown that the model has a globally asymptotically stable infection-free equilibrium whenever the reproduction is less than unity. Further, it has been demonstrated that whenever the model reproduction number is greater than unity then the model has a unique endemic equilibrium point which is globally asymptotically stable. Numerical simulations are carried out to illustrate the main results.

Although culling of symptomatic animals is a relatively easy strategy to implement, it is worth noting that some studies suggest that culling of both infected and susceptible animal may be more effective [24, 25]. The rationale is that by decreasing the host density, then the number of contacts per unit time between animals is low, thereby reducing disease transmission. In [25] it was demonstrated that culling of both susceptible and symptomatic animals only can be effective whenever the number of infected host is above a certain critical level [25]. We expect to improve this study in our future work by developing brucellosis model(s) with time delay that will enable us to compare aforementioned aspects.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

E. M. Potter, “Brucellosis,” in Foodborne Infectionsand Intoxication, J. G. Morris and E. M. Potter, Eds., Elsevier Inc., 4th, 2013.

E. Shevtsova, A. Shevtsov, K. Mukanov et al., “Epidemiology of brucellosis and genetic diversity of Brucella abortus in Kazakhstan,” PLoS ONE, vol. 11, no. 12, Article ID e0167496, 2016.

J. Zhang, G.-Q. Sun, X.-D. Sun et al., “Prediction and control of brucellosis transmission of dairy cattle in Zhejiang Province, China,” PLoS ONE, vol. 9, no. 11, Article ID e108592, 2014.

J. A. Al-Tawfiq and A. AbuKhamsin, “A 24-year study of the epidemiology of human brucellosis in a health-care system in Eastern Saudi Arabia,” Journal of Infection and Public Health, vol. 2, no. 2, pp. 81–85, 2009.

J. Godfroid, “Brucellosis in wildlife,” Revue Scientifique et Technique de l'OIE, vol. 21, no. 2, pp. 277–286, 2002.

Q. Hou and X.-D. Sun, “Modeling sheep brucellosis transmission with a multi-stage model in Changling COUnty of Jilin Province, China,” Applied Mathematics and Computation, vol. 51, no. 1-2, pp. 227–244, 2016.

P. O. Lolika, S. Mushayabasa, C. P. Bhunu, C. Modnak, and J. Wang, “Modeling and analyzing the effects of seasonality on brucellosis infection,” Chaos, Solitons & Fractals, vol. 104, pp. 338–349, 2017.

Q. Hou, X. D. Sun, J. Zhang, Y. J. Liu, Y. M. Wang, and Z. Jin, “Modeling the transmission dynamics of sheep brucellosis in Inner Mongolia Autonomous Region, China,” Mathematical Biosciences, vol. 242, no. 1, pp. 51–58, 2013.

P. Lou, L. Wang, X. Zhang, J. Xu, and K. Wang, “Modelling Seasonal Brucellosis Epidemics in Bayingolin Mongol Autonomous Prefecture of Xinjiang, China, 2010–2014,” BioMed Research International, vol. 2016, Article ID 5103718, 17 pages, 2016.

C. Yang, P. O. Lolika, S. Mushayabasa, and J. Wang, “Modeling the spatiotemporal variations in brucellosis transmission,” Nonlinear Analysis: Real World Applications, vol. 38, pp. 49–67, 2017.

V. Racloz, E. Schelling, N. Chitnis, F. Roth, and J. Zinsstag, “Persistence of brucellosis in pastoral systems,” Revue Scientifique et Technique de l'OIE, vol. 32, no. 1, pp. 61–70, 2013.

J. Tumwiine and G. Robert, “A mathematical model for treatment of bovine brucellosis in cattle population,” Journal of Mathematical Modelling, vol. 5, no. 2, pp. 137–152, 2017.

J. Zhang, Z. Jin, L. Li, and X.-D. Sun, “Cost assessment of control measure for brucellosis in Jilin province, China,” Chaos, Solitons & Fractals, vol. 104, pp. 798–805, 2018.

M.-T. Li, Z. Jin, G.-Q. Sun, and J. Zhang, “Modeling direct and indirect disease transmission using multi-group model,” Journal of Mathematical Analysis and Applications, vol. 446, no. 2, pp. 1292–1309, 2017.

M.-T. Li, G.-Q. Sun, Y.-F. Wu, J. Zhang, and Z. Jin, “Transmission dynamics of a multi-group brucellosis model with mixed cross infection in public farm,” Applied Mathematics and Computation, vol. 237, no. 11, pp. 582–594, 2014.

K. Aune, J. C. Rhyan, R. Russell, T. J. Roffe, and B. Corso, “Environmental persistence of Brucella abortus in the Greater Yellowstone Area,” The Journal of Wildlife Management, vol. 76, no. 2, pp. 253–261, 2012.

W. Beauvais, I. Musallam, and J. Guitian, “Vaccination control programs for multiple livestock host species: An age-stratified, seasonal transmission model for brucellosis control in endemic settings the LCNTDR Collection: Advances in scientific research for NTD control,” Parasites & Vectors, vol. 9, no. 1, p. 55, 2016.

J. K. Hale and S. M. Verduyn Lunel, Introduction to Functional-Differential Equations, Springer, Berlin, Germany, 1993.

P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.

J. P. LaSalle, The Stability of Dynamical Systems, SIAM, Philadelphia, Pa, USA, 1976.

H. Smith and X. Zhao, “Robust persistence for semi-dynamical systems,” Nonlinear Analysis: Theory, Methods & Applications, vol. 46, pp. 6169–6179, 2001.

L. Zhou, M. Fan, Q. Hou, Z. Jin, and X. Sun, “Transmission dynamics and optimal control of brucellosis in Inner Mongolia of China,” Mathematical Biosciences and Engineering, vol. 15, no. 2, pp. 543–567, 2018.

L. Arriola, “Sensitivity analysis for quantifying uncertainty in mathematical models,” in Mathematical and Statistical Estimation Approaches in Epidemiology, G. Chowell, J. M. Hyman, L. M. A. Bettencourt, and C. Castillo-Chavez, Eds., pp. 195–248, Springer, New York, NY, USA, 2009.

L. Bolzoni, V. Tessoni, M. Groppi, and G. A. De Leo, “React or wait: which optimal culling strategy to control infectious diseases in wildlife,” Journal of Mathematical Biology, vol. 69, no. 4, pp. 1–5, 2013.

W. Zhou, Y. Xiao, and R. A. Cheke, “A threshold policy to interrupt transmission of West Nile Virus to birds,” Applied Mathematical Modelling, vol. 40, no. 19-20, pp. 8794–8809, 2016.