DATA ANALYTICS AND SIR MODELING
OF COVID-19 IN BULGARIA
Miroslava Ivanova1, Lilko Dospatliev2 1Trakia University
Department of Informatics and Mathematics
Stara Zagora, 6000, BULGARIA 2Trakia University
Department of Pharmacology, Animal Physiology
and Physiological Chemistry
Stara Zagora, 6000, BULGARIA
The Novel Human Coronavirus (SARS-CoV-2) also well-known as
COVID-19 is the greatest public health challenge of the 21st
century. The aims of this article were to provide data analytic of
the COVID-19 cases in Bulgaria and to modeling their spread using
SIR (Susceptible-Infected-Recovered) model. We used the
covid19.analytics package from the statistical software R in
extracting time-series data of COVID-19 in Bulgaria (dating March
8, 2020 to November 24, 2020), in presenting the data analysis as
well as forecast the maximum number of infected people, peak time
and the basic reproduction number based on SIR model. As per SIR
model, the maximum number of infected people is reached after 61
days of starting of COVID-19 pandemic in Bulgaria (around May 13,
2020) with about 1.298446e^{+05} infections. The basic reproduction number R0 was found to 1.46, which means that on
average an infectious individual infects 1.46 susceptible
individuals during his infection period. We believe that performed
data analytics of COVID-19 cases in Bulgaria and the obtained
results of the SIR model will help government of Bulgaria when
restricting the spread of the virus.
You will need Adobe Acrobat reader. For more information and free download of the reader, please follow this link.
References
[1] L. Abadias, G. Estrada-Rodriguez and E. Estrada, Fractional-order
Susceptible-Infected model: definition and applications to the study of
COVID-19 main protease, Fractional Calculus and Applied Analysis, 23,
No 3 (2020), 635-655; doi: 10.1515/fca-2020-0033.
[2] A.M. Almeshal, A.I. Almazrouee, M.R. Alenizi et al., Forecasting the
spread of COVID-19 in Kuwait using compartmental and logistic regression
models, Applied Sciences, 10, No 10 (2020), # 3402; doi:
10.3390/app10103402.
[3] T.M. Alon, M. Doepke, J. Olmstead-Rumsey et al., The impact of COVID-
19 on gender equality, National Bureau of Economic Research (2020); doi:
10.3386/w26947.
[4] N. Anand, A. Sabarinath, S. Geetha et al., Predicting the spread of
COVID-19 using SIR model augmented to incorporate quarantine and testing,
Trans. Indian National Academy of Engin., 5 (2020), 141-148.
[5] L.M. Bettencourt and R.M. Ribeiro, Real time Bayesian estimation of the
epidemic potential of emerging infectious diseases, PLoS One, 3, No 5
(2008), Art. e2185.
[6] M. Blofield, B. Hoffmann and M. Llanos, Assessing the political
and social impact of the COVID-19 crisis in Latin America,
GIGA Focus Latein Amerika, 3 (2020); Available from URL: https:
//www.ssoar.info/ssoar/bitstream/handle/document/67260/ssoar-
2020-blofield_et_al-Assessing_the_Political_and_Social.pdf?
sequence=1&isAllowed=y&lnkname=ssoar-2020-blofield_et_al-
Assessing_the_Political_and_Social.pdf.
[7] Bulgaria, Ministry for Health, Order No RD-01-114/05.03.2020; Available
from URL: https://www.mh.government.bg/media/filer_public/
2020/03/05/rd-01-114_u7ofbPI.pdf (in Bulgarian).
[8] Bulgaria, Ministry for Health. Order No RD-01-183/06.04.2020; Available
from URL: https://www.mh.government.bg/media/filer_public/
2020/04/07/rd-01-183.pdf (in Bulgarian).
DATA ANALYTICS AND SIR MODELING... 1111
[9] Bulgaria, Moldova Report First Coronavirus Cases, Radio Free
Europe/Radio Liberty, 8 March 2020; Available from URL:
https://www.rferl.org/a/bulgaria-confirms-coronavirus-casesamid-
flu-outbreak/30475539.html.
[10] Y. Fang, Y. Nie and M. Penny, Transmission dynamics of the COVID-
19 outbreak and effectiveness of government interventions: A data-driven
analysis, J. of Medical Virology, 92, No 6 (2020), 645-659.
[11] R.T. Gandhi, J.B. Lynch and C. del Rio, Mild or moderate COVID-19,
The New England J. of Medicine, 383 (2020), 1757-1766.
[12] T. Harko, F.S.N. Lobo and M.K. Mak, Exact analytical solutions of the
Susceptible-Infected- Recovered (SIR) epidemic model and of the SIR
model with equal deaths and births, Applied Mathematics and Compu-
tation, 236 (2014), 184-194.
[13] https://botevgrad.com/news/78275/Prof-Todor-Kantardjiev-Ot-
Pleven-imame-dve-polojitelni-probi-za-koronavirus-ot-
Gabrovo-sashto/ (in Bulgarian).
[20] C. Huang, Y. Wang, X. Li et al., Clinical features of patients infected with
2019 novel coronavirus in Wuhan, China, Lancet, 395, No 10223 (2020),
497-506.
[21] N. Hussain and L. Baoming, Using R-studio to examine the COVID-19
patients in Pakistan implementation of SIR Model on cases, Intern. J. of
Scientific Research in Multidisciplinary Studies, 6, No 8 (2020), 54-59.
1112 M. Ivanova, L. Dospatliev
[22] J.H. Jones, Notes on R0. Califonia: Department of Anthropological Sciences.
Department of Anthropological Sciences Stanford University, (2007),
Available from URL: http://people.math.gatech.edu/~dmargalit7/
classes/math1553Spring2020/JonesR0.pdf.
[23] M. Keeling and P. Rohani, Modeling Infectious Diseases in Humans and
Animals, Princeton University Press, Princeton, NJ, USA (2008).
[24] W.O. Kermack and A.G. McKendrick, A contribution to the mathematical
theory of epidemics, Proc. Royal Soc. of London, Ser. A, 115 (1927), 700-
721.
[25] O. Kounchev, G. Simeonov and Z. Kuncheva, The TVBG-SEIR spline
model for analysis of COVID-19 spread, and a Tool for prediction scenarios,
arXiv Preprint: arXiv (2020), 2004.11338, 1-21.
[26] Y. Liu, A.A. Gayle, A. Wilder-Smith and J. Rockl¨ov, The reproductive
number of COVID-19 is higher compared to SARS coronavirus, J. of Travel
Medicine, 27, No 2 (2020), 1-4.
[27] B. Malavika, S. Marimuthu, M. Joy et al., Forecasting COVID-19
epidemic in India and high incidence states using SIR and logistic
growth models, Clinical Epidemiology and Global Health (2020); doi:
10.1016/j.cegh.2020.06.006.
[28] S. Margenov, N. Popivanov, I. Ugrinova et al., Mathematical and computer
modeling of COVID-19 transmission dynamics in Bulgaria by timedepended
inverse SEIR model, arXiv Preprint: arXiv (2020), 2008.10360,
1-16.
[29] T. Marinov and R. Marinova, Dynamics of COVID-19 using inverse problem
for coefficient identification in SIR epidemic models, Chaos, Solitons
& Fractals: X, 5 (2020); doi: 10.1016/j.csfx.2020.100041.
[30] W.J. McKibbin and R. Fernando, The global macroeconomic impacts of
COVID-19: Seven scenarios, CAMA (2020); doi: 10.2139/ssrn.3547729.
[31] K.R. Moran, G. Fairchild, N. Generous et al., Epidemic forecasting is
messier than weather forecasting: the role of human behavior and internet
data streams in epidemic forecast, J. Infectious Diseases, 214, No 4
(2016), S404-S408.
DATA ANALYTICS AND SIR MODELING... 1113
[32] U. Nguemdjo, F. Meno, A. Dongfack et al., Simulating the progression of
the COVID-19 disease in Cameroon using SIR models, Plos One (2020);
doi: 10.1371/journal.pone.0237832.
[33] M. Ponce and A. Sandhel, covid19.analytics: An R Package to Obtain,
Analyze and Visualize Data from the Corona Virus Disease Pandemic,
arXiv Preprint: arXiv (2020); 2009.01091, 1-47.
[34] A. Rajesh, H. Pai, V. Roy et al., CoVID-19 prediction for India from the
existing data and SIR(D) model study, arXiv Preprint: MedRxiv (2020);
doi: 10.1101/2020.05.05.20085902.
[35] W.C. Roda, M.B. Varughese, D. Han et al., Why is it difficult to accurately
predict the COVID-19 epidemic? Infectious Disease Modelling, 5 (2020),
271-281.
[36] P. Sun and K. Li, An SEIR model for assessment of current COVID-
19 pandemic situation in the UK, arXiv Preprint: MedRxiv (2020); doi:
10.1101/2020.04.12.20062588.
[37] A. Tob´ıas, Evaluation of the lockdowns for the SARS-CoV-2 epidemic in
Italy and Spain after one month follow up, Science of the Total Environ-
ment, 725, No 138539 (2020); doi: 10.1016/j.scitotenv.2020.138539.
[38] J. Tolles and T.B. Luong, Modeling epidemics with compartmental models,
JAMA Guide to Statistics and Methods, 323, No 24 (2020), 2515-2516.
[39] D. Wang, B. Hu, C. Hu et al., Clinical characteristics of 138 hospitalized
patients with 2019 novel coronavirus-infected pneumonia inWuhan, China,
JAMA, 323, No 11 (2020), 1061-1069.
[40] J. Wangping, H. Ke, S. Yang et al., Extended SIR prediction of the epidemics
trend of COVID-19 in Italy and compared with Hunan, China,
arXiv preprint: medRxiv (2020); doi: 10.1101/2020.03.18.20038570.
[41] WHO, World Health Organization Coronavirus disease (COVID-
19) advice for the public, (2020); Available from URL: https:
//www.who.int/emergencies/diseases/novel-coronavirus-
2019/advice-for-public.
[42] WHO, World Health Organization Managing epidemics: Key facts about
major deadly diseases, (2018); Available from URL: https://www.who.
int/emergencies/diseases/managingepidemics-interactive.pdf.
1114 M. Ivanova, L. Dospatliev
[43] WHO, World Health Organization Pneumonia of unknown cause -
China (2020); Available from URL: https://www.who.int/csr/don/05-
january-2020-pneumonia-of-unkown-cause-china/en/.
[44] S. Zhao, Q. Lin, J. Ran et al., Preliminary estimation of the basic reproduction
number of novel coronavirus (2019-nCoV) in China, from 2019 to
2020: A data-driven analysis in the early phase of the outbreak, Intern. J.
of Infectious Diseases, 92 (2020), 214-217.
[45] Z. Zhuang, S. Zhao, Q. Lin et al., Preliminary estimating the reproduction
number of the coronavirus disease (COVID-19) outbreak in Republic of
Korea and Italy by 5 March 2020, Intern. J. of Infectious Diseases, 95
(2020), 308-310.