COVID-19 Numeracy: More NYT

As I mentioned earlier today, The New York Times staff has published some excellent maps, charts and interactive simulations.

Their latest interactive simulation illustrates the effect of community mitigation (i.e., social distancing) for varying periods of time. I took snapshots at weekly intervals (social distancing for 7 days, 14 days, 21 days, etc.) [Click to enlarge.]

Effect of community mitigation for different time periods (NYT)

The curve in the lightest shade is the number of infected people. The middle shade curve is the number of people hospitalized and the darkest shade curve is the number of deaths. Each tick in the vertical scale is 10 million. Yes, millions. I recommend reading the article itself for assumptions and caveats on the mathematical model.

From this time series, it’s easy to see that two weeks (14 days) of community mitigation is wholly inadequate for suppressing the disease, let alone stopping it. To my eye, the shortest reasonable duration is 42 days (6 weeks).

Here is the latest NYT map (March 26, 2020) depicting the spread of COVID-19 across the United States. [Click to enlarge.]

Spread of COVID-19 (The New York Times, March 26, 2020)

Dedicated to Edward Tufte, who is still kicking around, thank goodness.

COVID-19 Numeracy: Geography

Updated March 26, 2020.

State-by-state

Governor Cuomo, in his briefing today, asked why New York State is leading the United States in positive COVID-19 cases. Let’s compare the top four states:

                     Total   Total   Percent    
State Tests Cases Positive Date
------------- ------ ------ -------- ----
New York 38,390 14,904 38.8% 3/24
Massachusetts 13,749 1,159 8.4% 3/24
California 27,650 2,102 7.6% 3/24
Washington 33,933 2,221 6.5% 3/23

Indeed, New York State has the highest percentage of positive tests — substantially higher. There must be a factor in play. Could there be a difference in testing (i.e., difference in the test kits? different test protocol?) Did SARS-CoV-2 hit New York State sooner? Is COVID-19 more prevalent due to the population density in the New York City region, as the Governor speculated?

Only a good epidemiological study could really measure the prevalence of SAR-CoV-2 in the community. As to timing, I don’t think we will ever know when SARS-CoV-2 hit our shores. New York City is a major global transportation hub. Even so, Washington and California have strong business and family ties across the Pacific and Pacific rim and have major transportation hubs (Los Angeles, San Francisco, San Diego, Seattle).

If I get a chance today, I’m going to look into our neighbor to the north, Vancouver. Vancouver is also quite connected with Asia and the Pacific rim.

The Washington Department of Health has had difficulties with its COVID-19 database. Although they posted a new number of positive cases, 2,469 as of March 24, the number of negative cases was not updated, showing the same total negative test count as of March 23.

The national picture

The New York Times has published many informative graphs and simulations during the COVID-19 crisis. The chart below shows the spread of SARS-CoV-2 (the virus) and COVID-19 (the disease caused by SARS-CoV-2) across the United States as of March 19, 2020.

Source: The New York Times

The chart illustrates how population density affects the number of COVID-19 cases. [Click to enlarge.] The biggest circles coincide with the major urban centers in the United States. On the west coast are Seattle-Tacoma-Washington, the San Francisco Bay area, and Los Angeles, for example.

The chart demonstrates why it isn’t appropriate to think numerically about the country as a whole. Sure, Americans want to know how the U.S. is doing overall. However, when it comes to actual COVID-19 containment and mitigation, it’s a local, regional problem. That’s why our county and municipal health departments are so important and local statistics are the most meaningful and practical.

All is not local, however. The urban centers are connected by highways and air routes, giving SARS-CoV-2 passage. As the disease plays out, the coastal urban centers with the greatest population and density will likely succumb first, followed by so-called second tier cities. When the disease subsides in major urban centers, the second tier will maintain active reservoirs of SARS-CoV-2. If major urban centers drop community mitigation (e.g., social distancing) too soon, the active pools in the second tier will cause COVID-19 flare-ups in any major urban center that drops its mitigration too early.

Major urban centers are more tightly connected to global business and commerce than the second tier. SARS-CoV-2 likely arrived much earlier than the second tier and SARS-CoV-2 has probably been circulating in the first tier population for a longer time.

Italy

Media outlets have drawn many comparisons between the Unitet States and Italy. Let’s consider a few facts about Italy from the CIA World Factbook.

Italy’s land area is 294,140 sq km (square kilometers), roughly the size of Arizona (294,207 sq km). Italy’s population is estimated to be 62,402,659 people (July 2020). California is the most populous state with 39,512,223 residents (July 2019). Overall, Italy has as many people as California and New York State combined in a smaller land mass (Arizona).

Italy’s Lombardy region is, arguably, the Italian region most severely affected by COVID-19. Lombardy is located in northern Italy with the major urban center Milan as its capital.

Lombardy has 10,078,012 people (August 2019). It is the most populated region in Italy. It’s population is roughly the size of Michigan (or North Carolina, or Georgia). Lombardy, by population, would be the tenth largest state in the U.S. Lombardy is 23,844 sq km, roughly the size of West Virginia. Lombardy has a population density of 420 people per sq km.

The city of Milan has 1.369 million people. (The World Fact Book claims 3.140 million people, but it depends on how one draws boundaries.) Milan has a population density of 7,572 people per sq km (2019). That’s roughly the same density as Huntington Park (Los Angeles metro) or Somerville (Boston metro). New York City has a population density of 10,431 people per sq km. [I taught in Somerville; it’s dense.]

Lombardy accounts for 68% of COVID-19 fatalities in Italy. 23% of Italians are 65 or older, making it the second oldest population in the world. (Japan is first.) Age and underlying health conditions are known COVID-19 vulnerabilities. High population density and more elderly people (with existing health issues) are a fatal combination of factors.

Temperature

One final observation because it involves maps. The Kinsa US Health Weather Map has received media attention including MSNBC’s The Rachel Maddow Show. (Hi, Rachel.)

