Vaccination Prioritization Strategies for COVID-19 in Korea: A Mathematical Modeling Approach

Coronavirus disease 2019 (COVID-19) vaccination has recently started worldwide. As the vaccine supply will be limited for a considerable period of time in many countries, it is important to devise the effective vaccination strategies that reduce the number of deaths and incidence of infection. One of the characteristics of COVID-19 is that the symptom, severity, and mortality of the disease differ by age. Thus, when the vaccination supply is limited, age-dependent vaccination priority strategy should be implemented to minimize the incidences and mortalities. In this study, we developed an age-structured model for describing the transmission dynamics of COVID-19, including vaccination. Using the model and actual epidemiological data in Korea, we estimated the infection probability for each age group under different levels of social distancing implemented in Korea and investigated the effective age-dependent vaccination strategies to reduce the confirmed cases and fatalities of COVID-19. We found that, in a lower level of social distancing, vaccination priority for the age groups with the highest transmission rates will reduce the incidence mostly, but, in higher levels of social distancing, prioritizing vaccination for the elderly age group reduces the infection incidences more effectively. To reduce mortalities, vaccination priority for the elderly age group is the best strategy in all scenarios of levels of social distancing. Furthermore, we investigated the effect of vaccine supply and efficacy on the reduction in incidence and mortality.


Introduction
Since the first case reported in Wuhan, China, in December 2019, Coronavirus Disease 2019 (COVID-19) has rapidly spread worldwide. COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The cases of the disease display similar symptoms to those of Middle East Respiratory Syndrome and Severe Acute Respiratory Syndrome, such as fever, dry cough, dyspnea, and diarrhea [1]. The disease has been characterized as a pandemic by the World Health Organization (WHO) on 11 March 2020 [2]. On 17 February 2021, it was reported that more than 108.8 million people were infected with COVID-19, and more than 2.4 million casualties were recorded across 223 countries [3]. To reduce the spread of the disease while COVID-19 vaccines were not available, countries initially used contact tracing [4] and later implemented non-pharmaceutical interventions (NPIs), which include social distancing (SD), wearing masks, epidemiological surveys, and work/school closure, following the guidelines made by their own government and by WHO.
On 16 February 2021, there have been 84,325 cases with 1534 casualties in Korea; 57.1% of these cases were from the densely populated metropolitan area-26,484 cases from Seoul and 21,648 cases from Gyeonggi Province [5]. As contact tracing was not sufficient to prevent the spread of the disease, NPIs have been implemented differently across regions in Korea, and owing to high floating population and population density,  [14] Recent studies have reported that the susceptibility [15,16], infectivity [16], severity [17], and fatality [18] of COVID-19 vary with age. Thus, the number of incidences, severe patients, and deaths will vary depending on how the vaccination priority is determined for each age group. Besides essential workers and patients in critical condition, Centers for Disease Control and Prevention recommends vaccination priority for the elderly people [19].
The Korean government plans to prioritize vaccination by dividing the entire population, excluding frontline essential workers and patients in critical condition, into different age groups [8]. To reflect this plan and simulate the effect of vaccination under various scenarios, we construct a mathematical model for describing the transmission of COVID-19 with the vaccination prioritization of four age groups: 0-19, 20-49, 50-64, and 65 or older. Several studies on mathematical modeling for COVID-19 vaccination have been conducted recently. In References [20,21], the critical vaccination coverage for various hypothetical vaccine efficacy scenarios is estimated in South Africa and Australia, respectively, to control the disease. In References [22][23][24][25], the effect of vaccination in terms of efficacy and coverage is investigated combined with other interventions, such as the rollout speed of COVID-19 vaccine, face mask usage, or SD. Age-structured models are considered in References [26][27][28][29], and optimal vaccine allocation strategies on age groups are identified under various control scenarios. A multi-region epidemic model is constructed in Reference [30], and the optimal control theory is applied to reduce the number of infectious individuals in the targeted domain with an optimal cost. The appropriate price for COVID 19 vaccine is suggested in Reference [31] using a mathematical model.
In this study, we aimed to investigate the effect of vaccination priority strategies for different age groups on infection incidence and mortality under different SD levels implemented in Korea. For this purpose, we first developed an age-structured mathematical model to describe the transmission dynamics of COVID-19 combined with vaccination. Using the age-structured mathematical model, we also compared the effect of age-dependent vaccination priority strategies on various levels of vaccine efficacy and supply.

