top of page

Our understanding of the dynamics of Covid is mostly based on models of the infection process. The basic SIR model is the base for most epidemic model.  It splits the population into 3 buckets; Susceptible (S), Symptomatic Infected (Is) and Recovered (R).

 

I have made a number of additions to the model. 

 ​

 1) The infections per person (Ipp), often called "Reproduction Rate (R(t)), changes over time as social distancing, and masking effects the transmission of the virus. There is a exponential time constant for the effect of behavior changes on Rn(t) .

 2)  There is a significant population who are Asymptomatic Infected (Ia) and develop some immunity without knowing it. The model shows 75% of the infected are in fact asymptomatic.   This has a major impact on the size of the susceptible population. 

 3) The vaccinated population decreases the Susceptible population but most people who been infected will still chose to get vaccinated. Therefore the additional protected  population from vaccination is reduced by the fraction of already protected. In the summer of '21, With 60% vaccinated and 40% infected, the total protected is 76% which means the US is very close to  herd immunity assuming that asymptomatic infection gives significant protection.  

4) In the fall of '21, all the US population are protected to some degree and secondary infections start to become significant. First there are  much less infectious breakthrough infections in the vaccinated. Then the Omicron variant appears that is mostly limited to the unvaccinated population.  

5) The treatments for Covid have also improved enormously, the doctors have learnt how keep really sick  people alive and eventually cure them. The deaths per day are calculated using a Case Fatality Rate and a time delay from infection to death. In the first wave in Texas, the CFR was 3% after 5 days, and has improved to 1.7% after 20 days. The Omicron variant is 5x less lethal than the initial Covid infection. 

The dynamics of Covid infections can be seen in the log graph of daily new cases for different countries, early in the pandemic. In the first days, there was a very rapid rise  of 10x in only 10 days. Almost immediately countries instituted lockdowns which stalled the infections and then saw them fall at different rates.  The straight line rise and fall in the log graph is characteristic of the exponential growth of an infection. 

 

 

The idea was to develop a version of the standard SIR
model of infections that accounts for  changes in social isolation behavior  and health care.

In the classic SIR model, the change in number of  infections over time dI/dt depends on the rate of new cases minus the rate of recovery. The new cases per day is given by the infection rate constant beta, multiplied by the number of Infected people I(t) and the fraction of people available to be infected S (Susceptibles). The rate of recovery dR/dt is given by the recovery time constant gamma multiplied by the concentration of Infections I(t).

The ratio of infection to recovery rate constant (beta/gamma) has useful meaning as the initial Infections per person or Ipp(0), often called R0. It is the gain in the chain reaction of infections. For Covid, the Ipp of 5 means that every infected person infects 5 more for the roughly 10 days they are infectious.  In thirty days there will be 5^3 =125 x increase in  infections !! That is why Covid is so dangerous.  

 

In reality, peoples social behavior changes to avoid infection so the rate of infection  changes over time beta(t) and so does Infections per person Ipp(t). 

It has become clear that there is a large number of "hidden" infected individuals who are  asymptomatic and are not caught by current US testing that is mostly limited to asymptomatic. At current testing levels, the number of asymptomatic people caught by random chance is very low. Estimates vary, but using a ratio of protected = 4 x symptomatic (Asymp : Symp = 3:1)  as noted below.

These differential equations were numerically integrated in an Excel spreadsheet, to model an infection. 
 
The model can be used to extract metrics of the
process. At the beginning of the infection, the susceptibles
can be assumed close to 1, so the infection grows exponentially and can visualized as a straight line in a log graph, with a slope of Ipp.

 

The daily infections d log(I)/dt, new cases d log(C)/dt, deaths d log(D)/dt, all have the form = f(I), so they will track each other. Hence the universal use of log plots to analysis pandemics.


Now we have a way to quantify the stages of the infection illustrated in the graph of daily new cases for different countries. In very the early stages, the daily cases have a steep log slope consistent with a initial Ipp of over 5. In response, various levels of lockdown were instituted, and the infection stalled and then started to drop with countries showing different declining log slopes with a Ipp of 0.7-0.9. The effect of the lockdown was to change behavior, which in turn changed in the infections per person. It appears that the evolution of Ipp over time can be generalized as a series of constant levels.

Obviously human behavior does not change as step function. Manually adjusting Ipp to fit gave a curve that looked like an exponential.  Using a exponential function, gave a reasonable fit to US, NY, TX  and WA data with a 1/e time constant of 7  days, and around 20 days for the change in infection per person to stabilize.  This leads to a generalized expression for Ipp(t) with a series of constant levels linked by exponential approach, that was applied to TX data. 

The idea that the transition between 2 different levels of social
behavior is an exponential seems inherently   reasonable, and  has been used to describe transitions in isolation for  other infections (Ref 1). A possible physical interpretation has been  proposed, borrowed from finance (Ref 2), which suggested that the transition might be  an "Ornstein-Uhlenbeck" process.

 