The Health Weather Map is a visualization of seasonal illness linked to fever as measured by Kinsa Smart Thermometers. The Kinsa Thermometer connects to your device (and the Kinsa App) via Bluetooth. The Kinsa Thermometer costs about $20-$25 USD retail and requires the Kinsa App for Web communication.

Here is today’s map (March 26). [Click to enlarge.]

Kinsa US Health Weather Map (March 26, 2020)

I look at the Health Weather Map with interest, but take it with a grain of salt. There are two factors to consider.

  • Although the Kinsa Thermometer is inexpensive, it needs an expensive smart device in order to communicate results via the Web.
  • It’s not clear how people use the thermometer, i.e., do they take their temperature every day? When they feel sick?

Given the need for a smart device and Internet connection, I’m not sure if Kinsa users represent the population at large. Use of the device may bias results, too. If users take their temperature only when they are sick, then the map will show more high temperatures. Ah, well, I need to dig into their study. [Thank you, Kinsa, for posting your methodology.]

Stay apart, mate, and we’ll both be healthy — P.J. Drongowski

COVID-19 Numeracy: New York City

Yesterday, Governor Jay Inslee issued a full, statewide “stay at home” order for Washington State. Up to this point, home-stay was recommended and did not have the force of law. Now it does.

Thank goodness. When I first glanced at the percentage of new positive cases per day (Washington State), I felt optimistic. After reflection, I decided. “Yes, it’s good that the percentage is not increasing rapidly, but even no change is still bad.” Even a constant flow of new patients into the health care system will eventually overwhelm it. A linear curve, as opposed to an exponential curve, merely postpones the day when we exceed capacity.

Someone might say, “Oh, linear isn’t bad.” Yes, it is and here’s why. It’s not like an acute care patient arrives, is given some magic cure and is sent home. Mean duration of hospitalization in Chine was 12.8 days. [Please see “Clinical Characteristics of Coronavirus Disease 2019 in China” in the New England Journal of Medicine.] Even with a so-called linear curve, acute care patients are admitted and stay. Most hospitals are “right sized” during normal times by administrators interested in efficency. So, it doesn’t take long before capacity and staff are exhausted.

China had hospitalization rates ranging from 15% to 20%. In the Chinese epidemic, 41.3% of patients received oxygen therapy and 6.1% received mechanical ventilation.

There is only one advantage to a linear curve vs. exponential — it gives hospitals and other facilities an opportunity to ramp up for the influx of acute care patients. The availability of respirators and ventilators is a critical concern as well as temporary beds and infrastructure for convalescence and recovery.

We need to push the number of new COVID-19 cases downward through social distancing and good hygiene. I’m glad that Governor Inslee has issued a full stay-at-home order. Now, we all need to listen to it and act.

New York City

Governor Cuomo of New York was visibly more alarmed today. I pulled together a quick table of new case and new testing statistics for New York City (NYC).

I’m alarmed, too.

New York City (NYC) 
Cum New Cum New Percent
Date Positive Positive Tests Tests Ratio Hospital
------ -------- -------- ------ ------ ----- --------
Mar 20 4408 1939 14386 4852 40.0% 18%
Mar 21 6211 1803 19463 5077 35.5% 15%
Mar 22 9045 2832 26389 6927 40.1% 13%
Mar 23 12305 3260 33003 6614 49.3% 13%
Mar 24 14904 2599 38390 5387 48.2% --
Mar 25 17856 2952 44076 5686 51.9% 12%
Mar 26 21393 3537 51031 6955 50.9% 14%
Mar 27 25398 4005 57954 6923 57.9% 15%
Mar 28 29766 4368 65902 7948 55.0% --
Mar 29 33768 4002 73104 7202 55.7% 14%
Mar 30 37453 3685 -- -- -- 14%

24% of hospitalizations are ICU patients (March 23)
23% of hospitalizations are ICU patients (March 24)
24% of hospitalizations are ICU patients (March 26)
24% of hospitalizations are ICU patients (March 27)
24% of hospitalizations are ICU patients (March 29)
25% of hospitalizations are ICU patients (March 30)

The percentage of new positive cases per day (which takes the ever-varying number of tests into account) for Washington State is around 6.5%. The ratio for New York City is 6 times higher than Washington. (Update: Washington State calls this ratio “Positivity Rate.”)

I fear that NYC is already in the weeds. Folks, NYC is the largest financial center in the United States. So, even if all you care about is money, this is the wrong time to lift restrictions and to restart the American economy. Plain and simple, that’s crazy talk. And, immoral.

Stay distant and stay healthy — P.J. Drongowski

Update quoting CNBC News:

New York and New Jersey are seeing coronavirus attack rates at least five times higher than other parts of the country, a U.S. official in charge of the White House’s pandemic response efforts said Monday.

“The New York metro area of New Jersey, New York City, and parts of Long Island have an attack rate close to one in 1,000,” Dr. Deborah Birx, a physician and the White House coronavirus response coordinator, said at a press briefing Monday evening. The attack rate is the percentage of a population that gets the disease.

She said roughly 28% of the specimens submitted in that region have tested positive for COVID-19, while less than 8% have tested positive for the disease in the rest of the country.

New York is currently the hardest-hit state in the country, ahead of New Jersey, California and Washington state. New York City, alone, accounts for 12,305 of the 20,875 confirmed infections in the state as of Monday morning.

COVID-19 Numeracy: Comparability

% Positive New Cases Per Day

It’s time to show some numbers of my own!

Thanks to The Seattle Times and the Washington State Department of Health (DOH), I have built a table that summarizes COVID-19 case statistics since late February. Data before March 10 is somewhat spotty and lacks the number of tests per day, etc. The raw numbers are shown in the table below.

            Cum      Cum    Cum 