Epidemiological Data
In this research, we used data of confirmed COVID-19 cases in Seoul and Gyeonggi Province between 1 February 2020 and 14 February 2021 [32,33]. Because the data for Seoul are provided in age groups of 10 years, we combine the data from both Seoul and Gyeonggi Province into 10-year age groups for consistent interpretation, namely 0-9, 10-19, 20-29, 30-39, 40-49, 50-59, 60-69, and 70 or older (70+). Figure 1 shows the epidemic curve of confirmed COVID-19 cases, and Table 2 shows the total number of confirmed cases for each age group in Seoul/Gyeonggi area.  In Seoul/Gyeonggi area, the number of confirmed cases per day remained below 100 before August 2020 but increased to approximately 200 after mid-August. Since mid-November, a large number of confirmed cases has emerged as winter approaches [34], which can be explained by the fact that the contact time in a confined space has increased, and indoor environments provide a stable condition for SARS-CoV-2 with less sunlight [35].
The infection incidences in Seoul/Gyeonggi area were heavily affected by different levels of SD implemented by the Korean government [36][37][38][39]. The definition of SD levels announced by the government varied in the course of our study. Initially, there were three levels of governmental SD, but, on 7 November 2020, the Korean government revised the governmental SD to five levels as shown in Table S1 in Supplementary Section A1 [40]. In addition, the government's decision and policy-making does not strictly follow the definitions, as the SD announced by the government is considered as a guideline. For example, elevation to governmental SD level 2 requires a monitoring duration of confirmed cases for at least 1 week at governmental SD level 1.5 [40], but governmental SD level 1.5 only lasted for 5 days [41]. Moreover, governmental SD level 2.5 was reinforced with the prohibition of gathering of five or more people on 24 December 2020 [42], despite the fact that governmental SD level 3 only prohibits the gathering of 10 or more people. For the above reasons, the definite classification of SD was difficult, so referring to the government policies that actually took effect, we set our own criteria for SD as presented in Table 3. In this research, the contact matrix, which characterizes the contact degree between age groups, is the linear combination of the location-specific matrices of workplace, school, household, and other locations [43], and SD levels directly impact the contact matrix by location and age as shown in Table 3; the details of the contact matrix will be explained in Section 2.2 and Supplementary Section B1. Contact in locations other than workplace, household, and school decreased by 70%. Contact in household increased by 50% for age less than 20 and 10% for age 20 and above [44].
* Governmental SD level 2.5 with reinforcements is assumed to be equivalent to governmental SD level 3 despite having slightly weaker standards [40] (Table S1 in Supplementary Section B1) due to stronger prohibition of gathering [40,42].

