A General Model of Zoonosis Where the Diseases Can Only Be Transmitted from Animals to Humans

Zoonosis is an important factor affecting human economic development and population mortality. This paper introduces a general model of zoonosis, in which the diseases can only be transmitted from animals to humans, such as rabies, brucellosis and so on. The basic reproduction number 0 R is derived. And then the global stability of the disease-free equilibrium and endemic equilibrium models is analyzed by using the method of comparison principle and Lyapunov function. Next, a numerical analysis is performed to elaborate the consistency of theoretical and numerical results and to prove the practical significance of zoonosis research. The numerical results show that our models are applicable to zoonosis with animal size larger than or smaller than population size. Finally, in order to see the most important factor for the epidemic of zoonosis a sensitive analysis is analyzed.


Introduction
Infectious diseases seriously endanger human health and the development of community economy all over the world [1] [2] [3], historically, infectious diseases have caused millions of deaths.SARS attracted worldwide attention in 2003, more than 8000 infected people have been found worldwide, resulting in more than 900 deaths [4]; Dengue fever in Thailand has reached 5284 people, and the disease is prevalent in Ecuador, Paraguay, Australia and other countries to varying degrees by 2006 [5].Thus, researching and controlling the dynamic process of infections is very significant to society.Among numerous infectious diseases, zoonosis is an important factor affecting human economic development and population mortality.There are about 200 confirmed zoonoses in the world, such as avian influenza, rabies, malaria, brucellosis and so on.From 1996 to 2015, 30,300 human rabies cases were reported [6].Every year, about 50,000 people are newly infected with brucellosis [7].What's more, there are about 300 to 500 million new cases because of malaria annual, among which maybe one million to three million deaths, most of them are children [8].Malaria kills a child every 40 seconds, costing the world more than 2000 young lives a day [8] [9].Many scholars have studied infectious disease models.At the beginning Kermack and Mckendrick opened the door of infectious disease research.In 1927 and 1932, they separately established the classical SIR [10] and SIS [11] model of plague.In 2002 van Driessche and Watmough [12] provided a method of establishing Lyapunov function which furnishes strong theoretical support for the latter study of the dynamic process of the infections.On this basis, more and more scholars established different models for different infectious diseases [13] [14] [15].However, a basic model on zoonosis has not been established at now.Many models consider the vaccination [12] [16], but there are only a few people get vaccinated in poor area, even there is not an effective vaccination for some newly discovered zoonoses.Some models consider human transmission to human or animal [17] [18], in fact there are some zoonosis that can only be transmitted from animals to humans such as rabies, brucellosis, tapeworm disease and so on.What's more, few animals are treated after they get sick in the real life.This paper studies on this kind of zoonosis.
In this article, a general model of the interaction between infected animals and humans are established.Referencing to classical model [10] [11] we generalize a SIS model which consider the cure rate those who are infected on human population and a SI model on animal population.The organizational structure of the article is shown below.In Section 2 a mathematical expression of the model is constructed.The third section studies the global dynamics of disease from the perspective of disease-free equilibrium and endemic equilibrium.Section 4 makes numerical simulation to explain the consistency of the theory and numerical results, after that the sensitive analysis is performed.Finally, the conclusion is shown in Section 5.

Model Formulation
We consider both animal and human, then classify each of them into tow subclasses: susceptible and infectious.Defining 1 1 , S I as the susceptible and infec- tious human, 2 2 , S I as the susceptible and infectious animals respectively.If a susceptible human individual is contact with infected animals, this human individual may be infected.If this one received effective treatment, he can change back to a susceptible individual, but if he isn't cure in time, he is in the risk of death from infection.The model proposed in this paper is a system composed of four ordinary differential equations: The dynamic transmissions for the model of infectious of animal and humans are demonstrated in Figure 1.factor with probability of 1 2 , α α .All the parameters with subscript 1 define the people population property while the parameters with subscript 2 define the animal population property.Based on the reality, we can know there is not a parameter is negative for human population.

Dynamic Analysis
Following, we express the derivative of t with 1 1 2 2 , , , S I S I ′ ′ ′ ′ .Adding the first two equations, we can get There is Accordingly, we can conclude an invariant set for the model

The Basic Reproduction Number
Obviously, setting the right-hand side of the four differential equations to zero we can find the equilibrium.When set The model ( 1) has a unique disease-free equilibrium 0 E : ( ) when 1 2 0 I I = ≠ , we can find a unique equilibrium * E : ( ) , , , where Based on this, the knowledge of next generation matrix and basic reproduction number are proposed [10], we now consider the follow auxiliary model It can be written as ) Following the recipe from [10] we have FV − is the next generation matrix and 0 R is the spectral of

Dynamic Analysis on Disease-Free Equilibrium
Theorem 1.If 0 1 R < , the disease-free equilibrium 0 E of the model is locally asymptotic stability, or else, 0 E is instability. Proof: The disease-free equilibrium Jacobian matrix of model ( 1) is ( ) The eigenvalues of  can be obtained from the above equation as follows: ( ) when 0 1 R < , we can know 4 0 λ < , then every eigenvalue of the model has the negative real parts, 0 E is locally asymptotic stability.Otherwise, 0 E is unstable and there will be exist diseases case.Theorem 2. If 0 1 R < , then the disease-free equilibrium 0 E is globally asymptotic stability.

Proof:
The Jacobian matrix of model ( 2) on the disease-free equilibrium 0 E is ( ) It is easy to see that when 0 1 R < , there is Since the invariant set of the model ( 1) is Ω , choose positive and small enough 1 ε , 2 ε , and 0 i t > , such for all i t t > , there is Consider the following auxiliary model ( ) ( ) Similarly, the Jacobian matrix of model ( 3) on the disease-free equilibrium 0 E is ( ) for positive and small enough 1 2 , ε ε .By the comparison principle [19], it can be concluded that ( ) . By the theory of [20], which defined that ( ) is called asymptotically autonomous-with limit equation ( ) Thus, there exist ( )

S t S →
as t → ∞ .Accordingly, when 0 1 R < , the disease-free equilibrium of the model 0 E is global asymptot-Advances in Pure Mathematics ic stability.

Dynamic Analysis on Endemic Equilibrium
We have calculated the endemic equilibrium * E in the above article.Next, the global exponentially stable of * E will be shown.Theorem 3. If 0 1 R > , then the endemic equilibrium * E of the model is globally asymptotic stability.

Proof:
On the endemic equilibrium * E , we derive the model ( 1) satisfied the follow equation: Thus, the (1) can be converted into It is easy to get The corresponding differential equation is We construct the Lyapunov function as ( ) Therefore, the derivative of L is as follow It is easy to see, 0 L′ ≤ for all 3 X .In addition, if and only if 3 0 Thus, the endemic equilibrium * E of the model is globally asymptotic stability in the interior of Ω .

Numerical Simulation
In this section, we discuss two situations in which animal population size is larger or smaller than population size.For the purpose of proving the universality of the model, we analyze the global stability by establishing the simulation data including animals and setting the parameters and initial values of the population size.After that, we give the sensitive analysis to see what is the most important factor to infect the existence or extinction of disease.

Numerical Results
Because it is not specific to infectious diseases, we use dimensionless parameter to conduct the numerical experiments.All the parameters values and initial values are displayed in Table 1, when people size is larger than animal size.
Since the transmission rate of animal to animal is bigger than animal to human and disease mortality rate is bigger than nature mortality rate, for 0 1 R < , we set We can conclude from Figure 2(a) that when 0 1 R < , the number of in- fected animal individuals is stable at 0, that is to say the disease-free equilibrium 0 E is globally stable.Otherwise, from Figure 2(b) when 0 1 R > , the result shows that the quantity of infected animal individuals is finally stable around 10.9, in the meantime, we calculate * 2 10.86957I = , it illustrates that the endemic equilibrium * E is globally stable when 0 1 R > .
Next, we analysis the situation of the people size is smaller than animal size.
Values of parameters and initial values are shown in Table 2.
For 0 1 R < , we set   Similarly, Figure 3(a) we know when 0 1 R < , the disease-free equilibrium 0 E is globally stable.When 0 1 R > , the number of infected animal individu- als is stable around 545.

Sensitive Analysis
With the view of investigating the role of each parameter on 0 R , we use the fol- lowing formula: ( ) where ρ is the parameters of 0 R , and ( ) SI ρ shows the sensitivity of 0 R to ρ .The greater the absolute value of ( ) SI ρ , the greater the effect of ρ on 0

Conclusion
In this article, we construct a general model of zoonoses transmitted from animals to humans.By dynamic analysis we can know that when 0 1 R < , the dis- ease-free equilibrium 0 E is globally stable.When 0 1 R > , * E is globally stable.Theory tells us that infectious diseases will not unrestrictedly grow, but if not take action to solve the epidemic of infectious diseases, it will exist all the time and then continuously make harm to human society.We can take measures to control the birth rate of animals to reduce the incidence of infectious diseases.

Figure 1
Figure 1 illustrate that the people and animal population import new susceptible individuals by the birth of baby with the number of 1 b and 2 b .Under the influence of infected animals, the susceptible individuals will transfer to infected one with the rate of 1 β and 2 β .In the people population, the infected individual can return back to a susceptible individual with the rate 1 σ by med- ical treatment.The infected individuals will die and then be removed from the population for the reason of nature factor with probability of 1 2 , d d and disease

Figure 1 .
Figure 1.Transmission diagram of zoonosis among animals and humans.

Figure 2 .
Figure 2. Numerical results of the infected animal individuals when people size is larger than animal size.

Figure 3 .
Figure 3. Numerical results of the infected animal individuals when people size is smaller than animal size.
follows that increase the mortal- ity of animal can promoting the extinction of disease.However, massive killing of animals is not advocated for those animals like dogs, sheep.So we can find ways to reduce the infection rate and birth rate.

Table 1 .
Values of parameters and initial values when people size is larger than animal size.

Table 2 .
Values of parameters and initial values when people size is smaller than animal size.