Estimating the Parameters of Covid-19 Cases in South Africa

in this paper we employ SiR model to study the Covid-19 data of South Africa for a chosen period. This model is solved using three numerical methods, namely, differential Transform Method (dTM), Multistage differential Transform Method (MsdTM), Repeated Multistage differential Transform Method (RMsdTM) to obtain approximations of the number of susceptible, active infected and recovered in South Africa for 60 days starting from June 1, 2021. The proximity of the solution of the RMsdTM to the actual data in comparison to solutions using the other two methods was observed. MsdTM is an improvement over dTM as it uses updated values of the variables as new initial conditions at each iteration of the method. RMsdTM, in which the values of parameters are also changed at suitable intervals of time, besides using updated values of variables is a further improvement over both these methods.

Time and again, many pandemics have wrecked havoc on mankind. Spanish flu, almost a century ago; H1N1 Swine flu pandemic, a decade back; and from recent past -the Ebola and Zika virus pandemics have affected mankind adversely. Coronavirus disease is a highly contagious disease caused by SARS-CoV-2 1 . Post its probable initiation in Wuhan, China, in December 2019, it spread rapidly across the globe, turning into one of the worst pandemics that human civilisation has ever seen. The Coronavirus was a novel virus, causing mild to moderate respiratory illness in most of the affected people. However, in some cases, the virus also led to serious diseases, sometimes even leading to death. There has been extensive research to understand the behaviour of this virus, in order to contain its spread and cure the infected. Many lives have been lost and a huge economic crisis has engulfed several countries due to the restrictions in travel and trade. Epidemics are region-specific and each country is equipped to handle epidemics affecting their territory. However, a pandemic spreads to several countries and sometimes affects the whole world. The response of each country to the same disease varies and is dependent on its socio-economic structure. A developed country's preparedness may also not be sufficient to handle such pandemics due to social behaviour, as could be seen in the first wave of CoViD-19 in the United States, United Kingdom and italy. on the other hand, many developing countries like Africa, have been able to contain its spread.
Several Mathematical models like SiS, SiR, SiRD, SEiRD etc., have been developed and employed to study the behaviour of epidemics and make predictions 2. Many such models have been utilized to study this new virus as well 3 in order to predict the number of infected, so that the medical systems are well equipped to handle the predicted number of cases at a given time.
SiR model is a classic compartmental model that was introduced in early 20th century 4. The three compartments of this model, as shown in Figure 1, consist of the susceptible, s(t), representative of that component of the population who have not been infected and are probable infectives; the infected, i(t), standing for individuals who have been infected and are capable of spreading the infection and the recovered/ removed, r(t), which include those who have recovered from the infection and developed immunity, as well as the deceased. in this model, the members of r compartment cannot re-enter the s compartment. The total population, n, is considered to be fixed i.e. n = s + i + r. The rate of transmission at which susceptible get infected is b and the rate at which infected get recovered is g.
The following system of differential equations is formed for the SiR model 5 in this paper we employ SiR model to study the Covid-19 data of South Africa for a period of 60 days 6. in South Africa, the first Covid-19 patient was confirmed on 5th March 2020 7 . The person who tested positive for the virus being a thirty-eight-year old man who had returned from a trip to northern italy 8. Since then, the country has seen several peaks in the number of daily infections, caused due to many variants. in fact, in August most mutated Covid variant C.1.2 was also detected in South Africa. The South African government has continuously been monitoring the pandemic situation and imposing or relaxing restrictions like lockdown, air-travel ban etc. from time to time 9 . in this paper, we have undertaken the study of Covid-19 data for a portion of the period of the third wave, starting from 1st June, 2021 to 30th July, 2021.
The paper is organised as follows: The Differential Transformation Method (DTM) is explained in Section 2. in this section, we also initiate a case study of the Covid-19 pandemic scenario for the above stated period using Differential Transform Method on SiR model. The case study is extended in Section 3 using the Multistage Differential Transform Method (MsDTM). We introduce Repeated Multistage Differential Transform Method (RMsDTM) in Section 4 and obtain the estimations of the undertaken case study using this method. A comparison of the outcomes of the three methods is attempted in Section 5. We have given the conclusion of the paper in Section 6.

Differential Transform Method (DTM)
Methodology. The Differential Transform Method dates back to 1986 when Zhou in his paper 10 on linear and nonlinear initial value problems in electric circuits introduced and applied the method to obtain an approximation of the solution. Due to better approximations of the solution in comparison to other methods for solving such systems, DTM has since been used extensively to study initial value problems of systems of linear and nonlinear ordinary differential equations and also systems of partial differential equations 11,12,13 Let f be a function defined on some open interval containing a point c. if f possesses derivatives of all orders at c, then the series is called the Taylor series for f about c. in particular, Taylor series for f about c = 0 is given by The function f is obtained as an inverse differential transform as follows Some algebraic and analytical properties of the transform function are listed in Table 1 in Differential Transform Method, the function f is approximated by finite degree polynomial of arbitrary degree K, i.e.
The remainder terms of the series in equation (1) represent the error in the above approximation and is negligibly small.

DTM for SIR model for Covid-19 in South Africa
in this section we use DTM to solve the SiR model to find susceptible, active infected and recovered of Covid-19 cases in South Africa, for a chosen time period. if s(t), i(t) and r(t) represent the number of susceptibles, active infections and recovered at time t and S, i and R are their respective differential transform then the transformed system for SiR model is 14 We consider a period of 60 days, starting from 1st June, 2021. The initial conditions on s, i and r, as obtained from actual data, and the parameters used 15 are as follows Using the inverse differential transform, the approximation of the polynomials for s, i and r upto degree 2 were found to be The graphs in Figure 2 show the infectives and recovered found using the above solution for the period of 60 days, along with the actual infectives and recovered of that period.
in DTM, the number of infected and recovered were obtained using the polynomials given in 2. These are quadratic polynomials, each with positive leading coefficient. So, these are increasing functions for x>0. This is the reason for the number of infected and recovered obtained using this method to be increasing. Whereas in actual data, the numbers are not showing such pattern.

Multistage Differential Transform Method (MsDTM) Methodology
Multistage Differential Transform Method (MsDTM) 16, 17 is an extension of DTM. in this method, the initial value of the function is dynamically updated at equal intervals of time (Tools of Mathematica 18 have been used for these calculations.). This leads to a better approximation of the We assume that the interval [0, T] is partitioned into p-sub-intervals [ti-1, ti], i = 1, … , p, where t 0 = 0 and tp = T. The sub-intervals are of equal length h, where, h = T/p. DTM is applied to the first sub-interval, i.e. [0, t 1 ] with the given initial conditions, to obtain the approximate value of the function at time t 1 . This value is now taken as the initial condition for the second interval, [t 1 , t 2 ] and DTM is applied. However, here and in subsequent intervals, the Taylor series for f about the point t i-1 is considered  Table 1 for obtaining the solution. in the i-th sub-interval, the solution is given by

MsDTM for SIR model for Covid-19 in South Africa
The case study initiated in Section is now attempted using MsDTM. The period of 60 days is sub-divided into 30 periods, each of length 2 days. So, p = 30 and h = 2. Using the polynomials obtained at each step and resetting the initial conditions as explained, the values of s, i and r were obtained for the period of 60 days (values for 2 days at each step).
The graphs in Figure 3 show the number of infectives and recovered obtained using MsDTM in comparison to the actual data.
in MsDTM at each step a different quadratic polynomial is obtained for the number of infected and recovered. However, as per the method, the polynomials obtained in successive iterations agree at the last values of these parameters obtained in the previous and the values taken as initial conditions for the succeeding iteration. So, the graph is a joined graph of different quadratic polynomials, with leading coefficient of each being positive. This justifies the increasing graphs, as seen in Figure 3.

Repeated Multistage Differential Transform Method (RMsDTM) Methodology
Multistage Differential Transform Method has an advantage over DTM in terms of dynamic       length, say sub-intervals of length each. MsDTM is applied to each subinterval, using the value of function obtained from the previous interval as the initial condition and also using new values of the parameters.

Repeated MsDTM for SIR model for Covid-19 in South Africa
The active infections of South Africa during the period of study indicated that the rate of transmission and rate of recovery were changing with time. Keeping this in view, we divided the interval of 60 days into four equal parts of a fortnight each. At each step, as explained in the previous sub-section, updated values of and were used as new initial conditions. The updated values, used as initial conditions at each step, are tabulated in Table 2.
The graphs of active infections and recovered, obtained for the period of 60 days using this method versus the actual data of both categories, are presented in Figure 4.
We observe from the graph of infected that the graph is increasing in some period and decreasing in some other, in the same way as the actual data. Such variation has been possible by choosing different values of the parameters and as per the prevailing conditions of that time. This leads to obtaining such polynomials which give number of infected close to the actual number. The graph of number of recovered also shows similar behaviour and can be seen almost overlapping the actual data throughout the chosen period.

ReSulTS AnD DISCuSSIon
in this study, we applied the three numerical methods, namely, DTM, MsDTM and RMsDTM to SiR model and obtained approximations of the susceptible, infected and recovered in South Africa, for a chosen time period. out of these three methods, Repeated MsDTM gives values which are converging to the actual values. The solutions of DTM and MsDTM are comparable but are in huge variance from the actual data. Figure 5 depicts the number of infected and recovered, obtained using all three methods, for the chosen period of 60 days. The proximity of solution of RMsDTM to the actual data in comparison to solutions using the other two methods can be observed. We can also see from the graphs that the solutions obtained using DTM and MsDTM are close to each other but far from the actual data.
We studied the actual data and inferred that rate of transmission and rate of recovery was showing variation at various time periods within the 60 days interval. This suggested that we choose appropriate values of the parameters and for each period as shown in Table 2. The recovery rate was increasing with time. The transmission rate, however, was showing a different pattern. After increasing in the second fortnight, it showed a steady decline. in RMsDTM, introduced in this paper, we take these variations into consideration and approximate these parameters at different time intervals, from the given data. These values of parameters are used at respective time intervals to obtain a solution of SiR model which is converging to the actual numbers.
The effectiveness of RMsDTM is also visible through the graphs in Figure 6, depicting errors in the number of infected and the number of recovered, obtained using these three methods. out of the three, the error can be seen to be least in RMsDTM's solution.
The error values are tabulated in Tables 3  and 4. it is clearly seen that the errors in solutions of RMsDTM are much lower than those of the other two methods. out of DTM and MsDTM, error in solution of the later is lesser.

ConCluSIon
MsDTM is an improvement over DTM as it uses updated values of the variables as new initial conditions at each iteration of the method. it is also verified in the case study undertaken in this paper. However, considering the need of updated parametric values as well, due to changing onground conditions during a pandemic from time to time, we have introduced Repeated Multistage Differential Transform Method (RMsDTM) in which the values of parameters are also changed at suitable intervals of time, besides using updated values of variables. in this paper, through the case study for Covid-19 scenario in South Africa, we have presented a comparison between number of infected and recovered obtained using DTM, MsDTM and RMsDTM. it may be concluded from the discussions that DTM and MsDTM solutions are convergent, if the time interval is divided into periods of short length (2 days in our study). However, the solutions of both the methods are differing from actual data by a considerable number. The reason behind introducing RMsDTM is justified when we see that the solutions obtained using this method are converging to the actual data.