^{1}

^{*}

^{1}

In the current work, we study two infectious disease models and we use nonlinear optimization and optimal control theory which helps to find strategies towards transmission control and to forecast the international spread of the infectious diseases. The relationship between epidemiology, mathematical modeling and computational tools lets us to build and test theories on the development and fighting with a disease. This study is motivated by the study of epidemiological models applied to infectious diseases in an optimal control perspective. We use the numerical methods to display the solutions of the optimal control problems to find the effect of vaccination on these models. Finally, global sensitivity analysis LHS Monte Carlo method using Partial Rank Correlation Coefficient (PRCC) has been performed to investigate the key parameters in model equations. This present work will advance the understanding about the spread of infectious diseases and lead to novel conceptual understanding for spread of them.

Recently, due to the fast spread of pandemic diseases, mathematical modeling in the field of epidemiology has attracted many scientists in different areas. Many mathematical models have been developed to describe the transmission of communicable diseases [

One of the new approaches in modeling dynamic systems is the theory of optimal control. For the first time, R. E. Bellman introduced a new method to solve dynamic systems by using the principle of optimality which reduces significantly the computation of the optimal controls [

To model a dynamic system, we usually use a set of ordinary differential equations. A system of ODEs for t 0 ≤ t ≤ t f can be described by

y ˙ = [ y ˙ 1 y ˙ 2 ⋮ y ˙ n ] = [ f 1 ( y 1 ( t ) , ⋯ , y n ( t ) , t ) f 2 ( y 1 ( t ) , ⋯ , y n ( t ) , t ) ⋮ f n ( y 1 ( t ) , ⋯ , y n ( t ) , t ) ]

Based on how the conditions at the endpoints of the domain are specified, we classify an ODE solving problem into initial value problems (IVP) and boundary value problems (BVP). For an initial-value problem, all the conditions are specified at the initial point. For a boundary-value problem, the conditions are needed for both initial and final points. There are many numerical methods to solve initial value problems, such as Euler, Runge-Kutta or adaptive methods and boundary value problems, such as shooting methods [

Euler method is the most common used single-step method. In this discretization technique, for differential equation x ˙ = f ( x ( t ) , t ) , we can make a convenient approximation of this:

x n + 1 ≃ x n + h f ( x ( t n ) , t n ) :

The approximation x n + 1 of x ( t ) at the point t n + 1 has an error of order h^{2}. There exists a trade-off between accuracy and complexity of calculation which depends heavily on the chosen value for h. As h is decreasing, the calculation would be longer however more exact. One of the disadvantages of this method is for many higher order systems. It is very difficult to have an effective Euler approximation. Thus, we need to use more accurate and elaborate methods and one of them is the Runge-Kutta method.

Runge-Kutta method is a multiple-step method. In this technique, we obtain the solution at time t k + 1 from the values t j − k , ⋯ , t k and j is the number of steps. To approximate a differential equation of the form x ˙ = f ( x ( t ) , t ) , we can use the second order Runge-Kutta method

x n + 1 ≃ x n + h 2 [ f ( x n ( t ) , t n ) + f ( x n + 1 ( t ) , t n + 1 ) ] ;

or the fourth order Runge-Kutta method

x n + 1 ≃ x n + h 6 ( k 1 + 2 k 2 + 2 k 3 + k 4 )

where

k 1 = f ( x ( t ) , t )

k 2 = f ( x ( t ) + h 2 k 1 , t + h 2 )

k 3 = f ( x ( t ) + h 2 k 2 , t + h 2 )

k 4 = f ( x ( t ) + h k 3 , t + h )

For the second and fourth order Runge-Kutta method, the approximation x n + 1 of x ( t ) at the point t n + 1 has an error of order h 3 and h 5 .

In this research, we study the most basic epidemiological models S-I-R model (composed by Susceptible-Infected-Recovered) and S-E-I-R model (Susceptible-Exposed-Infected-Recovered). For these models, we develop some analytical results that are useful in understanding of simple epidemic diseases. We continue this study by proposing the equivalent optimal control problems of the mentioned epidemic models and we numerically solve them using the backward-forward sweep method with fourth order Runge-Kutta. Finally, we perform global sensitivity analysis by LHS Monte Carlo method using PRCC to identify the key parameters that contribute most significantly to the spread or control of the infectious diseases.

Recently, due to the fast spread of pandemic diseases, mathematical modeling in the field of epidemiology has attracted many scientists in different areas. Many mathematical models have been developed to describe the transmission of communicable diseases and among these models, the classical Kermack-McKendrick SIR epidemic model builds the basic skeleton of all of them [

S := Susceptible (People who could potentially catch the disease)

I := Infective (People who currently have the disease)

R := Removed (People recovered or have died)

Assumptions:

1) Total population remains constant

2) Rate of increase in the infectives is proportional to the contact between susceptible and infective

3) Removal rate (death rate is constant)

Using these assumptions, the classical S-I-R model has the following form:

{ d S d t = − β I S + δ R d I d t = β I S − γ I d R d t = γ I − δ R (2.1)

where, β demonstrates rate of infection, γ implies to rate of recovery and δ represents rate of immunity loss. If δ = 0 , we assume a model without immunity loss. In the first equation of system (2.1), susceptible S decreases according to the number of contacts between infective I and susceptible S. Therefore, because of decreasing the rate of change of susceptible over time, in the first equation we get − β I S . The rate of change of infective I increases by IS and decreases by γ I . The term β I S has been added to the second equation of system (2.1) which is due to the increasing the contact between S and I. The negativity of γ I is showing decreasing the rate of change in infective I by moving to the next stage which is recovered or died. The term γ I has been added t the third equation which means that the rate of changing the recovered R is increasing by this factor. The time-evolution of system (2.1) over 300 days have been demonstrated in Figures 1-4.

Will the Disease Spread? What Is the Max Number of Infectives I_{max}? How Many People Catch the Disease?

To answer these questions consider the following general S-I-R model:

{ d S d t = − β I S d I d t = β I S − γ I d R d t = γ I (2.2)

At the start of outbreak we have S = S 0 , I = I 0 and R = 0 . Total population size remains constant during epidemic; therefore, the rate of change of S + I + R must be zero:

d d t ( S + I + R ) = 0 , S + I + R = S 0 + I 0 (2.3)

To find out if the disease will spread, we need to check that

d I d t = I ( β S − γ ) < 0 , S ≤ S 0

Therefore, if S 0 > γ β = 1 q , then disease will spread. Here, 1 q is the contact ration

which is the fraction of population that comes to contact with individual during the period of infectious. However, if the reproductive number or the ratio number

R 0 = β S 0 γ > 1 , we have epidemic. This ratio represents the number of secondary

infection in the population caused by initial primary infection, i.e. how many other people get the disease.

To find the maximum number of infectives or I max , we combine d S d t and d I d t :

d I d S = β I S − γ I − β I S = − 1 + γ β S = − 1 + 1 q S

Assuming

I + S − 1 q ln S = I 0 + S 0 − 1 q ln S 0 (2.4)

Then

I max = I 0 + S 0 − 1 q ( 1 + ln ( q S 0 ) )

Here, I max represents the maximum number of people who have the disease at a given time. For COVID-19, or similar worldwide diseases the value for q (contact parameter) is high since disease easily transmits. When q is large, it means that the number of people get infected is a lot.

To reduce the reproduction rate, one can reduce the number of susceptible, S 0 . One way to decrease the number of susceptible is using vaccination which is a common method to eradicate of infectious diseases. Vaccination can go further than being used for just individuals, but it can be beneficial in large scale communities by preserving the effective reproduction rate below the level which would allow an epidemic to spread. However, an epidemic can start and spread very quickly if the reproduction rate rises beyond the critical value for an epidemic [

To find out how many people catch the disease, based on the first assumption, the total population is constant and to end the disease, the number of infected need to go down to zero (end of out break):

S + I + R = S 0 + I 0

and

R ( e n d ) = − S ( e n d ) + I 0 + S 0

Here, S ( e n d ) is unknown. From (2.4), we have

S ( e n d ) − 1 q ln ( S ( e n d ) ) = I 0 + S 0 − 1 q ln S 0

The graph of S ( e n d ) is decreasing and shows at small value of S ( e n d ) and larger q, we have larger value for R ( e n d ) .

An SIR model is an epidemiological model that represents the number of people infected with a contagious illness in a closed population over time. In another word, there are some other important infections which include a significant latency or incubation period during which individuals have been infected but are not yet infectious themselves (for example this latency period is zero for cold). During incubation period the individual is exposed. See

Here, we write the total population as N = S + E + I + R . So, the S-E-I-R model has the form

{ d S d t = Λ − β I S − δ S d E d t = β I S − ε E − δ E d I d t = ε E − γ I − δ I d R d t = γ I − δ R (3.1)

where

S := Susceptible (People who could potentially catch the disease)

E := Exposed (People who are infected but are not yet infectious)

I := Infective (People who currently have the disease)

R := Removed (People recovered or have died)

δ := Constant death rate

Λ = μ × N := Constant influx of new susceptible ( μ Constant birth rate)

ε := Latency transfer rate to infectious

γ := Recovery rate of infectious

β I S := The bilinear (mass action) incidence.

For simplicity, we assume that the death rates are equal δ S = δ E = δ I = δ R .

If we have S > 0 and E = I = R = 0 , we have a disease free population or disease free equilibrium, which means that there is no disease. To find disease or endemic equilibrium point, we look for a feasible region Σ such that:

d N d t > 0 → d ( S + E + I + R ) d t = d S d t + d E d t + d I d t + d R d t > 0

Therefore, from (5) we have

( Λ − β I S − δ S ) + ( β I S − ε E − δ E ) + ( ε E − γ I − δ I ) + ( γ I − δ R ) = Λ − δ S − δ E − δ I − δ R = Λ − δ ( S + E + I + R ) = Λ − δ N > 0

Thus

d N d t = Λ − δ N ≥ 0 ⇒ Λ ≥ δ N ⇔ Λ δ ≥ N = S + E + I + R

Therefore, the feasible region Σ would be:

Σ = { ( S , E , I , R ) ∈ ℝ 4 | S + E + I + R ≤ Λ δ }

From equation three we have,

ε = ( δ + γ ) I ⇒ I * = ( ε δ + γ ) E

From equation four,

γ I − δ R = 0 ⇒ R * = γ δ ( ε δ + γ ) E

1) Case 1: If E * = 0 No Exposed. So, I * = 0 and R * = 0 .

From equation one,

0 = Λ − δ S * ⇒ S * = Λ δ

Therefore, the diseases free equilibrium would be:

P 0 = ( Λ δ , 0 , 0 , 0 )

2) Case 2: If E * ≠ 0 , then I * ≠ 0 and R * ≠ 0 .

If we add equation two and three, we get

( β I * S * − ε E * − δ E * ) + ( ε E * − γ I * − δ I * ) = 0 + 0 = 0

⇒ β I * S * − ( δ + γ ) I * − δ E * = 0

⇒ β I * S * − ( ε + δ ) E * = 0

⇒ β ( ε δ + γ ) E * S * − ( ε + δ ) E * = 0

⇒ [ β ( ε δ + γ ) S * − ( ε + δ ) ] E * = 0

⇒ ( ε β δ + γ ) S * − ( ε + δ ) = 0

⇒ S * = ( ε + δ ) ( δ + γ ) ε β

General replication number R 0 is the number of new cases any single infected individual is going to create and produce or infect susceptible. To find R 0 , at

equilibrium we have S * ≥ Λ δ . For S * = ( ε + δ ) ( δ + γ ) ε β , we have:

( ε + δ ) ( δ + γ ) ε β ≥ Λ δ ⇔ 1 ≥ Λ ε β δ ( ε + δ ) ( δ + γ ) = R 0

1) R 0 ≤ 1 ⇒ Disease Free equilibrium P 0 = ( Λ δ , 0 , 0 , 0 ) ; we can control and there is no disease.

2) R 0 > 1 ⇒ Endemic equilibrium P * = ( S * , E * , I * , R * ) > 0 .

The time-evolution of system (5) over 300 days have been demonstrated in

A general optimal control (OC) problem needs a cost functional ( J [ x ( t ) , u ( t ) ] ) , a set of state variables ( x ( t ) ∈ X ) , a set of control variables ( u ( t ) ∈ U ) in a time t, with t 0 ≤ t ≤ t f . The main goal is finding a piecewise continuous control u ( t ) and the associated state variable x ( t ) to maximize a given objective functional.

Definition 4.1 (Basic OC Problem in Lagrange formulation). An OC problem is in the form

max u = J [ x ( t ) , u ( t ) ] = ∫ t 0 t f f ( t , x ( t ) , u ( t ) ) d t

s .t x ˙ ( t ) = g ( t , x ( t ) , u ( t ) )

x ( t 0 ) = x 0

x ( t f ) could be free, which means that the value of x ( t f ) is unrestricted, or could be fixed, i.e., x ( t f ) = x [

We consider f and g to be continuously differentiable functions. We suppose that the control set U is a Lebesgue measurable function. Therefore, as long as the controls will always be piecewise continuous, the associated states will be piecewise differentiable.

We can change the maximization problem to a minimization problem by making the cost functional negative:

min { J } = − max { J }

Definition 4.2 (Bolza formulation). The Bolza formulation of the OC problem can be defined as

max u = J [ x ( t ) , u ( t ) ] = Φ ( t 0 , x ( t 0 ) , t f , x ( t f ) ) + ∫ t 0 t f f ( t , x ( t ) , u ( t ) ) d t

s .t x ˙ ( t ) = g ( t , x ( t ) , u ( t ) )

x ( t 0 ) = x 0

where Φ is a continuously differentiable function [

Definition 4.3 (Mayer formulation). [

max u = J [ x ( t ) , u ( t ) ] = Φ ( t 0 , x ( t 0 ) , t f , x ( t f ) )

s .t x ˙ ( t ) = g ( t , x ( t ) , u ( t ) )

x ( t 0 ) = x 0

Pontryagin proposed the idea of adjoint functions to append the differential equation to the objective functional which was one of the most important results of Mathematics in the 20th century and illustrates the necessary conditions to find the optimal control. Similar to Lagrange multipliers in multivariate calculus, Adjoint functions append constraints to the function of several variables to be maximized or minimized [

Definition 4.4 (Hamiltonian). Consider the OC problem in definition (4.1). The function

H ( t , x ( t ) , u ( t ) , λ ( t ) ) = f ( t , x ( t ) , u ( t ) ) + λ ( t ) g ( t , x ( t ) , u ( t ) )

is called Hamiltonian function and λ is the adjoint variable.

Theorem 4.5 (Pontryagin’s Maximum Principle). [

H ( t , x * ( t ) , u ( t ) , λ ( t ) ) ≤ H ( t , x * ( t ) , u * ( t ) , λ ( t ) )

for all controls u at each time t, where H is the Hamiltonian previously defined and

λ ′ ( t ) = ∂ H ( t , x * ( t ) , u * ( t ) , λ ( t ) ) ∂ x

λ ( t f ) = 0

The last condition, λ ( t f ) = 0 , called transversality condition, is only used when the OC problem does not have terminal value in the state variable, i.e., x ( t f ) is free.

This Pontryagin’s Maximum Principle converts the problem of finding a control which maximizes the objective functional subject to the state ODE and initial condition into the problem of optimizing the Hamiltonian pointwise. Therefore, with this adjoint equation and Hamiltonian, we have

∂ H ∂ u = 0

at u * for each t, meaning that the Hamiltonian has a critical point and we call this condition as optimality condition. Therefore, to find the necessary conditions, we do not need to calculate the integral in the objective functional and we only use the Hamiltonian.

Theorem 4.6. [

max u = J [ x ( t ) , u ( t ) ] = ∫ t 0 t f f ( t , x ( t ) , u ( t ) ) d t

s .t x ˙ ( t ) = g ( t , x ( t ) , u ( t ) )

x ( t 0 ) = x 0

Suppose that f ( t , x , u ) and g ( t , x , u ) are both continuously differentiable functions in their three arguments and concave in x and u. Suppose u * is a control with associated state x * , and λ a piecewise differentiable function, such that u * , x * and λ together satisfy on t 0 ≤ t ≤ t f :

f u + λ g u = 0 ,

λ ′ = − ( f x + λ g x ) ,

λ ( t f ) = 0,

λ ( t ) ≥ 0.

Then for all controls u, we have J ( u * ) ≥ J ( u )

Based on how the conditions at the endpoints of the domain are specified, we classify an ODE solving problem into initial value problems (IVP) and boundary value problems (BVP). For an initial-value problem, all the conditions are specified at the initial point. For a boundary-value problem the conditions are needed for both initial and final points. There are many numerical methods to solve initial value problems such as Euler, Runge-Kutta or adaptive methods and boundary value problems, such as shooting methods [

Numerical methods for solving OC problems started from the 1950s with the works of Bellman [

In this section, we present an optimal control (OC) problem to study the dynamics of S-I-R model, using a vaccination process (u) as a measure to control the disease. Let x 1 represents the susceptible population, x 2 represents the proportion of population that is infected and x 3 represents the proportion of population that is recovered or dead. The optimal control problem can be defined as:

min u = J [ x ( t ) , u ( t ) ] = ∫ t 0 t f ( x 2 + u 2 ) d t (5.1)

s .t d x 1 d t = − β x 1 x 2 + δ x 3 − u x 1 (5.2)

d x 2 d t = β x 1 x 2 − γ x 2 (5.3)

d x 3 d t = γ x 2 − δ x 3 (5.4)

x ( t 0 ) = ( x 1 ( 0 ) , x 2 ( 0 ) , x 3 ( 0 ) ) (5.5)

with x ( t ) = ( x 1 ( t ) , x 2 ( t ) , x 3 ( t ) ) and λ ( t ) = ( λ 1 ( t ) , λ 2 ( t ) , λ 3 ( t ) ) , with initial conditions x 1 ( 0 ) = 6 × 10 7 , x 2 ( 0 ) = 10 7 , x 3 ( 0 ) = 10 and the parameters β = 5 × 10 − 9 , γ = 0.12 , δ = 1 / 60 .

Let consider the problem (5.1) and constraints (5.2)-(5.4). With x ( t ) = ( x 1 ( t ) , x 2 ( t ) , x 3 ( t ) ) and λ ( t ) = ( λ 1 ( t ) , λ 2 ( t ) , λ 3 ( t ) ) , the Hamiltonian of this problem can be written as

H ( t , x ( t ) , u ( t ) , λ ( t ) ) = A x 2 + u 2 + λ 1 ( − β x 1 x 2 + δ x 3 − u x 1 ) + λ 2 ( β x 1 x 2 − γ x 2 ) + λ 3 ( γ x 2 − δ x 3 )

A is a weight parameter describing the comparative importance of the two terms in the functional. Using the PMP the optimal control problem can be studied with the state variables

x ˙ 1 = − β x 1 x 2 + δ x 3 − u x 1

x ˙ 2 = β x 1 x 2 − γ x 2

x ˙ 3 = γ x 2 − δ x 3

The adjoint variables are:

λ ˙ 1 = λ 1 ( u + β x 2 ) + λ 2 β x 2

λ ˙ 2 = − A + λ 1 β x 1 − λ 2 ( β x 1 − γ ) − λ 3 γ

λ ˙ 3 = λ 3 δ − λ 1 δ

with transversality conditions λ ( t f ) = 0 .

In this section, we present an optimal control (OC) problem to study the dynamics of S-E-I-R model, using a vaccination process (u) as a measure to control the disease. Let x 1 represents the susceptible population, x 2 represents the

proportion of population that is in the incubation period, x 3 represents the proportion of population that is infected and x 4 represents the proportion of population that is recovered or dead. The optimal control problem can be defined as:

min u = J [ x ( t ) , u ( t ) ] = ∫ t 0 t f ( x 3 + u 2 ) d t (6.1)

s .t d x 1 d t = Λ − β x 1 x 3 − δ x 1 − u x 1 (6.2)

d x 2 d t = β x 1 x 3 − ε x 2 − δ x 2 (6.3)

d x 3 d t = ε x 2 − γ x 3 − δ x 3 (6.4)

d x 4 d t = γ x 3 − δ x 4 (6.5)

x ( t 0 ) = ( x 1 ( 0 ) , x 2 ( 0 ) , x 3 ( 0 ) , x 4 ( 0 ) ) (6.6)

with initial conditions x 1 ( 0 ) = 6 × 10 7 , x 2 ( 0 ) = 10 7 , x 3 ( 0 ) = 10 , x 4 ( 0 ) = 1 and the parameters β = 5 × 10 − 9 , γ = 0.12 , δ = 1 / 60 .

Let consider the problem (6.1) and constraints (6.2)-(6.5). With x ( t ) = ( x 1 ( t ) , x 2 ( t ) , x 3 ( t ) , x 4 ( t ) ) and λ ( t ) = ( λ 1 ( t ) , λ 2 ( t ) , λ 3 ( t ) , λ 4 ( t ) ) , the Hamiltonian of this problem can be written as

H ( t , x ( t ) , u ( t ) , λ ( t ) ) = A x 3 + u 2 + λ 1 ( Λ − β x 1 x 3 − δ x 1 − u x 1 ) + λ 2 ( β x 1 x 3 − ε x 2 − δ x 2 ) + λ 3 ( ε x 2 − γ x 3 − δ x 3 ) + λ 4 ( γ x 3 − δ x 4 )

A is a weight parameter describing the comparative importance of the two terms in the functional. Using the PMP the optimal control problem can be studied with the state variables

x ˙ 1 = Λ − β x 1 x 3 − δ x 1 − u x 1

x ˙ 2 = β x 1 x 3 − ε x 2 − δ x 2

x ˙ 3 = ε x 2 − γ x 3 − δ x 3

x ˙ 4 = γ x 3 − δ x 4

The adjoint variables are:

λ ˙ 1 = λ 1 ( u + β x 3 + δ ) − λ 2 β x 3

λ ˙ 2 = λ 2 ( ε + δ ) − λ 3 ε

λ ˙ 3 = − A − λ 4 γ + λ 3 ( γ + δ ) − λ 2 β x 1 + λ 1 β x 1

λ ˙ 4 = λ 4 δ

with transversality conditions λ ( t f ) = 0 .

Global sensitivity analysis allows us to change all parameters simultaneously over the entire parameter interval. This is a way to evaluate the relative effects of each input parameter and also to identify the interactions between parameters to the model output. In global sensitivity analysis we determine that with variation of input parameters in a certain range, which parameters and interactions have the most influential impact on the overall behavior of our model [

There are several types of global sensitivity analyses, such as weighted average of local sensitivity analysis, partial rank correlation coefficient, multi parametric sensitivity analysis, Fourier amplitude sensitivity analysis (FAST) and Sobol’s method, which can be used for systems pharmacology models [

LHS method is a sampling method and requires fewer samples compare to simple random sampling to achieve the same accuracy [

In LHS method, sampling is independent for each parameter and can be done by randomly selecting values from each pdf. We may sample each interval once for each parameter without any replacement. The LHS matrix is consisting of N rows corresponding to the number of simulations or sample size and also it includes k columns corresponding to the number of varied parameters. Then, N model solutions may be simulated, using each combination of parameter values which they represent each row of the LHS matrix [

Here, a parameter sensitivity analysis has being conducted to identify the biological parameters that have the most significant effect on our model system by the LHS Monte Carlo method using PRCC with uniform distributions for the 95 percent confidence intervals. The global sensitivity results with p-values corresponding to S compartment, I compartment and R compartment have been demonstrated in

According to LHS, we simulated the responses of the model for each organ by randomly selecting values for the parameter set from the 95 percent confidence intervals. These analyses were done by developing a LHS/PRCC method with uniform distributions for the 95 percent confidence intervals. We found that some parameters illustrate significant performance in terms of sensitivity of the output to the variations of these parameters in some compartments while they do not have this effect for others. These results have been depicted in

Infectious diseases can be defined as diseases that can be transmitted from human to human, from human to animal, or from animal to animal. The mathematical modeling of infectious disease spread has been studied for many years

and recently it has been widely discussed due to the spread of the COVID-19 pandemic. To build up an appropriate infectious disease dynamic model, we may need to use a system of ordinary differential equations that cover the spread process, spread law, and spread trend of infectious diseases.

In this paper, we considered the S-I-R and S-E-I-R models and for these, we could develop some analytical results which can be useful in studying the simple epidemics. We displayed the evolution of these two compartmental models over time, Susceptible-Infected-Recovered and Susceptible-Exposed-Infected-Recovered for interesting values of parameters. We followed the optimal control perspective to study these models and because of the complexity of the presented optimal control problems, we could no longer solve them analytically and we ended up looking at the numerical solutions. The optimal curves for the states variables and optimal control were obtained and demonstrated for each control problem separately.

An uncertainty analysis can be applied on the epidemiological models to investigate the uncertainty in system output that is generated from uncertainty in parameter inputs. Sensitivity analysis assesses how variations in model outputs can be apportioned, qualitatively or quantitatively, to different inputs. The final objective of this study was to determine the key parameters in spread of infectious diseases using sampling-based method (Partial Rank Correlation Coefficient-PRCC). In this research, we applied LHS/PRCC method with uniform distributions for the 95 percent confidence intervals on the model Equation (2.1) and Equation (3.1). As we have seen, some parameters have positively and some others negatively affected the spread of disease.

The authors declare no conflicts of interest regarding the publication of this paper.

Azizi, T. and Alali, B. (2020) Studying the Impact of Vaccination Strategy and Key Parameters on Infectious Disease Models. Open Journal of Optimization, 9, 86-104. https://doi.org/10.4236/ojop.2020.93007