Date Positive Negative Tests New Cases New Tests Ratio Notes
------ -------- -------- ------- --------- --------- ----- --------
Feb 28 1
Feb 29 7
Mar 2 18
Mar 4 39 SnoEmerg
Mar 5 70 Pence
Mar 6 83
Mar 7 102
Mar 8 136
Mar 9 162 1149 1311
Mar 10 267 105
Mar 11 366 3037 3403 99 Ban>250
Mar 12 457 4350 4807 91 1404 6.48%
Mar 13 568 6001 6569 111 1762 6.30%
Mar 14 643 7122 7765 75 1196 6.27%
Mar 15 769 9451 10220 126 2455 5.13%
Mar 16 904 11582 12486 135 2266 5.96% HomeRec
Mar 17 1012 13117 14129 108 1643 6.57%
Mar 18 1187 15918 17105 175 2976 5.88%
Mar 19 1376 19336 20712 189 3607 5.24%
Mar 20 1524 21719 23243 148 2531 5.85%
Mar 21 1793 25328 27121 269 3878 6.94%
Mar 22 1996 28879 30875 203 3754 5.41% Everett
Mar 23 2221 31712 33933 225 3058 7.36%
Mar 24 2469 248 StayHome
Mar 25 2580 111
Mar 26 3207 43173 46380 627 12447 7.92%
Mar 27 3700 49015 52715 493 6335 7.78%
Mar 28 4300 54896 59196 600 6481 9.26%
Mar 29 4896 60566 65462 596 6266 9.51%
Mar 30 5515 DOHIssues
Mar 31 5984 68814 74798 1088 9336 11.65%
Apr 1 6585 72833 79418 601 4620 13.01%
Apr 2 6966 75633 82599 381 3181 11.98%
Apr 3 7591 80327 87918 625 5319 11.75%
Apr 4 7984 393 DOHIssues
Apr 5 8384 400
Apr 6 8682 298
Apr 7 9097 415
Apr 8 9609 512
Apr 9 9887 278
Apr 10 10224 337
Apr 11 10411 187
Apr 12 10538 127
Apr 13 10694 112160 122854 166
Apr 14 10783 113500 124283 89 1429 6.23%
Apr 15 11152 117748 128900 369 4617 7.99%
Apr 16 11445 120182 131627 293 2727 10.74%
Apr 17 11802 135706 357 4079 8.75%
Apr 18 11790 126852 138642 2936 DOHCorrection
Apr 19 12085 128926 141011 295 2369 12.45%
Apr 20 12282 132749 145031 197 4020 4.90%
Apr 21 12494 135459 147953 212 2922 7.26%
Apr 22 12753 140623 153376 259 5423 4.78%
Apr 23 12977 147347 160324 224 6948 3.22%
Apr 24 13319 157275 170594 342 10270 3.33%
Apr 25 13521 161956 175477 202 4883 4.14%
Apr 26 13686 165993 179679 165 4202 3.93%
Apr 27 13842 168673 182515 156 2836 5.50%
Apr 28 14070 173730 187800 228 5285 4.31%
Apr 29 14327 179654 193981 257 6181 4.16%
Apr 30 14637 184087 198724 310 4743 6.54%
May 1 15003 192312 207315 366 8591 4.26%
May 2 15185 196820 212005 182 4690 3.88%
May 3 15462 200858 216320 277 4315 6.42%
May 4 15594 203859 219453 132 3133 4.21%
May 5 15905 208908 224813 311 5360 5.80%
May 6 16231 214449 230680 326 5867 5.56%
May 7 16388 219447 235835 157 5155 3.05%
May 8 16674 226315 242989 286 7154 4.00%
May 9 16891 231984 248875 307 5886 5.22%
May 9 16891 231984 248875 307 5886 5.22%
May 10 17122 234986 252108 231 3233 7.15%
May 11 17330 238991 256321 208 4213 4.94%
May 12 17512 243568 261080 182 4759 3.82%
May 13 17773 250158 267931 261 6851 3.81%
May 14 17951 255352 273303 178 5372 3.31%
May 15 18288 262705 280993 337 7690 4.38%
May 16 18433 266810 285243 145 4250 3.41%
May 17 18611 270524 289135 178 3892 4.57%
May 18 18811 274309 293120 200 3985 5.02%
May 19 18971 278971 297942 160 4822 3.32%
May 20 19117 283810 302927 146 4985 2.93%
May 21 19265 289093 308358 148 5431 2.73%
May 22 19585 296691 316276 320 7918 4.04%
May 23 19828 306765 326593 243 10317 2.36%
May 24 20065 310533 330598 237 4005 5.92%
May 25 20181 312610 332791 116 2193 5.29% George Floyd
May 26 20406 315395 335801 225 3010 7.48%
May 27 20764 322327 343091 358 7290 4.91%
May 28 21071 327162 348233 307 5142 5.97%
May 29 21349 333005 354354 278 6121 4.54% SEA protests
May 30 21702 339197 360899 353 6545 5.39%
May 31 21977 343295 365272 275 4373 6.29%
Jun 1 22157 346642 368799 180 3527 5.10%
Jun 2 22484 354843 377327 327 8528 3.83%
Jun 3 22729 360858 383587 245 6260 3.91%
Jun 4 22993 367870 390863 264 7276 3.63%
Jun 5 23442 377146 400588 449 9725 4.62%
Jun 6 23729 381327 405056 287 4468 6.42%
Jun 7 24041 386249 410290 312 5234 5.96%
Jun 8 24354 390700 415054 313 4764 6.57%
Jun 9 24642 --- --- 288 --- ---
Jun 10 24779 400433 425212 137 10158 4.18%
Jun 11 25171 414691 439862 392 14550 2.69%
Jun 12 25538 424608 450146 367 10284 3.57%
Jun 13 25834 436768 462602 296 12456 2.38%
Jun 14 26158 445107 471265 324 8663 3.74%
Jun 15 26531 453495 480026 373 8761 4.26%
Jun 16 26784 408232 435016 253 --- --- DOH no sero
Jun 17 27192 419707 446899 408 11883 3.43%
Jun 18 27601 428340 455941 409 9042 4.52%
Jun 19 28225 437844 466069 624 10128 6.16%
Jun 20 28680 446258 474938 455 8869 5.13%
Jun 21 28870 448334 477204 190 2266 8.38%
Jun 22 29386 457673 487059 516 9855 5.24%
Jun 23 29869 465629 495498 483 8439 5.72%
Jun 24 30367 475428 505795 498 10297 4.84%
Jun 25 30855 483573 514428 488 8633 5.65%
Jun 26 31404 494398 525802 549 11374 4.83%
Jun 27 31752 502691 534443 348 8641 4.03%
Jun 28 32253 515967 548220 501 13777 3.64%
Jun 29 32824 524451 557275 571 9055 6.31%
Jun 30 33435 538529 571964 611 14689 4.16%
Jul 1 34151 550838 584989 716 13025 5.50%
Jul 2 34778 565197 599975 627 14986 4.18%
Jul 3 35247 572029 607276 469 7301 6.42%
Jul 4 35898 576808 612706 651 5430 11.99%
Jul 5 36985 592271 629256 1087 16550 6.57%
Jul 6 37420 598104 635524 435 6268 6.94%
Jul 7 37941 607131 645072 521 9548 5.46%
Jul 8 38581 621749 660330 640 15258 4.19%
Jul 9 39218 629248 668466 637 8136 7.83%
Jul 10 ---
Jul 11 40656 645349 686005 1438 17539 8.20%
Jul 12 41757 666517 708274 1101 22269 4.94%
Jul 13 42304 675930 718234 547 9960 5.49%
Jul 14 43046 690840 733886 742 15652 4.74%
Jul 15 44313 708861 753174 1267 19288 6.57% Yakima backlog
Jul 16 45067 722590 767657 754 14483 5.21%
Jul 17 46026 745760 791786 959 24129 3.97%
Jul 18 46946 762393 809339 920 17553 5.24%
Jul 19 47743 778611 826354 797 17015 4.68%
Jul 20 48575 792609 841184 832 14830 5.51%
Jul 21 49247 805905 855152 672 13968 4.81%
Jul 22 50009 820754 870763 762 15611 4.88%
Jul 23 50824 833158 883982 815 13219 6.17%
Jul 24 51849 851825 903674 1025 19692 5.21%
Jul 25 52635 866712 919347 786 15673 5.01%
Jul 26 53321 879983 933304 686 13957 4.92%
Jul 27 54205 891029 945234 884 11930 7.41%
Jul 28 54985 903322 958307 780 13073 5.97%
Jul 29 55803 917851 973654 818 15347 5.33%
Jul 30 --- DOH outage
Jul 31 57541 943987 1001528 1738 27874 6.24%
Aug 1 58173 950107 1008280 632 6752 9.36%
Aug 2 58715 --- 542 DOH issues
Aug 3 59379 --- 664
Aug 4 60084 --- 705
Aug 5 60917 --- 833

