{"id":1505,"date":"2020-03-17T12:29:15","date_gmt":"2020-03-17T16:29:15","guid":{"rendered":"http:\/\/josmfs.net\/?p=1505"},"modified":"2020-10-01T11:09:00","modified_gmt":"2020-10-01T15:09:00","slug":"coronavirus-mathematics","status":"publish","type":"post","link":"https:\/\/josmfs.net\/wordpress\/2020\/03\/17\/coronavirus-mathematics\/","title":{"rendered":"Coronavirus Mathematics"},"content":{"rendered":"<p><img loading=\"lazy\" decoding=\"async\" class=\"alignleft size-full wp-image-1504\" src=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/03\/Coronavirus-Math-Fig.jpg\" alt=\"\" width=\"200\" height=\"228\" \/>Given the mathematical nature of this website I feel reluctantly impelled to address the coronavirus pandemic. The mathematics behind the spread of infection is basically the same exponential growth that I discussed in the \u201c<a href=\"https:\/\/josmfs.net\/wordpress\/2019\/01\/10\/math-and-religion\/\">Math and Religion<\/a>\u201d post and has recently been explained by the ever-lucid Grant Sanderson at his <a href=\"https:\/\/www.youtube.com\/watch?v=Kas0tIxDvrg\">3Blue1Brown<\/a> website.<\/p>\n<p>What I wish to draw attention to is the series of posts on the coronavirus by Kevin Drum on his <a href=\"https:\/\/www.motherjones.com\/kevin-drum\/2020\/03\/coronavirus-growth-in-western-countries-march-16-update\/\">website at Mother Jones<\/a>. I have collected his recent posts comparing the spread of the virus in various countries and added some mathematical commentary of my own, which is the content of this post.<\/p>\n<p>But the bottom line seems to be that in virtually <em>all<\/em> the countries, including the US, the virus infection is spreading at the Italian rate of <em>doubling every 4 days<\/em>! The readers of this website are sufficiently numerate to realize the frightening import of that number. If that weren\u2019t enough, Kevin provides additional posts on the results of the <a href=\"https:\/\/www.motherjones.com\/kevin-drum\/2020\/03\/heres-the-punch-in-the-gut-version-of-the-imperial-college-coronavirus-study\/\">modeling at the Imperial College<\/a> that are truly nerve-wracking for someone such as myself in the most vulnerable cohort. The only blessing so far seems to be that, for once, the children are spared.<\/p>\n<p>At this time, I don\u2019t have the stomach to keep updating the post as new numbers come in. That may change. I could address the catastrophe of having ignorance and incompetence at the helm of the national ship of state, but it is too depressing.<\/p>\n<p>See the <a href=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/03\/Coronavirus-Drum-200317a.pdf\">Coronavirus Mathematics<\/a><\/p>\n<p><strong>(Updates 3\/17\/2020, 3\/21\/2020, 4\/17\/2020, 9\/20\/2020, 10\/1\/2020)\u00a0<\/strong> <!--more--><\/p>\n<p><strong>(Update 3\/17\/2020)\u00a0 <\/strong>That didn\u2019t take long. But Kevin Drum had an <a href=\"https:\/\/www.motherjones.com\/kevin-drum\/2020\/03\/a-chart-qa-where-is-the-coronavirus-pandemic-headed\/\">update<\/a> to his curves that I thought was relevant, so I have included it in a revised version of his posts. He is also saying his plots show a 5-day doubling rate instead of the 4-day that I had determined\u2014infinitesimal consolation.<\/p>\n<p><strong>(Update 3\/21\/2020) <\/strong>Kevin Drum had a significant update regarding the mathematical model for the growth of the infection. It is not ultimately exponential or even logistic. Contrary to religious conversion, which is essentially permanent, infections go away when the subject either is cured (likely) or dies (less likely). So the ultimate shape of the curve is closer to a Gaussian distribution.<\/p>\n<p style=\"padding-left: 40px; text-align: center;\"><a href=\"https:\/\/www.motherjones.com\/kevin-drum\/2020\/03\/coronavirus-growth-in-western-countries-march-18-update\/\">Coronavirus Growth in Western Countries: March 18 Update<\/a><br \/>Kevin Drum, March 19, 2020<\/p>\n<p style=\"padding-left: 40px;\">Here\u2019s the coronavirus growth rate through Wednesday. The US is now slightly above the Italian trendline, possibly because we\u2019ve picked up the pace of testing.<\/p>\n<p style=\"padding-left: 40px;\">A couple of notes. First, a friend emails to warn me against calling the spread of coronavirus \u201cexponential.\u201d Technically, that\u2019s true. At the beginning of a viral outbreak the growth rate is pretty much exponential, but it slows down later. This <a href=\"https:\/\/www.zdnet.com\/article\/graph-theory-suggests-covid-19-might-be-a-small-world-after-all\/\">article<\/a> suggests that it then follows a power law, and eventually it flattens out and then declines. When all\u2019s said and done, the outbreak will follow a Gaussian path that looks like a familiar bell curve. For example, here is the growth curve of the Spanish flu in 1918, skyrocketing from nothing to massive death tolls over the course of only four weeks:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-1513 aligncenter\" src=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/03\/Coronavirus-Math-Influenza-1918.jpg\" alt=\"\" width=\"500\" height=\"378\" srcset=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/03\/Coronavirus-Math-Influenza-1918.jpg 500w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/03\/Coronavirus-Math-Influenza-1918-300x227.jpg 300w\" sizes=\"auto, (max-width: 500px) 100vw, 500px\" \/><\/p>\n<p>(<strong>Update 4\/17\/2020<\/strong>) After my previous postings, I would hope you are continuing to monitor Kevin Drum\u2019s excellent updates and discussions. I would like to mention a different mathematical take on the situation, namely, simulations. Grant Sanderson at his excellent <a href=\"https:\/\/www.youtube.com\/channel\/UCYO_jab_esuFRV4b17AJtAw\/videos\">3Blue1Brown<\/a> Youtube blog (3\/27\/2020) has a marvelous <a href=\"https:\/\/www.youtube.com\/watch?v=gxAaO2rsdIs\">epidemic simulation<\/a> using a simple SIR model that illustrates the effects of<\/p>\n<ul>\n<li>Basic Setup (growth rate sensitive to # daily interactions, probability of infection, duration of illness)<\/li>\n<li>Identify and isolate<\/li>\n<li>Social distancing<\/li>\n<li>Travel restrictions<\/li>\n<li>R (effective reproductive number: #expected new infections by infected person during illness)<\/li>\n<li>Central hubs (periodically visit grocery store, etc.)<\/li>\n<\/ul>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignleft size-full wp-image-1564\" src=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/04\/Coronavirus-Math-Update-200417-Fig.jpg\" alt=\"\" width=\"942\" height=\"527\" srcset=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/04\/Coronavirus-Math-Update-200417-Fig.jpg 942w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/04\/Coronavirus-Math-Update-200417-Fig-300x168.jpg 300w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/04\/Coronavirus-Math-Update-200417-Fig-768x430.jpg 768w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/04\/Coronavirus-Math-Update-200417-Fig-624x349.jpg 624w\" sizes=\"auto, (max-width: 942px) 100vw, 942px\" \/><\/p>\n\n\n<p><strong>(Update 9\/20\/2020)<\/strong>&nbsp; Unusual Effectiveness of Masks. Henry Reich at Minute Physics has a great <a href=\"https:\/\/josmfs.net\/wordpress\/coronavirus-why-masks-work-better-than-youd-think-y47t9qlc9i4\/\">video<\/a> (at <a href=\"https:\/\/www.youtube.com\/watch?v=Y47t9qLc9I4\">Youtube<\/a> with more information and links) about how much more effective masks are when many people wear them.&nbsp; For example, if two people are wearing masks that are each 50% effective, then together they become 75% effective in reducing the spread of virus particles.&nbsp;<\/p>\n\n\n\n<div class=\"wp-block-image\"><figure class=\"aligncenter size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"400\" height=\"213\" src=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Two-Person-Mask-Effectiveness-100dpi.jpg\" alt=\"\" class=\"wp-image-1771\" srcset=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Two-Person-Mask-Effectiveness-100dpi.jpg 400w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Two-Person-Mask-Effectiveness-100dpi-300x160.jpg 300w\" sizes=\"auto, (max-width: 400px) 100vw, 400px\" \/><\/figure><\/div>\n\n\n\n<p>And in a large crowd if 50% of the people are wearing masks that are 50% effective, then the overall effectiveness is 44%.&nbsp; So there is effective protection even when half a population is not wearing masks.&nbsp;<\/p>\n\n\n\n<div class=\"wp-block-image\"><figure class=\"aligncenter size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"401\" height=\"228\" src=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Multi-Person-Mask-Effectiveness-100dpi.jpg\" alt=\"\" class=\"wp-image-1769\" srcset=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Multi-Person-Mask-Effectiveness-100dpi.jpg 401w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Multi-Person-Mask-Effectiveness-100dpi-300x171.jpg 300w\" sizes=\"auto, (max-width: 401px) 100vw, 401px\" \/><\/figure><\/div>\n\n\n\n<p>In fact, if 60% of the population would wear masks that are 60% effective, the infection rate R would drop below 1 and so the virus would stop spreading and be contained.<\/p>\n\n\n\n<div class=\"wp-block-image\"><figure class=\"aligncenter size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"400\" height=\"191\" src=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Stopping-Spread-100dpi.jpg\" alt=\"\" class=\"wp-image-1770\" srcset=\"https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Stopping-Spread-100dpi.jpg 400w, https:\/\/josmfs.net\/wordpress\/wp-content\/uploads\/2020\/09\/Coronavirus-Stopping-Spread-100dpi-300x143.jpg 300w\" sizes=\"auto, (max-width: 400px) 100vw, 400px\" \/><\/figure><\/div>\n\n\n\n<p>Henry Reich points to another <a href=\"https:\/\/www.youtube.com\/redirect?q=https%3A%2F%2Faatishb.com%2Fmaskmath&amp;redir_token=QUFFLUhqbEkyclhWdHU1Z2Zwd1kzT05IODdQZHluOVRvUXxBQ3Jtc0trS1lmNWliMC0ybS0yVU91MFdicDBEckRlZG9VWkdRcWN5RENNVllrVzh5RmxmYWRjeWQzYkViR01IVDhvZUZ1WGk3T19oNmx4LU04Zlo0aTZTLUdNT3NM\">website<\/a> by Aatish Bhatia that goes into more detail with animated graphics.<\/p>\n\n\n<p><strong>(Update 10\/1\/2020) <\/strong>Key to the Pandemic<\/p>\n<p>This article in the Atlantic on the pandemic is extremely significant and bears close reading.\u00a0 I have included some excerpts of main points, but the whole article is fascinating and helpful for understanding the pandemic.<\/p>\n<p>(<a href=\"https:\/\/www.theatlantic.com\/health\/archive\/2020\/09\/k-overlooked-variable-driving-pandemic\/616548\/\">https:\/\/www.theatlantic.com\/health\/archive\/2020\/09\/k-overlooked-variable-driving-pandemic\/616548\/<\/a>, retrieved 10\/1\/2020)<\/p>\n<h1 style=\"text-align: center;\">This Overlooked Variable Is the Key to the Pandemic<\/h1>\n<p style=\"text-align: center;\">It\u2019s not R.<\/p>\n<p style=\"text-align: center;\"><em>Zeynep Tufekci<\/em>, September 30, 2020<\/p>\n<p>\u2026 The definition of <em>k<\/em>, [the measure of its dispersion,] is a mouthful, but it\u2019s simply a way of asking whether a virus spreads in a steady manner or in big bursts, whereby one person infects many, all at once. After nine months of collecting epidemiological data, we know that this is an <em>overdispersed<\/em> pathogen, meaning that it tends to spread in clusters, but this knowledge has not yet fully entered our way of thinking about the pandemic\u2014or our preventive practices.<\/p>\n<p>The now-famed R0 (pronounced as \u201cr-naught\u201d) is an <em>average<\/em> measure of a pathogen\u2019s contagiousness, or the mean number of susceptible people expected to become infected after being exposed to a person with the disease. \u2026 Unfortunately, averages aren\u2019t always useful for understanding the distribution of a phenomenon, especially if it has widely varying behavior. \u2026 [Superspreading is] a pattern that\u2019s been observed many times since the pandemic begin, and that is similarly not captured by R. That\u2019s where the dispersion comes in.<\/p>\n<p>There are COVID-19 incidents in which a single person likely infected 80 percent or more of\u00a0 the people in the room in just a few hours. But, at other times, COVID-19 can be surprisingly much less contagious.\u00a0 Overdispersion and super-spreading of this virus are found in research across the globe. \u2026 Multiple studies from the beginning have suggested that as few as 10 to 20 percent of infected people may be responsible for as much as 80 to 90 percent of transmission, and that many people barely transmit it. \u2026<\/p>\n<p>We can think of disease patterns as leaning deterministic or stochastic: In the former, an outbreak\u2019s distribution is more linear and predictable; in the latter, randomness plays a much larger role and predictions are hard, if not impossible, to make. \u2026 Diseases like the flu are pretty nearly deterministic and R0 (while flawed) paints about the right picture (nearly impossible to stop until there\u2019s a vaccine). That\u2019s not necessarily the case with super-spreading diseases. \u2026<\/p>\n<p>In study after study, we see that super-spreading clusters of COVID-19 almost overwhelmingly occur in poorly ventilated, indoor environments where many people congregate over time\u2014weddings, churches, choirs, gyms, funerals, restaurants, and such\u2014especially when there is loud talking or singing without masks. For super-spreading events to occur, multiple things have to be happening at the same time, and the risk is not equal in every setting and activity, Muge Cevik, a clinical lecturer in infectious diseases and medical virology at the University of St. Andrews and a co-author of a recent extensive review of transmission conditions for COVID-19,\u00a0 told me.<\/p>\n<p>Cevik identifies \u201cprolonged contact, poor ventilation, [a] highly infectious person, [and] crowding\u201d as the key elements for a super-spreader event. Super-spreading can also occur indoors beyond the six-feet guideline, because SARS-CoV-2, the pathogen causing COVID-19, can travel through the air and accumulate, especially if ventilation is poor. Given that some people infect others before they show symptoms, or when they have very mild or even no symptoms, it\u2019s not always possible to know if we are highly infectious ourselves. We don\u2019t even know if there are more factors yet to be discovered that influence super-spreading. But we don\u2019t need to know all the sufficient factors that go into a super-spreading event to avoid what seems to be a necessary condition most of the time: many people, especially in a poorly ventilated indoor setting, and especially not wearing masks. As Natalie Dean, a biostatistician at the University of Florida, told me, given the huge numbers associated with these clusters, targeting them would be very effective in getting our transmission numbers down.<\/p>\n<p>Overdispersion should also inform our contact-tracing efforts. In fact, we may need to turn them upside down. Right now, many states and nations engage in what is called forward or prospective contact tracing. Once an infected person is identified, we try to find out with whom they interacted afterward so that we can warn, test, isolate, and quarantine these potential exposures. But that\u2019s not the only way to trace contacts. And, because of overdispersion, it\u2019s not necessarily where the most bang for the buck lies. Instead, in many cases, we should try to work backwards to see who first infected the subject.<\/p>\n<p>Because of overdispersion, most people will have been infected by someone who also infected other people, because only a small percentage of people infect many at a time, whereas most infect zero or maybe one person. \u2026<\/p>\n<p>Another significant consequence of overdispersion is that it highlights the importance of certain kinds of rapid, cheap tests. Consider the current dominant model of test and trace. In many places, health authorities try to trace and find forward contacts of an infected person: everyone they were in touch with since getting infected. They then try to test all of them with expensive, slow, but highly accurate PCR (polymerase chain reaction) tests. But that\u2019s not necessarily the best way when clusters are so important in spreading the disease. \u2026<\/p>\n<p>Meanwhile, researchers have shown\u00a0 that rapid tests that are very accurate for identifying people who do not have the disease, but not as good at identifying infected individuals, can help us contain this pandemic. \u2026 they are particularly valuable for cluster identification during an overdispersed one. \u2026<\/p>\n<p>In an overdispersed regime, identifying <em>transmission events<\/em> (someone infected someone else) is more important than identifying <em>infected<\/em> <em>individuals<\/em>.\u00a0 Consider an infected person and their 20 forward contacts\u2014people they met since they got infected. Let\u2019s say we test 10 of them with a cheap, rapid test and get our results back in an hour or two. \u2026 If everyone is negative, we can act as if nobody is infected, because the test is pretty good at finding negatives. However, the moment we find a few transmissions, we know we may have a super-spreader event, and we can tell all 20 people to assume they are positive and to self-isolate\u2014if there are one or two transmissions, there are likely more, exactly because of the clustering behavior. \u2026 [O]verdispersion also enhances the utility of other aggregate methods, such as wastewater testing, especially in congregate settings like dorms or nursing homes, allowing us to detect clusters without testing everyone. \u2026<\/p>\n<p>Unfortunately, until recently, many such cheap tests had been held up by regulatory agencies in the United States, partly because they were concerned with their relative lack of accuracy in identifying positive cases compared with PCR tests\u2014a worry that missed their population-level usefulness for this particular overdispersed pathogen.<\/p>\n\n\n<p><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Given the mathematical nature of this website I feel reluctantly impelled to address the coronavirus pandemic. The mathematics behind the spread of infection is basically the same exponential growth that I discussed in the \u201cMath and Religion\u201d post and has recently been explained by the ever-lucid Grant Sanderson at his 3Blue1Brown website. What I wish [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[65],"tags":[25,189,38,188,58],"class_list":["post-1505","post","type-post","status-publish","format-standard","hentry","category-other-voices","tag-3blue1brown","tag-coronavirus","tag-exponential-growth","tag-kevin-drum","tag-minute-physics"],"_links":{"self":[{"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/posts\/1505","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/comments?post=1505"}],"version-history":[{"count":13,"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/posts\/1505\/revisions"}],"predecessor-version":[{"id":1786,"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/posts\/1505\/revisions\/1786"}],"wp:attachment":[{"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/media?parent=1505"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/categories?post=1505"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/josmfs.net\/wordpress\/wp-json\/wp\/v2\/tags?post=1505"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}