 Research article
 Open access
 Published:
Priority and age specific vaccination algorithm for the pandemic diseases: a comprehensive parametric prediction model
BMC Medical Informatics and Decision Making volume 22, Article number: 4 (2022)
Abstract
Background
There have been several destructive pandemic diseases in the human history. Since these pandemic diseases spread through humantohuman infection, a number of nonpharmacological policies has been enforced until an effective vaccine has been developed. In addition, even though a vaccine has been developed, due to the challenges in the production and distribution of the vaccine, the authorities have to optimize the vaccination policies based on the priorities. Considering all these facts, a comprehensive but simple parametric model enriched with the pharmacological and nonpharmacological policies has been proposed in this study to analyse and predict the future pandemic casualties.
Method
This paper develops a priority and age specific vaccination policy and modifies the nonpharmacological policies including the curfews, lockdowns, and restrictions. These policies are incorporated with the susceptible, suspicious, infected, hospitalized, intensive care, intubated, recovered, and death submodels. The resulting model is parameterizable by the available data where a recursive least squares algorithm with the inequality constraints optimizes the unknown parameters. The inequality constraints ensure that the structural requirements are satisfied and the parameter weights are distributed proportionally.
Results
The results exhibit a distinctive third peak in the casualties occurring in 40 days and confirm that the intensive care, intubated, and death casualties converge to zero faster than the susceptible, suspicious, and infected casualties with the priority and age specific vaccination policy. The model also estimates that removing the curfews on the weekends and holidays cause more casualties than lifting the restrictions on the people with the chronic diseases and age over 65.
Conclusion
Sophisticated parametric models equipped with the pharmacological and nonpharmacological policies can predict the future pandemic casualties for various cases.
Background
The human history has witnessed a number of devastating pandemics such as the smallpox, cholera, plague, dengue, influenzas, Ebola, severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS), and novel coronavirus diseases 2019 (COVID19) [1]. Since the pandemic diseases mostly spread through humantohuman infection, the nonpharmacological policies including the restrictions, closures, and curfews have been imposed until an effective vaccine has been developed [2]. Even though a vaccine has been developed, problems in its production and distribution create some constraints on fighting the pandemic diseases. In this case, the state authorities seek policies that optimize the respective priorities such as reducing the deaths, easing the curfews, lifting the restrictions, and opening of the schools. This paper proposes a comprehensive parametric model with the priority and age specific vaccination policy which can be used for the prediction and analysis of the future casualties under the constructed policy. In this paper, the healthcare staff constitutes the highest priority group, and the elderly people are located into the risk groups based on their ages.
Accurate models can be a useful tool to understand the dynamics of the pandemic diseases and to identify the role of the internal (mutation) and external (pharmacological and nonpharmacological policies) impacts on the pandemic casualties [3]. Modelling of the pandemic diseases can be achieved by using the nonparametric (statistical and machine learning), and the parametric (mathematical) approaches. Statistical approaches usually reveal a mean and a standard deviation which can be used for characterizing the pandemic properties such as the incubation period and the infectious rate of the pandemic diseases. Overton et al. produced data for the incubation period of the COVID19 by assuming that it has the Gama distribution and fitted the data with the maximum likelihood estimator [4]. This research stated that the majority of the infected people develop symptoms in 14 days. Hong and Li estimated a timedependent reproduction number of disease with the Poisson model having a removal rate to account for the random uncertainties in the reported casualties [5]. It is concluded that China, Italy, Sweden, and the United States of America (USA) have high COVID19 reproduction numbers since they were unable to control the spread of the virus. Oehmke et al. determined the speed, acceleration, jerk, and 7daylag in the COVID19 transmission for the USA and determined the parameters with the ArellanoBond statistical estimator [6]. It is expressed that there were significant differences in the spread of the virus among the states of the USA due to lack of a national nonpharmacological policy.
In terms of the machine learning based modelling approaches, Pinter et al. proposed an adaptive networkbased fuzzy inference systems (ANFIS) and a multilayered perceptronimperialist competitive algorithm (MLPICA) to estimate the infected individuals and the mortality rate [7]. Tuli et al. considered the generalized Inverse Weibull distribution combined with the cloud computing to predict the growth of the epidemic and to design control strategies for the COVID19 spread [8]. Aydin and Yurdakul evaluated the policies of the 142 countries against fighting the COVID19 by using the kmeans clustering, decision trees, and random forest algorithms [9]. The research revealed that the economic welfare, smoking rates, and the diabetes rates are not directly related to the effectiveness level of the countries. Rustem et al. utilized the linear regression (LR), the least absolute shrinkage and selection operator (LASSO), the support vector machines (SVM), and the exponential smoothing (ES) to estimate the threating factors of the COVID19 [10]. The results confirmed that the ES outperforms the others while the SVM performs poorly. Bird et al. evaluated the countrylevel pandemic risks and classified the preparedness of the countries in terms of the transmission, the mortality, and the inability to test by applying the stack of gradient boosting, the decision trees, the stack of SVM, and the extra trees [11]. It is concluded that the geopolitics and the demographic attributes shape the risks caused by the COVID19.
Parametric modelling approaches suit their purpose and are also parameterizable by the available data, since they have a certain model structure representing the mathematical relationships as simple as possible [12]. Goel and Sharma proposed a mobilitybased susceptible, infected, recovered (SIR) model covering the population distribution and the lockdowns [13]. It is observed that the infected casualties are delayed and decreased in the presence of the lockdowns. Piovella provided a simplified analytical solution of the susceptible, exposed, infected, recovered (SEIR) model to predict the casualty peaks and asymptotic cases without iteratively solving the ordinary differential equations [14]. Even though the numeric and analytical solutions are close, there exist biases around the peak values. Piccolomini and Zama proposed a forced susceptible, exposed, infected, recovered, dead (fSEIRD) model with two different piecewise timedependent infection rates [15]. It is stated that the model fits the data and makes reliable predictions for Italy. However, even though the SIR, SEIR, fSEIRD models are simple and require few parameters, they do not consider the pharmacological and nonpharmacological policies which play important roles on the dynamics of the pandemic diseases. In addition, they do not include the hospitalized, intensive care, and intubated pandemic casualties. Lee et al. modelled an optimal age specific vaccination policy against the H1N1 pandemic influenza in Mexico [16]. The model suggested that the optimal vaccination can be achieved by allocating more vaccines for the young adults age between 20 and 39. Recently, we developed a suspicious, infected, recovered (SpID) model with the second order difference equations rather than the first order ordinary differential equations as in the SIR, the SEIR, and the SEIRD models [17]. The results confirmed that the SpID model can represent the higher order properties such as the peak in the COVID19 casualties. In our further research, we proposed a SpIDN model with the nonpharmacological policies (N) including the curfews, restrictions, and lockdowns [2]. The results highlighted the role of each nonpharmacological policy on the COVID19 casualties. In addition, recently we performed a research to analyses the linear and nonlinear dynamics of the COVID19 by only considering the pharmacological policies [18]. This paper developed three model structures from linear to strongly nonlinear and optimized their parameters with the mathematical optimization and machine learning approaches. As an alternative to the modelbased control of the pandemic casualties, an artificial intelligence approach, which is implicitly a model free approach, was constructed to generate the multidimensional nonpharmacological policies [19]. This artificial intelligence algorithm allowed to weight each nonpharmacological policy together with each pandemic casualty under a certain vaccination policy. It firstly aimed to stabilize the pandemic casualties and then minimize them in time. Zhao et al. recently built an agespecific transmission model to quantify the transmissibility in different age groups [20]. Matrajt et al. developed an optimal vaccine allocation algorithm aiming at reducing the deaths, infections, and hospitalizations [21].
This paper proposes a susceptible (S^{c}), suspicious (S^{p}), infected (I^{n}), hospitalized (H), intensive care (I^{t}), intubated (I^{b)}, recovered (R), death (D) with the priority and age specific vaccination (V) and nonpharmacological (N) policies (S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model). The key contributions of the paper can be summarized briefly as

