Loading [MathJax]/jax/output/SVG/jax.js
Research article

Kronecker product bases and their applications in approximation theory

  • The Kronecker product is widely utilized to construct higher-dimensional spaces from lower-dimensional ones, making it an indispensable tool for efficiently analyzing multi-dimensional systems across various fields. This paper investigates the representation of analytic functions within hyper-elliptical regions through infinite series expansions involving sequences of Kronecker product bases of polynomials. Additionally, we examine the growth order and type and Tρ-property of series composed of Kronecker product bases that represent entire functions. We also delve into the convergence properties of Kronecker product bases associated with special functions, including Bessel, Chebyshev, Bernoulli, Euler, and Gontcharoff polynomials. The obtained results extend and enhance the existing findings of such representations in hyper-spherical regions.

    Citation: Mohra Zayed, Gamal Hassan. Kronecker product bases and their applications in approximation theory[J]. Electronic Research Archive, 2025, 33(2): 1070-1092. doi: 10.3934/era.2025048

    Related Papers:

    [1] Saeedreza Tofighi, Farshad Merrikh-Bayat, Farhad Bayat . Designing and tuning MIMO feedforward controllers using iterated LMI restriction. Electronic Research Archive, 2022, 30(7): 2465-2486. doi: 10.3934/era.2022126
    [2] Qianqian Zhang, Mingye Mu, Heyuan Ji, Qiushi Wang, Xingyu Wang . An adaptive type-2 fuzzy sliding mode tracking controller for a robotic manipulator. Electronic Research Archive, 2023, 31(7): 3791-3813. doi: 10.3934/era.2023193
    [3] Wenxuan Li, Suli Liu . Dynamic analysis of a stochastic epidemic model incorporating the double epidemic hypothesis and Crowley-Martin incidence term. Electronic Research Archive, 2023, 31(10): 6134-6159. doi: 10.3934/era.2023312
    [4] Yang Sun, Ming Zhu, Yifei Zhang, Tian Chen . Adaptive incremental backstepping control of stratospheric airships using time-delay estimation. Electronic Research Archive, 2025, 33(5): 2925-2946. doi: 10.3934/era.2025128
    [5] Peng Yu, Shuping Tan, Jin Guo, Yong Song . Data-driven optimal controller design for sub-satellite deployment of tethered satellite system. Electronic Research Archive, 2024, 32(1): 505-522. doi: 10.3934/era.2024025
    [6] Lingling Zhang . Vibration analysis and multi-state feedback control of maglev vehicle-guideway coupling system. Electronic Research Archive, 2022, 30(10): 3887-3901. doi: 10.3934/era.2022198
    [7] Denggui Fan, Yingxin Wang, Jiang Wu, Songan Hou, Qingyun Wang . The preview control of a corticothalamic model with disturbance. Electronic Research Archive, 2024, 32(2): 812-835. doi: 10.3934/era.2024039
    [8] Krishna Vijayaraghavan . Robust-observer design for nonlinear systems with delayed measurements using time-averaged Lyapunov stability criterion. Electronic Research Archive, 2025, 33(6): 3857-3882. doi: 10.3934/era.2025171
    [9] Xuerong Shi, Zuolei Wang, Lizhou Zhuang . Spatiotemporal pattern in a neural network with non-smooth memristor. Electronic Research Archive, 2022, 30(2): 715-731. doi: 10.3934/era.2022038
    [10] Zilu Cao, Lin Du, Honghui Zhang, Lianghui Qu, Luyao Yan, Zichen Deng . Firing activities and magnetic stimulation effects in a Cortico-basal ganglia-thalamus neural network. Electronic Research Archive, 2022, 30(6): 2054-2074. doi: 10.3934/era.2022104
  • The Kronecker product is widely utilized to construct higher-dimensional spaces from lower-dimensional ones, making it an indispensable tool for efficiently analyzing multi-dimensional systems across various fields. This paper investigates the representation of analytic functions within hyper-elliptical regions through infinite series expansions involving sequences of Kronecker product bases of polynomials. Additionally, we examine the growth order and type and Tρ-property of series composed of Kronecker product bases that represent entire functions. We also delve into the convergence properties of Kronecker product bases associated with special functions, including Bessel, Chebyshev, Bernoulli, Euler, and Gontcharoff polynomials. The obtained results extend and enhance the existing findings of such representations in hyper-spherical regions.



    In the growing epidemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), effective control measures are required to reduce the final size of infection and mortality. While public health and social measures (PHSMs), especially controlling social mixing or lockdowns, are effective in reducing the final size by 12%-13% [1,2], they impose a considerable economic and social cost. The specific pharmaceutical intervention of mass vaccination is considered useful, especially if it helps to achieve herd immunity level in the population. Various types of vaccines against SARS-CoV-2 have been created in response to the spread of this virus [3]. SARS-CoV-2 has a large basic reproduction number (i.e., the average number of secondary infections due to a single infector in a totally susceptible population) and a relatively high infection fatality risk (IFR) compared with seasonal influenza. Therefore, in 2021, affected countries needed to urgently provide an effective vaccination program.

    Amid the growing pandemic worldwide, there was a limited number of vaccine stocks when the program had just begun. Therefore, effective vaccine prioritization became a critical public health issue and was widely debated [4,5]. Some studies show that prioritizing vaccines to people in old age [6,7,8], essential workers [7,9], or people with disabilities [9,10] would be crucial to minimize the spread of the viral infections. Other studies imply that demographic location [9] and the rollout speed [11,12,13] are also important factors to expedite the suppression of the epidemic. One study demonstrates that some US states do not prioritize or have no vaccination protocols to incarcerated people and points this out to be a potential risk of future spreading of the virus [14]. We raised an urgent question as to how to prioritize vaccines in order to reduce the cumulative number of coronavirus disease 2019 (COVID-19) cases and deaths by employing two distinct vaccination strategies in Japan more efficiently. One strategy is a severity reduction scheme, which mainly targets older people and people with underlying comorbidities. These hosts have been identified as having particularly high IFRs. Therefore, a reduction in the severity of COVID-19 aimed to protect these vulnerable people by a double dose. The other strategy is a transmission blocking strategy, which targets as many people as possible with only a single dose. The resulting efficacy with a single dose is smaller than that of a double dose. However, if a single dose offers a substantial protection from infection, elevating vaccination coverage would help reduce the reproduction number. Therefore, the herd immunity might be achieved more efficiently than the severity reduction scheme. A previous study showed that mortality and years of life lost were reduced most effectively if adults older than 60 years were prioritized [15]. However, other studies have shown that vaccinating younger essential workers or using transmission blocking strategies is more effective than simply covering older people [4,16].

    In February 2021, Japan was successful in importing vaccine stocks that were barely sufficient to cover all adults older than 65 years with the double-dose scheme. Nevertheless, because of forthcoming epidemic risks, policymakers were confronted with the choice to decide on providing a double dose versus a single dose. This decision can sometimes be controversial, especially in a situation where there is a limited number of vaccines available [17]. Therefore, theoretical support using mathematical models is required.

    In a series of two studies, we aimed to identify the optimal vaccine distribution strategy using a simple mathematical model. This part 1 study used the final size equation of a heterogeneous epidemic model and addressed the static situation in which an epidemic is expected only after completing the vaccination program.

    As of 8 March 2021, the Japanese government announced that 72 million vaccine shots had been imported to Japan, which could technically cover the population aged 65 years or older. As an alternative, although a single-dose program only provides a partial antibody response compared with that with a double dose, it can virtually double the vaccination coverage if all vaccines are widely distributed. To compare the herd immunity level for these different vaccination strategies, we computed the final size of infection (i.e. the total number of people who experiences infection by the end of epidemic) and cumulative risk of death for each strategy. The cumulative number of confirmed COVID-19 cases in Japan as of 20 April 2021 was 530,000, which is approximately 0.4% of the total population, while the final size (z) in a homogenous model was estimated to be 89.3% using the following equation:

    R0=ln(1z)z (1)

    where R0 is the basic reproduction number, which ranges from 1.5 to 3.5 with a representative value of 2.5 in COVID-19 [18]. We assumed that the initial size of infection was negligible and all people who were yet to be vaccinated remained susceptible. The final size equation is known to be extended to multi-host version, and if the type of host corresponds to age group, the age-specific cumulative incidence of age-structured susceptible-infectious-recovered (SIR) model can be conveniently calculated using a single equation in an iterative manner. It should be noted that compartmental extensions to susceptible-exposed-infectious-recovered (SEIR) model or SIR model with realistic infectious period distribution does not alter the final size.

    We categorized the Japanese population into 15 discrete age groups for every 5 years of age, with the oldest age group set at 70 years and older. The population data were retrieved from the census data, which were published by the Statistics Bureau of Japan as of February 2021 (Supplementary Table 1). We used the age-specific IFR that was calculated in a study conducted by Levin et al. [19].

    Table 1.  Contact frequencies between age groups as shown by a contact matrix.
    Age of contactees Age of contactors
    0-4 5-9 10-14 15-19 20-24 25-29 30-34 35-39 40-44 45-49 50-54 55-59 60-64 65-69 70-
    0-4 1.05 0.54 0.21 0.01 0.04 0.20 0.79 0.67 0.39 0.09 0.05 0.09 0.21 0.24 0.12
    5-9 0.56 1.03 0.59 0.21 0.01 0.07 0.24 0.49 0.57 0.36 0.07 0.03 0.08 0.14 0.05
    10-14 0.22 0.59 1.70 0.45 0.02 0.01 0.10 0.17 0.51 0.64 0.24 0.12 0.03 0.02 0.04
    15-19 0.01 0.24 0.50 0.92 0.34 0.08 0.04 0.08 0.30 0.57 0.48 0.35 0.08 0.03 0.06
    20-24 0.06 0.02 0.04 0.45 0.47 0.48 0.24 0.13 0.22 0.33 0.77 0.84 0.35 0.10 0.13
    25-29 0.30 0.10 0.01 0.11 0.48 0.80 0.65 0.23 0.26 0.25 0.48 0.86 0.71 0.37 0.13
    30-34 1.26 0.37 0.15 0.05 0.25 0.69 1.07 0.53 0.31 0.23 0.31 0.46 0.74 0.56 0.20
    35-39 1.17 0.84 0.29 0.12 0.15 0.27 0.58 0.88 0.57 0.33 0.22 0.28 0.56 0.48 0.26
    40-44 0.78 1.13 1.00 0.53 0.29 0.35 0.39 0.66 1.22 0.63 0.54 0.39 0.39 0.46 0.43
    45-49 0.20 0.76 1.34 1.10 0.49 0.37 0.31 0.41 0.68 0.98 0.68 0.59 0.31 0.23 0.41
    50-54 0.08 0.13 0.43 0.78 0.95 0.61 0.36 0.24 0.50 0.58 1.08 0.64 0.41 0.20 0.31
    55-59 0.14 0.04 0.17 0.46 0.86 0.89 0.44 0.24 0.30 0.41 0.53 1.10 0.49 0.22 0.25
    60-64 0.30 0.11 0.03 0.10 0.32 0.66 0.65 0.45 0.27 0.20 0.31 0.45 0.94 0.53 0.25
    65-69 0.41 0.23 0.04 0.05 0.11 0.42 0.59 0.47 0.39 0.18 0.18 0.24 0.65 0.62 0.43
    70- 0.57 0.25 0.21 0.26 0.42 0.41 0.60 0.71 1.02 0.90 0.79 0.77 0.87 1.22 1.46

     | Show Table
    DownLoad: CSV

    We used the age-dependent transmission model, which required an age-dependent contact matrix. Using a method that was described elsewhere [20,21], we quantified the contact matrix during the spring period of COVID-19 in 2000 before PHSM and used the contact matrix shown in Table 1 [22].

    We derived the next-generation matrix K by calculating the element-wise product of the normalized contact matrix (i.e., contact matrix divided by its eigenvalue) and the vaccination coverage matrix S, which was then scaled by the basic reproduction number R0:

    K=R0Cρ(C)S=(R00,00R00,70R70,00R70,70) (2)

    The subscript in the element of the matrix represents the age groups that the infectee and infector belong to (e.g., R25, 40 represents the average number of secondary cases in the age group of 25-29 years produced by a single primary case in the age group of 40-44 years). The matrix of susceptibility (S), including the vaccine effect, was a 15 × 15 matrix, which reflected the relative susceptibility of each age group as determined by vaccination coverage:

    (3)

    where pXs and pXd represent the fraction of people belonging to age group X who received a single and double dose, respectively (i.e., vaccination coverage), and vs and vd denote the relative susceptibility of vaccinated individuals who received a single and double dose, respectively. In the older group of people aged ≡ 60 years, a slightly reduced age-dependent vaccine efficacy compared with younger people was considered [15]. The age-dependent relative vaccine efficacy is denoted by w60, w65, and w70 for the subgroups of 60-64, 65-69, and ≡70 years, respectively. Thus, the relative susceptibility of among vaccinated individuals in these age subgroups, described by 1 - wi (i ∊ {60, 65, 70}), was multiplied to the vaccinated portions of 60-64, 65-69, and ≡ 70 years.

    The cumulative proportion of people in each age group who experience infection by the end of epidemic (i.e., the final size) is expressed in the vector form z:

    z=(z00z05z70)=(1exp[(R00,00z00+R00,05z05++R00,70z70)]1exp[(R05,00z00+R05,05z05++R05,70z70)]1exp[(R70,00z00+R70,05z05++R70,70z70)]) (4)

    Unfortunately, there is no analytical solution for this vector, but the recursive equation (4) can be iteratively solved. The final size of the infection in total (ZI) was calculated by taking a dot product of z and the population fraction vector p (i.e., the relative size of the discrete age-specific population):

    ZI=Nzp (5)

    where N is the total population size of Japan (i.e., 125,620,000 people). Similarly, the final size of death in total (ZD) was calculated by taking a dot product of z and the IFR vector f:

    ZD=Nz(pf) (6)

    As mentioned above, f was retrieved from the literature [19] and is available in the Supplementary material. Table 2 shows the parameter values and their references that we used. We assumed that the basic reproduction number (R0) was 2.5. On basis of a study conducted by Dagan et al. [23], the relative susceptibility after a single and double dose, denoted by vs and vd, were fixed at 50% and 5%, respectively. The age-dependent relative vaccine efficacy, w60, w65, and w70, were assumed to be 83.33%, 83.33%, and 66.67%, respectively [15]. To examine the sensitivity of our result (i.e., optimal vaccination target) regarding the age-dependent relative susceptibility, we also conducted the same computation with the relative susceptibility after a single-dose scheme (vs), which was set at 35% and 20% [21,23]. Lastly, we carried out a sensitivity analysis of the cumulative risks of infection and death, varying the relative risk of death to 1%, 5% and 10%, respectively.

    Table 2.  Parameter values used for calculating the final size of infection/death with comparative vaccination strategies against the COVID-19 epidemic in Japan.
    Parameter Value(s) Reference
    Basic reproduction number (R0) 2.5 [38,39,40]
    Relative susceptibility after one dose of vaccination (vs) 50%, 35%, and 20% Assumption and[21,23]
    Relative susceptibility after two doses of vaccination (vd) 1%, 5%, and 10% [23]
    Age-dependent relative vaccine efficacy in 60-64 (w60), 65-69 (w65), and 70- (w70) 83.33%, 83.33%, and 66.67% [15]

     | Show Table
    DownLoad: CSV

    All quantitative analyses as mentioned above were implemented using statistical package R version 4.1.0 (The Comprehensive R Archive Network; https://cran.r-project.org/).

    Figure 1 shows the expected cumulative risk of infection as a function of doses consumed using the following two different vaccination strategies: (i) all doses are randomly provided to 36 million people older than 65 years who receive the double-dose scheme; and (ii) 50% of the doses are randomly provided to all people older than 65 years, and the rest are randomly allocated to 69 million adults aged between 20 and 64 years who receive the single-dose scheme. Assuming that the relative susceptibility after a single dose was 50%, in the absence of vaccination, 89.1% of the population were expected to experience infection, leading to 2,652,000 deaths. The cumulative risk of infection gradually decreased as the vaccination coverage increased and changed to 78.2% and 78.3% when all 72 million doses were consumed when following strategies (i) and (ii), respectively (Figure 1A). The cumulative risk of death was expected to change more drastically in strategy (i) than in (ii). The cumulative mortality would decrease from 2,652,000 to 1,856,000 in strategy (i), while strategy (ii) would only decrease to 2,355,000 if all 72 million doses were administered (Figure 1D). However, if we assumed that single-dose vaccination would contribute to yielding a relative susceptibility at 20% (instead of 50%) and that all 72 million doses were consumed, the estimated cumulative risk of infection would be 78.2%, leading to 1856,000 deaths in strategy (i). Additionally, the estimated cumulative risk of infection would be 63.5%, leading to 1833,000 deaths in strategy (ii) (Figure 1F).

    Figure 1.  Cumulative risk of infection (percentage in relation to the total population of Japan; panels A, B, and C) and the cumulative number of deaths (panels D, E, and F). These data include the total number of doses randomly administered to people older than 65 years receiving the double-dose strategy (solid lines), or half administered to adults aged between 20 and 64 years and the other half administered to people older than 65 years taking the single-dose scheme (dotted lines). The relative susceptibility after a single dose is assumed to be 50% (panels A and D), 35% (panels B and E), and 20% (panels C and F).

    Figure 2 shows the cumulative risk of infection and death as a function of single-dose coverage among adults and different double-dose coverage among older people. When vs was assumed to be 50%, with double-dose coverage of 30% among older people combined with single-dose coverage of 100% among adults, the final size of infection would be 65.3% with 2,317,000 deaths (90.5 million doses are required to achieve these data; Figure 2A and 2D). With the same coverage, if vs was assumed to be 20% (instead of 50%), the final size of infection would be 23.9% with 1,253,000 deaths (Figure 2C and 2F). However, if double-dose coverage among older people was elevated to 90% (with vs assumed to be 50%), while single-dose coverage among adults was 100%, the final size of infection would change to 54.6% with 1,588,000 deaths (Figure 2A and 2D). If vs was assumed to be 20%, the final size of infection would be 12.0% with 421,000 deaths for the same of combination of parameters (134 million doses are required to achieve these data; Figure 2C and 2F). In this scenario, the effective reproduction number was reduced from 2.50 to 1.14, which is still above the herd immunity threshold.

    Figure 2.  Cumulative risk of infection (percentage in relation to the total population of Japan; panels A, B, and C) and the cumulative number of deaths (panels D, E, and F) with single-dose coverage among adults aged between 20 and 64 years, along with different types of double-dose coverage among people aged 65 years and older. The relative susceptibility after a single dose, vs, was assumed to be 50% (A and D), 35% (B and E), and 20% (C and F), while that after a double dose, vd, was fixed at 5%.

    Figure 3 shows the results from the sensitivity analysis with respect to vd. When vs was fixed at 35% and all 72 million doses were consumed for people aged 65 years and older following strategy (i) 77.2% and 79.5% of the population would eventually experience infection, leading to 1,770 thousand and 1,954 thousand deaths, if the relative susceptibility after a double dose, vd, was 1% and 10%, respectively.

    Figure 3.  Cumulative risk of infection (percentage in relation to the total population of Japan; panels A, B, and C) and the cumulative number of deaths (panels D, E, and F) with single-dose coverage among adults aged between 20 and 64 years, along with different types of double-dose coverage among people aged 65 years and older. The relative susceptibility after a double dose, vd, was assumed to be 1% (A and D), 5% (B and E), and 10% (C and F), while that after a single dose, vs, was fixed at 35%.

    To answer a policy question regarding the prioritization of COVID-19 vaccines, we used an age-dependent heterogeneous transmission model and derived a final size equation. We compared the final sizes in various possible scenarios of the total dose, using either a double-dose or a single-dose (partial) scheme. Our study showed that enlarging the coverage by using a single dose would be beneficial, especially when the reduced susceptibility with one shot was substantial. However, considering that a 50% reduction in susceptibility following one shot was plausible, the single-dose scheme did not outperform the conventional double-dose strategy.

    There are two important take-home messages from the present study. First, while the single-dose scheme is expected to contribute to reducing transmission and elevating herd immunity in the population, this strategy is not recommended as an actual policy. Our simulation showed that prioritizing people aged 65 years and older was vital for reducing the death toll. This finding is compatible with a previous study, which showed that vaccine priority for people older than 60 years would reduce mortality and years of life lost [15]. This finding can be explained by the high IFR among older people, especially among those aged 75 years and older [19]. Therefore, the decrease in the relative susceptibility among individuals vaccinated with two doses, especially in older age groups, directly leads to the decrease in the death size (Figure 3D to F). Reducing the incidence of infection was the primary goal of this study, and the single-dose scheme was not substantially inferior to the double-dose scheme. Furthermore, Japan has had intense anti-vaccine campaigns since investigations for adverse effects of human papillomavirus vaccine started in 2014 [24]. However, elevating the vaccination coverage among young adults by the single-dose scheme with insufficient efficacy could be challenging. Second, the vaccine stocks that only covered the older population were not sufficient to contain the epidemic, and a much greater amount was required to anticipate the herd immunity threshold. Because of the limited efficacy, our simulation showed that the minimum double-dose vaccine coverage needed to be 99.4% in those aged 20 years or older to reach the effective reproduction number to fall below 1.0 only using vaccination. Importing vaccines to cover virtually the whole nation has been a formidable task. Therefore, other PHSM, including social distancing and lockdowns, have been required as possible options to concurrently prevent the spread of the epidemic [25].

    Another concern of the ongoing epidemic of COVID-19 is that we have repeatedly faced the threat of variants of concern (VOCs), including Alpha, Delta, and Omicron variants. In addition to the elevated transmissibility of these VOCs compared with the wild type (which leads to a higher basic reproduction number), studies have shown slightly lower titers in vaccinated sera against a UK VOC and marked resistance against vaccines in other VOCs (e.g., the Beta variant), even after taking a double dose [3,26,27,28,29,30,31,32]. A previous study showed that increasing primary vaccination coverage against the Omicron variant would lead to a reduction in hospitalizations and deaths [33]. Nevertheless, as a practical recommendation to the present status in Japan, the partial effectiveness of vaccinations against these VOCs suggests that rapid implementation of mass vaccination is still urgent to reduce the overwhelming hospital burden. A reduction in the incidence of new COVID-19 cases is required to sustain the socioeconomic burden and to maintain human casualties at the lowest level before new types of vaccines that are effective against VOCs are available.

    In the present study, we addressed the vaccine prioritization issue in a static parameter setting; the cumulative number of cases and deaths can be calculated by simply solving a set of mathematical equations. Using this model, vaccination was assumed to be completed in advance of the epidemic, and key parameters, including the basic reproduction number, were not assumed to vary over the course of the epidemic. While the final size did not differ substantially with varying relative susceptibility after single- or double-doses of vaccines (Figure 1 and 3), the equation allowed us to solve the policy question in a simple manner. However, the use of a final size equation was based on the assumption that the dynamics of the entire epidemic were sufficiently captured by the assumed model. If the epidemic starts during the course of vaccination, and if several epidemiological parameters, including transmissibility, vary over the course of time, other approaches accounting for such dynamic process are required. Partly because of this technical issue, and also because modelling study was not necessarily considered during decision making process by all countries, policy decisions on this matter across the world were highly variable. From ethical point of view, it was reasonable to choose double dose strategy as a convention, but some dynamic modelling studies indicated that single dose would outperform [4,16]. As a consequence, it was evident in real time that modelling studies yield diverse policy recommendations depending on varying assumptions and methodologies, and in practice, it was difficult to make a proper conclusive judgement in real time. In fact, a published study emphasized an importance of sociocultural factors that change the transmission dynamics, such as cultural difference from region to region, attitude changes, educational backgrounds, and people*s adherence to prevention activity [34]. This point remained to be a future subject for policy science to resolve possibly in advance of pandemic event.

    There are several limitations to this study. First, the vaccine rollout rate was not taken into consideration. We only computed the cumulative risk of infection from the final size and the cumulative number of deaths as a weighted average, and these are what we expect to observe in a long span. The time gap between vaccination and the titers to reach a protective level was also disregarded, but they did not matter during the process of final size computation [35]. Second, the initial infection size was disregarded. As of 30 April 2021, Japan has experienced 590,000 reported cases with > 10,000 deaths, accounting for 0.5% and 0.01% of the whole population, respectively [36]. Even though these numbers are negligible, the computed infection size and death toll might differ if the initial infection size was considered, especially if using a similar approach at a later stage of the pandemic. Third, VOCs were ignored in this analysis, but Japan has continuously experienced a series of epidemics induced by VOCs [37]. As mentioned above, VOCs induced altered susceptibility against vaccines, and the cumulative risks could have been larger than what was computed in the present study. Fourth, heterogeneities other than age were not explicitly considered. Vaccine priority in Japan was given to medical professionals who are at a particularly elevated risk of infection.

    We believe that our simple method will be useful in public health decision-making in determining vaccine dosage and distribution. When the vaccine stock is limited, enlarging vaccination coverage by using a single-dose regimen would not substantially outperform a double-dose scheme in avoiding mortality.

    This work was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI (grant number to TK: 21K10467), the Health Care Science Institute (IKEN), and Fujiwara Foundation. HN received funding from Health and Labour Sciences Research Grants (20CA2024, 20HA2007, 21HB1002, and 21HA2016), the Japan Agency for Medical Research and Development (JP20fk0108140 and JP20fk0108535), the JSPS KAKENHI (21H03198), Environment Research and Technology Development Fund (JPMEERF20S11804) of the Environmental Restoration and Conservation Agency of Japan, the Japan Science and Technology Agency CREST program (JPMJCR1413), and the SICORP program (JPMJSC20U3 and JPMJSC2105). We thank local governments, public health centers, and institutes for surveillance, laboratory testing, epidemiological investigations, and data collection. We thank Ellen Knapp, PhD, from Edanz (https://jp.edanz.com/ac) for editing a draft of this manuscript. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

    The authors declare that there is no conflict of interests regarding the publication of this paper.

    Table S1.  Age-specific population fraction.
    Age group Fraction
    0-4 (p00) 0.037
    5-9 (p05) 0.040
    10-14 (p10) 0.043
    15-19 (p15) 0.045
    20-24 (p20) 0.051
    25-29 (p25) 0.050
    30-34 (p30) 0.052
    35-39 (p35) 0.059
    40-44 (p40) 0.066
    45-49 (p45) 0.078
    50-54 (p50) 0.070
    55-59 (p55) 0.063
    60-64 (p60) 0.059
    65-69 (p65) 0.064
    70- (p70) 0.224

     | Show Table
    DownLoad: CSV
    Table S2.  Age-specific IFR.
    Age group IFR
    0-4 (f00) 0.00004
    5-9 (f05) 0.00004
    10-14 (f10) 0.00004
    15-19 (f15) 0.00004
    20-24 (f20) 0.00004
    25-29 (f25) 0.00004
    30-34 (f30) 0.00004
    35-39 (f35) 0.00068
    40-44 (f40) 0.00068
    45-49 (f45) 0.0023
    50-54 (f50) 0.0023
    55-59 (f55) 0.0075
    60-64 (f60) 0.0075
    65-69 (f65) 0.025
    70- (f70) 0.085
    IFR, infection fatality risk. Retrieved from [19].

     | Show Table
    DownLoad: CSV


    [1] F. Ding, T. Chen, Iterative least-squares solutions of coupled Sylvester matrix equations, Syst. Control Lett., 54 (2005), 95–107. https://doi.org/10.1016/j.sysconle.2004.06.008 doi: 10.1016/j.sysconle.2004.06.008
    [2] F. Ding, T. Chen, On iterative solutions of general coupled matrix equations, SIAM J. Control Optim., 44 (2006), 2269–2284. https://doi.org/10.1137/S036301290444135 doi: 10.1137/S036301290444135
    [3] P. M. Bentler, S. Y. Lee, Matrix derivatives with chain rule and rules for simple, Hadamard, and Kronecker products, J. Math. Psychol., 17 (1978), 255–262. https://doi.org/10.1016/0022-2496(78)90020-2 doi: 10.1016/0022-2496(78)90020-2
    [4] A. Graham, Kronecker Products and Matrix Calculus with Applications, John Wiley & Sons, 1982.
    [5] J. R. Magnus, H. Neudecker, Matrix differential calculus with applications to simple, Hadamard, and Kronecker products, J. Math. Psychol., 29 (1985), 474–492. https://doi.org/10.1016/0022-2496(85)90006-9 doi: 10.1016/0022-2496(85)90006-9
    [6] W. H. Steeb, Y. Hardy, Matrix Calculus and Kronecker Product: A Practical Approach to Linear and Multilinear Algebra, World Scientific Publishing Company, 2011. https://doi.org/10.1142/8030
    [7] F. Ding, Transformations between some special matrices, Comput. Math. Appl., 59 (2010), 2676–2695. https://doi.org/10.1016/j.camwa.2010.01.036 doi: 10.1016/j.camwa.2010.01.036
    [8] Y. Shi, H. Fang, M. Yan, Kalman filter-based adaptive control for networked systems with unknown parameters and randomly missing outputs, Int. J. Robust Nonlinear Control, 19 (2009), 1976–1992. https://doi.org/10.1002/rnc.1390 doi: 10.1002/rnc.1390
    [9] Y. Shi, B. Yu, Output feedback stabilization of networked control systems with random delays modeled by Markov chains, IEEE Trans. Autom. Control, 54 (2009), 1668–1674. https://doi.org/10.1109/TAC.2009.2020638 doi: 10.1109/TAC.2009.2020638
    [10] X. L. Xiong, W. Fan, R. Ding, Least-squares parameter estimation algorithm for a class of input nonlinear systems, J. Appl. Math., 2012, 684074. https://doi.org/10.1155/2012/684074
    [11] M. Dehghan, M. Hajarian, An iterative algorithm for the reflexive solutions of the generalized coupled Sylvester matrix equations and its optimal approximation, Appl. Math. Comput., 202 (2008), 571–588. https://doi.org/10.1016/j.amc.2008.02.035 doi: 10.1016/j.amc.2008.02.035
    [12] L. Jodar, H. Abou-Kandil, Kronecker products and coupled matrix Riccati differential systems, Linear Algebra Appl., 121 (1989), 39–51. https://doi.org/10.1016/0024-3795(89)90690-3 doi: 10.1016/0024-3795(89)90690-3
    [13] P. A. Regaliat, S. K. Mitra, Kronecker products, unitary matrices and signal processing applications, SIAM Rev. 31 (1989), 586–613. https://doi.org/10.1137/1031127
    [14] J. G. Nagy, M. K. Ng, L. Perrone, Kronecker product approxmations for image restoration with reflexive boundary conditions, SIAM J. Matrix Anal. Appl., 25 (2004), 829–841. https://doi.org/10.1137/S0895479802419580 doi: 10.1137/S0895479802419580
    [15] F. M. Fernández, The Kronecker product and some of its physical applications, Eur. J. Phys., 37 (2016), 065403. https://doi.org10.1088/0143-0807/37/6/065403 doi: 10.1088/0143-0807/37/6/065403
    [16] J. R. Magnus, H. Neudeckery, Matrix Differential Calculus with Applications in Statistics and Econometrics, John Wiley & Sons, 2019. https://doi.org/10.1002/9781119541219
    [17] C. C. Santos, C. C. Dias, Matrix Multiplication without size restrictions: An approach to the Kronecker product for applications in statistics, J. Math. Stat. Res., 4 (2022), 156. https://doi.org/10.36266/JMSR/156 doi: 10.36266/JMSR/156
    [18] E. Tyrtyshnikov, Kronecker-product approximations for some function-related matrices, Linear Algebra Appl., 379 (2004), 423–437. https://doi.org/10.1016/j.laa.2003.08.013 doi: 10.1016/j.laa.2003.08.013
    [19] H. Zhang, F. Ding, On the Kronecker products and their applications, J. Appl. Math., 2013 (2013), 296185. https://doi.org/10.1155/2013/296185 doi: 10.1155/2013/296185
    [20] B. Cannon, On the convergence of series of polynomials, Proc. London Math. Soc., 43 (1937), 348–365. https://doi.org/10.1112/plms/s2-43.5.348 doi: 10.1112/plms/s2-43.5.348
    [21] B. Cannon, On the representation of integral functions by general basic series, Math. Z., 45 (1939), 185–208. https://doi.org/10.1007/BF01580282 doi: 10.1007/BF01580282
    [22] J. M. Whittaker, On Series of poynomials, Q. J. Math., 5 (1934), 224–239. https://doi.org/10.1093/qmath/os-5.1.224
    [23] J. M. Whittaker, Interpolatory Function Theory, Cambridge Tracts in Math, 1935.
    [24] J. M. Whittaker, Sur Les Se Ries De Base De Polynomes Quelconques, AGauthier-Villars, 1949.
    [25] A. El-Sayed Ahmed, On Some Classes and Spaces of Holomorphic and Hyperholomorphic Functions, Ph.D thesis, Bauhaus University Weimar-Germany, 2003.
    [26] A. El-Sayed Ahmed, Z. M. G. Kishka, On the effectiveness of basic sets of polynomials of several complex variables in elliptical regions, in Progress in Analysis, (2003), 265–278. https://doi.org/10.1142/9789812794253_0032
    [27] Z. M. G. Kishka, A. El-Sayed Ahmed, On the order and type of basic and composite sets of polynomials in complete reinhardt domains, Period. Math. Hung. 46 (2003), 67–79. https://doi.org/10.1023/A: 1025705824816
    [28] Z. G. Kishka, M. A. Saleem, M. A. Abul-Dahab, On Simple Exponential Sets of Polynomials, Mediterr. J. Math., 11 (2014), 337–347. https://doi.org/10.1007/s00009-013-0296-7 doi: 10.1007/s00009-013-0296-7
    [29] M. Nassif, Composite sets of polynomials of several complex variables, Publicationes Math. Debrecen, 18 (1971), 43–52. https://doi.org/10.5486/PMD.1971.18.1-4.05 doi: 10.5486/PMD.1971.18.1-4.05
    [30] M. A. Abul-Ez, D. Constales, Basic sets of polynomials in Clifford analysis, Complex Variables Theory Appl. Int. J., 14 (1990), 177–185. https://doi.org/10.1080/17476939008814416 doi: 10.1080/17476939008814416
    [31] M. A. Abul-Ez, D. Constales, On the order of basic series representing Clifford valued functions, Appl. Math. Comput., 142 (2003), 575–584. https://doi.org/10.1016/S0096-3003(02)00350-8 doi: 10.1016/S0096-3003(02)00350-8
    [32] G. F. Hassan, M. Zayed, Approximation of monogenic functions by hypercomplex Ruscheweyh derivative bases, Complex Var. Elliptic Equations, 68 (2022), 2073–2092. https://doi.org/10.1080/17476933.2022.2098279 doi: 10.1080/17476933.2022.2098279
    [33] G. F. Hassan, M. Zayed, Expansions of generalized bases constructed via Hasse derivative operator in Clifford analysis, AIMS Math., 8 (2023), 26115–26133. https://doi.org/10.3934/math.20231331 doi: 10.3934/math.20231331
    [34] G. F. Hassan, E. Abdel-salam, R. Rashwan, Approximation of functions by complex conformable derivative bases in Fréchet spaces, Math. Method. Appl. Sci., 46 (2022), 2636–2650. https://doi.org/10.1002/mma.8664 doi: 10.1002/mma.8664
    [35] M. Zayed, G. F. Hassan, On the approximation of analytic functions by infinite series of fractional Ruscheweyh derivatives bases, AIMS Math., 9 (2024), 8712–8731. https://doi.org/10.3934/math.2024422 doi: 10.3934/math.2024422
    [36] M. Zayed, G. F. Hassan, E. A. B. Abdel-Salam, On the convergence of series of fractional Hasse derivative bases in Fréchet spaces, Math. Meth. Appl. Sci., 47 (2024), 8366–8384. https://doi.org/10.1002/mma.10018 doi: 10.1002/mma.10018
    [37] A. El-Sayed Ahmed, Hadamard product of simple sets of polynomials in Cn, Theory Appl. Math. Comput. Sci., 4 (2014), 26–30.
    [38] M. A. Abel Dahab, A. M. Saleem, Z. M. Kishka, Effectiveness of Hadamard product of basic sets of polynomials of several complex variables in hyperelliptical regions, Electrumic J. Math. Anal. Appl., 3 (2015), 52–65. https://doi.org/10.21608/ejmaa.2015.310695 doi: 10.21608/ejmaa.2015.310695
    [39] H. V. Henderson, F. Pukelsheim, S. R. Searle, On the History of the Kronecker product, Linear Multilinear Algebra, 14 (1983), 113–120. https://doi.org/10.1080/03081088308817548 doi: 10.1080/03081088308817548
    [40] B. Holmquist, The direct product permuting matrices, Linear Multilinear Algebra, 17 (1985), 117–141. https://doi.org/10.1080/03081088508817648 doi: 10.1080/03081088508817648
    [41] R. A. Horn, C. R. Johnson, Topics in Matrix Analysis, Cambridge University Press, 1991. https://doi.org/10.1017/CBO9780511840371
    [42] A. A. Abd-El-Monem, M. Nassif, On the convergence of certain classes of Gontcharoff polynomials, Indagationes Math., 68 (1965), 615–621. https://doi.org/1010.1016/s1385-7258(65)50062-7
    [43] M. Abdalla, M. A. Abul-Ez, J. Morais, On the construction of generalized monogenic Bessel polynomials, Math. Methods Appl. Sci., 41 (2018), 9335–9348. https://doi.org/10.1002/mma.5274 doi: 10.1002/mma.5274
    [44] M. A. Abul-Ez, Bessel polynomial expansions in spaces of holomorphic functions, J. Math. Anal. Appl., 221 (1998), 177–190. https://doi.org/10.1006/jmaa.1997.5840 doi: 10.1006/jmaa.1997.5840
    [45] M. A. Abul-Ez, M. Zayed, Criteria in Nuclear Fréchet spaces and Silva spaces with refinement of the Cannon-Whittaker theory, J. Funct. Spaces, 2020 (2020), 15. https://doi.org/10.1155/2020/8817877 doi: 10.1155/2020/8817877
    [46] G. F. Hassan, L. Aloui, Bernoulli and Euler polynomials in Clifford analysis, Adv. Appl. Clifford Algebras, 25 (2015), 351–376. https://doi.org/10.1007/s00006-014-0511-z doi: 10.1007/s00006-014-0511-z
    [47] M. Nassif, A class of a Gontcharoff polynomials, Assiut Univ. Bull. Sci. Technol., 1 (1958), 1–8.
    [48] K. A. M. Sayyed, Basic Sets of Polynomials of Two Complex Variables and Their Convergence Properties, Ph.D thesis, Assiut University, 1975.
  • This article has been cited by:

    1. Ali Kandil, Lei Hou, Mohamed Sharaf, Ayman A. Arafa, Configuration angle effect on the control process of an oscillatory rotor in 8-pole active magnetic bearings, 2024, 9, 2473-6988, 12928, 10.3934/math.2024631
    2. Jun Liu, Jinxiang Zhou, Xue Han, Shiqiang Zheng, High-speed magnetically levitated centrifugal hydrogen recirculation pump in proton exchange membrane fuel cells, 2024, 87, 03603199, 268, 10.1016/j.ijhydene.2024.09.002
    3. Zhihang Huang, Changhe Li, Zongming Zhou, Bo Liu, Yanbin Zhang, Min Yang, Teng Gao, Mingzheng Liu, Naiqing Zhang, Shubham Sharma, Yusuf Suleiman Dambatta, Yongsheng Li, Magnetic bearing: structure, model, and control strategy, 2024, 131, 0268-3768, 3287, 10.1007/s00170-023-12389-8
    4. Yi Yang, Xin Cheng, Rougang Zhou, Position Servo Control of Electromotive Valve Driven by Centralized Winding LATM Using a Kalman Filter Based Load Observer, 2024, 17, 1996-1073, 4515, 10.3390/en17174515
    5. Jiyuan Sun, Gengyun Tian, Pin Li, Chunlin Tian, Zhenxiong Zhou, Levitating Control System of Maglev Ruler Based on Active Disturbance Rejection Controller, 2024, 14, 2076-3417, 8069, 10.3390/app14178069
    6. Lixin Liu, Cheng Liu, Changjin Che, Yunbo Wu, Qing Zhao, Research on the Coordinated Control of Mining Multi-PMSM Systems Based on an Improved Active Disturbance Rejection Controller, 2025, 14, 2079-9292, 477, 10.3390/electronics14030477
    7. Hongkui Zhang, Qinwei Zhang, Hang Shen, Yipeng Lan, Jiaqi Wen, Chuan Zhao, Improve LADRC Strategy for Variable Air Gap Permanent Magnetic Levitation System, 2025, 13, 2169-3536, 61641, 10.1109/ACCESS.2025.3549146
    8. Gengyun Tian, Chunlin Tian, Jiyuan Sun, Shusen Diao, Horizontal Control System for Maglev Ruler Based on Improved Active Disturbance Rejection Controller, 2025, 15, 2076-3417, 4938, 10.3390/app15094938
  • Reader Comments
  • © 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(374) PDF downloads(19) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog