Prediction
of Third Wave in India
Dr
Himanshu Shekhar
I
dedicate this blog to my Math Professor Late Shri SB Singh of MIT,
Muzaffarpur, Bihar who nurtured and honed my mathematical skills.
A
predictive approach using sinusoidal variation is attempted in this post of the
blog. This post was created due to absolute faith of my Friend TK Sadhana Singh
in me and my calculations. The encouragement from each reader is my asset and I
look forward for more input to give me correct direction for prediction.
Introduction
England has claimed to have achieved
very nice progress in vaccination, which resulted in poor penetration of virus at
the general health of people. Although England has attributed control, due to
lockdown, but at the same time, now vaccination is stated to be the main reason
for non-fatal next variant or wave in England. India is not able to progress rapid
vaccination due to vast population and scarcity of vaccines. However, based on
the available data till date, an attempt is made to predict possibility of
third wave, using various mathematical technics. As number of data-points specifying
waves is insufficient, some assumptions are made for the prediction. First an
attempt is made to model both the waves in India, using single equation of
damped vibration.
Nature of Waves
India has seen two waves, each
attained a peak, followed by a decline. The reasons for advent of second peak
was neither known nor investigated, properly. In absence of any correct social
engineering parameters known, it is purely mathematical exercise, which is
repeated. At first nature of the variation of daily confirmed cases are
observed. Although normal distribution curve fits well, but this is not
sufficient to predict the re-occurrences and wave nature. Normal distribution
has only one peak at mean position with proper standard deviation representing span
of each wave. To get rid of this, damped vibration equations are fitted, after
assessing the nature of both the waves. First concern is peak confirmed cases
in both the waves, which has a multiplication factor of 4. First wave has peak
cases of 1 lakh, while second wave has 4 lakh peak daily cases. Second feature
is span of the wave on time scale. Time span of both the waves has a factor of
4.
For modelling this, curves are
generated incorporating these features, where peak value of each wave changes
and so is the time-scale. First curve represents equation of the type y = a.[1+Sin(2pt/T)], where a is amplitude or
peak value of the sinusoidal wave, T is time period of wave and t is time
lapsed. The addition of unity is to make all the terms of the curve positive. The
plotted curve is shown below.
The indicated curve has same amplitude
and time-period for each oscillation. Exponential decay is introduced in peak
by incorporating exponential term. The governing equation becomes y = a.e(-bt).
[1+Sin(2pt/T)],
where b is exponential decay factor. A typical representative curve is
indicated below for decaying peak of sinusoidal wave.
However, the actual curve is
not decaying in peak but each subsequent wave has higher peak than the previous
wave in current situation. This can be handled by taking a negative value of
exponential decay factor. For time compression, the time-period is also assumed
to be varying with time, linearly. The same is incorporated in equation and T
is replaced by (T-td), where d is time-period decay factor. The produced curve
is shown below.
The below represented
exponential decay and varying time step. If now the mirror image of the curve
is taken in y-axis, the designed curve can represent the actual situation of
confirmed cases. The curve has reduced peak for each wave and the time period
is high for first wave followed by sharp decline for the second wave and so on.
This type of equation is used for simulating the curve for daily confirmed
cases.
Integrated Equation for Waves
The process of equation development is followed for the
prediction and the variation of daily confirmed cases is simulated using
equation mentioned in previous section. The peak cases is amplified by a factor
4 and the time period is reduced by a factor of 4, as obvious from the used
equation. The curve generated is superimposed over the actual number of daily
confirmed cases.
It is clear that curve is
representing to a fair degree of accuracy, the daily confirmed cases. As per
this, the third wave of curve will follow immediately after end of second wave
and probably it rise will be much faster. The peak may be of the order of
around 4.74 lakh and the phase may start by 07 June 2021 and end by 29 June
2021. However, if somehow, the lockdown is extended up to 15 June 2021, all
over India, the occurrences can be avoided and probably then no third wave will
come.
Conclusion
As concluding remarks, I would
like to reiterate that the daily confirmed cases may be simulated as sinusoidal
variation with varying peak and time-period. If mathematics is to be believed,
the third wave will be faster and not more severe than second wave. If lockdown
is extended till 15 June 2021, the third wave can be avoided. This is a purely
mathematical exercise and the assessment needs validation in course of time. I
personally don’t see any chances of third wave in India, whether vaccination is
complete or incomplete.
Stay safe, stay healthy,
follow social distancing, and remain connected through social media. All the
best.
Regards,
Dr Himanshu Shekhar


