A comprehensive S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model has been constructed by referring the known relationships among the COVID19 casualties illustrated in Fig. 1.

Priority and age specific vaccination policy has been formulated and incorporated into the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model together with the nonpharmacological policies.

Constrained recursive least squares (RLS) optimizer has been modified to learn the unknown parameters of the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model by satisfying the structural and proportional contribution requirements of the design.

An extensive analysis has been performed to assess the role of the priority and age specific vaccination policy and the nonpharmacological policies.
It is important to note that even though this paper mostly refers the COVID19 pandemic, it can be implemented to all the pandemic diseases having the architecture shown in Fig. 1, which is constructed based on the epidemiological facts. In the rest of the paper, the proposed model structures, the proposed S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model, the constrained RLS for the multidimensional models, and the analysis of the model have been provided.
Methods
The proposed model architecture
Individuals in the susceptible \(\left( {S_{k}^{c} } \right)\) group are vulnerable to the pandemic diseases where the suspicious \(\left( {S_{k}^{p} } \right)\) ones leave the group (Fig. 1, number 1) and the noninfected \(\left( {I_{k}^{nn} } \right)\) ones rejoin the susceptible group (number 5). The vaccinated \(\left( {V_{k} } \right)\)(number 2), the recovered \(\left( {R_{k} } \right)\) (number 22), and the death \(\left( {D_{k} } \right)\)(number 21) become nonsusceptible \(\left( {S_{k}^{nc} } \right)\)(number 3) and leave the susceptible \(\left( {S_{k}^{c} } \right)\) group where the remained ones constitute the current susceptible \(\left( {S_{k + 1}^{c} } \right)\)(number 4) group. The individuals in the suspicious \(\left( {S_{k}^{p} } \right)\) group, who are tested and/or quarantined, either move to the infected \(\left( {I_{k}^{n} } \right)\)(number 6) group or the noninfected \(\left( {I_{k}^{nn} } \right)\)(number 7) group where some individuals in the infected \(\left( {I_{k}^{n} } \right)\) group can return the suspicious \(\left( {S_{k}^{p} } \right)\) group again. Also, since the infected \(\left( {I_{k}^{n} } \right)\) individuals spread the virus until they are isolated, they act like as an excitation signal \(\left( {I_{k  1}^{na} } \right)\) (number 8) on the suspicious casualties. Individuals in the infected \(\left( {I_{k}^{n} } \right)\) group can be in the hospitalized \(\left( {H_{k} } \right)\)(number 9) group or in the nonhospitalized \(\left( {H_{k}^{n} } \right)\)(number 10) group where the nonhospitalized \(\left( {H_{k} } \right)\) individuals join the recovered \(\left( {R_{k} } \right)\)(number 15) group after a quarantine period. The individuals in the hospitalized \(\left( {H_{k} } \right)\) group can union with the intensive care \(\left( {I_{k}^{t} } \right)\)(number 11), the intubated \(\left( {I_{k}^{b} } \right)\)(number 12), the death \(\left( {D_{k} } \right)\)(number 13), or the recovered \(\left( {R_{k} } \right)\)(number 14) groups. The individuals in the intensive care \(\left( {I_{k}^{t} } \right)\) group can move to the intubated \(\left( {I_{k}^{b} } \right)\)(number 16), the death \(\left( {D_{k} } \right)\)(number 17), or the hospitalized \(\left( {H_{k} } \right)\)(number 18) groups. Similarly, the individuals in the intubated \(\left( {I_{k}^{b} } \right)\) group can join either the intensive care \(\left( {I_{k}^{t} } \right)\)(number 19), or the death \(\left( {D_{k} } \right)\)(number 20) groups. The nonpharmacological policies \(\left( {u_{k} } \right)\)(number 23) and priority and age specific vaccination policy \(\left( {V_{k}^{*} } \right)\)(number 24) act like an external inhibitor on all the casualties at varying rates.
The ScSpInHItIbRDVN model
This section initially formulates the parametric submodels, and then the vaccination and the nonpharmacological policies of the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model.
The S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN submodels
This subsection constructs the parametric models of each submodel illustrated in Fig. 1.
The susceptible \(S_{k}^{c}\) submodel
Considering the connections coming in and leaving out the susceptible \(S_{k}^{c}\) group in Fig. 1, one can formulate the \(S_{k}^{c}\) submodel with a difference equation. We can initially write the difference equation of the nonsusceptible \(S_{k}^{nc}\) group shown in Fig. 2 as
where \(S_{k}^{nc}\) represents the nonsusceptible individuals who have gained immunity and also the individuals who lost their lives, \(R_{k}\) represents the recovered individuals, \(D_{k}\) represents the dead individuals, \(V_{k}\) represents the vaccinated individuals, \(a_{14}\),\(a_{15}\),\(c_{1}\) are the unknown parameters.
The representation of the susceptible \(S_{k + 1}^{c}\) group in Fig. 2 is
where \(S_{k}^{c}\) represents the individuals who may be infected and have a lack of immunity, \(S_{k}^{p}\) represents the suspicious individuals, \(I_{k}^{nn}\) represents the noninfected individuals, \(a_{11}\),\(a_{12}\),\(a_{13}\) are the unknown parameters.Substituting Eq. (1) in Eq. (2) yields
All the parameters in Eq. (3) are unknown and will be learned from the available data with the RLS algorithm subject to the inequality constraints in the next section.
The next subsection provides the modelling steps of the suspicious \(S_{k}^{p}\) submodel.
The suspicious \(S_{k}^{p}\) submodel
Some of the susceptible \(S_{k}^{c}\) individuals become suspicious \(S_{k}^{p}\) as they exhibit symptoms or contact an infected individual, or return from the regions where the pandemic disease is a threat. These individuals are either tested or quarantined for a time duration. In this paper, we define the suspicious \(S_{k}^{p}\) individuals as the number of the people tested daily. Therefore, the model can predict the number of the required tests in the future. We can represent the \(S_{k}^{p}\) submodel shown in Fig. 3 as
where \(I_{k}^{na}\) represents the individuals who can become suspicious again and excitation effect of the infected individuals on the suspicious casualties (related to filiation time), \(u_{k}^{{}}\) is the nonpharmacological policy, \(V_{k}^{{S^{p} }}\) is the vaccination policy, \(a_{21}\), \(a_{22}\), \(a_{23}\), \(a_{24}\), \(a_{25}\) \(b_{2}\), \(c_{2}\) are the parameters.
The next subsection presents the modelling steps of the infected \(I_{k}^{n}\) submodel.
The infected \(I_{k}^{n}\) submodel
Some of the suspicious \(S_{k}^{p}\) individuals becomes infected \(I_{k}^{n}\) where they either become hospitalized \(H_{k}\) or nonhospitalized \(H_{k}^{n}\), who are quarantined for a period of time, as illustrated in Fig. 1. We can formulate its model by considering the corresponding connections in Fig. 4 as
where \(V_{k}^{{I^{n} }}\) is the vaccination policy, \(a_{31}\), \(a_{32}\), \(a_{33}\), \(a_{34}\), \(a_{35}\) \(b_{3}\), \(c_{3}\) are the parameters.
The next subsection introduces the hospitalized \(H_{k}\) submodel.
The hospitalized \(H_{k}\) submodel
Some of the infected \(I_{k}^{n}\) individuals requiring standard treatments join the hospitalized \(H_{k}\) group. The hospitalized \(H_{k}\) individuals can join the intensive care \(I_{k}^{t}\), the intubated \(I_{k}^{b}\), the recovered \(R_{k}\), or the death \(D_{k}\) groups as shown in Fig. 5. We can formulate the hospitalized model as
where \(V_{k}^{H}\) is the vaccination policy, \(a_{41}\), \(a_{42}\), \(a_{43}\), \(a_{44}\), \(a_{45}\), \(a_{46}\), \(b_{4}\), \(c_{4}\) are the parameters,
The next subsection presents the formulation of the intensive care \(I_{k}^{t}\) submodel.
The intensive care \(I_{k}^{t}\) submodel
Some of the hospitalized \(H_{k}\) individuals move to the intensive care \(I_{k}^{t}\) group where some of them move back to the hospitalized \(H_{k}\) group as shown in Fig. 6. Similarly, some of the intensive care \(I_{k}^{t}\) patients become intubated \(I_{k}^{b}\) where some of them rejoin the intensive care \(I_{k}^{t}\) group, and the rest join the death \(D_{k}\) group. We can construct the intensive care \(I_{k}^{t}\) model as
where \(V_{k}^{{I^{t} }}\) is the vaccination policy, \(a_{51}\), \(a_{52}\), \(a_{53}\), \(b_{5}\), \(c_{5}\) are the parameters.
The next subsection provides the intubated \(I_{k}^{b}\) submodel.
The intubated \(I_{k}^{b}\) submodel
Some of the hospitalized \(H_{k}\) individuals and the intensive care \(I_{k}^{t}\) patients become intubated \(I_{k}^{b}\) as shown in Fig. 7. A number of the intubated \(I_{k}^{b}\) patients move back to the intensive care \(I_{k}^{t}\) unit while the rest join the death \(D_{k}\) group. We can construct the intubated model as
where \(V_{k}^{{I^{b} }}\) is the vaccination policy, \(a_{61}\), \(a_{62}\), \(a_{63}\), \(a_{64}\) \(b_{6}\), \(c_{6}\) are the parameters.
The next subsection formulates the recovered \(R_{k}\) submodel.
The recovered \(R_{k}\) submodel
A number of the hospitalized \(H_{k}\) and the nonhospitalized \(H_{k}^{n}\) individuals join the recovered \(R_{k}\) group who become nonsusceptible \(S_{k}^{nc}\) as illustrated in Fig. 8. We can formulate the recovered \(R_{k}\) submodel as
where \(V_{k}^{R}\) is the vaccination policy, \(a_{71}\), \(a_{72}\), \(a_{73}\), \(a_{74}\) \(b_{7}\), \(c_{7}\) are the parameters.
The next subsection expresses the death submodel.
The death \(D_{k}\) submodel
Some of the hospitalized \(H_{k}\), the intensive care \(I_{k}^{t}\), and the intubated \(I_{k}^{b}\) individuals join the death \(D_{k}\) group and become nonsusceptible \(S_{k}^{nc}\) as illustrated in Fig. 9. We can form the death \(D_{k}\) model as
where \(V_{k}^{D}\) is the vaccination policy, \(a_{81}\), \(a_{82}\), \(a_{83}\), \(a_{84}\), \(a_{85}\) \(b_{8}\), \(c_{8}\) are the parameters.
The next subsection formulates the vaccination policy \(V_{k}^{*}\) and reviews the nonpharmacological \(u_{k}\) policies.
The vaccination \(V_{k}^{*}\) and nonpharmacological \(u_{k}\) policies
This section firstly introduces the priority and age specific vaccination policy \(V_{k}^{*}\) and reviews the nonpharmacological policies \(u_{k}\) that we have developed recently for the first time in the literature [2].
The priority and age specific vaccination policies \(V_{k}^{*}\)
The \(*\) in the priority and age specific vaccination policy \(V_{k}^{*}\) represents the \(S^{c}\),\(S^{p}\),\(I^{n}\),\(H\),\(I^{t}\),\(R\), and \(D\) in the submodels given by Eqs. from (4) to (10). The priority and age specific vaccination policy basis \(V_{k}^{b}\) is defined in terms of the number of the daily vaccinated people in each group as
where \(H_{k}^{s}\) is the healthcare staff, \(A_{k}^{80 + }\) is the people age 80 and over, \(A_{k}^{65  79}\) is the people age between 65 and 79, \(A_{k}^{50  64}\) is the people age between 50 and 64, \(A_{k}^{25  49}\) is the people age between 25 and 49, \(A_{k}^{15  24}\) is the people age between 15 and 24.
Since the people age under 15 are not considered for the vaccination, they are not included in the basis \(V_{k}^{b}\). The corresponding weight parameter vector \(w_{k}^{*}\) scales the contribution of the vaccination policy for each submodel. For example, the weight parameter vector for the hospitalized \(w_{k}^{H}\) is
where the parameters of the \(w_{k}^{H}\) are \(w_{k}^{s}\) is the percentage of the hospitalized \(H_{k}^{s}\), \(w_{k}^{80 + }\) is the percentage of the hospitalized \(A_{k}^{80 + }\), \(w_{k}^{65  79}\) is the percentage of the hospitalized \(A_{k}^{65  79}\), \(w_{k}^{50  64}\) is the percentage of the hospitalized \(A_{k}^{50  64}\), \(w_{k}^{25  49}\) is the percentage of the hospitalized \(A_{k}^{25  49}\), \(w_{k}^{15  24}\) is the percentage of the hospitalized \(A_{k}^{15  24}\).
Now we can formulate the priority and age specific vaccination policy for the hospitalized \(V_{k}^{H}\) in Eq. (6) as
Similarly, we can construct the priority and age specific vaccination policy \(V_{k}^{*}\) for the other submodels by following the same steps introduced in this section. The next subsection provides the revised nonpharmacological policies \(u_{k}\).
The nonpharmacological policies \(u_{k}\)
The authorities impose various curfews and restrictions to confine the spread of the virus. The most common ones are the curfews on the people age over 65, age under 20, and people with the chronic diseases which have been parametrized in [2] (since there is no available data) as
where \(u_{k}^{s}\) is the response of the curfew (in closed form solution), \(n^{s}\) is the number of the people under the curfew, \(k\) is the number of the days and \(k_{i}\) is the start day of the curfew, \(\alpha\) is the discount factor of the response, where \(\alpha^{k} \approx 0\) for \(\alpha = 0.71\) and \(k = 14\) (quarantine duration), \(\sigma_{k}^{s}\) is the random nonparametric uncertainty in the response.
The other common precaution is the curfews on the weekends and holidays, which has a transient ascent part as
where \(u_{i,k}^{wh}\) is the response of the curfews on the weekends and holidays, \(n^{wh}\) is the number of the people under the curfews on the weekends and holidays, \(\sigma_{i,k}^{wh}\) is the random uncertainty in the response.
Its transient descent part is modelled as
The overall response \(u_{k}^{wh}\) is
In terms of the closure of the schools and universities, it is not a curfew as it only hinders mass gatherings of the students; hence, they can come together in smaller groups. Therefore, the response has a transient ascent part as in Eq. (15) and transient descent part as in Eq. (16). These parts are essentially for removing the negative impacts of the schools being open. Then an uncertain saturated part \(u_{sat}\) represents the small gatherings after the closure of the schools. After the transient ascent and descent parts, the saturated part can be represented as
where \(n^{su}\) is the number of the students, \(\sigma_{k}^{su}\) is the random uncertainty in the response, \(k_{n} = k_{i}^{su} + k_{n}^{su}\) where \(k_{i}^{su}\) is the start day and \(k_{n}^{su}\) is the duration of the closure.
Comparison of the prediction models
One can summarize the main advantages of the constructed S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model over the wellknown models such as the SIR, SEIR models in terms of the solution and analysis as

