Image default
Home » How to Think Like an Epidemiologist

How to Think Like an Epidemiologist

There is a statistician’s rejoinder — sometimes offered as wry criticism, sometimes as honest advice — that could hardly be a better motto for our times: “Update your priors!”

In stats lingo, “priors” are your prior knowledge and beliefs, inevitably fuzzy and uncertain, before seeing evidence. Evidence prompts an updating; and then more evidence prompts further updating, so forth and so on. This iterative process hones greater certainty and generates a coherent accumulation of knowledge.

In the early pandemic era, for instance, airborne transmission of Covid-19 was not considered likely, but in early July the World Health Organization, with mounting scientific evidence, conceded that it is a factor, especially indoors. The W.H.O. updated its priors, and changed its advice.

This is the heart of Bayesian analysis, named after Thomas Bayes, an 18th-century Presbyterian minister who did math on the side. It captures uncertainty in terms of probability: Bayes’s theorem, or rule, is a device for rationally updating your prior beliefs and uncertainties based on observed evidence.

Reverend Bayes set out his ideas in “An Essay Toward Solving a Problem in the Doctrine of Chances,” published posthumously in 1763; it was refined by the preacher and mathematician Richard Price and included Bayes’s theorem. A couple of centuries later, Bayesian frameworks and methods, powered by computation, are at the heart of various models in epidemiology and other scientific fields

As Marc Lipsitch, an infectious disease epidemiologist at Harvard, noted on Twitter, Bayesian reasoning comes awfully close to his working definition of rationality. “As we learn more, our beliefs should change,” Dr. Lipsitch said in an interview. “One extreme is to decide what you think and be impervious to new information. Another extreme is to over-privilege the last thing you learned. In rough terms, Bayesian reasoning is a principled way to integrate what you previously thought with what you have learned and come to a conclusion that incorporates them both, giving them appropriate weights.”

With a new disease like Covid-19 and all the uncertainties it brings, there is intense interest in nailing down the parameters for models: What is the basic reproduction number, the rate at which new cases arise? How deadly is it? What is the infection fatality rate, the proportion of people with the virus that it kills?

But there is little point in trying to establish fixed numbers, said Natalie Dean, an assistant professor of biostatistics at the University of Florida.

“We should be less focused on finding the single ‘truth’ and more focused on establishing a reasonable range, recognizing that the true value may vary across populations,” Dr. Dean said. “Bayesian analyses allow us to include this variability in a clear way, and then propagate this uncertainty through the model.”

A textbook application of Bayes’s theorem is serology testing for Covid-19, which looks for the presence of antibodies to the virus. All tests are imperfect, and the accuracy of an antibody test turns on many factors including, critically, the rarity or prevalence of the disease.