I have tried to make this table as accurate as possible. Please forgive any transcription errors. The Washington Department of Health has suffered technical issues, so there are gaps in the data. Their database has been overwhelmed and they have had trouble tallying negative test results and the number of tests.

Don’t forget about the University of Washington Virology Laboratory COVID-19 dashboard as another reliable data source.

The State DOH reports the cumulative number of positive cases (tests) negative test results daily — the second and third columns in the table. The fourth columb, cumulative number of tests, is the sum of the cumulative positive and negative results each day.

From these numbers, we can compute the number of new positive cases and the number of new tests for each 24 hour period:

    Tests.Today = Positive.Today + Negative.Today 
New cases = Positive.Today - Positive.Yesterday
New tests = Tests.Today - Tests.Yesterday
Ratio = (New cases / New Tests) * 100.0

The seventh column, ratio, is the percentage of positive new cases per day. I believe that this percentage is an appropriate measure for tracking COVID-19 over time. The ratio normalizes the number of positive cases for the number of tests performed. Thus, it makes daily results comparable even as our test capacity and number of daily tests grow.

Update: Washington State calls this ratio/percentage “Positivity Rate.”

The bar chart below is a quick and dirty plot of the daily percentage of positive cases over the last ten days. Be sure to check the vertical scale when interpreting the chart!

Percentage of positive COVID-19 tests per day (Washington State)

I hope this trend declines. Please remember that any measurements taken today show what COVID-19 was doing 10 to 14 days ago. The delay is due to the time required for COVID-19 incubation, worsening of symptoms, clinical presentation and testing.

The last column in the table, Notes, are my mental notes marking significant events along the COVID-19 timeline in our area. The events are:

  • Snohomish County (my home) declares a state of emergency.
  • Vice President Pence visits the State of Washington.
  • Governor Inslee bans groups of more than 250 people.
  • The State advises people to stay at home, especially the elderly.
  • The City of Everett (a stone’s throw away from home) imposes a stay at home order.

Up to today, Governor Inslee has not imposed a blanket “stay at home” order. Suggestions made by the state are similar to a stay at home order, but do not have the force of law.

Don’t do this

Here is the egregious chart of the day by way of Cleveland.com. I still check in with my original hometown. Hope the Cavs can hit the court again.

To be fair, the author did put a few caveats about number of tests, etc. in the article. The caveats are lost in the TL:DR age. Please don’t publish meaningless charts and comparisons.

If I were teaching a math or science class today, I think I would use this chart in class about critical thinking and interpetation of the graphs thrown at us by the media.

No one (CNN, etc.) should compare state-by-state number at this point in time. There isn’t a standard for reporting and test capacity is, literally, all over the map.

Stay distant and stay healthy — pj

COVID-19 Numeracy: Charts

New York, New York

Life has dropped into a bit of a routine. Get up, start the coffee pot and watch the briefing by Andrew Cuomo, Governor of New York. The Guv can be a little tedious and pedantic (well-known characteristics to New Yorkers), but he is a wonderful communicator. First, he speaks slowly and in simple, clear sentences — good for people for which English is not their first language.

He also has been presenting clear statistics on testing and the number of new COVID-19 cases. His interpetation and remarks are correct. The more you test, the more you will find the disease. One of his charts today had three columns: locality (county and entire state), cumulative COVID-19 cases to date and the number of new cases identified for the preceding 24 hour period. At the very least, the Governor needs to add a fourth column, the number of tests performed during the preceding 24 hour period. The fourth column would show the concurrent growth of new cases and tests.

