1 Introduction

The novel coronavirus disease 2019 (COVID-19) is an emerging epidemic responsible for infecting lots of people and killing lakhs globally since the first report till today, according to the World Health Organization (WHO). The WHO China Country Office has informed of cases of pneumonia unknown aetiology detected in Wuhan, Hubei Province of China, on 31 December 2019. A total of 44 case patients with pneumonia of unknown aetiology have reported to WHO by the national authorities in China since 31 December 2019 through 3 January 2020. Further, WHO received the information from the National Health Commission China that the outbreak is associated with exposures in one seafood market in Wuhan on 11 and 12 January 2020. The Chinese authorities identified a new type of coronavirus, which has been isolated on 7 January 2020. China shared the genetic sequence of the novel coronavirus for countries to use in developing specific diagnostic kits. During the period 13–20 January 2020, the Ministry of Public Health, Thailand, reported the first imported case of laboratory-confirmed novel COVID-19 from Wuhan, China, and the Ministry of Health, Labour and Welfare, Japan (MHLW), reported an imported case of laboratory-confirmed COVID-19 from Wuhan, China. Also, National IHR Focal Point (NFP) for Republic of Korea reported the first case of novel coronavirus in the Republic of Korea. According to WHO, globally 5,213,483 cases confirmed and 338,225 total deaths until 23 May 2020 in which the USA holding the 1st position in the world ranking in the number of confirmed cases of COVID-19 as in total cases 1,555,537. Across the world, countries have implemented a lot of control procedures to respond to COVID-19, with the aim of reducing communication and dropping mortality. Particularly, some countries are in the process of measures public health and social processes based on local epidemiology, while others are considering shrink these interventions. Besides, countries have employed public health and social measures based on local risk assessments, such as full lockdown, restricted local movement, and educational institutes are implementing the virtual classes and industries are encouraging work from home, and regulated universal travel measures.

Numerous papers have appeared modelling the forecasting the trends of epidemiology of the disease [1,2,3,4,5], and many analyses have been done for the currently available epidemic data [6,7,8]. The analysis made by using practical data shows that the power law kinetics with fractal exponent provides a better fit to the current data for the number of deaths and spared rate compared than the classical epidemiological approach that assumes an exponential growth of the disease [9,10,11,12]. The idea of fractal dimension is useful in the categorization and quantification of shape and texture which is used to describe systems with simple self-similarity and regular fractals [13]. However, it is not enough to characterize sets having heterogeneous scaling properties like growth of epidemic disease. For such a characterization Grassberger has generalized the analysis with the notion of multifractals which implies a continuous spectrum of exponents for the characterization of the system [14]. In this oversimplification, an inhomogeneous fractal is considered to be associated with countably many subfractal sets of dissimilar dimensions [13,14,15,16,17,18]. We analyse the COVID-19 infection data from the top 15 affected countries in which we observed the growth and death patterns of the COVID-19 are heterogeneous. Hence, this paper applies the multifractal formalism on COVID-19 data with the notion that country-specific transmission rates follow a power law growth nature.

Herd immunity is related to a scenario under which people develop immunity against contagious infectious disease that materializes when the people is immune either through vaccination or immunity developed through former infection and they become immune to that disease. This provides indirect protection or herd immunity to those people who are not immune to the disease, see for instance Fig. 1.

Fig. 1
figure 1

The principle of herd immunity

As on 25 May 2020, India is under the top 10 most infected countries with over 6000 per day recorded cases for four consecutive days. Data from Johns Hopkins University [6] show that the total number of cases in India has surpassed that of Iran. After a sharp spikes of cases since the second week of May, it can be predicted that the scenario will continue till it can attain a maximum peak, may be in the month of June. After this intense period, herd immunity can occur by a natural course of infection. In this report, using a brief comparative analysis, we have shown that India can develop a significant mutation against the disease. We have restricted our data analysis till 22nd May.

The work is organized as: In Sect. 2, necessary materials and method of multifractal formalism are presented. In Sect. 3, the COVID-19 infection data from the top 15 affected countries are described and discussed the heterogeneous nature of COVID-19 growth patterns. Further, multifractal formalism is applied on COVID-19 data with the notion of power law growth behaviour of virus. The generalized fractal dimension curves are estimated, and using these results comparison is made between India with other countries. Conclusions of the work are provided in Sect. 4.