Mathematical Model
We developed an age-structured mathematical model to describe the transmission dynamics of COVID-19 with vaccination. In this model, the population is separated into compartments based on their characteristics for each age group i: S i = susceptible, V i = vaccinated, E i = exposed, P i = pre − symptomatically infectious, A i = asymptomatically infectious, I i = asymptomatically infectious, H M i = hospitalized with mild symptoms, H S i = hospitalized with severe symptoms, R i = recovered and D i = dead. The age classes i = 1, 2, . . . , 8 represent individuals aged 0-9, 10-19, 20-29, 30-39, 40-49, 50-59, 60-69, and 70+, respectively. The schematic diagram of the model is shown in Figure 2. The susceptible population can be vaccinated or become infected. Infection can occur for both the susceptible and the vaccinated if pre-symptomatically, asymptomatically, or symptomatically infectious population makes contact. Once considered infected, this population is exposed and becomes pre-symptomatically infectious after the latent period. The pre-symptomatically infectious population does not show any symptom and moves on to either asymptomatically infectious or symptomatically infectious. The asymptomatically infectious population continues to show no symptoms; thus, no isolation is made until they become recovered. On the other hand, the symptomatically infectious population will be hospitalized based on the severity of their symptoms-either mild or severe-once they are confirmed with COVID-19. The hospitalized population is completely prohibited from making contacts with others. Both hospitalized groups recover, except for those who die from a group with severe symptoms. .
where the infection force Λ i for each age group i = 1, 2, . . . , The model parameters in Equation (1) are described in Table 4, and N j is the total contact possible population of age group j, equivalently N j = S j + E j + P j + A j + I j + R j . Note that m ij is an entry of the contact matrix,C M , which is an 8 × 8 matrix estimated from [36,43] reflecting the actual transmission rates between individuals of different ages in the target area. The detailed computation for contact matrix is given in Supplementary Section B1.   In this study, we find the effective reproduction number R t which measures the mean number of secondary cases infected by an infectious individual. Mathematically, it is computed as R t = ρ(G), where ρ is the spectral radius of the next generation matrix G. For model (1), R t is derived as follows: Here, M A is the matrix computed as follows: where S i is the susceptible population of age group i, V i is the vaccinated population of age group i, C M is the 8 × 8 contact matrix, and diag{ } n denotes the diagonal matrix with n diagonal entries. The detailed derivation of the above formula for R t is given in Supplementary Section B2.

Vaccination Strategies
The Korean governmental vaccination plan on 28 January 2021 states that [8], excluding frontline essential workers and patients in critical condition, vaccines will be administered based on age groups, such as 18-49, 50-64, and 65 or older, starting with the elderly, and 70% of the population will be vaccinated in 6 months from its start date [45]. Though the vaccination of people under the age of 18 has not been decided, the government is planning to administer the vaccine depending on the clinical outcomes [8]. To emphasize the effects of vaccine administration scenarios based on age groups, and because of nonexact vaccination administration schedule based on the producers, the overall vaccination efficacy was used in this research, which was calculated as 88% by taking the average of efficacy of different vaccines in Table 1 with doses as weights. The details of the calculation are explained in Supplementary A2. As the demographic in our study considered 10-year age groups, we investigated four different vaccination priority strategies, which initially distributed vaccination to target age groups of 0-19, 20-49, 50-64, and 65 or older (65+) by 70% of the target population, and then vaccination was administered to other age groups proportional to their population until 70% of the total population was vaccinated. We will refer this vaccination priority strategies for different age groups 0-19, 20-49, 50-64, and 65+ by "0-19 first", "20-49 first", "50-65 first", and "65+ first", respectively. Note that our simulation was based on 10-year age groups, but we implemented two vaccination strategies "50-64 first" and "65+ first", which are not separable by 10-year intervals. To resolve the discrepancy in the age intervals, vaccines distributed to age groups 60-64 and 65-69 were based on population ratio [48]. Four of the five vaccines planned to be introduced in Korea require two vaccine doses per person as shown in Table 1. In this study, we considered the date of the last vaccination as the vaccination date. To compare the case with no prioritization, an additional strategy was considered where vaccination was administered to 70% of all age groups with proportion of population from the start, which will be referred to herein as the POP strategy. Figure 3 illustrates the vaccination allocation for age groups in the five vaccination strategies. The slight difference of distribution of vaccines between age groups within each strategy is caused by the population difference of age groups. The general population ratio between age groups can be observed in the distribution of vaccines in the POP strategy where all age groups were administered with equal coverage rate.
Since the SD can significantly affect the transmission dynamics of COVID-19 [36] and it is very likely that vaccination will be administered while appropriate SD policies are implemented in Korea, we investigated the effect of the five vaccination priority strategies on the number of confirmed cases and mortalities under the assumption of the different SD levels. Furthermore, as the Korean government has criteria for changing the intensity of SD according to the number of confirmed cases [33], we reflected these criteria in the simulation to investigate the effect of the age-dependent vaccination priority strategies on the number of confirmed cases and mortalities.