#styln-briefing-block { font-family: nyt-franklin,helvetica,arial,sans-serif; background-color: #ffffff; color: #121212; box-sizing: border-box; margin: 30px auto; max-width: 510px; width: calc(100% – 40px); border-top: 5px solid #121212; border-bottom: 2px solid #121212; padding: 5px 0 10px 0; } @media only screen and (min-width: 600px) { #styln-briefing-block { margin: 40px auto; } } #styln-briefing-block a { color: #121212; } #styln-briefing-block ul { margin-left: 15px; } #styln-briefing-block a.briefing-block-link { color: #121212; border-bottom: 1px solid #cccccc; font-size: 0.9375rem; line-height: 1.375rem; } #styln-briefing-block a.briefing-block-link:hover { border-bottom: none; } #styln-briefing-block .briefing-block-bullet::before { content: ‘•’; margin-right: 7px; color: #333; font-size: 12px; margin-left: -13px; top: -2px; position: relative; } #styln-briefing-block .briefing-block-bullet:not(:last-child) { margin-bottom: 0.75em; } #styln-briefing-block .briefing-block-header-section { margin-bottom: 16px; } #styln-briefing-block .briefing-block-header { font-weight: 700; font-size: 1.125rem; line-height: 1.375rem; display: inline-block; margin-bottom: 5px; } @media only screen and (min-width: 600px) { #styln-briefing-block .briefing-block-header { font-size: 1.25rem; line-height: 1.5625rem; } } #styln-briefing-block .briefing-block-header a { text-decoration: none; color: #333; } #styln-briefing-block .briefing-block-header a::after { content: ‘›’; position: relative; font-weight: 500; margin-left: 5px; } #styln-briefing-block .briefing-block-footer { font-size: 14px; margin-top: 1.25em; /* padding-top: 1.25em; border-top: 1px solid #e2e2e2; */ } #styln-briefing-block .briefing-block-briefinglinks a { font-weight: bold; margin-right: 6px; } #styln-briefing-block .briefing-block-footer a { border-bottom: 1px solid #ccc; } #styln-briefing-block .briefing-block-footer a:hover { border-bottom: 1px solid transparent; } #styln-briefing-block .briefing-block-header { border-bottom: none; } #styln-briefing-block .briefing-block-lb-items { display: grid; grid-template-columns: auto 1fr; grid-column-gap: 20px; grid-row-gap: 15px; line-height: 1.2; } #styln-briefing-block .briefing-block-update-time a { color: #999; font-size: 12px; } #styln-briefing-block a { color: #D0021B; } #styln-briefing-block .briefing-block-footer-meta { display: none; justify-content: space-between; align-items: center; } #styln-briefing-block .briefing-block-ts { color: #D0021B; font-size: 12px; display: block; } @media only screen and (min-width: 600px) { #styln-briefing-block a.briefing-block-link { font-size: 1.0625rem; line-height: 1.5rem; } #styln-briefing-block .briefing-block-bullet::before { content: ‘•’; margin-right: 10px; color: #333; font-size: 12px; margin-left: -15px; top: -2px; position: relative; } #styln-briefing-block .briefing-block-update-time a { font-size: 13px; } } @media only screen and (min-width: 1024px) { #styln-briefing-block { width: 100%; } }

Updated 2020-08-05T05:01:36.114Z

The first SARS-CoV-2 antibody test approved by the F.D.A., in April, seemed to be wrong as often as it was right. With Bayes’s theorem, you can calculate what you really want to know: the probability that the test result is correct. As one commenter on Twitter put it: “Understanding Bayes’s theorem is a matter of life and death right now.”

Joseph Blitzstein, a statistician at Harvard, delves into the utility of Bayesian analysis in his popular course “Statistics 110: Probability.” For a primer, in lecture one, he says: “Math is the logic of certainty, and statistics is the logic of uncertainty. Everyone has uncertainty. If you have 100 percent certainty about everything, there is something wrong with you.”

By the end of lecture four, he arrives at Bayes’s theorem — his favorite theorem because it is mathematically simple yet conceptually powerful.

“Literally, the proof is just one line of algebra,” Dr. Blitzstein said. The theorem essentially reduces to a fraction; it expresses the probability P of some event A happening given the occurrence of another event B.


“Naïvely, you would think, How much could you get from that?” Dr. Blitzstein said. “It turns out to have incredibly deep consequences and to be applicable to just about every field of inquiry” — from finance and genetics to political science and historical studies. The Bayesian approach is applied in analyzing racial disparities in policing (in the assessment of officer decisions to search drivers during a traffic stop) and search-and-rescue operations (the search area narrows as new data is added). Cognitive scientists ask, ‘Is the brain Bayesian?’ Philosophers of science posit that science as a whole is a Bayesian process — as is common sense.

Take diagnostic testing. In this scenario, the setup of Bayes’s theorem might use events labeled “T” for a positive test result — and “C” for the presence of Covid-19 antibodies:

Now suppose the prevalence of cases is 10 percent (that was so in New York City in the spring), and you have a positive result from a test with accuracy of 87.5 percent sensitivity and 97.5 percent specificity. Running numbers through the Bayesian gears, the probability that the result is correct, and that you do indeed have antibodies is 79.5%. Decent odds, all things considered. If you want more certainty, get a second opinion. And continue to be cautious.

An international collaboration of researchers, doctors and developers created another Bayesian strategy, pairing the test result with a questionnaire to produce a better estimate of whether the result might be a false negative or a false positive. The tool, which has won two hackathons, collects contextual information: Did you go to work during lockdown? What did you do to avoid catching Covid-19? Has anyone in your household had Covid-19?

“It’s a little akin to having two ‘medical experts,’” said Claire Donnat, who recently finished her Ph.D. in statistics at Stanford and was part of the team. One expert has access to the patient’s symptoms and background, the other to the test; the two diagnoses are combined to produce a more precise score, and more reliable immunity estimates. The priors are updated with an aggregation of information.