Finally, published data on mobility based on cell phone data also shows an exponential response to the lockdowns in March.

 

Vaccinations started in late December. Vaccinations reduce the Susceptible population, and hence the new infections. However, there is a random chance that people getting vaccinated were already protected by being infected, so the contribution must be reduced by the fraction of already protected. Therefore the new vaccinations per day multiplied by (1-cume cases*asymp ratio) were added to the Susceptible every day.  The vaccinated were lagged by 10 days to account for jab to work. 

 

It looks like the average distancing behavior has remained unchanged from early November '20 to April '21 which seems consistent with observations on the ground. 

 

Finally in April '21 and June '21, distancing and mask use effectively end, causing increases in Ipp and  bumps in cases through June '21 when Ipp has returned to the levels at the beginning of the pandemic. 

The wave in August '21 finally rolled off in Sept '21. Using a Asym : Sym ratio of 3 resulted in the model for the waves in April and August to roll off  with constant distancing (Ipp). This suggests that the model is on the right track. In Oct, the remaining unprotected from either vaccine or infection is down to 7% as shown in the graph of Susceptible fraction over time.

At a susceptible fraction of 7% the primary pandemic is over. At that point breakthrough infections with a low  Ipp of 1.3 start to dominate. Shortly after the Omicorn variant produces a rapid increase and then decrease in cases.  The wave is best fit with a Ipp infection of 8 on a 30% susceptible population. This is almost certainly the unvaccinated  as we know that vaccination reduces the risk of infection by at least 5x. 

The waves of infection in NY mimic the same trends as Tx. 

Quantifying the effect of social distancing 

 

The Infections per person (Ipp) are driven by social distancing and the infectiousness of the dominant variant.  Goggle created a mobility index based on individuals volunteering to make  their cell phone location data available, which can be used as a proxy for social distancing.  The retail, grocery, transit and workplace mobility metrics showed the greatest variation  during the pandemic. Assuming that the mobility at the peak in the early pandemic was closet to zero, then at a mobility of -20% of normal the Ipp becomes<1 and the infection stops growing. The Delta variant was responsible for the peak in the fall of 2021. 

There is a roughly 1 month lag between a change in mobility and a change in cases. 

https://covid.cdc.gov/covid-data-tracker/#variant-proportions

https://www.google.com/covid19/mobility/

 

 


Modelling the time line for fatality rate. 

 

For the unlucky few who die, there is a time lag L between infection and  death, and there is evidence (shown below) that this time has lengthened as the doctors have learnt how to manage the disease. Therefore, the number of deaths per day is given by the death rate times the number of infected people at the time of infection  I(t-L). Similarly, the Case Fatality Rate can be followed over time as deaths per day / cases per day lagged by time to death L.  

The number of  deaths per day are obviously to the number of cases per day. The log plot of daily  deaths and cases shows the log distance between deaths and cases varying, which  means that the ratio of deaths to cases (CFR) is changing. Also the time lag between the peaks in March are smaller then than the lag in August, which means that the time lag from infection to death is also changing. 

TX had a relatively mild infection in March, so it avoided any effects of patient load on outcomes. A series of CFR's at different lags for TX are shown. Assuming that CFR is likely to be relatively constant, it appears that prior to June, the lag was 7 days with a CFR of 3%, and post June the  lag was 20 days with a CFR of 2%.  The changing CFR and lag presumably a result of improvements in medical care. 

Based on these observations the death rate was modelled as a series of constant steps. 

 

In January, vaccinations have started to take hold with about 8% of the population. This is focused on the >65 group who constitute 15% of the population and are where almost all of the deaths occur. About 40% of the at risk community have been vaccinated. There is a 20 day lag from case to death and 10 days to generate immunity. In 1 months time the death rate should be halved. The graph for TX  shows the projected death account after adjusting for vaccination. 

Conclusions 

The modelling analysis led to a number of key characteristics of Covid infection. 

The Infections per person (Ipp) vary as behavior changes.

 

Under "normal life" behavior, Ipp is around 5. 

Under lockdown conditions Ipp follows an exponential approach to a new constant value that takes about 20 days. 

The death rate and the time from infection to death both change, probably due to improvements in medical care. 

References


Thanks to Andrew Zachary for pointing out the link to financial
modelling.

Ref 1 https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5348083/

Ref 2 https://www.linkedin.com/pulse/modeling-term-structure-
pandemic-negative-interest-peter-cotton-phd/?
trackingId=K6StB6cgwAAA8dSumxqOfQ%3D%3D

Ref 3 https://en.wikipedia.org/wiki/Ornstein%E2%80%
93Uhlenbeck_process

Ref 4 "Modelling to Inform Infectious Disease Control" by Neils
Becker (2013)

Ref 5 "Repeat contact and the spread of disease" by Peter Cotton
(2020)

Trends_in_cases_.jpg
Eq5.jpg
TX.jpg
Susceptible.jpg
NY.jpg
Mobility1.jpg.png
Lag2.jpg
Lag1.jpg
bottom of page