It has difference equations rather than the differential equations. Therefore, it can be solved iteratively without requiring an ordinary differential equation solver.

It has coupled and linear dynamics instead of the slightly coupled nonlinear dynamics. Thus, the mathematical analysis of the parametric model is straightforward.

Its unknown parameters are learned from the reported data by using the wellknown multidimensional optimization approaches rather than the single dimensional statistical approaches.
The next section forms the RLS approach with the inequality constraints to learn the unknown parameters of the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model.
The constrained RLS algorithm
In this paper, the constrained optimization is considered for two reasons: The first one is that the submodels have certain parameter structures together with the corresponding parameter signs and the second reason is to reflect the contributions of the data having huge magnitude differences (for example, while the susceptible \(S_{k}^{c}\) group covers millions of the individuals, the hospitalized \(H_{k}\) group covers only thousands of them). In this section, initially we will divide the optimization problem in terms of the estimated submodel casualties (outputs) and the real casualties. Then, the RLS algorithm with the inequality constraints are modified to learn the unknown parameters.
The estimated submodels
We can represent the estimated submodels \(\hat{y}_{k}^{*}\) in terms of the known basis vector \(b_{k}^{*}\) and the unknown parameter vector \(w_{k}^{*}\), where the \(*\) is denoted for the \(S^{c}\),\(S^{p}\),\(I^{n}\),\(H\),\(I^{t}\),\(R\), and \(D\) in the submodels given by Eqs. from (3) to (10) as
For example, the basis \(b_{k}^{{S^{c} }}\) of the estimated susceptible \(\hat{y}_{k}^{{S^{c} }}\) submodel is formed with respect to the left hand side of Eq. (3) as
And the corresponding unknown parameter vector \(w_{k}^{{S^{c} }}\) of the estimated susceptible \(\hat{y}_{k}^{{S^{c} }}\) submodel with respect to the right hand side of Eq. (3) is
The other estimated submodels, their bases and parameter vectors are formed by following the same procedures as in Eqs. (19), (20), and (21), respectively. The next subsection introduces the modified RLS algorithm with the inequality constraints to learn the unknown parameter vectors \(w_{k}^{*}\).
Learning the Unknown Parameters with the Constrained RLS
The reported casualties are the outputs of the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN submodels and we call them as the real outputs \(y_{k}^{*}\). For example, the real output of the susceptible submodel is the left hand side of Eq. (3), which is \(S_{k + 1}^{c}\). The objective function is constructed with the 2norm of the instant estimation error defined as
where \(\alpha\) is the inequality constraints which are the lower bound of the parameters. We can construct the Lagrange multipliers used for solving the optimization problems as
Getting partial derivative of \(L\left( {w_{k}^{*} ,\lambda } \right)\) with respect to the \(w_{k}^{*}\) yields
Getting partial derivative of \(L\left( {w_{k}^{*} ,\lambda } \right)\) with respect to the \(\lambda\) gives
Reorganizing Eq. (24) as \(w_{k}^{*}\) is on the left and the rest are on the right, and then substituting it in Eq. (25) yields
The Lagrange multiplier \(\lambda\) from Eq. (26) is obtained as
Then by reinserting Eq. (27) into Eq. (24), the unknown parameter vector \(w_{k}^{*}\) can be attained. The next section extensively analyses the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model.
Results
This section initially presents the parameters of the proposed S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model and then analyses the training and prediction results.
Parameters of the model
Table 1 provides the parameters of the model.
The next subsection compares the real and estimated COVID19 casualties with the constrained RLS algorithm.
Real and estimated casualties
Figure 10 shows the real (reported) and estimated casualties for Turkey.
As can be seen from Fig. 10, the estimated casualties with the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model closely follow the real casualties. The model can track the steep peaks and also the daily variations in the casualties even though the constructed parameter spaces are limited (in machine learning approaches, we randomly manipulate the parameter spaces until we have close estimations). The casualties in Fig. 10 have two distinctive peaks and estimated future casualties in Fig. 13 shows the third peak, which will occur in 40 days. It is clear that the susceptible \(S_{k}^{c}\) casualties have noticeable reduction with the initiation of the vaccination process. It seems that this vaccination process has affected the other casualties since they sharply decrease as well. The decrease in the casualties has also been supported with the nonpharmacological policies \(u_{k}\). Figure 11 shows the mean errors and the corresponding standard deviations in the estimates.
As illustrated by Fig. 11, even though all the mean errors are small, the standard deviations are quite large. This is due to existence of the steep peaks shown in Fig. 10. These peaks occurred in December when there were not any active pharmacological and nonpharmacological policies. This implies that the character of the casualties is largely shaped based on the external impacts such as the pharmacological and nonpharmacological policies. Our recent work highlighted that the pharmacological and nonpharmacological policies have damping impact on the casualties whereas they also have natural frequency determined by the internal and coupling dynamics. The next subsection presents the priority and age specific vaccination policy results.
Priority and age specific vaccination policy
Figure 12 shows the priority and age specification vaccination policy for the hospitalized \(V_{k}^{H} = w_{k}^{{H^{T} }} V_{k}^{b}\) and death \(V_{k}^{D} = w_{k}^{{D^{T} }} V_{k}^{b}\) submodels.
Table 2 provides the background parameters in Fig. 12.
Since there is no reported data for the healthcare staff \(H^{s}\)(‘Staff’ in Table 1), the largest values are assigned for them as they are in the highest risk group. Therefore, the vaccination of the healthcare staff has the largest hospitalized \(w_{k}^{H}\) and death \(w_{k}^{D}\) policy values (Table 2, number 1). People age 80 and over \(\left( {A^{80 + } } \right)\) has 28% death percentage and 2.4% hospitalization percentage. Henceforth, the corresponding \(V_{k}^{D}\) is larger than the \(V_{k}^{H}\)(number 2). Same comments are valid for the people age between 65 and 79 \(\left( {A^{65  79} } \right)\) (number 3) as they share the similar percentages. However, with respect to the people having age 50 and 64, they have close \(w_{k}^{H}\) and \(w_{k}^{B}\) parameters; hence, the \(V_{k}^{H}\) and \(V_{k}^{D}\) policy values are close to each other (number 4). For the people age between 25 and 49 \(\left( {A^{25  49} } \right)\), \(w_{k}^{H}\) is 51.7 and \(w_{k}^{D}\) is 5.5; thus, the corresponding \(V_{k}^{H}\) is larger than the \(V_{k}^{D}\)(number 5). Similar comments can be made for the people age between 15 and 24 (number 6). Lastly, since the healthcare staff \(H^{s}\) and people age 80 and over \(\left( {A^{80 + } } \right)\) have the smallest population among the all age groups, they have the smallest regions on the horizontal axis representing the number of days \(k\).
The estimated future casualties
Figure 13 shows the future casualties estimated by the S^{c}S^{p}I^{n}HI^{t}I^{b}DVN model.
The future estimates are obtained under the assumption that 100.000 peoople are vaccinated daily. It is clear that the number of the susceptible \(S_{k}^{c}\) individuals reduces almost linearly since the vaccinated \(V_{k}\) individuals leave the group. Since the susceptible \(S_{k}^{c}\) and the suspicious \(S_{k}^{p}\) groups are strongly coupled, reduction in the susceptible \(S_{k}^{c}\) group is reflected onto the suspicious \(S_{k}^{p}\) group as well. Figure 13 also clearly shows the third peak in the COVID19 casualties. In addition, it is noticaeble that even though the suspicious \(S_{k}^{p}\), the infected \(I_{k}^{n}\), and the hospitalized \(H_{k}\) converge to zero in 500 days, the intensive care \(I_{k}^{t}\), the intubated \(I_{k}^{b}\), and the death \(D_{k}\) converge to zero around 120 days. This fast convergence is due to priority and age specific vaccination policy which focus on vaccination of the people in the high risk groups.
Analysis of the vaccination policy
Figure 14 shows the number of the daily vaccinations and the corresponding average future casualties.
As can be clearly seen from Fig. 14, all the COVID19 casualties reduce depending on the number of the daily vaccinations. The largest reductions occur in the number of the intensive care \(I_{k}^{t}\), the intubated \(I_{k}^{b}\), and the death \(D_{k}\) when the number of the vaccination rises from 50.000 to 100.000. The further noticeable reduction occurs in the number of the suspicious \(S_{k}^{p}\) and the infected \(I_{k}^{n}\) when the daily vaccination number rises from 100.000 to 200.000.
Analysis of the nonpharmacological policies
Figure 15 shows the role of the nonpharmacological policies on the casualties.
As can be seen from Fig. 15 when all the nonpharmacological policies are in place, all the casualties are small and they increase when the curfews are lifted (blue bar). Removing the restrictions on the people age over 65 and people with the chronic diseases has limited effects as they are in the priority group and most of them have been already vaccinated (orange bar). With respect to the partial opening of the schools, since the majority of the students are not attending the schools, their impact is bounded as well (yellow bar). However, curfews on the weekends, holidays, and nights cover the whole population; henceforth, their roles on the casualties are distinctive (purple bar).
Conclusions
This paper developed a comprehensive parametric S^{p}S^{c}I^{n}HI^{t}I^{b}RDVN model to analyse and estimate the role of the priority and age specific vaccination policy and the nonpharmacological policies. The model has a structure constructed by using the key insights about the pandemic diseases. To satisfy the model structural requirements and avoid dominant effects of the large susceptible and suspicious data, a constrained RLS algorithm has been formed. The results clearly show the importance of the priority and age specific vaccination policy on all the casualties. The future hospitalized, intensive care, intubated, and death casualties converge to zero before the other casualties since they have larger importance. However, the future susceptible, suspicious, infected, and recovered casualties are large due to the people in the lower risk groups are not vaccinated yet. The paper also addresses the relationships among the various daily vaccinations, nonpharmacological polices, and the corresponding COVID19 casualties. The results confirm that the curfews on the weekends and holidays has an overwhelming role on reducing the casualties.
Limitations of the work
Effects of the nonpharmacological policies on each age and chronic diseases group are not weighted. Moreover, vaccine effectiveness for each age group has not been added the model. Besides, climate and environmental effects are not considered.
4.2 Future works
The nonpharmacological policies of the S^{c}S^{p}I^{n}HI^{t}I^{b}RDVN model should be also modified to consider the priority and age specific impacts on each casualty. In addition, effectiveness of the different brand of vaccines on each age group should be considered in the model. Moreover, the model can be expanded by considering the unknown uncertainties. Finally, a toolbox should be constructed and provided for free to help the researchers when applying the proposed model.
Availability of data and materials
The data are available from: https://corona.cbddo.gov.tr/.
Abbreviations
 SARS:

Severe acute respiratory syndrome
 MERS:

Middle East respiratory syndrome
 COVID19:

Novel coronavirus diseases
 ANFIS:

Adaptive networkbased fuzzy inference systems
 MLPICA:

Multilayered perceptronimperialist competitive algorithm
 LR:

Linear regression
 LASSO:

Least absolute shrinkage and selection operator
 SVM:

Support vector machines
 ES:

Exponential smoothing
References
Qiu W, Rutherford S, Mao A, Chu C. The pandemic and its impacts. Health Cult Soc. 2017;9:1–11.
Tutsoy O, Polat A, Colak S, Balikci K. Development of a multidimensional parametric model with nonpharmacological policies for predicting the COVID19 pandemic casualties. IEEE Access. 2020;8:225272–83.
Walters CE, Meslé MMI, Hall IM. Modelling the global spread of diseases: a review of current practice and capability. Epidemics. 2018;25:1–8.
Overton CE, Stage HB, Ahmad S, CurranSebastian J, Dark P, Das R, Fearon E, Felton T, Fyles M, Gent N, Hall I, House T, Lewkowicz H, Pang X, Pellis L, Sawko R, Ustianowski A, Vekaria B, Webb L. Using statistics and mathematical modelling to understand infectious disease outbreaks: COVID19 as an example. Infect Dis Model. 2020;5:409–41.
Hong HG, Li Y. Estimation of timevarying reproduction numbers underlying epidemiological processes: a new statistical tool for the COVID19 pandemic. PLoS ONE. 2020;15(7):e0236464.
Oehmke JF, Moss CB, Singh LN, Oehmke TB, Post LA. Dynamic panel surveillance of COVID19 transmission in the United States to inform health policy: observational statistical study. J Med Internet Res. 2020;22(10):e21955.
Pinter G, Felde I, Mosavi A, Ghamisi P, Gloaguen R. COVID19 pandemic prediction for Hungary; a hybrid machine learning approach. Mathematics. 2020;8(6):890.
Tuli S, Tuli S, Tuli R, Gill SS. Predicting the growth and trend of COVID19 pandemic using machine learning and cloud computing. Internet Things. 2020;11:100222.
Aydin N, Yurdakul G. Assessing countries’ performances against COVID19 via WSIDEA and machine learning algorithms. Appl Soft Comput. 2020;97:106792.
Rustam F, Reshi AA, Mehmood A, Ullah S, On BW, Aslam W, Choi GS. COVID19 future forecasting using supervised machine learning models. IEEE Access. 2020;8:101489–99.
Bird JJ, Barnes CM, Premebida C, Ekárt A, Faria DR. Countrylevel pandemic risk and preparedness classification based on COVID19 data: a machine learning approach. PLoS ONE. 2020;15(10):e0241332.
Keeling MJ, Rohani P. Modeling infectious diseases in humans and animals. Princeton: Princeton University Press; 2011.
Goel R, Sharma R. Mobility based SIR model for pandemics—with case study Of COVID19 (2020). arXiv:2004.13015.
Piovella N. Analytical solution of SEIR model describing the free spread of the COVID19 pandemic. Chaos Solitons Fractals. 2020;140:110243.
Piccolomini EL, Zama F. Monitoring Italian COVID19 spread by a forced SEIRD model. PloS ONE. 2020;15(8):e0237417.
Lee S, Golinski M, Chowell G. Modeling optimal agespecific vaccination strategies against pandemic influenza. Bull Math Biol. 2012;74(4):958–80.
Tutsoy O, Colak S, Polat A, Balikci K. A novel parametric model for the prediction and analysis of the COVID19 casualties. IEEE Access. 2020;8:193898–906.
Tutsoy O, Polat A. Linear and nonlinear dynamics of the epidemics: system identification based parametric prediction models for the pandemic outbreaks. ISA Trans. 2021. https://doi.org/10.1016/j.isatra.2021.08.008.
Tutsoy O. Pharmacological, nonpharmacological policies and mutation: an artificial intelligence based multidimensional policy making algorithm for controlling the casualties of the pandemic diseases. IEEE Trans Pattern Anal Mach Intell. 2021. https://doi.org/10.1109/TPAMI.2021.3127674.
Zhao Z, et al. The optimal vaccination strategy to control COVID19: a modelling study based on the transmission scenario in Wuhan city, China. (2020). https://papers.ssrn.com/sol3/papers.cfm?abstract_id=3719045.
Matrajt L, Eaton J, Leung T. Vaccine optimization for COVID19: who to vaccine first? Sci Adv. 2021;7(6):1–11.
Acknowledgements
The author thanks A. Polat, S. Colak, K. Balikci, N.F. Ozdil, and A. Inan for helpful discussions.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
OT was responsible to develop and analyse the proposed algorithm. MYT collected the data, analysed them and revised the paper. Both authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Tutsoy, O., Tanrikulu, M.Y. Priority and age specific vaccination algorithm for the pandemic diseases: a comprehensive parametric prediction model. BMC Med Inform Decis Mak 22, 4 (2022). https://doi.org/10.1186/s12911021017206
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12911021017206