2 Material and method of multifractality

This section presents the explanation of the material and method description applied in this paper.

Prediction of the number of infected case for different diseases is the foremost task in the healthcare system. In general, epidemiological analysis models assume that the number of confirmed positive cases of COVID-19 grows exponentially with respect to fixed reproduction rate n. Each person infected by the COVID-19 transmits to n number of new persons when \(n>1\), and then the total number of COVID-19 cases grows as

$$\begin{aligned} n^{\frac{t}{\tau }}=e^{at} \end{aligned}$$
(1)

here \(\tau \) is the incubation time which depends on the nature of COVID-19. The limitation of exponential approaches is that it assumes each person transmits to n number of new person without considering the inhibition due to the interaction with already infected people such as lockdown, quarantine, restricted local movement, distance-learning, teleworking, regulated travel, and other prophylactic measures. Therefore, in order to create short-term predictions to prepare for the extent of the global pandemic infection data from the top 25 affected countries have investigated in [9]. Also, all country-specific infection rates follow a power law growth behaviour and the scaling exponents per country are calculated in [9]. In this paper, the partition function of COVID 19 data for top 15 affected countries has been estimated to apply analogy of multifractal formalism with COVID-19 transition, since the partition function exhibits the power law distribution.

A set S of points describing an object can be divided into boxes labelled by an index i such that the ith box has \(N_i\) points of the total N points of the set. These points are sample points describing the content of the underlying measure. Let us use the mass or probability \(\mu _i=N_i/N\) in the ith cell to construct weighted d-measure which we write

$$\begin{aligned} M_d(q,\delta )=\sum _{i=1}^{N}\mu _i^q\delta ^d=Z(q,\delta )\delta ^d \end{aligned}$$
(2)

The mass exponent \(\tau (q)\) for the set depends on the moment

$$\begin{aligned} Z(q,\delta )=\sum _i\mu _i^q, \end{aligned}$$
(3)

of order q chosen. Like \(N(\delta )\) of bare d-measure the partition function (if we may call) \(Z(q,\delta )\) also often exhibits power law distribution

$$\begin{aligned} Z(q,\delta )\sim \delta ^{-\tau (q)}, \end{aligned}$$
(4)

where the exponent \(\tau (q)\) is called the mass exponent not the fractal dimension if it is nonlinear in q. The mass exponent is more revealing than the simple fractal dimension as we shall soon find it out. It can equivalently be defined as

$$\begin{aligned} \tau (q)=-\lim _{\delta \rightarrow 0}{{\ln Z(q,\delta )}\over {\ln \delta }}. \end{aligned}$$
(5)

The generalized fractal dimension (GFD) \(D_q\) can be defined as

$$\begin{aligned} D_q=\frac{1}{q-1}\tau (q)=\lim _{\delta \rightarrow 0}\frac{1}{q-1}{{\ln Z(q,\delta )}\over {\ln \delta }}. \end{aligned}$$
(6)

As \(q \longrightarrow 1\), \(D_q\) converges to \(D_1\), which is given by

$$\begin{aligned} D_1= \lim _{\delta \rightarrow 0} \frac{\sum ^{N}_{i=1} \mu _i \ln \mu _i}{\ln \delta }, \end{aligned}$$
(7)

where \(D_1\) is the information dimension and \(D_q\) is a monotonically decreasing function of q such that \(D_0 \ge D_1 \ge D_2\). Here \(D_0\) and \(D_2\) denote the fractal dimension and correlation dimension, respectively.

The set S is a country and it can be divided into boxes labelled by an index i such that the ith box has \(N_i\) number of person affected by COVID-19 of the total number of persons infected by every positive case N of the country S. Let us use the mass \(\mu _i=N_i/N\) in the ith cell to estimate the generalized fractal dimension by using Eq. (6).

3 Results and discussion

COVID-19 data have been downloaded from the publicly available JHU-CSSE (2020) data repository provided by the Johns Hopkins University Center for Systems Science and Engineering (JHUCSSE) [6] and Statistics and Research Coronavirus Pandemic (COVID-19) available in the website of our world in data [7]. We have selected the data sets with the 15 most infections as presented in [6] on 22-05-2020 11:00 p.m. IST. The countries and their infection counts and death counts are given in Table 1.

