Are fast test results preferable to high test sensitivity in contact-tracing strategies? ======================================================================================== * Jonas L. Juul * Kaare Græsbøll ## Abstract Across the world, countries are fighting to reduce the spread of COVID-19. The backbone of the response is a test-trace-isolate strategy, where suspected infected get tested and isolated and possible secondary cases get traced, tested and isolated. Because more accurate tests often take longer to analyze, and the benefits of contact tracing are strengthened by rapid diagnosis, there exists a trade-off in test sensitivity and test waiting time in test-trace-isolate strategies. Here we ask: How many false negatives can be tolerated in a rapid test so that it reduces transmission better than a slower, more accurate test? How does this change with contact tracing efficiency and test waiting time? We find that a rapid, less sensitive test performs best for many test-parameter choices and that this is true even for modest contact tracing efficiency. For COVID-19-like viral parameters, a test with 40% false negatives and immediate result might reduce transmission as well as a test with no false negatives and a 3-day waiting time. Our analysis suggests employing rapid tests to reduce test waiting times as a viable strategy to reduce transmission when testing infrastructure is under stress. ## I. INTRODUCTION Rapid diagnosis and isolation of COVID-19 cases is critical in reducing further transmission of the virus. Since a large fraction of infections take place before the infected develops symptoms,1–3 isolation following the onset of symptoms is not sufficient to control the pandemic.4,5 To ensure rapid diagnosis and reduce transmission before symptom onset, the World Health Organization therefore recommends tracing, testing and isolating close contacts of COVID-19 positives.6 The success of the test-trace-isolate strategy depends on many factors. In addition to social factors such as the population being informed and complying with guidelines, successful rapid diagnosis and isolation depends heavily on the following 3 factors. *i*. Test waiting time: How long it takes from a person wants to get tested to the result of the test arrives and contact tracing begins. *ii*. Test sensitivity: How often a test returns a false negative. *iii*. Tracing efficiency: The fraction of secondary cases that are successfully found through contact tracing confirmed positives. These 3 strategy parameters – test waiting time, test sensitivity, and tracing efficiency – are not necessarily fixed. A surge in cases might cause an increase in test demands. Higher test demand could in turn cause longer turnaround times as seen in some countries during the fall and winter of 2020.7,8 Longer turnaround times lead to slower tracing of secondary infections and thus a vicious cycle begins. The 3 strategy parameters are also not independent. Faster test results often come at the cost of lower sensitivity. PCR tests are highly accurate, but take significant time to analyze. Faster results can be achieved by pooling tests9–11 or using rapid antigen tests,12 lateral flow devices13 or saliva tests.14 Recent analyses concluded that test sensitivity is secondary to frequency and turnaround times in population-screening strategies15,16 as implemented by Slovakia in November 2020.17 In context of the more-widely applied test-trace-isolate strategies, rapid tests increase the effect of contact tracing. This raises the natural question: When is it better to get a fast result than an accurate result? How does this change with the efficiency of contact tracing efforts, and with turnaround times? Here we confront these questions. ## II. MODEL DESCRIPTION To quantify the impact of test-trace-isolate strategies on growing epidemics, we simulate the branching structure of the chains of infections, also referred to as the “epidemic tree”.18,19 Our model lets us keep track of who infected whom, which is essential when simulating contact tracing and quarantining infectious people. In the model (Figure 1A), infected individuals give rise to new cases, unless quarantined following testing and tracing efforts. ![FIG. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/19/2021.02.17.21251921/F1.medium.gif) [FIG. 1.](http://medrxiv.org/content/early/2021/02/19/2021.02.17.21251921/F1) FIG. 1. **A** Illustration of our model. Infectious people transmit the disease unless they are placed in quarantine following testing and tracing efforts. The success of these efforts are contingent on the three key parameters of the model: the test waiting time *δ*, test sensitivity 1 *− p*false, and tracing efficiency *p*trace. **B** Infectiousness over time, *P*time(*t*), in our model. Background colors specify whether the infected experiences symptoms and is infectious at the given time (see legend in Figure 1A) The value on the vertical axis for day *x* after infection is the probability that a given secondary case gets infected on day *x*. For asymptomatic cases, *P*time(*t*) is as depicted – the only difference being that the final, symptomatic, phase is replaced with an asymptomatic phase. We initiate a simulation with some number of newly infected people, *N*seed *∈* N. The simulation progresses in discrete timesteps, corresponding to days, and we make the simplifying assumption that every infected person, goes through the same phases before recovering: 3 days of being presymptomatic and noninfectious followed by 8 infectious days. As for COVID-19,20 some fraction of cases, *p*asymp, remain asymptomatic for the entire infectious period; all other infected cases experience symptoms starting on day 7 after infection. In the absence of testing, tracing and isolation, an infectious person would give rise to *k* secondary cases. For each infected, we assume that *k* is drawn from the probability distribution *P* (*k*). When each of these *k* secondary cases is infected, we determine the time of infection by drawing an integer from the probability distribution *P*time(*t*). *P*time(*t*) takes positive values on the days where the infected is infectious and, mimicking COVID-19,1,15 peaks around symptom onset (Figure 1B). These infections take place unless the infected is in quarantine at the time the infection would occur. Infectious people quarantine only when waiting for a test to be taken, receiving a test result, or after testing positive. In our model, an infectious person orders a test if either of two things happens: 1) The person is traced; 2) The person develops symptoms. In either case, the person orders a test immediately and then waits *δ* = *δ*test +*δ*result days for the result. The test waiting time is divided into *δ*test days waiting for the test to be taken followed by *δ*result days to receive the result. The test correctly identifies the case with probability equal to its sensitivity, 1 *− p*false, where *p*false is the false negative rate. If the test comes back positive, each of the person’s secondary cases is traced with independent probability *p*trace. We obtain indistinguishable results when simulating the same *δ* with varying values of *δ*test and *δ*result. Thus, the key parameters of the model are the test waiting time *δ*, test sensitivity 1 *− p*false, and tracing efficiency *p*trace. Following *t*max timesteps, we count the number of nodes that completed their whole infectious period, *n*parents. We also count the number of people these nodes infected, *n*children. The output of the simulation is the effective reproduction number of simulated disease: *R*eff = *n*children*/n*parents. ## III. RESULTS We use the model to examine the trade-off between test waiting time and test accuracy depending on the tracing efficiency. In our simulations, we therefore vary the parameters *δ*, 1 *− p*false, and *p*trace and fix all other parameters (for *δ ≥* 1, we set *δ*result = 1). For *P* (*k*) we choose a Poisson distribution with mean *R* = 2 (slightly lower than estimated in early stages of the pandemic21,22) and for *P*time(*t*) we choose the right-skewed distribution depicted in Fig. 1B. Finally, we choose *t*max = 50. To develop some intuition, let us first introduce the results we obtain when fixing the tracing efficiency, *p*trace = 0.80. This constraint leaves 2 free parameters: The test sensitivity, 1 *− p*false, and the test delay, *δ*. We now compare the effective reproduction number ![Graphic][1] obtained by using a slow, but accurate test (parameters: ![Graphic][2] and some *δ*slow *≥* 1 day) to the reproduction number, ![Graphic][3], obtained with a less accurate, but rapid test (parameters: some ![Graphic][4] and *δ*rapid = 0 days). To evaluate whether speed or accuracy is to be preferred, we compute the difference in obtained effective reproduction numbers of the virus under the different choices of tests, ![Graphic][5]. Let us choose some test delay, e.g. *δ*slow = 2 days (perhaps corresponding to a PCR test with a 1-day waiting time to get tested and a subsequent 1-day waiting time to get the result). In this case, how will Δ*R*eff depend on the risk of getting a false negative test result? For very high sensitivity (low *p*false), this faster test will be almost as accurate as the slower test it is being compared to. For this reason, the fast test will be preferable to the slower one (Fig. 2A top colorbar). If we now imagine slowly decreasing the test sensitivity, Δ*R*eff will gradually increase until it reaches a breaking point where the slow and rapid tests reduce the effective reproduction number equally well: Δ*R*eff = 0. Decreasing the sensitivity even further makes Δ*R*eff positive, meaning that for this high probabilities of false negatives, the accurate test is to be preferred. ![FIG. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/19/2021.02.17.21251921/F2.medium.gif) [FIG. 2.](http://medrxiv.org/content/early/2021/02/19/2021.02.17.21251921/F2) FIG. 2. **A** Comparison of effective reproduction number, ![Graphic][6], when using a slow but accurate test with sensitivity 1 (no false negatives), and a test waiting time of *δ*slow days and using a rapid but less accurate test with a given sensitivity and test waiting time of *δ*rapid = 0 days. In this inset, we assume that 80% of secondary cases are successfully traced following a positive test. The colorbars show results obtained for different choices of *δ*slow (top colorbar: *δ*slow = 2 days, middle: *δ*slow = 3 days, bottom: *δ*slow = 4 days), and the rapid-test sensitivity (horizontal axis in each colorbar). In each colorbar a breaking point separates sensitivity values favoring the slower, accurate test and values favoring the rapid test. **B** Comparison of effective reproduction numbers when the result of the slower test arrives after *δ*slow = 2 days, as a function of the sensitivity and tracing efficiency. **C** Same as in B but with *δ*slow = 3 days. **D** Same as in B, C but with the slower result arriving after *δ*slow = 4 days. **E** Heatmap of the effective reproduction number obtained using the rapid, less accurate test, ![Graphic][7]. The breaking-point lines of B, C, D are plotted in black. If we now make the same plot but with a higher waiting time to receive the slow test, the breaking point will move to a lower test sensitivity. The second colorbar in Fig. 2A shows the result for *δ*slow = 3 days. The breaking point moved from sensitivity *≈* 0.91 to *≈* 0.63 with this single day increase in waiting time – corresponding to a fourfold increase in *p*false! Increasing the waiting time once more, setting *δ*slow = 4, moves the breaking point further to the left. In this latter case, a 54% chance of a false negative is better than a 4 day wait for accuracy. Having established some intuition for the simulations, we proceed to varying the third parameter: the tracing efficiency. By varying the tracing efficiency, for each choice of *δ*slow we get 2-dimensional heatmaps instead of the one-dimensional colorbars presented in the previous paragraphs (Fig. 2B-D). In these heatmaps, the breaking points become white curves. Every point to the right of the breaking-point curve is a parameter combination where a faster test is preferable. Every point to the left of the breaking-point curve is a parameter combination that favors an accurate test. Notice how all the breaking-point curves start in the lower-right corner (where tests are completely accurate but no contact tracing is done) and how quickly they move to the left with increasing tracing efficiency. Figure 2E plots the obtained ![Graphic][8]. For each simulated choice of parameters, each computed *R*eff is averaged over 10 simulations. For clarity, the heatmaps in Fig. 2 have been smoothed with a Gaussian filter. ## IV. DISCUSSION Testing, tracing and isolating positive cases is central in many countries’ strategy to fight the current COVID-19 pandemic.23–25 We have demonstrated that there is a sizeable trade-off between test sensitivity and test waiting times in such strategies, and that it is often beneficial to prioritize test speed over test sensitivity. Moreover, we find that this benefit of rapid tests increases quickly with increases in test waiting times, and that even modest tracing efficiency unlocks the advantages of rapid tests. This indicates that additional waiting time for test results must be avoided and that it often makes sense to reduce test sensitivity in order to do so. It is to be expected that testing systems will occasionally get under stress during a pandemic, and having a way to avoid build-up of waiting times in this scenario is crucial. Designing such stress-relieve strategies presents an interesting direction for future research. Some of the assumptions we have made can be questioned. Three such assumptions are: That quarantine hinters any transmission; that symptomatic individuals do not quarantine after testing (false) negative; and that the probability of getting a false negative result does not depend on the infected’s infectiousness at the time the test was taken. We note that all of these assumptions will favor reducing transmission by slower, more accurate tests: That people do not break isolation benefits the tests with long waiting times; false negative results leading to completely normal behavior damages only the tests that allow for false negative results; lastly, making false negatives less likely at high viral load would make the rapid tests more reliable early in the course of disease, when many new secondary cases could be avoided following diagnosis. That we have chosen our assumptions as to disfavor rapid low-sensitivity tests means that our results can be interpreted as conservative estimates of the benefits of reducing test waiting time with less sensitive tests. Our choices of the probability distributions *P* (*k*) and *P*time(*t*) can also be questioned. A better choice for *P* (*k*) might be a heavy-tailed distribution that could account for superspreading behavior.26,27 Choosing a geometric distribution for *P* (*k*) with the same mean yields indistinguishable results (not shown). On the other hand, the choice of *P*time(*t*) influences results: Higher infectivity early in the course of disease makes the benefit of rapid tests greater (simulations not shown). Two other limitations arise not from parameter choices, but the model itself. First, our model does not take depletion of susceptibles into account. Such effects are unimportant to what we were interested in: the effect of testing and tracing on the effective reproduction number at a given stage of the pandemic. Secondly, our implementation of contact tracing allows only for descendants of the confirmed positive to be traced. Not including two-way contact tracing or the possibility of identifying individuals through non-parents could cause us to underestimate the effect of contact tracing.28–31 This, and not taking into account possible beneficial effects on contact tracing that might be caused by infected being present at the test site when receiving the diagnosis, again biases our results against the rapid tests that so rely on tracing efforts. Overall, our analysis suggests employing rapid tests to reduce test waiting times as a viable strategy to reduce transmission even at modest waiting times and contact tracing efficiency. ## Data Availability The code and data necessary to reproduce the content of this paper is available on [www.github.org/jonassjuul/TestTraceTradeoff](https://www.github.org/jonassjuul/TestTraceTradeoff) [https://www.github.org/jonassjuul/TestTraceTradeoff](https://www.github.org/jonassjuul/TestTraceTradeoff) ## AUTHOR CONTRIBUTIONS KG and JLJ designed the research. JLJ designed the model, carried out the simulations, and wrote first draft of the manuscript. KG contributed to designing the model and made significant contributions to writing the manuscript. ## COMPETING INTERESTS STATEMENT The authors declare that they have no competing interests. ## CODE AVAILABILITY The code and data necessary to reproduce the content of this paper is available on [www.github.org/jonassjuul/TestTraceTradeoff](http://www.github.org/jonassjuul/TestTraceTradeoff) ## ACKNOWLEDGMENTS We thank SSI for contributing to funding this research. JLJ thanks Steven H. Strogatz for helpful discussions and Gorm G. Jensen for comments on an early version of the manuscript. * Received February 17, 2021. * Revision received February 17, 2021. * Accepted February 19, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. X. He, E. H. Lau, P. Wu, X. Deng, J. Wang, X. Hao, Y. C. Lau, J. Y. Wong, Y. Guan, X. Tan, et al., Nature medicine 26, 672 (2020). 2. D. P. Oran and E. J. Topol, Annals of Internal Medicine (2020). 3. M. M. Arons, K. M. Hatfield, S. C. Reddy, A. Kimball, A. James, J. R. Jacobs, J. Taylor, K. Spicer, A. C. Bardossy, L. P. Oakley, et al., New England journal of medicine (2020). 4. S. M. Moghadas, M. C. Fitzpatrick, P. Sah, A. Pandey, A. Shoukat, B. H. Singer, and A. P. Galvani, Proceedings of the National Academy of Sciences 117, 17513 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzMwLzE3NTEzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDIvMTkvMjAyMS4wMi4xNy4yMTI1MTkyMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. N. C. Grassly, M. Pons-Salort, E. P. Parker, P. J. White, N. M. Ferguson, K. Ainslie, M. Baguelin, S. Bhatt, A. Boonyasiri, N. Brazeau, et al., The Lancet Infectious Diseases 20, 1381 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30630&link_type=DOI) 6. WHO-China Joint Mission, “Report of the WHO-China Joint Mission on Coronavirus Dis-ease 2019 (COVID-19).” (2020), [https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf](https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf) (Accessed Dec 10, 2020). 7. “Lange ventetider p∼a test i hele landet: Flere steder i hovedstaden er der ingen tider inden jul,” [https://www.dr.dk/nyheder/indland/lange-ventetider-paa-test-i-hele-landet-flere-steder-i-hovedstaden-er-der-ingen](https://www.dr.dk/nyheder/indland/lange-ventetider-paa-test-i-hele-landet-flere-steder-i-hovedstaden-er-der-ingen), accessed: 2021-01-14. 8. “Heading into holidays, us covid-19 testing strained again,” [https://apnews.com/article/](https://apnews.com/article/) us-covid-19-testing-strained-holidays-db20ebbcc1fa8a411be8f9ebc241af3b, accessed: 2021-01-14. 9. N. Barak, R. Ben-Ami, T. Sido, A. Perri, A. Shtoyer, M. Rivkin, T. Licht, A. Peretz, J. Magenheim, I. Fogel, et al., MedRxiv (2020). 10. N. Thierry-Mieg, BMC bioinformatics 7, 28 (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-7-28&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16423300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) 11. M. Gopalkrishnan and S. Krishna, Journal of the Indian Institute of Science, 1 (2020). 12. J. Dinnes, J. J. Deeks, A. Adriano, S. Berhane, C. Davenport, S. Dittrich, D. Emperador, Y. Takwoingi, J. Cunningham, S. Beese, et al., Cochrane Database of Systematic Reviews (2020). 13. L. Y. Lee, S. Rozmanowski, M. Pang, A. Charlett, C. Anderson, G. J. Hughes, M. Barnard, L. Peto, A. S. Richard Vipond, S. Hopkins, D. W. C. John Bel and, N. Gent, A. S. Walker, D. W. Eyre, and T. E. Peto, “An observational study of sars-cov-2 infectivity by viral load and demographic factors and the utility lateral flow devices to prevent transmission.” [http://modmedmicro.nsms.ox.ac.uk/wp-content/uploads/2021/01/infectivity\_manuscript\_20210119\_merged.pdf](http://modmedmicro.nsms.ox.ac.uk/wp-content/uploads/2021/01/infectivity_manuscript_20210119_merged.pdf) (2021). 14. A. L. Wyllie, J. Fournier, A. Casanovas-Massana, M. Campbell, M. Tokuyama, P. Vijayakumar, J. L. Warren, B. Geng, M. C. Muenker, A. J. Moore, C. B. Vogels, M. E. Petrone, I. M. Ott, P. Lu, A. Venkataraman, A. Lu-Culligan, J. Klein, R. Earnest, M. Simonov, R. Datta, R. Han-doko, N. Naushad, L. R. Sewanan, J. Valdez, E. B. White, S. Lapidus, C. C. Kalinich, X. Jiang, D. J. Kim, E. Kudo, M. Linehan, T. Mao, M. Moriyama, J. E. Oh, A. Park, J. Silva, E. Song, T. Takahashi, M. Taura, O.-E. Weizman, P. Wong, Y. Yang, S. Bermejo, C. D. Odio, S. B. Omer, C. S. Dela Cruz, S. Farhadian, R. A. Martinello, A. Iwasaki, N. D. Grubaugh, and A. I. Ko, New England Journal of Medicine 383, 1283 (2020), [https://doi.org/10.1056/NEJMc2016359](https://doi.org/10.1056/NEJMc2016359). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) 15. M. J. Mina, R. Parker, and D. B. Larremore, New England Journal of Medicine 383, e120 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) 16. D. B. Larremore, B. Wilder, E. Lester, S. Shehata, J. M. Burke, J. A. Hay, M. Tambe, M. J. Mina, and R. Parker, MedRxiv (2020). 17. E. Mahase, “Covid-19: Mass testing in slovakia may have helped cut infections,” (2020). 18. J. S. Juul and S. H. Strogatz, Physical Review Research 2, 033005 (2020). 19. D. T. Haydon, M. Chase-Topping, D. Shaw, L. Matthews, J. Friar, J. Wilesmith, and M. Wool-house, Proceedings of the Royal Society of London. Series B: Biological Sciences 270, 121 (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2002.2191&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12590749&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180748500002&link_type=ISI) 20. K. Mizumoto, K. Kagaya, A. Zarebski, and G. Chowell, Eurosurveillance 25, 2000180 (2020). 21. S. Zhao, Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, Y. Lou, D. Gao, L. Yang, D. He, et al., International journal of infectious diseases 92, 214 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.01.050&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32007643&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) 22. Y. Liu, A. A. Gayle, A. Wilder-Smith, and J. Rocklöv, Journal of travel medicine (2020). 23. J. Panovska-Griffiths, C. C. Kerr, R. M. Stuart, D. Mistry, D. J. Klein, R. M. Viner, and C. Bonell, The Lancet Child & Adolescent Health 4, 817 (2020). 24. E. Jeong, M. Hagose, H. Jung, M. Ki, and A. Flahault, International journal of environmental research and public health 17, 9571 (2020). 25. “Covid-19 strategy update,” [https://www.who.int/publications/m/item/covid-19-strategy-update](https://www.who.int/publications/m/item/covid-19-strategy-update), accessed:2021-02-16. 26. B. M. Althouse, E. A. Wenger, J. C. Miller, S. V. Scarpino, A. Allard, L. H’ebert-Dufresne, and H. Hu, arXiv preprint arXiv:2005.13689 (2020). 27. J. O. Lloyd-Smith, S. J. Schreiber, P. E. Kopp, and W. M. Getz, Nature 438, 355 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature04153&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16292310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233300200048&link_type=ISI) 28. K. T. Eames and M. J. Keeling, Proceedings of the Royal Society of London. Series B: Biological Sciences 270, 2565 (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2003.2554&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14728778&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187989900005&link_type=ISI) 29. I. Z. Kiss, D. M. Green, and R. R. Kao, Mathematical biosciences 203, 124 (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mbs.2006.03.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16620875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F19%2F2021.02.17.21251921.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241255700007&link_type=ISI) 30. S. Kojaku, L. H’ebert-Dufresne, and Y.-Y. Ahn, arXiv preprint arXiv:2005.02362 (2020). 31. W. J. Bradshaw, E. C. Alley, J. H. Huggins, A. L. Lloyd, and K. M. Esvelt, Nature communications 12, 1 (2021). [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: F2/embed/inline-graphic-6.gif [7]: F2/embed/inline-graphic-7.gif [8]: /embed/inline-graphic-8.gif