## 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: , , , , , ,

## 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.