Testing positive for the Coronavirus doesn’t mean one actually has a SARS-CoV-2 infection. Like statistical models, one might say that all medical tests are flawed*, but some (hopefully all that we make decisions based off of) are useful. While creating these tests, scientists balance cost (in **many** ways), speed, generalizability, and of course, accuracy. Most of the time, in the case of a binary test (ex. you are pregnant or you are not, etc.) ‘Accuracy’ can be broken down into two components - sensitivity and specificity.

Test **sensitivity** is the probability of a test correctly identifying those with the disease as having the disease.

Test **specificity** is the probability of a test correctly identify those without the disease as not having the disease.

We can imagine creating a test with 100% sensitivity by always providing a positive reading at the complete expense of specificity and vice versa. Depending on what is more important (not always clear), scientists can improve one at the cost of its counterpart.

Understanding these concepts and using them in practice can lead to some interesting and sometimes counterintuitive test result interpretation when we start looking at health from a *population* level. What is the probability of having the disease given I tested positive? What is the probability that my neighbor Dave actually does have the disease given he told me that he tested negative for COVID-19 last week as he continues to try to talk to me while coughing at extremely close distances?!?! You get the idea…

To figure out (read: estimate) the probability of having the virus given testing positive (or even negative) we can use Bayes Theorem - a simple derivation of the definition of conditional probability:

\[P(A|B) = \frac{P(B|A)P(A)}{P(B)}\]

equivalently:

\[P(A|B) = \frac{P(B|A)P(A)}{P(B|A)P(A) + P(B|A^{-})P(A^{-})}\]

Rephrasing in the context of COVID-19:

The probability of having COVID-19 given you have tested positive for it (\(P(D^{+}|T^{+})\)) is equal to the probability of testing positive given you have the disease (\(P(T^{+}|D^{+})\)) scaled by the proportion of the population with the disease (\(P(D^{+})\)), divided by a member of your population’s probability of testing positive for the disease (\(P(T^{+})\)). If we knew (for a fact) each piece of this puzzle, we could correctly state the probability of an individual having COVID19 given they tested positive.

\[P(D^{+}|T^{+}) = \frac{P(T^{+}|D^{+})P(D^{+})}{P(T^{+})} = \frac{P(T^{+}|D^{+})P(D^{+})}{P(T^{+}|D^{+})P(D^{+}) + P(T^{+}|D^{-})P(D^{-})}\]

## Sensitivity (AKA True Positive Rate or Recall): \(P(T^{+}|D^{+})\)

Measures of sensitivity will differ greatly based on the tests conducted. It is reported as being as high as 97% for some chest CT scans and as low as 30-60% for certain RT-PCR tests.

Getting specific, Cellex Inc.’s Rapid diagnostic test (RDT) which can determine if a sample of blood serum contains antibodies against the virus in less than a half hour has a reported sensitivity of 93.8%. Meaning that 93.8% of those with COVID-19 will be receive a positive result from the test, while the remaining 6.8% will receive a false negative.

## Disease Prevalence: \(P(D^{+})\)

To get disease prevalence we can simply divide total number of cases by population size. Except we need to: 1.) Define the population. and 2.) Estimate the prevalence since we cannot know the true number.

As of 2020-04-18 there are 17,550 confirmed cases in my home state of Connecticut - with an estimated population of 3.5 million, ~0.5% of the population has at one point had a positive test result for the disease. Of course we know that actual prevalence is likely higher than that with. The very controversial preprint from Stanford University was released the other day, suggesting a prevalence in Santa Clara County under certain assumptions to be as high as ~4.19%.

While these estimates are fine to get some sort of ballpark estimate. The *population* when considering *population prevelance* consists of more than just geographic bounds. For example, in Connecticut, African American’s tested positive for COVID-19 at a rate 20% higher than Whites. The population of certain nursing homes also have a much higher prevalence than populations of neighboring communities. We could also imagine that simply displaying a fever or a cough could place you in a high risk ‘group’ where people with those symptoms have a different prevalence than the general population.

For this example, we will go with a population prevalence of 4%.

## False Positive Rate (AKA Fallout AKA 1 - specificity): \(P(T^{+}|D^{-})\)

The Cellex RDT’s sensitivity is reported to be 95.6%. In other words, its false positive rate is about 4.4%.

## Putting this all together now.

So we know have some ballpark estimates for the sensitivity, specificity, and prevalence, what does that get us?

\[ \frac{P(T^{+}|D^{+})P(D^{+})}{P(T^{+}|D^{+})P(D^{+}) + P(T^{+}|D^{-})P(D^{-})}\]

\(P(T^{+}|D^{+}) = 0.938\)

\(P(D^{+}) = 0.04\)

\(P(T^{+}|D^{-}) = 0.044\)

`p = (0.938*0.04)/(0.938*0.04 + 0.044*0.96)`

With these estimates for sensitivity, specificity, and prevalence, a positive test result indicates that you are not actually very likely (probability of 47%) to have the disease.

This estimate is, however, extremely sensitive to the three inputs. For example, imagine you have all of the symptoms of COVID-19 - you are now no longer in the general population with \(p_{g}(D^+)=4\%%\). Instead, the prevalence *in this high risk group* may be much higher, perhaps 50% (I don’t really know what it is).

`p_hr = (0.938*0.5)/(0.938*0.5 + 0.044*0.5)`

If that was the case, with the same test, a positive test would indicate a COVID-19 infection 95.5% of the time.

This is precisely the reason why we commonly don’t test for diseases (especially rare diseases) randomly. A test with very high sensitivity and specificity is almost useless if the baseline prevalence is near zero as most positives will actually be false positives.

## Extra work

While this is one of the most basic applications of Bayes theorem that is always covered in undergraduate statistics courses as soon as the theorem is brought up, there are some extensions we can use to get a more credible estimate and (perhaps more importantly) quantify our uncertainty about the estimate. All three parameters we utilized have some uncertainty surrounding them, we could capture these uncertainties in a MCMC to produce estimates and credible intervals around these estimates based on the data we examine.

* I am not a doctor and this is not medical advice :)