Got a Positive Covid Test Result? It might not mean what you think it does.
Ed note 8/12/20 11:30am: Based on feedback from my readers here is an important clarification: This blog post refers to a situation of an individual Covid-19 test case and what the inferences are regarding a positive result when the prevalence of infection is relatively small.
Infection rates at the population level (e.g. state or nation level reporting) are:
- Based on an entirely different set of assumptions and are unaffected by the phenomenon reported herein.
- Press reports of Covid infections are likely under-reported by a factor of at least 6X - meaning if the press reports 5 million infections in the US the real number is probably greater than 30 million.
- In addition, I'd like to re-emphasize that if you get a positive test result you need to behave as if there is 100% chance that you have Covid-19 - and isolate for a 14 day period.
You are asymptomatic and just got tested for Covid-19. It came back positive. What are your chances of actually having Covid?
38.5%
i.e. there is a 61.5% chance you don't actually have it.
This is assuming the Covid test is 98% accurate (regarding its specificity).
How could a positive result from a 98% accurate test give a 38.5% probability of being infected? This goes against the grain of common sense.
There is actually nothing wrong with the test or the methodology. It also has nothing to do with being asymptomatic and actually having the disease (this is a separate issue and has been widely reported).
This is just simple math and this is what we will explore in today's blog post.
Important note: Despite the lower chance of having the infection that is described here, all positive test results should always be followed by at least 14 days of isolation and quarantine per the CDC's recommendations.
The Bayesian Approach
Thomas Bayes (/beɪz/; c. 1701 – 7 April 1761) was an English statistician, philosopher and Presbyterian minister who is known for formulating a specific case of the theorem that bears his name: Bayes' theorem (source: Wikipedia).
Here is a quick review of the three steps:
- The prior: What kind of evidence do we currently have about Covid-19 (such as how common it is in the population)?
- What new evidence is being presented (= the Covid test results) and how strong is that evidence (e.g. how accurate is the test)?
- After factoring in both the strength of the prior and the confidence we have in the new evidence, what is the new probability (i.e. how probable is it that you have Covid given the prior and the new evidence)?
Defining the prior
- Covid-19 prevalence in Maryland is estimated to be 1.37% (for August 11, 2020).
New Evidence
- False negative rates for PCR tests range between 2% and 29%. This implies a sensitivity of between 71% and 98%. We'll give the test the benefit of the doubt and use 90% as a sensitivity.
- False positive rates for PCR tests range between 0 and 4%. This implies a specificity of between 96% and 100%. We will take the average of 98%.
Calculating the Posterior Probability
- The test will catch 90% of the 1.37% infected (=1.23% of the population) - these are the true positives.
- But the test will also report 2% of the remaining 98.63% (=1.97% of the population) that actually doesn't have Covid - these are the false positives.
Now what?
The Second Test
Appendix
The First Test
- P(cv | +) = the probability of having Covid given a positive test result (this is the answer we are looking for)
- P(+ | cv) = the probability of getting a positive test result given you have Covid (this is the sensitivity of the test) = 90%
- P(cv) = the prior probability of having Covid (the prevalence of the infection) = 1.37%
- P(+) = the probability of a positive test result (see below)
- P(+ | no cv) = probability of getting a positive test result where there is no Covid (this is the false positive rate) = 1 - Specificity = 1 - 98% = 0.02
- P(no cv) = probability of not having Covid = 1 - 1.37% = 0.9863
The Second Test (assume positive test result)
- P(cv | +) = the probability of having Covid given a second positive test result (this is the answer we are looking for)
- P(+ | cv) = the probability of getting a positive test result given you have Covid (this is the sensitivity of the test) = 90%
- P(cv) = the prior probability of having Covid (the posterior probability from the first test) = 37.5%
- P(+) = the probability of a second positive test result (see below)
- P(+ | no cv) = probability of getting a positive test result where there is no Covid (this is the false positive rate) = 1 - Specificity = 1 - 98% = 0.02
- P(no cv) = probability of not having Covid = 1 - 37.5% = 0.615
The Second Test (assume negative test result)
- P(cv | -) = the probability of having Covid given a second negative test result (this is the answer we are looking for)
- P(- | cv) = the probability of getting a negative test result given you have Covid (this is the false negative rate of the test) = 10%
- P(cv) = the prior probability of having Covid (the posterior probability from the first test) = 37.5%
- P(-) = the probability of a negative test result (see below)
- P(- | no cv) = probability of getting a negative test result where there is no Covid (this is the true negative rate) = Specificity = 98%
- P(no cv) = probability of not having Covid = 1 - 37.5% = 0.615
Negative Results for the First Test
- P(cv | -) = the probability of having Covid given a negative test result (this is the answer we are looking for)
- P(- | cv) = the probability of getting a negative test result given you have Covid (this is the false negative rate of the test) = 10%
- P(cv) = the prior probability of having Covid = 1.37%
- P(-) = the probability of a negative test result (see below)
- P(- | no cv) = probability of getting a negative test result where there is no Covid (this is the true negative rate) = Specificity = 98%
- P(no cv) = probability of not having Covid = 1 - 1.37% = 0.9863
Comments
Post a Comment