I also heard the Governor struggling, just a bit, to put the number of new cases in context with the growth in testing. [New York has done an excellent job of ramping up its testing, BTW.] That’s where the ratio of new cases to tests comes into play:

Percentage of new positive cases = ((New cases) / (New tests)) * 100 

I get that most people do not think in terms of ratios between 0.0 and 1.0, so this measure can and should be expressed as a percentage.

The percentage (ratio) accounts for the growth in testing. I believe that this percentage (ratio) is a better tracking metric than the raw number of new cases alone. It is certainly better than the cumulative number of COVID-19 cases since the beginning of the epidemic. The cumulative number expresses the depth of suffering, but it does not adequately convey the growth or decline of the disease. Given that we will experience recurring waves of COVID-19 (if we suppress this disease successfully), the cumulative number of cases alone will obscure the rise and ebb of COVID-19.

I was encouraged to see New York using a metric similar to my COVID-19 Key Performance Indicator (KPI). New York officials are monitoring the number of all acute care patients. This is the most relevant number with respect to the capacity of the health care system. Please note that acute care includes all of the “every day” acute care patients (heart attacks, strokes, etc.) and COVID-19 patients. All of these patients need to be hospitalized and all of them put a load on the health care system.

If I have to add a public service announcement (PSA) of my own, please don’t drive under the influence and please don’t drive distracted (i.e., texting). Hospital emergency rooms don’t need additional patients. And be careful with that axe, Eugene.

The Seattle Times

I’d like to give a shout-out to my (new) hometown newspaper, The Seattle Times. They publish a concise daily summary of new and cumulative test results (March 19, 2020 below; click to enlarge).

The summary has the daily number of tests performed, positive cases identified and negative cases. Charts break out positive cases and deaths by age. Tests per day has its own chart and shows the growing number of tests. There is a chart showing the cumulative number of positive cases and deaths over time since February 28.

All good work. The raw figures come from the Washington State Department of Health, who should be commended for their transparency.

My one suggestion is to add a chart showing the test results over time, like the graph published by the University of Washington Virology Laboratory (below).

UW Virology

The UW Virology Laboratory upped their game and publish daily results in the UW Virology COVID-19 Dashboard. [Click images to enlarge.]

University of Washing Virology Laboratory COVID-19 Dashboard

This is a screenshot of the dashboard for March 19, 2020. It’s easy to see the number of tests (“Sample Count”) performed each day, the number of positive and negative cases identified, and the number of inconclusive results. You can hover the cursor over a bar segment and get the raw count for the chosen segment.

Very nice work. The trend of positive cases is easy to see. The chart shows the raw number of positive cases, so interpretation is still biased by the number of tests performed. A separate chart is needed which plots the percentage of new positive cases over time.

One additional suggestion for UW Medicine — Please explain “Inconclusive”.

As I mentioned yesterday, each day’s results are a glimpse into the past. COVID-19 has a 2 to 14 day incubation period according to the CDC. Other medical experts estimate the infection to test delay at 10 to 14 days (infection, incubation, symptoms, clinical presentation, test). Even though the trend of new COVID-19 cases is relatively flat, we won’t really know for another 14 to 21 days if we successfully flattened the curve and delayed onset of the disease in the community.

Keep the faith and keep healthy — P.J. Drongowski

COVID-19 Numeracy: Time

On March 16, President Trump warned that the COVID-19 outbreak could last until July or August 2020. This was the first official recognition that the outbreak could take months to quell, not a few weeks. On March 18, The New York Times and other media outlets published the U.S. Government COVID-19 Response Plan. The plan states that a moderately severe outbreak could last 12 to 18 months with multiple “waves.”

The duration and potential number of waves may shock some. However, the timing and waves of resurgence/subsidence are fairly typical for epidemics and other phenomena in nature. Consider, for example, predator-prey population models. Predators thrive and reproduce when prey are plentiful. When the food supply dwindles, predators die off and the predator population dwindles. With fewer predators, prey survive and multiply. The cycle of growth and decline repeats.

Viruses need susceptible, infectable hosts in order to reproduce. Thanks to bodily immunilogical response, human hosts may develop an immunity as a result of infection and recovery. Immunity may be limited to the specific infecting strain of the virus. Unfortunately, viruses are adaptable and a strain may mutate into a new different strain which is not suppressed by immunity to the original strain.

The specific behavior of SARS-CoV-2 — the virus which causes the disease COVID-10 — is not fully known, including human immunilogical response to SARS-CoV-2. Citing Dr. Anthony Fauci, no scientific study to date has demonstrated acquired immunity to SARS-CoV-2. Dr. Fauci and other experts believe that acquired immunity is likely.

Timeline

The U.S. Government COVID-19 Response Plan states, “Late December 2019, authorities from the People’s Republic of China (PRC) announced a possible epidemic of pneumonia of unknown etiology centralized on a local large seafood and live animal market in Wuhan, China. Estimated case onset was early December.”

Given the lack of Chinese transparency, slow governmental response and the speed/extent of modern jet travel, the SARS-CoV-2 virus has probably been circulating in Snohomish and King County Washington for quite some time. The report states, “The first U.S. case of COVID-19 was confirmed in Washington State on January 20 and was travel-related.” This patient landed from Wuhan, China on January 15 without symptoms and eventually presented on January 19 with symptoms of pneumonia.

The first case of community spread in the U.S. (California) was confirmed by the CDC on February 26, 2020. The first case of community spread in Snohomish County, WA was confirmed on February 28. Thus, SARS-CoV-2 circulated in the Snohomish and King County area for nearly six weeks. Detection of community spread is a red flag for a developing epidemic.

On March 9, the Washington Department of Health confirmed 162 total cases including 22 deaths. Most of these deaths occurred at Life Care Center of Kirkland, a nursing home linked to at least 54 confirmed cases as of March 9. King and Snohomish Counties had 116 and 37 total cases, respectively, as of March 9. If you are interested in the details, the March 18, 2020 CDC MMWR is short and to the point.