Parameter Estimation
We estimated the infection probability b i by fitting the actual confirmed cases for each age group by using a MATLAB-embedded function, lsqcurvefit, which is a nonlinear solver that finds the coefficient b i that minimizes gap of confirmed cases between the actual data and the simulated results from Equation (1) in the least-squares sense. The detailed explanation of the function is given in Supplementary Section B3. In the estimation, we used the contact matrix for each age group in Korea [36,43]. For SD level 0, 1, 2, and 3, to reflect the impacts of SD levels, we used different contact matrices M 0 , M 1 , M 2 , and M 3 , respectively. As shown by the parameter values in Table 4, the average time from the presymptomatic stage to confirmation is approximately 6 days. Thus, under the assumption that SD affects pre-symptomatic infection immediately, it is estimated that the number of confirmed cases is influenced approximately 6 days after the implementation of new SD. The details on the SD affecting confirmed cases and the infection probability b i are shown in Supplementary Section B4. Figure 4 and Table 5 show the results of the estimation of infection probability from the actual data of confirmed cases in the four distinct periods of SD levels 0, 1, 2, and 3. As SD level increases, the effective reproduction number R t decreases as shown in Table 5. The results of the confirmed cases fitting for all ages for each SD level are shown in Figure S2 in Supplementary Section C. The results of the (cumulative) confirmed cases fitting for each age group in different SD levels are shown in (Figures S7-S10) Figures S3-S6 in Supplementary Section C.

Ethical Considerations and Data Sharing Policy
COVID-19 data for Seoul and Gyeonggi Province are accessible in Reference [32,33], respectively. We used the fully anonymized data available in Reference [32,33], so there is no ethical issue to consider in this work.

