Mathematical Modeling to Predict COVID-19 Infection and Vaccination Trends
by
, , , and
Bogdan Doroftei
1
,
Ovidiu-Dumitru Ilie
2,*
,
Nicoleta Anton
1
,
Sergiu-Ioan Timofte
2 and
Ciprian Ilea
1 1
Faculty of Medicine, University of Medicine and Pharmacy “Grigore T. Popa”, University Street, No. 16, 700115 Iasi, Romania
2
Department of Biology, Faculty of Biology, “Alexandru Ioan Cuza” University, Carol I Avenue, No. 20A, 700505 Iasi, Romania
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2022, 11(6), 1737; https://doi.org/10.3390/jcm11061737
Submission received: 25 January 2022
/
Revised: 15 March 2022
/
Accepted: 18 March 2022
/
Published: 21 March 2022
(This article belongs to the Collection Coronavirus Disease 2019: Clinical Presentation, Pathogenesis and Treatment)
Abstract
:Background: COVID-19 caused by the Severe Acute Respiratory Syndrome Coronavirus 2 placed the health systems around the entire world in a battle against the clock. While most of the existing studies aimed at forecasting the infections trends, our study focuses on vaccination trend(s). Material and methods: Based on these considerations, we used standard analyses and ARIMA modeling to predict possible scenarios in Romania, the second-lowest country regarding vaccinations from the entire European Union. Results: With approximately 16 million doses of vaccine against COVID-19 administered, 7,791,250 individuals had completed the vaccination scheme. From the total, 5,058,908 choose Pfizer–BioNTech, 399,327 Moderna, 419,037 AstraZeneca, and 1,913,978 Johnson & Johnson. With a cumulative 2147 local and 17,542 general adverse reactions, the most numerous were reported in recipients of Pfizer–BioNTech (1581 vs. 8451), followed by AstraZeneca (138 vs. 6033), Moderna (332 vs. 1936), and Johnson & Johnson (96 vs. 1122). On three distinct occasions have been reported >50,000 individuals who received the first or second dose of a vaccine and >30,000 of a booster dose in a single day. Due to high reactogenicity in case of AZD1222, and time of launching between the Pfizer–BioNTech and Moderna vaccine could be explained differences in terms doses administered. Furthermore, ARIMA(1,1,0), ARIMA(1,1,1), ARIMA(0,2,0), ARIMA(2,1,0), ARIMA(1,2,2), ARI-MA(2,2,2), ARIMA(0,2,2), ARIMA(2,2,2), ARIMA(1,1,2), ARIMA(2,2,2), ARIMA(2,1,1), ARIMA(2,2,1), and ARIMA (2,0,2) for all twelve months and in total fitted the best models. These were regarded according to the lowest MAPE, p-value (p < 0.05, p < 0.01, and p < 0.001) and through the Ljung–Box test (p < 0.05, p < 0.01, and p < 0.001) for autocorrelations. Conclusions: Statistical modeling and mathematical analyses are suitable not only for forecasting the infection trends but the course of a vaccination rate as well.
1. Introduction
Even though the number of human coronavirus infections per year is low [1], the current health crisis caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has placed each country worldwide in an unprecedented state of emergency [2]. Originally was theorized to be another low-level respiratory outbreak similar to the previous health crises [3,4,5,6] and pathophysiology analogies with SARS-CoV-2 [7].
It led initially to a relatively insignificant number of pneumonia cases [2] with an unknown etiology. Additional experiments brought new data and insight regarding the emergence of a novel beta-strain that belongs to the zoononic coronavirus (CoV) that is particularly virulent to humans [8]. Until the genome of this pathogen was fully sequenced, the number of cases confirmed already reached 15 with 1 fatality [9].
With fifty-nine confirmed cases at the beginning of 2020 [10], it became the main priority reflected by the associated mortality rate and high tropism towards the respiratory system [11]. COVID-19 has caused >280 million infections, and >5 million individuals passed away, from which more than 1.8 million infections and >58.00 deaths in Romania. The peak since the pandemic started in terms of infections and deaths was reached on 20th October with 18,863 cases of infections, and on 3 November with 591 deaths.
The need to model daily cases is essential to predict future directions. It is even more imperative for Eastern European countries such as Romania [12,13]. There are multiple reasons and limitations compared to Westernized countries [14,15]. Therefore, mathematical and statistical models became an integrated component within the current methodologies, dedicated to understanding the infectious and/or vaccination trend over a dedicated time interval [16]. Distinct approaches with relatively high accuracy in predicting trends failed because of randomness tendency of epidemics [17].
Systematically designed, approved, and distributed vaccines manufactured based on distinct procedures are already extensively discussed elsewhere [18,19,20,21] as a countermeasure to prophylactic strategies [22]. The first person vaccinated in Romania was a medical doctor with the candidate vaccine BNT162b1/BNT162b2 manufactured by Pfizer–BioNTech on 27 December 2020. On 4th February, the Romanian government approved and began using mRNA-1273 by Moderna vaccine, whereas AZD1222 (ChAdOx1) by AstraZeneca received approval 11 days later on 15 February. However, only on 4 May Ad26.COV2-S from Johnson & Johnson was available for the general population.
Such studies are crucial because of the low vaccination rate registered in middle-class countries. Presently, Romania has the second-lowest rate of vaccination among all European Union countries after Bulgaria. This subject caused a lot of controversies around individuals. Among the main factors that contribute to this generalized skepticism are of personal nature or beliefs.
Based on all these considerations, this article focuses on the number of doses administered depending on the manufacturer, those who completed or not the vaccination scheme, adverse reactions displayed by the recipients, and forecasts.
2. Materials and Methods
2.1. Data Collection and Parameters Analyzed
Data related to (1) the number of doses administered (dose 1, dose 2, and booster), (2) those that had or not completed the vaccination scheme, and (3) adverse reactions of local or general type were taken from the official website of the Romanian Government (https://vaccinare-covid.gov.ro/ (accessed on 28 December 2021)). We centralized all data by creating a time-series database using MS Excel for the following predetermined interval (27 December 2020–27 December 2021).
2.2. ARIMA
2.3. Statistical Analyses
Data analysis for three out of four parameters was carried out using Microsoft Excel 2010, whereas the software used for ARIMA modeling is STATGRAPHICS Centurion (18.1.14).
3. Results
We observe that a significant percentage of the Romanian citizens were willing to receive the serum from Pfizer–BioNTech. Following the centralization and analysis of data, we noted a fluctuating trend that lasted several months before reaching the peak on 21st April 2021. With 55,643 (SE-724.01, SD-13,851.18, CI95%-1423.76) individuals immunized in a single day, we expected an increase since the second dose became available starting with 17 January 2021. Precisely, 11 May 2021 was the day in which was registered the second most numerous groups with 54,757 (SE-732.40, SD-13,603.80, CI95%-1440.55). Unfortunately, this trend gradually decreased again from May until mid-October. After this decline, the phase normalized with 2 successful waves of recipients and a maximum of 35,425 (SE-691.87, SD-6600.03, CI95%-1374.52) recipients of booster on 29 September 2021. From that point onwards, all three trends significantly decrease, this possibly having to do with the relaxation of restrictions in our country. Even though the technology behind manufacturing is the same as BNT162b1/BNT162b2, the Moderna vaccine did not benefit from the same interest among citizens. With figures suggesting that approximately 15.27% of recipients have chosen mRNA-1273 by comparison with Pfizer–BioNTech, on 5 February 2021, 8501 (SE-107.58, SD-1930.60, CI95%-211.66) people received the first dose. The situation remained the same also towards the introduction of the second dose on 12 February 2021. The highest number of individuals that attended a vaccination center was 8225 (SE-112.24, SD-1944.12, CI95%-220.88) on 5 March 2021. On 9 October 2021, a total of 21,951 (SE-232.45, SD-2217.51, CI95%-461.81) people receive a booster dose. Furthermore, the highest number of individuals vaccinated with AZD1222 on one day was on 26 February and 23 April last year. Despite the relatively high figures in both cases, 11,284 (SE-189.48, SD-3067.14, CI95%-373.12 and 10.685 (SE-203.60, SD-3067.57, CI95%-401.20), trends declined almost entirely in short-time. Two hundred and forty-five (0.05%) individuals decided to get the second dose of AstraZeneca vaccine and fewer (n = 4) for booster only in December if we refer to the overall situation. Not only the number of doses administered did not exceed 100, but there is also a gap between 2 November and 14 December 2021. As opposed to the rest by design and technology, Johnson & Johnson is the only available vaccine that requires a single dose. Similar to other trends, this was linear most of the time over the analyzed interval, with one exception on 27 October 2021 with 57,359 (SE-552.61, SD-8525.30, CI95%-1088.66). On 8 November 2021 started to be administered the second dose of Ad26.COV2-S according to WHO, CDC, and EMA directives. Thus, a maximum of 568 (SE-25.24, SD-178.48, CI95%-50.72) booster doses on 21 December 2021 were administered (Figure 1).
At the end of the analyzed interval, a total of 7,791,250 individuals were completely immunized. From this, 64.93%, n = 5.058.908 with BNT162b1/BNT162b2, and 24.59%, n = 1,913,978 with Ad26.COV2-S. Only a small percentage selected mRNA-1273 or ChAdOx1 (10.49% n = 818,364). As observed, both reached a plateau phase, this highlighting the human reluctance as a consequence of high reactogenicity. In this context, only 5.12%, n = 399,327 and 5.37%, n = 419,037 were attributed to the last two COVID-19 vaccines (Figure 2).
There was a cumulative total of 2147 local adverse reactions and 17,542 general. Individually, 73.63%, n = 1581 were local (max = 56 on 5 January 2021) and 48.17%, n = 8451 (max = 170 on 6 February 2021) generalized in the case of BNT162b1/BNT162b2 followed by ChAdOx1 with 6.42%, n = 138 (max = 10 on 16 March 2021) and 34.39%, n = 6033 (max = 232 on 13th March 2021). Despite this, the Romanian citizens still opted for Pfizer–BioNTech. Through the prism of one critical argument might be explained this situation: the number of cases ≥100 on twenty-two different occasions, hence the lack of confidence in the serum from AstraZeneca. Continuing with this concept, 15.46%, n = 332 (max = 11 on 6 April 2021), whereas 11.03%, n = 1936 (max = 53 on 22 February 2021) of all adverse reactions were caused by Moderna. In more than eight months since its release, only 4.47% of the cases, n = 96 (max = 5 on 5 June 2021) led to local adverse reactions and 6.39%, n = 1122 (max = 45 on 26 October 2021) to those of generalized reported for Johnson & Johnson (Figure 3).
According to Elevli et al. [25], to successfully create any ARIMA model, it must first meet four conditions and evaluate if a series of values are constant throughout the analyzed period. Autocorrelation Function (ACF), and Partial Autocorrelation Function (PACF) (Figure 4), are the time-series plots developed to evaluate the seasonality and stationarity. ACF is a metric that describes if the previous values are related to the next ones, whereas PACF determines the degree of correlation coefficient between variable and lag [26]. The performance of the model and misspecification detection is measured through the Bayesian information criterion of Schwarz (BIC), and Akaike information criteria expression (AIC) [27]. Straight lines points to the limit of two standard deviations, while the bars that cross the lines indicate statistically meaningful autocorrelations (Figure 4).
Performances of multiple models were generated and interpreted. MAPE with the lowest value per statistical analysis was regarded as the best model. Among all models, ARIMA(1,1,0), ARIMA(1,1,1), ARIMA(0,2,0), ARIMA (2,1,0), ARIMA(1,2,2), ARIMA(2,2,2), ARIMA(0,2,2), ARIMA(2,2,2), ARIMA(1,1,2), ARIMA(2,2,2), ARIMA(2,1,1), ARIMA(2,2,1), and ARIMA(2,0,2) were chosen for all twelve months and total. The fitted models are presented in Figure 4 and Table 1 and Table 2 with a minimum MAPEJanuary = 6.81315, MAPEFebruary = 0.214755, MAPEMarch = 0.474466, MAPEApril = 0.117225, MAPEMay = 0.105518, MAPEJune = 0.0354902, MAPEJuly = 0.0251916, MAPEAugust = 0.0234608, MAPESeptember = 0.0352714, MAPEOctober = 0.0774352, MAPENovember = 0.0570842, MAPEDecember = 0.0187244, and MAPETotal = 0.733775.
Table 2 highlights the parameter estimates for the best models with CI95% and p-values < 0.05, which were further confirmed through the Ljung–Box test. The forecasted and fitted values are detailed in Table 3 and Figure 5 for the next 7 days. Thus, the forecasts for the next week may be between 538,694–681,206. in January, 919,209–991,628. in February, 1.94174 × 106–2.10285 × 106 in March, 3.20358 × 106–3.55743 × 106 in April, 4.25759 × 106–4.44011 × 106 in May, 4.70459 × 106–4.75999 × 106 in June, 4.97002 × 106–5.02345 × 106 in July, 5.23552 × 106–5.28803 × 106 in August, 5.49969 × 106–5.60622 × 106 in September, 6.8339 × 106–7.46267 × 106 in October, 7.72849 × 106–7.77887 × 106 in November, 7.92705 × 106–7.94037 × 106 in December, and between 7.92819 × 106–7.96656 × 106 in Total.
4. Discussion
Ad26.COV2-S from Johnson & Johnson and BNT162b1/BNT162b2 from Pfizer–BioNTech were the top two opted vaccines against COVID-19 in Romania with over >50,000 individuals immunized (57,359 unique dose vs. 55,643 with dose 1 and 54,757 with dose 2), excepting the third dose (35,425 with the booster). mRNA-1273 by Moderna and AZD1222 by AstraZeneca were the least favorite (8501 with dose 1, 8225 with dose 2, 21,951 booster vs. 11,284 with dose 1, 10,685 with dose 2, and 9 with the booster). These results are further consolidated by the number of individuals who had a complete vaccination scheme and adverse reactions, both local and general; BNT162b1/BNT162b2 = 5,058,908 (1581 vs. 8451), Ad26.COV2-S = 1,913,978 (96 vs. 1122), mRNA-1273 = 399,327 (332 vs. 1936), and AZD1222 = 419.037 (138 vs. 6033).
As per ARIMA generated, the optimal models through the Ljung–Box Test are: ARIMA(1,1,0), ARIMA(1,1,1), ARIMA(0,2,0), ARIMA(2,1,0), ARIMA(1,2,2), ARIMA(2,2,2), ARIMA(0,2,2), ARIMA(2,2,2), ARIMA(1,1,2), ARIMA(2,2,2), ARIMA(2,1,1), ARIMA(2,2,1), and ARIMA (2,0,2) for all twelve months and in total.
The most recent article with a similar design is that of Cihan [28]. With over four hundred million people forecasted to be fully vaccinated by June 2021, the figures per hundred million were as follows: 147 in Europe (17%), 139 in Asia (2.3%), 130 in South America (8.8%), 129 in the US (41.8%), and 109 in Africa (0.6%) and 5.6% of the World. The optimal ARIMA models through the Ljung–Box test are ARIMA (5,2,2) in the US, ARIMA (1,2,3) in Asia, ARIMA (5,2,0) in Europe, ARIMA (2,2,1) in Africa, ARIMA (1,2,1) in South America, and ARIMA (5,2,1) in the World.
There will be an increased interest in vaccines between 5–10% in the first quarter of 2022. The short-term forecasts apply for influenza, HPV, pneumococcal, and polio vaccines, without being detected a decline in the overall interest for the COVID-19 vaccine [29]. However, confinement is still one of the most suitable prevention measures. The number of pediatric consultations/antenatal visits decreased by 52%/45% in April and 34%/34% in May 2020 compared to the same periods of 2019 (p = 0.0001), and demand for immunization significantly decreased as well [30].
Another nationwide article conducted is by Lumbreras-Marquez et al. [31], in which they briefly discuss the risk of morbidity and mortality among Mexican pregnant women. With 934 deaths in 2020, the maternal mortality ratio (MMR) was 46.6 per 100.00 live births, with 202 attributed to COVID-19. Around 31% (286/934) was associated with respiratory failure in contrast to 5% between 2011–2019 since the Mexican government launched the vaccination program on 11 May 2021. Assuming 100% vaccination among women, the authors forecasted weekly maternal deaths that might occur and obtained 993 deaths with an MMR of 46.5; RMSE (0,1,1) was 5.57 and 6.15 in 2021 compared to 2020 (21.6%). The overall figures would decrease to 885 and to an MMR of 41.5.
Distinct authors employed other mathematical models to perform predictive analyses in various countries [32,33]. Hwang [32] adopted a heterogeneous autoregression (HAR) model due to the long-memory feature of COVID-19 based on the growth and vaccination rates. Three novel perspectives derive from this protocol. The first refers to the combination of both growth and vaccination rates, construction and comparison of three types of predictions and coverage probability improvement, and mean interval score of prediction periods via bootstrap technique.
The Susceptible–Infected–Recovered (SIR) model [34] fits within the scope, with an expected effective reproduction term (tR) less than 1. According to a recent article, the tR reduces at fast rates when the values of tR are high, the slope being dependent on the promptness response parameter.
A multinomial autoregressive model for time series of counts was introduced with the aim of analyzing the finite-range integer-valued data. For this, the estimations of the parameters were calculated using conditional least squares (CLS), weighted conditional least squares (WCLS), and conditional maximum likelihood (CML). The authors established the asymptotic properties of the estimators and performed simulation studies to certify the current procedure [35]. Bartolucci et al. [36] and proposed multinomial Bayesian and Dirichlet auto-regressive models for series of time-dependent data points centered on counting patients exclusive and exhaustive categorized on predefined groups. Specifically, they were allocated based on the severity and required treatments in either regular wards or intensive care units, along with individuals that passed away and went through the disease. Not only were formulated assumptions on the transition probabilities between categories over a specific period that previously had a normal distribution allowed, but the accessibility to incorporate hypotheses was offered. Markov chain Monte Carlo (MCMC) was employed to estimate the posterior distribution and transition matrices, also allowing to make predictions and compute the reproduction number (Rt), accuracy measured through Bayesian inference. In this way, the authors offer insight regarding data collection during the first wave in Lombardia, Italy, and the effect of non-pharmaceutical interventions. Furthermore, the Dirichlet-multinomial model is adequate in fitting/providing predictive performance for patients admitted in regular and intensive care units.
5. Conclusions
Our results emphasize the willingness of Romanian residents to get the vaccine. On the opposed pole, there is a significant discrepancy between the internal administration of Westernized countries with >50% of the overall population vaccinated and post-communist middle-class country such as Romania. Approximately 16 million doses have been administered since inception on 27 December 2020. BNT162b1/BNT162b2 and Ad26.COV2-S were the top two choices among the Romanian citizens, with figures comparable in contrast with mRNA-1273 and AZD1222 among all analyzed parameters. Statistical models still play a crucial role in making different predictions. As opposed to the existing literature, our study was focused on forecasting the vaccination rate, and not for establishing the infections trends. Following the analyses performed, ARIMA(1,1,0), ARIMA(1,1,1), ARIMA(0,2,0), ARIMA(2,1,0), ARIMA(1,2,2), ARI-MA(2,2,2), ARIMA(0,2,2), ARIMA(2,2,2), ARIMA(1,1,2), ARIMA(2,2,2), ARIMA(2,1,1), ARIMA(2,2,1), and ARIMA (2,0,2) are the best models that fit within the current situation from all scenarios generated based on their MAPE, p-value and through the Ljung–Box test. Conclusively, mathematical and statistical algorithms proved efficient in providing forecasts of either infectious or, in our case, vaccination trends in a country.
Author Contributions
B.D., O.-D.I. and S.-I.T. Conceptualization, Methodology, Data curation, Investigation, Formal analysis, Validation, Writing—original draft, N.A. and C.I. Conceptualization, Methodology, Writing—review and editing. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Van der Hoek, L. Human coronaviruses: What do they cause? Antivir. Ther. 2007, 12, 651–658. [Google Scholar] [PubMed]
- Harapan, H.; Itoh, N.; Yufika, A.; Winardi, W.; Keam, S.; Te, H.; Megawati, D.; Hayati, Z.; Wagner, A.L.; Mudatsir, M. Coronavirus disease 2019 (COVID-19): A literature review. J. Infect. Public Health 2020, 13, 667–673. [Google Scholar] [CrossRef] [PubMed]
- Zhou, N.; He, S.; Zhang, T.; Zou, W.; Shu, L.; Sharp, G.B.; Webster, R.G. Influenza infection in humans and pigs in southeastern China. Arch. Virol. 1996, 141, 649–661. [Google Scholar] [CrossRef] [PubMed]
- Yang, Y.; Peng, F.; Wang, R.; Yange, M.; Guan, K.; Jiang, T.; Xu, G.; Sun, J.; Chang, C. The deadly coronaviruses: The 2003 SARS pandemic and the 2020 novel coronavirus epidemic in China. J. Autoimmun. 2020, 109, 102434. [Google Scholar] [CrossRef]
- Ji, J.S. Origins of MERS-CoV, and lessons for 2019-nCoV. Lancet Planet. Health 2020, 4, e93. [Google Scholar] [CrossRef] [Green Version]
- Su, S.; Wong, G.; Liu, Y.; Gao, G.F.; Li, S.; Bi, Y. MERS in South Korea and China: A potential outbreak threat? Lancet 2015, 385, 2349–2350. [Google Scholar] [CrossRef]
- Yuki, K.; Fujiogi, M.; Koutsogiannaki, S. COVID-19 pathophysiology: A review. Clin. Immunol. 2020, 215, 108427. [Google Scholar] [CrossRef]
- Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef]
- Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
- Chen, Z.-L.; Zhang, Q.; Lu, Y.; Guo, Z.-M.; Zhang, X.; Zhang, W.-J.; Guo, C.; Liao, C.-H.; Li, Q.-L.; Han, X.-H.; et al. Distribution of the COVID-19 epidemic and correlation with population emigration from Wuhan, China. Chin. Med. J. 2020, 133, 1044–1050. [Google Scholar] [CrossRef]
- Kim, H.-J.; Hwang, H.; Hong, H.; Yim, J.-J.; Lee, J. A systematic review and meta-analysis of regional risk factors for critical outcomes of COVID-19 during early phase of the pandemic. Sci. Rep. 2021, 11, 9784. [Google Scholar] [CrossRef] [PubMed]
- Dascalu, S. The Successes and Failures of the Initial COVID-19 Pandemic Response in Romania. Front. Public Health 2020, 8, 344. [Google Scholar] [CrossRef] [PubMed]
- Al-Salem, W.; Moraga, P.; Ghazi, H.; Madad, S.; Hotez, P.J. The emergence and transmission of COVID-19 in European countries, 2019–2020: A comprehensive review of timelines, cases and containment. Int. Health 2021, 13, 383–398. [Google Scholar] [CrossRef] [PubMed]
- Vlădescu, C.; Scîntee, S.G.; Olsavszky, V.; Hernández-Quevedo, C.; Sagan, A. Romania: Health system review. Health Syst. Transit. 2016, 10, 1–172. [Google Scholar]
- Pullano, G.; Pinotti, F.; Valdano, E.; Boëlle, P.-Y.; Poletto, C.; Colizza, V. Novel coronavirus (2019-nCoV) early-stage importation risk to Europe, January 2020. Euro Surveill. 2020, 25, 2000057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Wang, L.; Li, J.; Guo, S.; Xie, N.; Yao, L.; Cao, Y.; Day, S.W.; Howard, S.C.; Graff, J.C.; Gu, T.; et al. Real-time estimation and prediction of mortality caused by COVID-19 with patient information based algorithm. Sci. Total Environ. 2020, 727, 138394. [Google Scholar] [CrossRef]
- Smith, G.D. Epidemiology, epigenetics and the ‘Gloomy Prospect’: Embracing randomness in population health research and practice. Int. J. Epidemiol. 2011, 40, 537–562. [Google Scholar] [CrossRef] [Green Version]
- Forni, G.; Mantovani, A.; Forni, G.; Mantovani, A.; Moretta, L.; Rappuoli, R.; Rezza, G.; Bagnasco, A.; Barsacchi, G.; Bussolati, G.; et al. COVID-19 vaccines: Where we stand and challenges ahead. Cell Death Differ. 2021, 28, 626–639. [Google Scholar] [CrossRef]
- Yan, Z.-P.; Yang, M.; Lai, C.-L. COVID-19 Vaccines: A Review of the Safety and Efficacy of Current Clinical Trials. Pharmaceuticals 2021, 14, 406. [Google Scholar] [CrossRef]
- Wu, Q.; Dudley, M.Z.; Chen, X.; Bai, X.; Dong, K.; Zhuang, T.; Salmon, D.; Yu, H. Evaluation of the safety profile of COVID-19 vaccines: A rapid review. BMC Med. 2021, 19, 173. [Google Scholar] [CrossRef]
- Francis, A.I.; Ghany, S.; Gilkes, T.; Umakanthan, S. Review of COVID-19 vaccine subtypes, efficacy and geographical distributions. Postgrad. Med. J. 2021. ahead-of-print. [Google Scholar] [CrossRef] [PubMed]
- Smit, M.; Marinosci, A.; Agoritsas, T.; Calmy, A. Prophylaxis for COVID-19: A systematic review. Clin. Microbiol. Infect. 2021, 27, 532–537. [Google Scholar] [CrossRef] [PubMed]
- Ilie, O.-D.; Ciobica, A.; Doroftei, B. Testing the Accuracy of the ARIMA Models in Forecasting the Spreading of COVID-19 and the Associated Mortality Rate. Medicina 2020, 56, 566. [Google Scholar] [CrossRef] [PubMed]
- Ilie, O.-D.; Cojocariu, R.-O.; Ciobica, A.; Timofte, S.-I.; Mavroudis, I.; Doroftei, B. Forecasting the spreading of COVID-19 across nine countries from Europe, Asia, and the American continents using the arima models. Microorganisms 2020, 8, 1158. [Google Scholar] [CrossRef] [PubMed]
- Elevli, S.; Uzgören, N.; Bingöl, D.; Elevli, B. Drinking water quality control: Control charts for turbidity and pH. J. Water Sanit. Hyg. Dev. 2016, 6, 511–518. [Google Scholar] [CrossRef]
- He, Z.; Tao, H. Epidemiology and ARIMA model of positive-rate of influenza viruses among children in Wuhan, China: A nine-year retrospective study. Int. J. Infect. Dis. 2018, 74, 61–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Chen, P.; Niu, A.; Liu, D.; Jiang, W.; Ma, B. Time Series Forecasting of Temperatures using SARIMA: An Example from Nanjing. IOP Conf. Ser. Mater. Sci. Eng. 2018, 394, 52024. [Google Scholar] [CrossRef]
- Cihan, P. Forecasting fully vaccinated people against COVID-19 and examining future vaccination rate for herd immunity in the US, Asia, Europe, Africa, South America, and the World. Appl. Soft Comput. 2021, 111, 107708. [Google Scholar] [CrossRef]
- Sycinska-Dziarnowska, M.; Paradowska-Stankiewicz, I.; Woźniak, K. The Global Interest in Vaccines and Its Prediction and Perspectives in the Era of COVID-19. Real-Time Surveillance Using Google Trends. Int. J. Environ. Res. Public Health 2021, 18, 7841. [Google Scholar] [CrossRef]
- Chelo, D.; Nguefack, F.; Enyama, D.; Nansseu, R.; Feudjo Tefoueyet, G.; Mbassi Awa, H.D.; Mekone Nkwelle, I.; Nguefack-Tsague, G.; Ndenbe, P.; Koki Ndombo, P.O. Impact and projections of the COVID-19 epidemic on attendance and routine vaccinations at a pediatric referral hospital in Cameroon. Arch. Pédiatr. 2021, 28, 441–450. [Google Scholar] [CrossRef]
- Lumbreras-Marquez, M.I.; Fields, K.G.; Campos-Zamora, M.; Rodriguez-Bosch, M.R.; Rodriguez-Sibaja, M.J.; Copado-Mendoza, D.Y.; Acevedo-Gallegos, S.; Farber, M.K. A forecast of maternal deaths with and without vaccination of pregnant women against COVID-19 in Mexico. Int. J. Gynecol. Obstet. 2021, 154, 566–567. [Google Scholar] [CrossRef] [PubMed]
- Hwang, E. Prediction intervals of the COVID-19 cases by HAR models with growth rates and vaccination rates in top eight affected countries: Bootstrap improvement. Chaos Solitons Fractals 2022, 155, 111789. [Google Scholar] [CrossRef] [PubMed]
- Chaturvedi, D.; Chakravarty, U. Predictive analysis of COVID-19 eradication with vaccination in India, Brazil, and U.S.A. Infect. Genet. Evol. 2021, 92, 104834. [Google Scholar] [CrossRef]
- Kermack, W.O.; McKendrick, A.G.; Walker, G.T. A contribution to the mathematical theory of epidemics. Proc. R. Soc. London. Ser. A Contain. Pap. Math. Phys. Character 1927, 115, 700–721. [Google Scholar] [CrossRef] [Green Version]
- Zhang, J.; Wang, D.; Yang, K.; Xu, Y. A multinomial autoregressive model for finite-range time series of counts. J. Stat. Plan. Inference 2020, 207, 320–343. [Google Scholar] [CrossRef]
- Bartolucci, F.; Pennoni, F.; Mira, A. A multivariate statistical approach to predict COVID-19 count data with epidemiological interpretation and uncertainty quantification. Stat. Med. 2021, 40, 5351–5372. [Google Scholar] [CrossRef]
Figure 1.
Graphic revealing trends for each dose depending on the manufacturer among Romanian individuals.
Figure 1.
Graphic revealing trends for each dose depending on the manufacturer among Romanian individuals.
Figure 2.
Graphic showing the total number of Romanian individuals who completed the vaccination scheme.
Figure 2.
Graphic showing the total number of Romanian individuals who completed the vaccination scheme.
Figure 3.
Graphic highlighting the total number of adverse reactions (local vs. general) reported in Romanian individuals.
Figure 3.
Graphic highlighting the total number of adverse reactions (local vs. general) reported in Romanian individuals.
Figure 4.
The estimated ACF and PACF graphs used to predict the vaccination trend against COVID-19 in Romania per each month and cumulative. In this figure (left and right) are displayed the associated plots for the estimated partial and autocorrelations between residuals at distinct lags. Specifically, the lag x, partial and autocorrelation, coefficients evaluate the affinity between the residuals at a time t and time (x − t) (for autocorrelation)/(x + t) (partial autocorrelations) at 95.0% probability to be close to 0. Distinct from autocorrelation is the condition that t + x accounts for the correlations at all lower lags, observation used to appreciate the order of autoregressive model where needed to fit the data. Valid for both functions if the probability at a specific for autocorrelation/particular for partial autocorrelation lag do not contain the estimated coefficient, indeed exists a statistically significant correlation at that lag at CI 95.0%.
Figure 4.
The estimated ACF and PACF graphs used to predict the vaccination trend against COVID-19 in Romania per each month and cumulative. In this figure (left and right) are displayed the associated plots for the estimated partial and autocorrelations between residuals at distinct lags. Specifically, the lag x, partial and autocorrelation, coefficients evaluate the affinity between the residuals at a time t and time (x − t) (for autocorrelation)/(x + t) (partial autocorrelations) at 95.0% probability to be close to 0. Distinct from autocorrelation is the condition that t + x accounts for the correlations at all lower lags, observation used to appreciate the order of autoregressive model where needed to fit the data. Valid for both functions if the probability at a specific for autocorrelation/particular for partial autocorrelation lag do not contain the estimated coefficient, indeed exists a statistically significant correlation at that lag at CI 95.0%.
Figure 5.
Forecast plots for ARIMA models in the next week per month and cumulative. In Figure 5 (left and right) are displayed the plots of the observed and forecasted maximized values for all twelve months and in total at 95.0% prediction limit for each forecast. Limits presented indicate the true value of each month and total at any point in future with likely to be with 95.0% confidence.
Figure 5.
Forecast plots for ARIMA models in the next week per month and cumulative. In Figure 5 (left and right) are displayed the plots of the observed and forecasted maximized values for all twelve months and in total at 95.0% prediction limit for each forecast. Limits presented indicate the true value of each month and total at any point in future with likely to be with 95.0% confidence.
Table 1.
ARIMA models comparison.
| Month | Model | RMSE | MAE | MAPE |
|---|---|---|---|---|
| January | ARIMA(1,1,0) | 5983.61 | 3469.22 | 6.81315 |
| ARIMA(2,1,1) | 6193.02 | 3478.39 | 6.81429 | |
| ARIMA(2,1,0) | 6074.41 | 3529.04 | 6.8465 | |
| ARIMA(0,2,0) | 5993.82 | 3589.63 | 6.98607 | |
| February | ARIMA(1,1,1) | 2443.14 | 1496.05 | 0.214755 |
| ARIMA(1,1,2) | 2404.11 | 1496.8 | 0.212518 | |
| ARIMA(2,1,2) | 2449.65 | 1497.25 | 0.212571 | |
| ARIMA(2,1,0) | 2366.97 | 1498.86 | 0.212441 | |
| ARIMA(2,1,1) | 2409.78 | 1500.66 | 0.212628 | |
| March | ARIMA(0,2,0) | 10,125.3 | 6020.7 | 0.474466 |
| ARIMA(1,1,0) | 10,069.7 | 6038.5 | 0.479113 | |
| ARIMA(1,1,1) | 9725.72 | 6091.12 | 0.500062 | |
| ARIMA(1,2,0) | 9870.77 | 6138.27 | 0.490985 | |
| ARIMA(2,1,0) | 9901.43 | 6183.01 | 0.50242 | |
| April | ARIMA(2,1,0) | 4506.49 | 2856.4 | 0.117225 |
| ARIMA(1,1,0) | 4429.84 | 2860.85 | 0.117355 | |
| ARIMA(1,1,1) | 4506.19 | 2864.04 | 0.117446 | |
| ARIMA(2,1,1) | 4582.38 | 2872.9 | 0.117986 | |
| ARIMA(2,1,2) | 4645.91 | 2896.7 | 0.119043 | |
| May | ARIMA(1,2,2) | 5701.78 | 3844.24 | 0.105518 |
| ARIMA(2,2,2) | 5817.73 | 3852.44 | 0.105606 | |
| ARIMA(2,2,0) | 6058.59 | 4005.54 | 0.110237 | |
| June | ARIMA(2,2,2) | 2058.41 | 1593.3 | 0.0354902 |
| ARIMA(2,2,1) | 2268.92 | 1694.63 | 0.0377449 | |
| ARIMA(1,2,2) | 2406.98 | 1741.78 | 0.0388305 | |
| ARIMA(0,2,2) | 2401.52 | 1784.57 | 0.0397096 | |
| July | ARIMA(0,2,2) | 1610.49 | 1222.21 | 0.0251916 |
| ARIMA(2,2,2) | 1633 | 1233.64 | 0.0254273 | |
| ARIMA(1,2,2) | 1634.18 | 1236.74 | 0.0254953 | |
| ARIMA(2,1,2) | 1986.95 | 1476.09 | 0.0304824 | |
| ARIMA(1,1,2) | 1949.75 | 1497.82 | 0.0309203 | |
| August | ARIMA(2,2,2) | 1760.73 | 1196.04 | 0.0234608 |
| ARIMA(1,2,2) | 1750.49 | 1217.35 | 0.023894 | |
| ARIMA(0,2,2) | 1719.34 | 1220.36 | 0.0239556 | |
| ARIMA(1,2,1) | 2142.33 | 1372.74 | 0.0269142 | |
| ARIMA(2,1,2) | 2216.94 | 1594.14 | 0.0313006 | |
| September | ARIMA(1,1,2) | 2641.24 | 1891.2 | 0.0352714 |
| ARIMA(2,1,2) | 2695.21 | 1929.22 | 0.0359849 | |
| ARIMA(2,2,2) | 2674.02 | 1981.65 | 0.0369697 | |
| ARIMA(0,2,2) | 2746.87 | 2089.46 | 0.0389583 | |
| ARIMA(1,2,2) | 2747.8 | 2096.01 | 0.0390938 | |
| October | ARIMA(2,2,2) | 6287.29 | 4643.42 | 0.0774352 |
| ARIMA(2,2,1) | 6897.6 | 5108.31 | 0.0847775 | |
| ARIMA(1,2,2) | 7032.4 | 5193.96 | 0.085918 | |
| ARIMA(1,2,1) | 6908.67 | 5262.54 | 0.0870042 | |
| November | ARIMA(2,1,1) | 6121.27 | 4144.73 | 0.0570842 |
| ARIMA(2,1,0) | 5986.3 | 4187.46 | 0.0577338 | |
| ARIMA(2,1,2) | 6207.1 | 4200.65 | 0.0579711 | |
| December | ARIMA(2,2,1) | 1891.64 | 1470.03 | 0.0187244 |
| ARIMA(2,2,0) | 1856.3 | 1473.05 | 0.0187625 | |
| ARIMA(2,2,2) | 1916.62 | 1499.57 | 0.0191095 | |
| ARIMA(0,2,2) | 2016.04 | 1643.22 | 0.0209753 | |
| Total | ARIMA(2,0,2) | 5360.97 | 3259.29 | 0.733775 |
| ARIMA(2,2,0) | 5394.91 | 3268.17 | 0.695145 | |
| ARIMA(2,2,1) | 5399.27 | 3268.9 | 0.695924 | |
| ARIMA(2,2,2) | 5406.99 | 3271.66 | 0.696387 | |
| ARIMA(1,1,1) | 5393.81 | 3274.41 | 0.688878 |
Table 2.
ARIMA models parameters.
| Month | Parameter | Estimate | Standard Error | t-Statistic | p-Value | Ljung–Box Test |
|---|---|---|---|---|---|---|
| January | AR(1) | 0.982382 | 0.0548607 | 17.9069 | 0 | 0.102632 |
| February | AR(1) | 0.9564 | 0.0333551 | 28.6733 | 0 | 0.864548 |
| MA(1) | −0.168594 | 0.118083 | −1.42776 | 0.164042 | ||
| March | no parameter(s) | 0.477973 | ||||
| April | AR(1) | 1.03694 | 0.18565 | 5.58548 | 0.000005 | 0.248501 |
| AR(2) | −0.0224815 | 0.189402 | −0.118698 | 0.906333 | ||
| May | AR(1) | 0.759005 | 0.158601 | 4.78564 | 0.000059 | 0.986002 |
| MA(1) | 0.431691 | 0.165075 | 2.61512 | 0.01465 | ||
| MA(2) | 0.575053 | 0.156229 | 3.68084 | 0.001069 | ||
| June | AR(1) | 1.18927 | 0.0786618 | 15.1188 | 0 | 0.0169788 |
| AR(2) | −0.975677 | 0.0728746 | −13.3884 | 0 | ||
| MA(1) | 1.26288 | 0.159909 | 7.89746 | 0 | ||
| MA(2) | −0.808264 | 0.123585 | −6.54012 | 0.000001 | ||
| July | MA(1) | 0.0470354 | 0.0932973 | 0.504145 | 0.618249 | 0.0043751 |
| MA(2) | 0.868584 | 0.0876427 | 9.91051 | 0 | ||
| August | AR(1) | 0.0639607 | 0.207405 | 0.308385 | 0.760246 | 0.000105 |
| AR(2) | −0.186658 | 0.196426 | −0.950269 | 0.350726 | ||
| MA(1) | 0.0050157 | 0.0645373 | 0.0777177 | 0.938648 | ||
| MA(2) | 0.95 | 0.053573 | 17.7328 | 0 | ||
| September | AR(1) | 1.03825 | 0.0200273 | 51.8419 | 0 | 0.0136246 |
| MA(1) | 0.29945 | 0.190727 | 1.57004 | 0.127639 | ||
| MA(2) | 0.371627 | 0.175248 | 2.12057 | 0.042949 | ||
| October | AR(1) | −0.211891 | 0.182648 | −1.16011 | 0.256963 | 0.690025 |
| AR(2) | −0.702692 | 0.168935 | −4.15955 | 0.000329 | ||
| MA(1) | −1.00778 | 0.157958 | −6.38005 | 0.000001 | ||
| MA(2) | −0.780661 | 0.166924 | −4.67675 | 0.000086 | ||
| November | AR(1) | 1.49501 | 0.212479 | 7.03601 | 0 | 0.0059234 |
| AR(2) | −0.530011 | 0.211406 | −2.50708 | 0.018252 | ||
| MA(1) | 0.267886 | 0.318328 | 0.841539 | 0.407176 | ||
| December | AR(1) | −0.0584976 | 0.202084 | −0.289471 | 0.774516 | 0.175271 |
| AR(2) | −0.745858 | 0.135203 | −5.51656 | 0.000009 | ||
| MA(1) | −0.0445405 | 0.265186 | −0.167959 | 0.867915 | ||
| Total | AR(1) | 1.97226 | 0.0144948 | 136.067 | 0 | 1.11 × 10−16 |
| AR(2) | −0.972185 | 0.0145503 | −66.8156 | 0 | ||
| MA(1) | −0.147522 | 0.0546424 | −2.69977 | 0.007264 | ||
| MA(2) | 0.103778 | 0.0542596 | 1.91261 | 0.056587 | ||
Table 3.
Prediction of vaccinated individuals against COVID-19 per month and total for the next week according to our ARIMA with CI95%.
Table 3.
Prediction of vaccinated individuals against COVID-19 per month and total for the next week according to our ARIMA with CI95%.
| Lower 95% | Upper 95% | ||
|---|---|---|---|
| Period | Forecast | Limit | Limit |
| January | |||
| 28 January 2021 | 538,694 | 526,474 | 550,914 |
| 29 January 2021 | 563,514 | 536,381 | 590,647 |
| 30 January 2021 | 587,896 | 542,802 | 632,991 |
| 31 January 2021 | 611,849 | 546,277 | 677,421 |
| 1 February 2021 | 635,380 | 547,182 | 723,579 |
| 2 February 2021 | 658,497 | 545,792 | 771,201 |
| 3 February 2021 | 681,206 | 542,325 | 820,086 |
| February | |||
| 28 February 2021 | 919,209 | 914,179 | 924,239 |
| 1 March 2021 | 932,663 | 920,849 | 944,477 |
| 2 March 2021 | 945,530 | 925,560 | 965,501 |
| 3 March 2021 | 957,836 | 928,655 | 987,018 |
| 4 March 2021 | 969,606 | 930,370 | 1.00884 × 106 |
| 5 March 2021 | 980,862 | 930,880 | 1.03084 × 106 |
| 6 March 2021 | 991,628 | 930,327 | 1.05293 × 106 |
| March | |||
| 28 March 2021 | 1.94174 × 106 | 1.92096 × 106 | 1.96251 × 106 |
| 29 March 2021 | 1.96859 × 106 | 1.92213 × 106 | 2.01504 × 106 |
| 30 March 2021 | 1.99544 × 106 | 1.91771 × 106 | 2.07317 × 106 |
| 31 March 2021 | 2.02229 × 106 | 1.9085 × 106 | 2.13608 × 106 |
| 1 April 2021 | 2.04914 × 106 | 1.89507 × 106 | 2.20322 × 106 |
| 2 April 2021 | 2.076 × 106 | 1.87781 × 106 | 2.27418 × 106 |
| 3 April 2021 | 2.10285 × 106 | 1.85703 × 106 | 2.34867 × 106 |
| April | |||
| 28 April 2021 | 3.20358 × 106 | 3.19435 × 106 | 3.2128 × 106 |
| 29 April 2021 | 3.26041 × 106 | 3.23948 × 106 | 3.28134 × 106 |
| 30 April 2021 | 3.31808 × 106 | 3.28272 × 106 | 3.35345 × 106 |
| 1 May 2021 | 3.37661 × 106 | 3.32444 × 106 | 3.42878 × 106 |
| 2 May 2021 | 3.436 × 106 | 3.36486 × 106 | 3.50714 × 106 |
| 3 May 2021 | 3.49627 × 106 | 3.40416 × 106 | 3.58838 × 106 |
| 4 May 2021 | 3.55743 × 106 | 3.44246 × 106 | 3.6724 × 106 |
| May | |||
| 28 May 2021 | 4.25759 × 106 | 4.24549 × 106 | 4.26969 × 106 |
| 29 May 2021 | 4.28472 × 106 | 4.25408 × 106 | 4.31536 × 106 |
| 30 May 2021 | 4.31365 × 106 | 4.26305 × 106 | 4.36424 × 106 |
| 31 May 2021 | 4.34394 × 106 | 4.27326 × 106 | 4.41462 × 106 |
| 1 June 2021 | 4.37527 × 106 | 4.28497 × 106 | 4.46557 × 106 |
| 2 June 2021 | 4.40739 × 106 | 4.29825 × 106 | 4.51653 × 106 |
| 3 June 2021 | 4.44011 × 106 | 4.31303 × 106 | 4.56718 × 106 |
| June | |||
| 28 June 2021 | 4.70459 × 106 | 4.70025 × 106 | 4.70894 × 106 |
| 29 June 2021 | 4.71453 × 106 | 4.70509 × 106 | 4.72396 × 106 |
| 30 June 2021 | 4.72637 × 106 | 4.71166 × 106 | 4.74108 × 106 |
| 1 July 2021 | 4.73785 × 106 | 4.71808 × 106 | 4.75761 × 106 |
| 2 July 2021 | 4.74701 × 106 | 4.7222 × 106 | 4.77182 × 106 |
| 3 July 2021 | 4.75379 × 106 | 4.72335 × 106 | 4.78424 × 106 |
| 4 July 2021 | 4.75999 × 106 | 4.72278 × 106 | 4.79721 × 106 |
| July | |||
| 28 July 2021 | 4.97002 × 106 | 4.9667 × 106 | 4.97334 × 106 |
| 29 July 2021 | 4.97892 × 106 | 4.97164 × 106 | 4.98621 × 106 |
| 30 July 2021 | 4.98783 × 106 | 4.97788 × 106 | 4.99777 × 106 |
| 31 July 2021 | 4.99673 × 106 | 4.98455 × 106 | 5.00892 × 106 |
| 1 August 2021 | 5.00564 × 106 | 4.99142 × 106 | 5.01986 × 106 |
| 2 August 2021 | 5.01454 × 106 | 4.99842 × 106 | 5.03067 × 106 |
| 3 August 2021 | 5.02345 × 106 | 5.0055 × 106 | 5.0414 × 106 |
| August | |||
| 28 August 2021 | 5.23552 × 106 | 5.2318 × 106 | 5.23923 × 106 |
| 29 August 2021 | 5.24423 × 106 | 5.23572 × 106 | 5.25274 × 106 |
| 30 August 2021 | 5.25285 × 106 | 5.24159 × 106 | 5.26412 × 106 |
| 31 August 2021 | 5.26164 × 106 | 5.24849 × 106 | 5.2748 × 106 |
| 1 September 2021 | 5.27046 × 106 | 5.25557 × 106 | 5.28535 × 106 |
| 2 September 2021 | 5.27925 × 106 | 5.26269 × 106 | 5.29582 × 106 |
| 3 September 2021 | 5.28803 × 106 | 5.26988 × 106 | 5.30618 × 106 |
| September | |||
| 28 September 2021 | 5.49969 × 106 | 5.49415 × 106 | 5.50524 × 106 |
| 29 September 2021 | 5.51582 × 106 | 5.5047 × 106 | 5.52695 × 106 |
| 30 September 2021 | 5.53257 × 106 | 5.51633 × 106 | 5.54881 × 106 |
| 1 October 2021 | 5.54996 × 106 | 5.52844 × 106 | 5.57148 × 106 |
| 2 October 2021 | 5.56801 × 106 | 5.54092 × 106 | 5.59511 × 106 |
| 3 October 2021 | 5.58676 × 106 | 5.5537 × 106 | 5.61982 × 106 |
| 4 October 2021 | 5.60622 × 106 | 5.5668 × 106 | 5.64564 × 106 |
| October | |||
| 28 October 2021 | 6.8339 × 106 | 6.82067 × 106 | 6.84713 × 106 |
| 29 October 2021 | 6.93473 × 106 | 6.89546 × 106 | 6.97401 × 106 |
| 30 October 2021 | 7.04299 × 106 | 6.97166 × 106 | 7.11432 × 106 |
| 31 October 2021 | 7.14982 × 106 | 7.04635 × 106 | 7.25329 × 106 |
| 1 November 2021 | 7.25173 × 106 | 7.1128 × 106 | 7.39067 × 106 |
| 2 November 2021 | 7.35569 × 106 | 7.17533 × 106 | 7.53606 × 106 |
| 3 November 2021 | 7.46267 × 106 | 7.23757 × 106 | 7.68778 × 106 |
| November | |||
| 28 November 2021 | 7.72849 × 106 | 7.71583 × 106 | 7.74114 × 106 |
| 29 November 2021 | 7.73911 × 106 | 7.70822 × 106 | 7.77 × 106 |
| 30 November 2021 | 7.74866 × 106 | 7.69434 × 106 | 7.80297 × 106 |
| 1 December 2021 | 7.7573 × 106 | 7.67552 × 106 | 7.83907 × 106 |
| 2 December 2021 | 7.76515 × 106 | 7.65286 × 106 | 7.87745 × 106 |
| 3 December 2021 | 7.77232 × 106 | 7.6272 × 106 | 7.91745 × 106 |
| 4 December 2021 | 7.77887 × 106 | 7.59921 × 106 | 7.95853 × 106 |
| December | |||
| 28 December 2021 | 7.92705 × 106 | 7.92313 × 106 | 7.93097 × 106 |
| 29 December 2021 | 7.92826 × 106 | 7.91954 × 106 | 7.93697 × 106 |
| 30 December 2021 | 7.92976 × 106 | 7.91742 × 106 | 7.9421 × 106 |
| 31 December 2021 | 7.93335 × 106 | 7.91754 × 106 | 7.94916 × 106 |
| 1 January 2022 | 7.9366 × 106 | 7.91599 × 106 | 7.9572 × 106 |
| 2 January 2022 | 7.93831 × 106 | 7.91207 × 106 | 7.96455 × 106 |
| 3 January 2022 | 7.94037 × 106 | 7.90872 × 106 | 7.97201 × 106 |
| Total | |||
| 28 December 2021 | 7.92819 × 106 | 7.91765 × 106 | 7.93873 × 106 |
| 29 December 2021 | 7.9335 × 106 | 7.90879 × 106 | 7.95822 × 106 |
| 30 December 2021 | 7.93926 × 106 | 7.89825 × 106 | 7.98028 × 106 |
| 31 December 2021 | 7.94546 × 106 | 7.88616 × 106 | 8.00477 × 106 |
| 1 January 2022 | 7.95208 × 106 | 7.87278 × 106 | 8.03139 × 106 |
| 2 January 2022 | 7.95912 × 106 | 7.85832 × 106 | 8.05991 × 106 |
| 3 January 2022 | 7.96656 × 106 | 7.84298 × 106 | 8.09014 × 106 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Doroftei, B.; Ilie, O.-D.; Anton, N.; Timofte, S.-I.; Ilea, C. Mathematical Modeling to Predict COVID-19 Infection and Vaccination Trends. J. Clin. Med. 2022, 11, 1737. https://doi.org/10.3390/jcm11061737
AMA Style
Doroftei B, Ilie O-D, Anton N, Timofte S-I, Ilea C. Mathematical Modeling to Predict COVID-19 Infection and Vaccination Trends. Journal of Clinical Medicine. 2022; 11(6):1737. https://doi.org/10.3390/jcm11061737
Chicago/Turabian StyleDoroftei, Bogdan, Ovidiu-Dumitru Ilie, Nicoleta Anton, Sergiu-Ioan Timofte, and Ciprian Ilea. 2022. "Mathematical Modeling to Predict COVID-19 Infection and Vaccination Trends" Journal of Clinical Medicine 11, no. 6: 1737. https://doi.org/10.3390/jcm11061737
APA StyleDoroftei, B., Ilie, O. -D., Anton, N., Timofte, S. -I., & Ilea, C. (2022). Mathematical Modeling to Predict COVID-19 Infection and Vaccination Trends. Journal of Clinical Medicine, 11(6), 1737. https://doi.org/10.3390/jcm11061737
Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.