Washington State immediately announced and slowly ramped up community mitigation. Governor Inslee issued the most stringent order on March 16. As of today, March 20, the state government has not issued a “stay/shelter at home” order like California.

State of Washington Social Distancing Summary

Washington’s new cases

The timeline is important in understanding current new case data and in knowing what to look for in the coming two to three weeks.

The University of Washington Virology Laboratory published a graph of its testing results as of Friday, March 20, 2020 (below).

UW Virology COVID-19 test results (March 20, 2020)

On a positive note (literally), the number of new COVID-19 cases detected by UW Virology is relatively flat. We are now 5 days into the most stringent social distancing rules and practices (March 20).

Although the results are promising, this is not the time to let up.

This graph is a snapshot of the past. According to the CDC, “The following symptoms may appear 2 to 14 days after exposure: fever, cough, shortness of breath.” COVID-19 tests are (effectively) rationed. A patient must meet testing criteria such as COVID-19 symptoms, contact with a known COVID-19 case, or travel from a COVID-19 hot-spot. Thus, the graph shows cases that are 2 to 14 days after infection.

Social distancing in Washington State did not “get serious” until March 16. The next three weeks (21 days) will tell us if social distancing in Washington State has really been effective in slowing the spread of COVID-19.

The curves (again)

Quoting the San Francisco Chronicle, “More than half of Californians could be infected with the new coronavirus over an eight-week period if nothing is done to stop it, according to a projection from Gov. Gavin Newsom’s office.”

Let’s hope this statement motivates California’s to adopt social distancing, hygiene and other measures suggested by health autorities.

Gov. Newsome is referring, of course to the orange curve in the now well-known CDC “flattening the curve” chart. If Californians (or the citizens of any other region or city) do not practice social distancing, then new cases will grow exponentially and hospitals, doctors and nurses will be overwhelmed.


Source: CDC, Drew Harris (Connie Hanzhang Jin/NPR)

Why “eight weeks” and why “more than half?”

Epidemiologists have studied disease for a very long time. And, like any other science, they have built mathematical models to predict the growth and decline of epidemics. The SIR model (Susceptible, Infectious, Recovered model) is one of the simplest — and classic — models. Typically, these models are driven by parameters such as the reproduction number, R0, pronounced “R naught”. R0 is the average number of people who catch the disease from one contagious person.

  • If R0 is less than one, then the disease will die out.
  • If R0 is one, the disease remains alive, but does not break out into an epidemic.
  • If R0 is greater than one, then the disease will spread, possibly growing into an epidemic.

The epidemiologists advising Gov. Newsome have plugged different values of parameters like R0 into their models, run the simulation and assessed the resulting statistics.

“Flattening the curve” essentially means reducing R0 as much as possible.

Here’s the rub which may have escaped some folks. The same number of people will catch the disease (COVID-19) in both the orange unmitigated scenario and the blue mitigated (social distancing) scenario. For the mathematically inclined, the area under the orange curve is equal to the area under the blue curve. The onset of the disease is delayed for some people in order to reduce the strain on the health care system.

The long haul

I started this blog post with statements that more realistically indicate the duration of our battle with COVID-19.

Ideally, we could like to control the growth of COVID-19 (the ramp up) and then stomp it out of existence (the ramp down). Unfortunately, real-world systems tend to behave more like the predator/prey model.

Viruses are opportunistic and adaptible. As long as SARS-CoV-2 circulates in the community, it will find susceptible individuals and infect them. Therefore, as soon as we reduce or abandon social distancing, COVID-19 will have a resurgence. The community will always be in danger of a resurgence as long as SARS-CoV-2 can find enough susceptible people to infect and thrive.

In the opinion of many health experts, our bodies likely acquire an immunity to SARS-CoV-2 as we fight and recover from the disease. At some point, when the so-called “herd immunity” is high enough, the disease will have trouble sustaining itself. We also anticipate a vaccination program against SARS-CoV-2 which will add to the herd immunity. [Yeah, no one likes to be thought of as a “herd”.]

Viruses are composed of proteins and nucleic acid. Without deep-diving molecular biology, viruses contain a form of genetic material: RNA, DNA or possibly both. Viruses usually depend upon a host cell for replication. Short story: viruses can steal genetic material anywhere in this process and mutate. Some mutations are bad and the mutants die out. Some mutations are good and the mutants thrive. That’s adaptation, that’s evolution.

SARS-CoV-2 could mutate. The mutated strain could thrive and grow because acquired immunity and/or the vaccine are specific to the original SARS-CoV-2 strain. Mutation could bring about a resurgence of COVID-19 although the underlying pathogen would be genetically different.

Realistically, it will be difficult to maintain stringent social distancing. We want to be together socially and we will definitely want to get our economy back to full strength, too. In the long haul, we should expect to see cycles of COVID-19 resurgence and decline. We may need to re-impose social distancing to keep COVID-19 in check.

Stay healthy and stay distant — P.J. Drongowski

COVID-19 Numeracy: Testing

Current situation: 18 March 2020

The United States continues to struggle with COVID-19 testing. Communities do not have sufficient COVID-19 test kits and they are, effectively, rationing tests. Generally, a patient must satisfy several criteria before a COVID-19 test is authorized and performed. Example criteria include cough, fever (high temperature), shortness of breadth, exposure to an infected individual, or travel to a known COVID-19 hot-spot. A patient may be tested first for influenza, receiving a COVID-19 test after influenza is ruled out, i.e., the patient fails the flu test.

Even in the face of shortage, the number of COVID-19 tests performed each day is steadily increasing albeit slowly.

With the current situation in mind, let’s look at few metrics.

Media reports

The most frequently reported statistics are the cumulative number of confirmed COVID-19 cases and the number of deaths due to coronavirus. A few news organizations report the number of new confirmed cases for the preceeding 24 hour period (AKA “per day”).

Most numerate people realize that “If you test, you will find the disease.” If you test more, then the raw number of confirmed cases will also rise.