Results
In this section, we investigated the effect of the age-dependent vaccination priority strategies under different levels of SD as described in Section 2.3. We assumed that 70% of the total population was vaccinated over 6 months, and the number of vaccinations per day was equal.
As shown in Figure 5, we compared the effects of the vaccination priority strategies on the number of confirmed cases and deaths under different levels of SD. The exact values of the results displayed in Figure 5 are shown in Table 6.  3. Labels "0-19", "20-49", "50-64", "65+" denote the vaccination priority strategies "0-19 first", "20-49 first", "50-64 first", and "65+ first", respectively.  Figure 5 shows that any vaccination strategies significantly reduced the number of incidences and mortalities, compared with no vaccination. Furthermore, Table 6 showed that any vaccination strategies under any SD level can reduce the number of incidences and fatalities by more than 10 times, compared with no SD at all. In particular, the highest SD level reduced the number by more than 4 times, compared with other SD levels. However, if the highest SD level persists for months, economic and social damage must be considered. Under SD level 0 or 1, when SD was minimally or not at all implemented, vaccination priority for the 20-49 age group reduced the cumulative incidence the most, and prioritizing the 65 or older group had the greatest reduction in the number of deaths. This result is consistent with the results of previous studies that did not assume SD [26]. However, under SD level 2 and 3, vaccination priority for 50-64 and 65 or older age groups reduced the cumulative confirmed cases the most, respectively.
We also investigated the effect of vaccination priority strategies on the number of confirmed cases and deaths when the SD level changed adaptively according to the Korean governmental strategies [40]. The SD level changing criteria are presented in Table 3. Figure 6 illustrates the transition of SD level according to the incidence, and Figure 7 shows the reduction rates in confirmed cases and deaths for each vaccination strategies, compared with the case of no vaccination. The exact values of cumulative confirmed cases and deaths of each scenario in Figure 7 are shown in Table 7. The changes of SD level according to confirmed cases and SD level criteria when no vaccination was given is presented in Figure  S11 in Supplementary Section C. Both in Figures 6 and 7, the colored circles for strategies "0-19 first", "20-49 first", "50-64 first", and "65+ first" denote the point when vaccination of the priority target population has reached 70%. Beyond this colored circle, vaccination was administered to other age groups proportional to their population until 70% of the total population was vaccinated.   Figure 6 shows that SD level was higher in the case of no vaccination than in other vaccination cases for most of the time, because, as shown in Figure 7, the number of confirmed cases in the case of no vaccination was much larger than that in other cases, making the SD level higher than that in other cases for most of the time. Figure 6 shows that without vaccination, SD level alternates between 2 and 3, while any vaccination strategy will lower the required SD level to 1. Figure 7 shows that reduction in death is substantially effective with strategy "65+ first", compared with other strategies since vaccination targets the age group with the highest death rate, and any vaccination strategy successfully reduces the cumulative confirmed cases, compared with no vaccination scenario.  Figure 8 shows the effect of vaccine supply on the reduction in cumulative confirmed cases and mortalities for different SD levels. Instead of focusing on the vaccine strategy introduced previously at 70% coverage of the total population, we explored the effects for different vaccination coverages. As vaccine supply varied, we assumed that, in each vaccination priority strategy, vaccine was initially administered to the target age groups by x% of the target population and then administered to the other age groups proportionally to their population until x% of the total population was vaccinated, where x% ranged from 0% to 100%. Similarly, the POP strategy also considered variation in vaccine supply using the same total coverage as that in the other strategies.  Figure 8 shows that the reduction rate for the number of confirmed cases and deaths increased as vaccine supply rate increased in any SD levels. In particular, in the absence of SD, the reduction rate for the number of confirmed cases and deaths increased rapidly as the rate of vaccine supply increased. However, under SD, the reduction rate gradually increased with increasing vaccine supply. For the reduction in deaths, "65+ first" strategy was the most effective for any vaccine supply in all SD levels. However, the most effective strategy for reducing cumulative confirmed cases varied depending on SD level or supply. In SD level 0, when vaccine supply was relatively limited, "50-64 first" was the most effective strategy in reducing cumulative confirmed cases, but as supply increases, "20-49 first" was the best strategy. Strategies "20-49 first", "50-64 first", and "65+ first" were the most effective for any vaccine supply in SD level 1, 2, and 3, respectively, although the difference of reduction rate is minimal, compared with other strategies. Figure 9 shows the effect of vaccine efficacy on the reduction in cumulative confirmed cases and mortalities.  Figure 9 illustrates that an increase in vaccine efficacy resulted in an increase in reduction rate under any SD level. Similar to the results about vaccine supply, as the vaccine efficiency increased, the reduction rate of the number of confirmed cases and deaths increased much more rapidly in the absence of SD than in the case of SD. Moreover, "65+ first" strategy was the most effective for any vaccine efficacy in all SD levels for reducing deaths. The most effective strategy for reducing cumulative confirmed cases varied depending on SD level or vaccine efficacy. At SD level 0, when vaccine efficacy was relatively less effective, "50-64 first" was the most effective strategy in reducing cumulative confirmed cases, but with higher efficacy, "20-49 first" became the best strategy. Strategies "20-49 first","50-64 first", and "65+ first" were the most effective for any vaccine supply in SD level 1, 2, and 3, respectively, although the reduction rate difference is minimal, compared with other strategies. Figures 10 and 11 show the combined effect of vaccine supply and efficacy on the reduction in cumulative confirmed cases and mortalities, respectively. Given the vaccine supply and efficacy, the best vaccination priority strategy for reducing cumulative confirmed cases and deaths was found for each SD level.   Figure 10 shows that, in the case of no SD, "20-49 first" was the best strategy in terms of reducing the number of confirmed cases when both vaccine supply and efficacy were sufficiently large, but "50-64 first" gave the best reduction when vaccine supply and efficacy were relatively small. When SD level 1 was implemented, "20-49 first" resulted in the best reduction. However, in the case of SD level 2 and 3, "50-64 first" and "over 65 first" were the best strategy, respectively.
Regarding mortalities, Figure 11 illustrates that "over 65 first" was the best strategy for reducing the number of mortalities under any levels of SD, except for the case that vaccine supply was relatively low in SD levels 0 or 1.
Additionally, we investigated the effects of vaccine allocation order based on the age groups 0-19, 20-49, 50-64, and 65+. Unlike the previous strategies, in which vaccine was distributed to the rest of the population after covering the target coverage rate of the initial age group, the next age group was selected for vaccination after the first group, where all four groups were vaccinated eventually. We assumed that only one group was vaccinated at a time with a coverage of 70% until the next vaccination group started vaccinating. Since there are four groups, there were 4! = 24 different scenarios to vaccinate all the age groups in order. Like previous results, the effects of vaccine allocation order were examined based on reduction of the number of cumulative confirmed cases and deaths, compared with when vaccine was not given for each SD level. Figure 12 shows the details and orders of the vaccine administration for each scenario. Figure 13 shows the effects of each scenario at different SD levels. The optimal vaccine allocation order was determined for reducing cumulative confirmed cases and deaths for each SD level; the optimal results were colored red in Figure 13.