Table 1 Countries with the most COVID-19 infections as on 22.05.2020
Fig. 2
figure 2

Generalized fractal dimensions of confirmed case and confirmed death of COVID-19

Fig. 3
figure 3

Daily confirmed case and confirmed death of COVID-19 and associated generalized fractal dimensions

Fig. 4
figure 4

Daily confirmed case and confirmed death of COVID-19 and associated generalized fractal dimensions

Fig. 5
figure 5

Daily confirmed case and confirmed death of COVID-19 and associated generalized fractal dimensions

Figure 2a shows generalized fractal dimension result for the confirmed cases of 15 countries. For each country day zero is considered to be the day of the first confirmed case of COVID-19. In each day the confirmed cases are normalized by dividing them by the total number of case and obtained the partition function \( Z(q,\delta )\) by Eq. (3). Also, Fig. 2b reveals generalized fractal dimensions of confirmed deaths corresponding to the confirmed case of top 15 countries.

In Fig. 2a, a comparison of the 80 days of 15 different countries’ case counts shows that India has ranked 7th spread rate higher growth than UK, Canada, France, Germany, Iran. Mexico, Russia and China, but lower growth than USA, Brazil, Turkey, Spain, Peru and Italy. It is important to note, however, that due to the discrepancies of COVID-19 testing, these countries reported numbers given by the World Health Organization are dependent on the amount of conducted tests. Further, daily confirmed cases are subject to change depending on a variety of factors.

However, in Fig. 2b, a comparison of confirmed death counts between top 15 COVID-19-infected countries from Fig. 2b elucidates that India has lower than all other 14 countries. Preventing people from interacting, each infected person is probable to spread the COVID-19 virus to less people, slowing the speed of the virus spread and permitting healthcare systems to handle the influx of patients. This shows the effectiveness of India’s continuous lockdowns as follows: first 21 days full lockdowns 25 March 2020–14 April 2020, second lockdown 19 days from 15 April 2020–3 May 2020, third lockdown 14 days 4 May 2020–17 May 2020 and finally forth ongoing lockdown during 18–31 May 2020.

Figure 3 narrates the daily confirmed cases and confirmed deaths of USA, UK, Spain, Brazil, France, Canada, and their corresponding generalized fractal dimensions begin with 3.89, 3.657, 3.557, 3.547, 3.155 and 2.72, respectively, which are shown in subplot. We observed that the transmit rate of USA, UK, Spain, Brazil, Canada seems to follow the same patterns, whereas France transmit rate is irregular not gradually increasing. Figure 4 describes the daily confirmed cases and confirmed deaths of Peru, Mexico, Germany, China, Russia, Iran and associated generalized fractal dimensions curve staring with 2.707, 2.963, 2.946, 2.642, 2.725, 2.732, respectively, which are presented in subplot. Further, transmit patterns of Peru, Mexico, Russia, Iran are similar and expect China remaining countries have mortality rate 2.785 in average.

Figure 5 shows the daily confirmed cases and confirmed deaths of Italy, Turkey, India, and their generalized fractal dimensions are revealed in subplot. Figure 5c clearly shows that generalized fractal dimension of death rate of India is monotonically decreasing from 2.51, whereas Italy and Turkey have initial value 3.016, 3.002, respectively. It should also be noted that remaining all 14 countries initial mortality rate starts above the value 2.51; hence, death rate of India is lesser than any other countries. This discrepancy might eventually get decreased or increased when more data on the outbreak in India will become available. However, the empirical data available as on 22.05.2020 are evident that India has more immunity against COVID-19. This is likely to be connected with many culture-related behaviours and to the presumed diversity in health conditions of the people.

4 Conclusion

The analyses presented in this study demonstrated that India has less mortality rate of COVID-19 as its total infected population is concern. Across the world, countries have executed a number of control procedures to respond to COVID-19, with the aim of slowing down transmission and dropping mortality. Particularly, India is in the process of scaling up public health and social measures based on local epidemiology and on local risk assessments, such as continuously four full lockdown, restricted local movement, and educational institutes are implementing the distance-learning and businesses to teleworking, and regulated international and national travel measures. Finally, the death rate is likely to be connected with many culture-related behaviours and to the presumed diversity in health conditions of the people. Notice that mortality and transmit rate might eventually get decreased or increased when more data on the outbreak of all countries will become available.