Very few news organizations report the total number of tests performed per day as well as the number of new confirmed cases per day. Lacking the number of tests performed for each day, one cannot put the number of new confirmed cases in context. Nor can an investigator determine a trend over time, i.e., , is community mitigation (e.g., social distancing, etc.) controlling the spread of COVID-19?

We should insist on receiving three per-day metrics: the number of tests performed that day, the number of new confirmed cases and the number of negative cases. The number of tests performed will tell us if COVID-19 test capacity has improved or not. Then, I recommend computing, tracking and comparing the ratio of new per-day COVID-19 cases (confirmed positives) divided by the number of tests performed that day. The ratio better indicates the growth (or decline) of the disease in the community.

The Washington State Department of Health does report a daily break down. Their page breaks out confirmed cases by county, by age and by sex. The Department of Health also provides testing information. Unfortunately, the number of individuals tested, positive and negative are cumulative. As of March 17, the test figures are:

  • Positive: 1,012
  • Negative: 13,117
  • Individuals tested: 14,129 (Positive+Negative)
  • Percent positive: 1,012 / 14,129 = 7.7%

Please give us nonaggregated results for each day!

Beware sampling bias

One might be inclined to interpret the daily percentage of new confirmed cases as the rate of incidence of the disease in the community (the population at large). Unfortunately, the pre-testing criteria guarantee a biased sample. We are testing patients for COVID-19 only if they have symptoms, have a certain risk factor and/or do not have the flu. The ratio may overstate the incidence of disease or, worse, understate the incidence by missing asymptomatic, infected individuals.

The lack of COVID-19 test capacity leaves public health officials in a terrible bind. One would like to randomly sample (test) the population at large in order to estimate the number of infected people (symptomatic and asymptomatic) in the community. This is akin to taking a political poll. We need to select randomly (say, 1,600 people) from the community and test the selected individuals for COVID-19. We estimate the total number of infected people by multiplying the proportion of infected people as measured by sampling (“the poll”) times the overall population size.

One needs the right kind of test, of course. We want to measure current, active infections, not past infections that have resolved (recovered). This leads us (almost) to the classic SIR model (susceptible, infected, recovered) for epidemics — yet another topic! Also, time to ask the Google about “real-time reverse transcription – polymerase chain reaction (RT-PCR)” and COVID-19 testing. Who said self-isolation was uninteresting? 🙂

From statistics theory, 1,600 tests has a margin of error of plus/minus 2.45%. 1,600 tests is a modest number and, unfortunately, we can barely perform medical diagnostic testing let alone conduct a necessary statistical study at this time. Ideally, we would have sufficient capacity to conduct “tracking polls” to determine the overall trend, i.e., are there more (or less) infected people today than last week? Our officials are truly flying blind at the worst possible time.

Recent studies suggest that asymptomatic people are a significant source of new infections. That’s why public health officials are saying, “Act and behave as though you are infected.” We need tests. We need tools.

Stay apart and say healthy — P.J. Drongowski

COVID-19 Numeracy: NYT curve flattening by country

The New York Times published a collection of graphs comparing new COVID-19 cases over time by country. The article asks the question “Which Country Has Flattened the Curve for the Coronavirus?”

If you are interested in COVID-19 numeracy, I suggest reading the comments section. There are good insights there.

Here is my own response.

Thank you for posting your study. I started blogging about “COVID-19 numeracy” in response to the rather poor way media outlets have portrayed the disease numerically. We are in this fight for the long haul (12 to 18 months) and need metrics that guide our actions.

Please listen to the recommendations in this comment section. Many of us have spent years in measurement and statistics.

Raw numbers, e.g., the number of (new) confirmed cases, are not always meaningful or useful. We know that the number of new confirmed cases per-day will rise dramatically as testing increases. The number of new confirmed cases needs to be “normalized” against the number of tests performed. I suggest tracking the ratio of new confirmed cases divided by the number of tests.

I strongly agree that we need to track the progress of COVID-19 on a daily basis. (A seven day moving average is a good idea.) The total (cumulative) number of cases/deaths — as tracked and reported by most media outlets — will not be useful 2, 3, 4 months into the crisis, especially when there will be “waves” of resurgence and subsidence. We need to understand the dynamics of the pandemic.

Health authorities need to report the number of tests per day, the number of positive cases for that day and the number of negative cases. The raw number of tests per day will tell us if authorities are meeting their commitment to increase the number of tests and allow us to compute ratios, etc.

Should we ever get to the point of spare testing capacity beyond diagnosis, we need to conduct periodic community studies, something akin to a political poll. Take a random sample of the community and determine the number of symptomatic and asymptomatic cases (by age, by sex, etc.) Such polling will allow us to track the actual infection rate in the population at large.

COVID-19: KPI

First off, I would like to say how sad I am for those who have lost their lives to COVID-19, offering my support and empathy to their loved ones. We live in Sonomish County, Washington, not very far from King County and Kirkland. We’ve already seen the devastation which COVID-19 can wreak on care facilities for elderly people. This disease is all too real.

By now, you’ve heard the phrase “flattening the curve” and have probably seen the graph from the CDC (below). Our global goal is to slow the spread of COVID-19 through the population such that severe cases do not overwhelm the health care system. In a nutshell, the health care system has a fixed number of beds, doctors, nurses, caretakers, respirator, ventilators, etc. When capacity is exhausted, the system cannot treat all incoming patients (COVID-19 plus the regular, on-going stream of emergency situations like heart attacks, strokes, etc.) and, quite frankly, people will die.


Source: CDC, Drew Harris (Connie Hanzhang Jin/NPR)

The graph has done a good job of educating us as to the need for social distancing, hygiene, and other measures which slow the spread of disease.

It also suggests a metric — a key performance indicator (KPI) — which can tell us how well we are doing. A KPI is a measurable value that demonstrates how effectively an organization is achieving its objective(s). It’s important to note that a KPI can and should be applied at different levels in the organization. We also need to know how the KPI changes over time.