“As new information comes in, we update our priors all the time,” said Susan Holmes, a Stanford statistician, via unstable internet from rural Portugal, where she unexpectedly pandemicked for 105 days, while visiting her mother.

.css-1wxds7f{margin-bottom:10px;font-family:nyt-franklin,helvetica,arial,sans-serif;font-weight:700;font-size:0.875rem;line-height:1.25rem;color:#333 !important;}.css-2al2sh{font-family:nyt-franklin,helvetica,arial,sans-serif;font-size:0.9375rem;line-height:1.25rem;color:#333;margin-bottom:0.78125rem;margin-top:20px;margin-bottom:5px;font-weight:700;}@media (min-width:740px){.css-2al2sh{font-size:1.0625rem;line-height:1.5rem;margin-bottom:0.9375rem;}}@media (min-width:740px){.css-2al2sh{margin-bottom:10px;}}.css-1yyoic1{font-family:nyt-franklin,helvetica,arial,sans-serif;font-size:0.9375rem;line-height:1.25rem;color:#333;margin-bottom:0.78125rem;}@media (min-width:740px){.css-1yyoic1{font-size:1.0625rem;line-height:1.5rem;margin-bottom:0.9375rem;}}.css-zkk2wn{margin-bottom:20px;font-family:nyt-franklin,helvetica,arial,sans-serif;font-size:0.875rem;line-height:1.5625rem;color:#333;}.css-1dvfdxo{margin:10px auto 0px;font-family:nyt-franklin,helvetica,arial,sans-serif;font-weight:700;font-size:1.125rem;line-height:1.5625rem;color:#121212;}@media (min-width:740px){.css-1dvfdxo{font-size:1.25rem;line-height:1.875rem;}}.css-16ed7iq{width:100%;display:-webkit-box;display:-webkit-flex;display:-ms-flexbox;display:flex;-webkit-align-items:center;-webkit-box-align:center;-ms-flex-align:center;align-items:center;-webkit-box-pack:center;-webkit-justify-content:center;-ms-flex-pack:center;justify-content:center;padding:10px 0;background-color:white;}.css-pmm6ed{display:-webkit-box;display:-webkit-flex;display:-ms-flexbox;display:flex;-webkit-align-items:center;-webkit-box-align:center;-ms-flex-align:center;align-items:center;}.css-pmm6ed > :not(:first-child){margin-left:5px;}.css-5gimkt{font-family:nyt-franklin,helvetica,arial,sans-serif;font-size:0.8125rem;font-weight:700;-webkit-letter-spacing:0.03em;-moz-letter-spacing:0.03em;-ms-letter-spacing:0.03em;letter-spacing:0.03em;text-transform:uppercase;color:#333;}.css-5gimkt:after{content:’Collapse’;}.css-rdoyk0{-webkit-transition:all 0.5s ease;transition:all 0.5s ease;-webkit-transform:rotate(180deg);-ms-transform:rotate(180deg);transform:rotate(180deg);}.css-eb027h{max-height:5000px;-webkit-transition:max-height 0.5s ease;transition:max-height 0.5s ease;}.css-6mllg9{-webkit-transition:all 0.5s ease;transition:all 0.5s ease;position:relative;opacity:0;}.css-6mllg9:before{content:”;background-image:linear-gradient(180deg,transparent,#ffffff);background-image:-webkit-linear-gradient(270deg,rgba(255,255,255,0),#ffffff);height:80px;width:100%;position:absolute;bottom:0px;pointer-events:none;}#masthead-bar-one{display:none;}#masthead-bar-one{display:none;}.css-19mumt8{background-color:white;margin:30px 0;padding:0 20px;max-width:510px;}@media (min-width:740px){.css-19mumt8{margin:40px auto;}}.css-19mumt8:focus{outline:1px solid #e2e2e2;}.css-19mumt8 a{color:#326891;-webkit-text-decoration:none;text-decoration:none;border-bottom:2px solid #ccd9e3;}.css-19mumt8 a:visited{color:#333;-webkit-text-decoration:none;text-decoration:none;border-bottom:2px solid #ddd;}.css-19mumt8 a:hover{border-bottom:none;}.css-19mumt8[data-truncated] .css-rdoyk0{-webkit-transform:rotate(0deg);-ms-transform:rotate(0deg);transform:rotate(0deg);}.css-19mumt8[data-truncated] .css-eb027h{max-height:300px;overflow:hidden;-webkit-transition:none;transition:none;}.css-19mumt8[data-truncated] .css-5gimkt:after{content:’See more’;}.css-19mumt8[data-truncated] .css-6mllg9{opacity:1;}.css-a8d9oz{border-top:5px solid #121212;border-bottom:2px solid #121212;margin:0 auto;padding:5px 0 0;overflow:hidden;}

That was the base from which Dr. Holmes refined a preprint paper, co-authored with Dr. Donnat, that provides another example of Bayesian analysis, broadly speaking. Observing early research in March about how the pandemic might evolve, they noticed that classic epidemiological models tend to use fixed parameters, or constants, for the reproduction number — for instance, with an R0 of 2.0.

But in reality, the reproduction number depends on random, uncertain factors: viral loads and susceptibility, behavior and social networks, culture and socioeconomic class, weather, air conditioning and unknowns.

With a Bayesian perspective, the uncertainty is encoded into randomness. The researchers began by supposing that the reproductive number had various distributions (the priors). Then they modeled the uncertainty using a random variable that fluctuates, taking on a range of values as small as 0.6 and as large as 2.2 or 3.5. In something of a nesting process, the random variable itself has parameters that fluctuate randomly; and those parameters, too, have random parameters (hyper-parameters), etcetera. The effects accumulate into a “Bayesian hierarchy” — “turtles all the way down,” Dr. Holmes said.

The effects of all these up-and-down random fluctuations multiply, like compound interest. As a result, the study found that using random variables for reproductive numbers more realistically predicts the risky tail events, the rarer but more significant superspreader events.

Humans on their own, however, without a Bayesian model for a compass, are notoriously bad at fathoming individual risk.

“People, including very young children, can and do use Bayesian inference unconsciously,” said Alison Gopnik, a psychologist at the University of California, Berkeley. “But they need direct evidence about the frequency of events to do so.”

Much of the information that guides our behavior in the context of Covid-19 is probabilistic. For example, by some estimates, if you get infected with the coronavirus, there is a 1 percent chance you will die; but in reality an individual’s odds can vary by a thousandfold or more, depending on age and other factors. “For something like an illness, most of the evidence is usually indirect, and people are very bad at dealing with explicit probabilistic information,” Dr. Gopnik said.

Even with evidence, revising beliefs isn’t easy. The scientific community struggled to update its priors about the asymptomatic transmission of Covid-19, even when evidence emerged that it is a factor and that masks are a helpful preventive measure. This arguably contributed to the world’s sluggish response to the virus.

“The problems come when we don’t update,” said David Spiegelhalter, a statistician and chair of the Winton Centre for Risk and Evidence Communication at the University of Cambridge. “You can interpret confirmation bias, and so many of the ways in which we react badly, by being too slow to revise our beliefs.”

There are techniques that compensate for Bayesian shortcomings. Dr. Spiegelhalter is fond of an approach called Cromwell’s law. “It’s heaven,” he said. In 1650, Oliver Cromwell, Lord Protector of the Commonwealth of England, wrote in a letter to the Church of Scotland: “I beseech you, in the bowels of Christ, think it possible you may be mistaken.”

In the Bayesian world, Cromwell’s law means you should always “keep a bit back — with a little bit of probability, a little tiny bit — for the fact that you may be wrong,” Dr. Spiegelhalter said. “Then if new evidence comes along that totally contradicts your main prior belief, you can quickly ditch what you thought before and lurch over to that new way of thinking.”

“In other words, keep an open mind,” said Dr. Spiegelhalter. “That’s a very powerful idea. And it doesn’t necessarily have to be done technically or formally; it can just be in the back of your mind as an idea. Call it ‘modeling humility.’ You may be wrong.”

[Like the Science Times page on Facebook. | Sign up for the Science Times newsletter.]


Related posts

Nevada voters cope with ballot issues amid count


CNN anchor Brooke Baldwin details ‘very dark’ battle with COVID-19


UK ‘allowed thousands of infected people into country’ in early days of pandemic


Leave a Comment