Discussion
In this study, we developed an age-structured mathematical model for describing the transmission of COVID-19, including vaccination. Using the model, we investigated the effect of vaccination priority strategies for different age groups on the transmission dynamics of COVID-19 under various scenarios of SD in Seoul/Gyeonggi area in Korea.
We estimated the infection probability of each age group under different levels of SD in the focus area by fitting the actual data of the confirmed cases by the least squares method. SD is one of the main NPIs implemented in Korea, and it played an important role in preventing a massive increase in the number of incidences and mortalities [36][37][38][39]. The effect of SD was reflected in terms of the contact matrix in our model, and the change in the number of confirmed patients due to direct/indirect effects of SD was used in the estimation of the infection probability. The change in the intensity of SD influenced the magnitude of the infection force, which was proportional to the infection probability and the entries of the contact matrix in the model Equation (1).
We observed that the reproduction number R t value decreased as the intensity of SD increased (Table 5), which shows the effect of NPIs implemented by the government in the course of our study.
To investigate the effect of vaccination in the target area, we reflected the vaccination plan prepared by the Korean government. The vaccination efficacy was calculated as 88%, based on the actual vaccine supply plan [8] in Korea and the efficacy of each vaccine as shown in Table 1. In the same manner as in the Korean government's vaccination plan, we divided the total population into four age groups, 0-19, 20-49, 50-64, and 65 or older, and investigated the effect of various vaccination priorities for the age groups on the reduction in the number of confirmed cases and mortalities under different levels of SD implemented in Korea. We found that prioritization of COVID-19 vaccination for the eldest age group of 65 or older resulted in the greatest reduction in the total mortalities under any levels of SD. This result is consistent with WHO recommendations to prioritize COVID-19 vaccination for the elder age groups [49]. Recent studies suggested prioritizing vaccine allocation for younger age groups with higher contact rate to effectively reduce the incidence of infection [26], which is consistent with our results for the case that no or a low-level SD was implemented during the vaccination period. However, when a higher level of SD was implemented, our results showed that prioritizing COVID-19 vaccination for the elder age groups was the better strategy to reduce the number of confirmed cases. These results can be explained by the characteristics of the elderly life in Korea: when the intensity of SD is high, the elderly people who have close contacts in confined spaces, such as nursing homes and religious facilities, are more vulnerable to infection.
The Korean government has criteria for changing the intensity of SD according to the size of the number of confirm case [40]. By reflecting these criteria in the simulation, we investigated the effect of age-dependent vaccination priority strategies. The simulation results showed that vaccine priority for 65 or older group minimized the number of mortalities (Figure 7). One interesting thing is that, in the case of "20-49 first", the reduction rate in mortality was lower than even that in the case of no vaccination. This can be explained by the fact that it takes more time to complete the vaccination at 70% for the "20-49" age group than for the other age group, as the former had the highest population ratio among all age groups [48], and thus vaccination for the elder age group with the highest mortality rate was delayed.
We investigated the effect of various vaccine supply rates on the reduction in the number of confirmed cases and mortalities (Figure 8). The results showed that the number of confirmed cases and deaths decreased very rapidly with increasing vaccine supply in the case of no SD, but, as the level of SD increased, the reduction rate in the number of confirmed cases and deaths decreased when vaccine supply increased. In other words, as the SD intensity decreased, the effectiveness of the vaccine supply in reducing the number of confirmed patients and deaths increased. Similar results were obtained qualitatively in the simulation about vaccine efficacy ( Figure 9). Thus, in the case of no or low-level SD, a sufficient supply of vaccines with high efficacy is very important in reducing the number of confirmed cases and deaths.
We investigated the effect on the reduction in cumulative confirmed cases and deaths under different SD levels when vaccine supply and efficacy were considered together. Considering the currently expected vaccine supply and vaccine efficiency in Korea, "over 65 first" was the best strategy for reducing mortalities under all SD levels. Concerning the cumulative confirmed cases, "50-64 first" and "over 65 first" were the best strategy under SD level 2 and 3, respectively, but in SD level 1, prioritizing age group 20-49 gave the best result. In SD level 0, "50-64 first" was the best strategy in reducing the cumulative confirmed cases, except for the case that both of vaccine efficacy and supply are sufficiently large, for which "20-49 first" gave the best result.
Lastly, we examined the effects of vaccine allocation order on age groups 0-19, 20-49, 50-64, and 65+, where all four groups were vaccinated in succession with 70% coverage. The reduction in cumulative confirmed cases and deaths were tested for 24 different scenarios in each SD level. Depending on SD level, the optimal results varied, but for all SD levels, the optimal vaccination order for reducing deaths started with vaccinating the eldest group 65+.
The present study had some limitations. Priority vaccination for essential workers has been planned by the Korean government, but it is not considered in our simulation. Currently, the contact matrix is not available for essential workers related to COVID-19 in Korea; however, the total number of essential workers is only approximately 1% of the population in Korea. As we focused on the effect of vaccination priority strategy on different age groups, essential workers were not included in our model. Moreover, the infection probability in our model is affected by the contact matrix and incidence data, which differ by countries and regions. Thus, if the approach used in our study is applied to other countries or regions, the results would be different from ours.
Nevertheless, despite these limitations, we successfully analyzed the effect of the vaccination priority policy for different age groups on the reduction in the number of confirmed cases and mortalities by using our newly developed mathematical model, which reflected the actual SD policy and vaccination plan implemented in Korea. In addition, we investigated the effective age-dependent vaccination priority strategies to minimize the number of confirmed cases and deaths for various vaccine supplies and vaccine efficiencies. We believe that the modeling approach in this study can be used to investigate the potential effect of age-dependent vaccination priority strategies with various vaccine supplies and efficacies on the reduction in confirmed cases and mortalities in other target areas.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/ijerph18084240/s1, Section A1: Social Distancing (SD) in Korea (Table S1: Description of social distancing determined by Korean government); Section A2: Overall Vaccine Efficacy Calculation; Section B1: Contact Matrix ( Figure S1: Contact matrix for different SD levels; SD level (a) 0, M 0 , (b) 1, M 1 , (c) 2, M 2 , and (d) 3, M 3 ); Section B2: The Effective Reproduction Number R t ; Section B3: lsqcurvefit; Section B4: Infection Probability b i ; Section C: additional result figures ( Figure S2: Estimation of transmission rate: confirmed cases of each SD level; Figure S3: Estimation of transmission rate: confirmed cases of each age group for SD level 0; Figure S4: Estimation of transmission rate: confirmed cases of each age group for SD level 1; Figure S5: Estimation of transmission rate: confirmed cases of each age group for SD level 2; Figure S6: Estimation of transmission rate: confirmed cases of each age group for SD level 3; Figure S7: Estimation of transmission rate: cumulative confirmed cases of each age group for SD level 0; Figure S8: Estimation of transmission rate: cumulative confirmed cases of each age group for SD level 1; Figure S9: Estimation of transmission rate: cumulative confirmed cases of each age group for SD level 2; Figure S10: Estimation of transmission rate: cumulative confirmed cases of each age group for SD level 3); Figure S11. Effects of vaccination priority strategies when SD level changes adaptively according to confirmed cases. Time series of (a) cumulative confirmed cases and (b) deaths.