Let’s consider a KPI which measures the number of all critical care patients above or below the capacity of the health care system (below). This delta tells us if we are successfully suppressing the spread of COVID-19 or not. We must measure all critical care patients because they all are vying for the same medical personnel, beds and equipment. If the KPI is positive, i.e., the number of critical care cases exceeds capacity, then we are failing. If the KPI is negative, then we are succeeding.

COVID-19 Key Performance Indicator (KPI)

A national KPI value is only somewhat useful. The KPI needs to be measured at the state level and regional (metro) level. State-level values are useful only for small states such as Rhode Island with one major population center. Regional-level values are more useful in large states like California or Washington with two or more major population (and health care) centers. The U.S. is a very large country and health care capacity in West Virginia, for example, is not available to residents in San Francisco. Thus, regional numerical break out is required.

We need to track the KPI over time. The trend will tell us if we are successfully supressing the spread of COVID-19 or not. Tracking by region over time will tell us if hot spots are cooling off or if new hot spots are developing.

Why am I proposing this KPI, especially now? As a people, we need to walk a fine line between actionable concern and fearful panic.

I see tables and graphics in the media which tally the total number of confirmed cases and the total number of deaths. Yes, we mourn the loss of our neighbors. If we have even an ounce of humanity, can we not? I agree that such figures convey the sense of urgency needed to motivate new behavior that slows the spread of disease.

However, the total number of confirmed cases and fatalities alone are flawed measures for decision making. With respect to confirmed cases, the number can only go up — drastically — as we ramp up testing. Media need to at least report the number of people tested each day along with the number of negative results as well as the number of new confirmed cases. Even then, perhaps only the ratio of new confirmed cases to the number of new tests is truly meaningful.

Just to be clear, “each day” means the number of people tested, confirmed positive and negative in a specific 24 hour period, not an accumulated tally to date.

Aggregated and accumulated measurements are not practical and may be unnecessarily misleading. News media please take notice.

We also need to break out and track new critical COVID-19 cases which require hospitalization. These are the people who need the health care system. By tracking this number, we will see if the load on the health care system is trending upward toward catastrophe or trending downward successfully. As much as I love my younger and/or healthier brothers and sisters, if they are successfully recovering in self-isolation, then they are not loading the health care system with negative implications, and possibly death, for critical care patients.

I apologize to anyone who may feel offended by my frank discusion. I’m trying to come to grips with all of the information thrown at me by media outlets. Hopefully, you will find this approach to be useful. Even if you do not adopt my proposal, I hope to have started a practical discussion.

Wishing you safe passage — P.J. Drongowski

Back In The Day

A Keyboard Corner member asked what people did for keyboard amplification before PA. Man, that question really kicked off some memories.

Back in the day (1966), I played a Farfisa Mini Compact Deluxe through an Ampeg SB-12 bass amp. It was all I could afford. Mom and Dad lent me the money and I mopped floors at the local donut shop to pay them back. The shop had a wooden floor that was impregnated with grease. I still can’t face donuts to this day. 🙂

Farfisa Mini Compact Deluxe

While packing for the house move, I found an original Farfisa brochure from the 60s era.

The Portaflex was pretty cool with its flip top. The amp was mounted to a covered board which acted as a base for the head and, when flipped over, it became the cabinet cover. The clamps held the base/cover board in place and did double-duty as the speaker connections to the cabinet.

Ampeg SB-12 bass amplifier

The SB-12 had a 12″ Jensen speaker powered by a 25 Watt tube amp. It weighed 47 pounds — the first of a long line of heavy schleps.

Being the 1960s, of course, that wasn’t enough. Since I started playing with electronics and DIY at an early age, I tried my hand at an extension cabinet. I somehow came into a 15″ JBL speaker with a small tear in the cone. Impedance be damned, I just hooked it up in parallel with the Jensen via the clamps. Before building a cabinet, I would carry the JBL around in a suitcase which doubled as a “cabinet.” (!) The tear in the cone lent more bad-itude. [Why fizzy digital distortion doesn’t cut the mustard.]

In the R&B band, the guitars and vocals went through matching Ampeg guitar amps (probably Gemini’s). Only the top local bands could really afford PA for vocals (typically Fenders). Nobody put instruments through PA. The bass player had a Guild ThunderBass amp with that funky head. The bass player was quite good and laid down decent grooves. Can’t remember too much about the psychedelic band…

My failed experiments at extension and PA speaker cabs wound up as end-of-gig props. When we saw The Who trash their gear, we thought “What the heck!” I’d pull one of the legs off of the Farf and ram it through one of the prop cabinets.

My dream rig would have been a Vox Continental through a Fender Twin Reverb or Fender Super Reverb. I copped the Fender tilt-back idea and built a tilt-back stand for the SB-12. That got the speaker pointing up toward my ears.

Both the Connie and the Fenders were out of my financial reach. It took me three years to pay back my folks. By then, I had to sell the whole rig in order to make the college tuition nut. Given the rigors of college math, physics and computer science, it was the end of playing for quite a while. I can’t believe how much a vintage SB-12 fetches on the market these days!

The Farfisa Mini Compact Deluxe left me with no delusions about 1960s electronics. I tried tuning the F# oscillator and bunged the tuning coil. That was an unnecessary repair expense. That’s why I’m happy as a clam to play the Yamaha Reface YC today. The YC does a good job nailing the Farfisa and Vox.

A few other memories stick out like playing music fairs on stage/demo gear with the psychedelic band. One stage was incredibly small and I had a horn driver literally right in my ear. We played Doors, Steppenwolf, Vanilla Fudge, etc. at phenomenally loud volume, attracting every biker within earshot. They loved us. I think I still suffer hearing loss from those jobs.

Before signing off, I want to plug “Classic Keys: Keyboard Sounds That Launched Rock Music” by Alan Lenhoff and David Robertson. I received a copy yesterday and gave it a quick browse. The photography is excellent and the example gear is in tip-top shape. The book is long on history — less on playing technique and artistis, so some may be disappointed. Discount and used copies are coming onto the market and you may be able to save a few bucks if you can’t pony up the full $60USD. Recommended.

Copyright © 2020 Paul J. Drongowski