## The shutdown, while painful, appears to have saved a LOT of lives so far

If you recall, the growth of the new corona virus disease in the US (and many other countries) at first looked to be exponential, meaning that the number of cases (and deaths) were rising at an alarming, fixed percent each and every single day.

Even if you slept through your high school or middle school math lessons on exponential growth, the story of the Shah and the chessboard filled with rice may have told you that the equation 2^x gets very, very hairy after a while. Pyramid schemes eventually run out of suckers people. Or perhaps you have seen a relatively modest credit-card bill get way out of hand as the bank applies 8 percent interest PER MONTH, which ends up multiplying your debt by a factor of 6 after just 2 years!

(If the total number of deaths were still increasing by 25 percent per day, as they were during the middle of March, and if that trend somehow continued without slowing down, then every single person residing inside America’s borders would be dead before the end of May. Not kidding! But it’s also not happening.)

However, judging by numbers released by the CDC and reported by my former colleague Ron Jenkins, I am quite confident that THE NUMBER OF CASES AND DEATHS FROM COVID-19 ARE NO LONGER following a fixed exponential curve. Or at least, the daily rate of increase has been going down. Which is good. But it’s still not zero.

Let me show you the data and fitted curves in a number of graphs, which often make complex things easier to visualize and understand.

My first graph is the total reported number of deaths so far in the US, compared to a best-fit exponential graph:

During the first part of this pandemic, during the first 40 or so days, the data actually fit an exponential graph pretty well – that is, the red dotted line (the exponential curve of best fit) fit the actual cumulative number of deaths (in blue). And that’s not good. However, since about day 50 (last week) the data is WAY UNDER the red dots. To give you an idea of how much of a victory that is: find day 70, which is May 9, and follow the vertical line up until it meets the red dotted line. I’ll wait.

Did you find it? If this pandemic were still following exponential growth, now and into the future, at the same rate, we would have roughly a MILLION PEOPLE DEAD BY JUNE 9 in just the US, just from this disease, and 2 million the week after that, and 4 million the next week, then 8 million, then 16 million, and so on.

THAT AIN’T HAPPENIN’! YAY! HUZZAH!

As you can see — the blue and red graphs have diverged. Ignore the relatively high correlation value of 0.935 – it just ain’t so.

But what IS the curve of best fit? I don’t know, so I’ll let you look for yourself.

Is it linear?

This particular line of best doesn’t fit the data very well; however, if we start at day 36 or thereabouts, we could get a line that fits the data from there on pretty well, like so:

The purple line fits the blue dots quite well after about day 37 (about April 6), and the statistics algorithms quite agree. However, it still calls for over 80,000 Americans dead by May 8. I do not want the slope of that line to be positive! I want it to turn to the right and remain horizontal – meaning NOBODY ELSE DIES ANY MORE FROM THIS DISEASE.

Perhaps it’s not linear? Perhaps it’s one of those other types of equations you might remember from some algebra class, like a parabola, a cubic, or a quartic? Let’s take a look:

This is a parabolic function, or a quadratic. The red dots do fit the data pretty well. Unfortunately, we want the blue dots NOT to fit that graph, because that would, once again, mean about a hundred thousand people dead by May 8. That’s better than a million, but I want the deaths to stop increasing at all. Like this piecewise function (which some of you studied). Note that the purple line cannot go back downwards, because generally speaking, dead people cannot be brought back to life.

Well, does the data fit a cubic?

Unfortunately, this also fits pretty well. If it continues, we would still have about a hundred thousand dead by May 8, and the number would increase without limit (which, fortunately, is impossible).

How about a quartic (fourth-degree polynomial)? Let’s see:

I admit that the actual data, in blue, fit the red calculated quartic red curve quite well, in fact, the best so far, and the number of deaths by Day 70 is the lowest so far. But it’s impossible: for the curve to go downwards like that would mean that you had ten thousand people who died, and who later came back to life. Nah, not happening.

What about logarithmic growth? That would actually be sweet – it’s a situation where a number rises quickly at first, but over time rises more and more slowly. Like this, in red:

I wish this described the real situation, but clearly, it does not.

One last option – a ‘power law’ where there is some fixed power of the date (in this case, the computer calculated it to be the date raised to the 5.377 power) which explains all of the deaths, like so:

I don’t think this fits the data very well, either. Fortunately. It’s too low from about day 38 to day 29, and is much too high from day 50 onwards. Otherwise we would be looking at about 230,000 dead by day 70 (May 8).

But saying that the entire number of deaths in the US is no longer following a single exponential curve doesn’t quite do the subject justice. Exponential growth (or decay) simply means that in any given time period, the quantity you are measuring is increasing (or decreasing) by a fixed percentage (or fraction). That’s all. And, as you can see, for the past week, the daily percentage of increase in the total number of deaths has been in the range of three to seven percent. However, during the first part of March, the rate of increase in deaths was enormous: 20 to 40 percent PER DAY. And the daily percent of increase in the number of cases was at times over A HUNDRED PERCENT!!! – which is off the chart below.

The situation is still not good! If we are stuck at a daily increase in the number of deaths as low as a 3%/day increase, then we are all dead within a year. Obviously, and fortunately, that’s probably not going to happen, but it’s a bit difficult to believe that the math works out that way.

But it does. Let me show you, using logs.

For simple round numbers, let’s say we have 50,000 poor souls who have died so far from this coronavirus in the USA right now, and that number of deaths is increasing at a rate of 3 percent per day. Let’s also say that the US has a population of about 330 million. The question is, when will we all be dead if that exponential growth keeps going on somehow? (Fortunately, it won’t.*) Here is the first equation, and then the steps I went through. Keep in mind that a growth of 3% per day means that you can multiply any day’s value by 1.03, or 103%, to get the next day’s value. Here goes:

Sound unbelievable? To check that, let us take almost any calculator and try raising the expression 1.03 to the 300th power. I think you’ll get about 7098. Now take that and multiply it by the approximate number of people dead so far in the US, namely 50,000. You’ll get about 355,000,000 – well more than the total number of Americans.

So we still need to get that rate of increase in fatalities down, to basically zero. We are not there yet. With our current highly-incompetent national leadership, we might not.

===================================================================

* what happens in cases like this is you get sort of an s-shaped curve, called the Logistic or logit curve, in which the total number levels off after a while. That’s shown below. Still not pleasant.

I have no idea how to model this sort of problem with a logistic curve; for one thing, one would need to know what the total ‘carrying capacity’ – or total number of dead — would be if current trends continue and we are unsuccessful at stopping this virus. The epidemiologists and statisticians who make models for this sort of thing know a lot more math, stats, biology, and so on than I do, but even they are working with a whole lot of unknowns, including the rate of infectiousness, what fraction of the people feel really sick, what fraction die, whether you get immunity if you are exposed, what is the effect of different viral loads, and much more. This virus has only been out for a few months…

What’s the best approach – should we lock down harder, or let people start to go back to work? Some countries have had lockdowns, others have not. How will the future play out? I don’t know. I do know that before we can decide, we need to have fast, plentiful, and accurate tests, so we can quarantine just the people who are infected or are carriers, and let everybody else get back on with their lives. We are doing this lockdown simply because we have no other choice.

Published in: on April 27, 2020 at 12:33 am  Comments (1)
Tags: , , , , , ,

## COVID-19 Numbers in the US do not seem to be growing exponentially

Looking at the past month of CDC-reported infections and deaths from the new corona virus, I conclude that there has been some good news: the total number of infections and deaths are no longer following an exponential growth curve.

The numbers are indeed growing, by either a quadratic (that is, x^2) or a quartic (x^4) curve, which is not good, and there is no sign of numbers decreasing.

BUT it looks as though the physical-social distancing and self-quarantining that I see going on around me is actually having an effect.

Yippee!

Here is my evidence: the actual numbers of infected people are in blue, and the best-fit exponential-growth equation is in red. You can see that they do not match well at all.

If they did match, and if this were in fact exponential growth, we would have just about the entire US population infected by the end of just this month of April – over 300 million! That no longer seems likely. Take a look at the next graph instead, which uses the same data, but polynomial growth:

Just by eyeballing this, you can see that the red dots and blue dots match really, really well. When I extend the graph until the end of April, I get a predicted number of ‘only’ 1.5 million infected. Not good, but a whole lot better than the entire US population!

Also, let’s look at total cumulative reported deaths so far. Here are the CDC-reported numbers plotted against a best-fit exponential curve:

Up until just a few days ago, this graph was conforming pretty well to exponential growth. However, since about April 8, that seems to be no longer the case. If the total numbers of deaths were in fact growing at the same percentage rate each day, which is the definition of exponential growth, then by the end of April we would have 1.5 million DEAD. That’s THIS MONTH. Continued exponential growth would have 1.2 BILLION dead in this country alone by the end of May.

Fortunately, that is of course impossible.

Unfortunately all that means is that the virus would run out of people to infect and kill, and we would get logistic growth (which is the very last graph, at the bottom).

This fourth-degree mathematical model seems to me to work much better at describing the numbers of deaths so far, and has a fairly good chance of predicting what may be coming up in the near future. It’s still not a good situation, but it shows to me that the social and physical distancing we are doing is having a positive effect.

But let’s not get complacent: if this model correctly predicts the next month or two, then by the end of April, we would have about 60 thousand dead, and by the end of May we would have 180 thousand dead.

But both of those grim numbers are much, much lower than we would have if we were not doing this self-isolation, and if the numbers continued to grow exponentially.

====================================================================

FYI, a logistic curve is shown below. Bacteria or fungi growing in a broth will grow exponentially at first, but after a while, they not only run out of fresh broth to eat, but they also start fouling their own environment with their own wastes. WE DO NOT WANT THIS SITUATION TO HAPPEN WITH US, NAMELY, THAT WE ALL GET INFECTED!!!

## Various graphs for deaths from COVID-19, so far

I wrote that I would show you what various graphs of various types of simple models look like for deaths so far due to the current corona virus: linear, exponential, polynomial, and so on. I think that a fourth-degree (not third-degree, like I wrote earlier) seems to fit the data best so far, and that’s better than exponential growth.

First, let’s look at a straight-line best-fit model, superimposed by Excel on the data. (Note: deaths are on the Y, or vertical, axis; the X-axis represents days since the beginning of March. So today, the 6th of April, is day 37 (31 + 6). The dotted red line represents the line of best fit, and the blue dots are the CDC-announced numbers of deaths so far.

As you can see, the straight dotted line doesn’t fit the data very well at all. R-squared, known as the correlation coefficient, tells us numerically how well it fits. If R or R^2 equals 1.000, then you have absolutely perfect correlation of the data to your model. Which we do NOT have here. By the way, in that model, then by mid-June we would have about 22,000 dead from this disease.

OK, let’s look at an exponential curve-of-best fit next:

As you can see, this red curve fits the data a LOT better, and R-squared is a lot higher.

Unfortunately.

We do NOT WANT EXPONENTIAL GROWTH OF THIS OR ANY OTHER DISEASE, BECAUSE IT MEANS WE ALL GET IT! In fact, if this model is accurate and isn’t slowed down, then by mid-June, just plugging in the numbers, we would have 3.3 BILLION (not million) people dead in the US alone. Fortunately, that won’t happen.

BUT there are some parts of the data where the curve doesn’t fit perfectly — let me point them out:

At the upper right-hand end, the red dotted line is quite a bit higher than the blue dots. Fortunately. And near the middle of the graph, the blue dots of death are higher than the red line.

OK, let’s look at some polynomial models instead:

This is a fancy version of the simple y=x^2 parabolas you may have graphed in Algebra 1. Once again, this doesn’t do a terrific job of conforming to the actual data. At the right-hand end, the blue dots of death are higher than the curve. In addition, if we continued the red curve to the left, we would find that something like two thousand people had already died in the US, and presumably came back to life. Which is ridiculous.

However, if this model were to hold true until mid-June, we would have 127 thousand dead. Not good.

Let’s try a third-degree polynomial (a cubic):

That’s pretty remarkable agreement between the data and the equation! That’s the equation I was using in my earlier post. The R-squared correlation is amazing. Unfortunately, if this continues to hold, then we would have about 468 thousand dead in the US.

Let’s continue by looking at a fourth-degree polynomial curve fitted to the data:

That is an amazingly good fit to the data! Unfortunately, let’s hope that it won’t continue to fit the data, because if it does, then we are looking at a little over a MILLION dead.

Let’s hope we can get these totals to level off by physically distancing ourselves from other households, washing our hands, and getting proper protective garments and testing technology to our medical personnel.

=============

Here’s another model that unfortunately does NOT work: logarithmic growth. If it were the case, then we would have about 10,700 deaths by mid-June.