Math Modeling: Day 31

Monday and Wednesday will be summarizing what we've seen (and hopefully learned) over the course of the last month and a half.

So I asked you to do three things by Monday, 4/13:

1. I would like you to read starting at section 3 "Analysis": pp. 8-12 (through to the conclusion). There are several tables and figures (figures are at the end, unfortunately, rather than embedded in the text). If you have questions, please ask.

2. Now here is a special development, that we want to investigate (mentioned in those sections of the paper):
   • There is a fair amount of evidence that the diurnal temperature range (DTR) -- the difference between maximum and minimum temperatures -- has actually decreased over much of the land surface of the globe.

     This is because minima are increasing faster than maxima, so the two are getting closer together.

     In some cases, it appears that maxima are not increasing much at all, perhaps because of cloudier conditions.

     Minima may be increasing faster in part because of night-time warming (when it is typically cooler).

   • So this changes our initial sense about what we would expect under climate change.

   I let you know about some new data -- modern temperature data -- from a small town (Custar) with a good weather station not far from Bowling Green (county seat of Wood County, and the station from which most of our modern Fletcher data derived).

   • The Custar data (csv format -- years 1982--2020 -- more than a climate period's worth!)
   • Bowling Green Climate Normals (csv format)

   I asked you to find the DTR for the Custar data (you have columns for the MaxAirTemp and MinAirTemp -- the DTR would be the difference of those two variables -- DTR=max-min).

   Regress the DTR, modeling it as a linear function: i.e. \[ DTR(t)=\beta t + \alpha \]

   The question is whether the parameter $\beta$ is significantly different from 0, and whether in fact it's negative -- so that the DTR in Custar has been decreasing over the 38 years or so of data.

3. Next, I asked you to look at climate normals of Bowling Green, OH (Wood County seat, and my hometown!:)
   1. to compute one of those sinusoidal models to the climate normals (max and min temperatures) of Bowling Green: e.g. \[ Max(t) = \alpha + A\sin{(t - t_0)} \] It's not hard, and you've certainly done a fair number of these. You'll perhaps be interested and surprised to see how sinusoidal the normals actually are.

      Here are some Mathematica commands to get started with the normals (creating the important variables).

   2. But the interesting twist is to fit a model of the form \[ Max(t) = \beta t + \alpha + A\sin{(t - t_0)} \] (for example), and to see what is happening in $\beta$ -- is it significant?

      If there were no climate change, then we would expect the climate normals to be periodic functions -- December wrapping right back into January. However, if, as we expect, there's global warming going on, then we'll see that normals at the end of the year are actually a little higher than the normals at the beginning of the year.

      And, if our hypothesis is right that the minima are changing more rapidly than the maxima, then the coefficient of the minimum will be larger than the coefficient of the maximum.

      See what you discover!

This is an interesting one. When I started looking for interesting Covid-19 models, I found this one InsightMaker, and it incorporated real data.

As applied math modelers, we need to make sense of real data. We're going to look at how some modelers tried to make sense of real data; how we need to change their model to reflect the real world, which intruded upon their theories....

Some good links that I might recommend (a few of which we'll focus on):

1. A great intro to many of the most important question: "When a new virus emerges, no one is immune. A highly transmissible virus, like the coronavirus behind the current pandemic, can spread like wildfire, quickly burning through the dry kindling of a totally naive population. But once enough people are immune, the virus runs into walls of immunity, and the pandemic peters out instead of raging ahead. Scientists call that the herd immunity threshold."

2. How epidemics like covid-19 end (and how to end them faster): A coronavirus causing a disease called covid-19 has infected more than 70,000 people since it was first reported in late 2019. To predict how big the epidemic could get, researchers are working to determine how contagious the virus is.

3. The SIR Model for Spread of Disease - Introduction (from the MAA).

   In particular, we will implement The SIR Model for Spread of Disease - The Differential Equation Model in InsightMaker.

4. Mathematics of the Corona outbreak, with Professor Tim Britton

   (An excellent introduction to SIR models, from both the infectious disease and mathematical sides)

   Questions:

   1. Where is the "calculus moment" in this video?

   2. How does Britton suggest breaking $R_0$ down into "actionable" pieces?

5. An on-line model developed by Ashleigh Tuite and David Fisman, Dalla Lana School of Public Health, University of Toronto

6. This one incorporates space into the model. These models are called "agent-based models".

7. Could Coronavirus Cause as Many Deaths as Cancer in the U.S.? Putting Estimates in Context

   This on-line estimator (i.e., a model!) allows one to estimate deaths, as well as death by age-category.

8. A nice R-based introduction to infectious diseases and nonlinear differential equations.

9. Use of a log-scale (which they take pains to explain) to illustrate deaths by country, updated daily. There is even a separate page at the New York Times to explain log plots: "A Different Way to Chart the Spread of Coronavirus: Those skyrocketing curves tell an alarming story. But logarithmic graphs can help reveal when the pandemic begins to slow."