STeCC: Smart Testing with Contact Counting Enhances Covid-19 Mitigation by Bluetooth App Based Contact Tracing ============================================================================================================== * Hossein Gorji * Markus Arnoldini * David F. Jenny * Wolf-Dietrich Hardt * Patrick Jenny ## Abstract SARS-CoV2 spread is hard to control, as asymptomatic people contribute to transmission. Currently, Covid-19 mitigation imposes social distancing and isolates the diseased. This slows down virus spread, eases stress on health care systems and thereby reduces the death toll. However, this strategy takes a high economic toll, and virus transmission will surge again if measures are lifted. App-based contact tracing of symptomatic cases and isolating their contacts has been proposed as an alternative, but may not suffice for mitigation, as asymptomatic infections remain unidentified. Here, we evaluate complementary mitigation strategies relying on virus-RNA testing to detect and quarantine both, symptomatic and asymptomatic cases. Epidemic dynamics modeling shows that stopping the pandemic by mass testing alone is unrealistic, as we lack enough tests. However, realistic numbers of tests may suffice in a smart-testing strategy, e.g. when biasing tests towards people with exceptionally high numbers of contacts. These people are at particularly high risk to become infected (with or without symptoms) and transmit the virus. A mitigation strategy combining smart testing with contact counting (STeCC) and contact tracing in one app would reduce R by 2.4-fold (e.g. from R=2.4 to R=1) with realistic test numbers (≈166 per 100’000 people per day) when a realistic fraction of smartphone owners use the app (≈72%, i.e. ≈50% in total population). Thereby, STeCC expands the portfolio of mitigation strategies and may help easing social distancing without compromising public health. ## Main The Covid-19 pandemic has evaded initial containment measures. Current responses have therefore shifted towards mitigating the effects. However, proven vaccines and therapies are lacking and the current capacity for detecting the virus via its genomic RNA is limited (*1*). Thus, mitigation in many countries relies on extreme social distancing and diagnosing virus (SARS-CoV2) infections in infected people showing mild to severe symptoms. In combination, these two measures slow down pandemic spread and avoid overburdening healthcare systems, thus easing the demand for intensive care. However, this strategy has three major shortcomings. First, the limitation of testing to diseased people handicaps the efficiency assessment of current mitigation measures (i.e., reduce the basic virus reproduction number R<1). Second, it leaves many infected people with mild or no symptoms undetected (*2*), and therefore renders them more likely to infect others. Third, as social distancing measures limit non-essential business, it imposes severe economic consequences, which worsen over time. A broad social-distancing-based approach is therefore not sustainable, but mitigation measures need to stay in place until effective therapies or vaccines become available to avoid a second wave of virus spread. As these options are still months away, we need to consider alternative mitigation strategies. Smartphone apps for Bluetooth-based contact tracing, such as the European PEPP-PT project, the British NHSX, or the Swiss D3PT, can help to identify individuals that have recently had an infection-relevant contact (i.e. one that confers a risk of transmission) with known Covid-19 cases and might therefore have benn infected. This can help to mitigate the pandemic (*3*, *4*). However, quantitative modeling suggests that this may not suffice to replace social distancing completely, i.e. if some of the assumptions about the disease dynamics turn out to be unrealistic. Specifically, this approach may fail if <56% of the average population use the app, if the epidemic doubling time is <3.5 days, if >35% of the infected people are asymptomatic, or if they were more than 0.3-fold as infectious as symptomatic cases ((*3*); **Fig. S5**). Infectiousness of asymptomatic cases has been estimated to be between 10% and 100% of symptomatic cases *(4–6).* The asymptomatic cases are a particular problem, as app-based contact tracing detects such cases insufficiently. Therefore, we have looked for an alternative approach, based on testing for virus RNA, a marker for active infections with or without symptoms (**Fig. 1**, highlighted in dark or light blue). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F1) Figure 1. Graphical illustration of the modeling approach showing the dependencies within the system describing the dynamics of the susceptible, undetected and detected infected populations. It is crucial that the model distinguishes between individuals detected by symptoms (with or without testing; light blue), and those detected by virus testing (dark blue). First, we asked how well mass testing random samples of the population could complement or replace social distancing. In theory, this can be achieved by testing large fractions of the population at regular intervals, and isolating people who test positive. Using a deterministic modeling approach that explicitly considers infected cases with and without symptoms, we can estimate how many tests per 100’000 people per day would suffice to achieve the same effect as current social distancing measures which have reduced R from 2.4–3 to less than 1, e.g. in Germany or Switzerland ((*7*, *8*); for technical details, see **supplementary text**, **table S2**). We will first present our findings and then discuss how testing could be realistically implemented, e.g. by combining it with serological testing or app-based technologies. To estimate how mitigation strategies affect R, we use a mathematical model that employs a set of ordinary differential equations to describe the dynamics of the infection of a susceptible population (**Fig. 1**, **supplementary text** section 1, **Fig. S1**, **table S1**). To obtain realistic simulations, we parametrize our model using published data ((*5*, *9*–*14***)**; **table S2; supplementary text** section 2). Infection of susceptible people will lead to a latency phase (**Fig. 1**, exposed). The exposed will later become transmissive, and either remain asymptomatic and recover or become pre-symptomatic and later develop mild symptoms (**Fig. 1**, dark blue; see parameters in **table S2**). These two transmissive groups do not know that they are infected, and remain unidentified in current mitigation approaches. Some people with mild symptoms of disease will self-isolate (**Fig. 1**, light blue). Due to self-isolation, they will have a reduced probability to infect others (we assume a 90% reduction). Infected people with severe symptoms are hospitalized, immediately isolated under strict quarantine and do not infect others (**Fig. 1**, light blue). The same applies for anyone else who is tested virus-positive. In addition, we consider the effect of overloading the health system. After all intensive care units (ICUs) are occupied, additional cases requiring intensive care will suffer an elevated death rate (**Fig. 1; supplementary text**, section 2). As our model consists of a set of ordinary differential equations, which become linear in the early stages of the pandemic (when nearly the whole population still is susceptible), we can now analytically test how particular mitigation strategies affect Ro (**supplementary text**, sections 3, 4 and 5). In order to compare different mitigation strategies to a baseline, we first describe what the model predicts if we omit any mitigation. Here, we assume a value for R that is in the mid-range of published estimates, i.e. that one infected person infects 2.4 others in average (R=2.4; **Fig. 2A;** alternative R values are tested in **Fig. S2A,B**). Under these circumstances, 87% will either recover or die from the disease within ≈250 days, which compares well to the 81% predicted by (*5*) for UK and US populations in the absence of mitigation plans. We also predict more than 4% of the population to be killed by the virus, and ICUs to be at maximum capacity for ≈150 days of that year. For the unmitigated case these numbers are going to change only by 2-fold or less, if alternative plausible values are assumed for R (1.9 or 2.9; **Fig. S2A,B**) or for the availability of ICU beds (**Fig. S2C**). Thus, our basic predictions are robust, even if input parameters for the virus infection dynamics (**table S2**) might be subject to change when more precise parameters became available. Please note that **Fig. 2A** depicts the worst-case scenario, which can be improved by mitigation. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F2) Figure 2. Social distancing and extensive mass testing alone have qualitatively equal mitigating effects. We used the model described in the **supplementary text**. (A) Model outcome if no mitigation strategies are in place (R=2.4; see **table S3** for additional input parameters). (B) Model outcome if extreme social distancing (71% lower infection rate, leading to R=0.7; **table S2**) is in place between days 50 and 200 of the outbreak. After day 200, social distancing is discontinued. (C) Model outcome if mass testing with isolation of detected cases is applied between days 50 and 200 of the outbreak (as shown in A). The shown effect is achieved if 37’300 people per 100’000 are tested every day (5% false negative rate; test speed = 1 day). **Tables S2 and S3** show all other model parameters. Testing is discontinued after day 200. Dashed blue lines represent recovered plus infected plus deceased, dashed-dotted green lines infected, and solid red lines deceased people. Lower panels show the same data as upper panels, but with a log10-scale for Y-axes. We then compare the effects of extreme social distancing with those of extensive mass testing. Extreme social distancing is imposed for a period of 150 days (day 50–200 of the pandemic), reducing R to 0.7 (**table S1**). This is achieved if the rate of infection of susceptible people is reduced by 71% (**Fig. 2B**; data for R=1, as recently achieved in Switzerland (*8*), is shown in **Fig. S4A**). In line with current observations, extreme social distancing dramatically decreases the fraction of people experiencing an infection and reduces the number of deaths (< 0.0001% in **Fig. 2B**), at least if social distancing is initiated early on (e.g. by day 50) and strictly adhered to for months. The fraction of the infected population declines during extreme social distancing. This is in line with other modeling studies (*3*). However, once these measures are abandoned, the infection starts spreading again, leading to a similar death toll and ICU overloading as in the case without mitigation (**Fig. 2B**, days 200–600). Next, we analyze a mitigation strategy that is based on mass testing alone (**supplementary text** section 3). We would test random samples of the population with concomitant isolation of detected cases. The test is assumed to take one day to process, to yield 5% false negative results (discussed below), which is technically reasonable, and to result in immediate quarantine of the virus-positive cases. **Figure S3A** shows the factor by which R changes as a function of number of tests and for different test processing times; the same is shown in **Fig. S3B** for a higher false negative rate. Again, we assume that the mitigation strategy is implemented during days 50–200 of the pandemic. Applying 37’300 tests per 100’000 people per day yields the same outcome as extreme social distancing (R=0.7; **Fig. 2C**, compare to **Fig. 2B; supplementary text** section 3 for details). Reducing the mass testing to 12’600 tests per 100’000 people per day is sufficient to keep the number of infections constant (R= 1, **Fig. S4B**) and yields equivalent results as a moderate form of social distancing (**Fig. S4A**). Both interventions can dramatically slow the epidemiological dynamics, but similarly to the previously discussed social distancing, the number of infections and deaths will start to rise again after the testing regime is abandoned (**Fig. 2C**, days 200-600). It is easy to see that both, mass testing and social distancing, can have qualitatively identical effects on the epidemiological dynamic, but they work differently. Social distancing decreases the overall infection rate by reducing the freedom of movement for all, whereas mass testing will allow us to limit the isolation to the infected fraction of the population that would transmit the virus. The latter would have the important advantage of inflicting much smaller burdens on the economy than social distancing, as far smaller fractions of the population need to be isolated. However, the number of tests we predict to be necessary to reach R≤1 (12’600 tests per 100’000 people per day; **Fig. S4B**) is clearly unrealistic for now, and >50-fold above the current testing capacities even of countries with highly developed and well-funded healthcare systems. Optimizing the test specificity or the test speed changes the required number of tests by merely <2-fold (**Fig. S3**). Thus, mass testing alone is not an option. We therefore analyzed mitigation strategies combining virus testing with other measures. This should reduce the number of tests needed to achieve R ≤1. First, we tested a mitigation strategy combining mild social distancing, which reduces the infection rate by 33% (this would yield R=1.6, if applied alone; **table S1**), and mass testing (assumed test speed = 1 day, 5% false negatives). This may approximate a scenario where mass testing is applied to balance the effects of easing current social distancing measures. Using this combined strategy, we would require 4’500 tests per 100’000 people per day to reach R=1 (**Fig. 3A**). However, this still exceeds the current test capacities. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F3) Figure 3. Effective combinations of mass testing with other mitigation strategies can be achieved with realistic numbers of tests. **(A)** A combination of mild social distancing (33% reduction in infection rate; R=1.6 if it were applied alone) and mass testing (4’500 tests per 100’000 per day; false negative rate = 5%; test speed = 1 day) during days 50-300 can reduce R to 1. **(B)** Smart selection of the tested subpopulation reduces the required test numbers, if additional information allows identifying a subpopulation with higher than average infection prevalence (details, below). The Y-axis shows the factor by which the infection prevalence increases in the tested subpopulation. This will reduce the required number of tests. Blue dashed line: smart testing alone is used to achieve R=1 (if R=2.4 for the case without mitigation); red solid line: smart testing combined with mild social distancing (33% reduction in infection rate by social distancing). Green and red points illustrate examples discussed, below. **(C)** The same effect as in A can be achieved with smart testing, assuming that we can find a subpopulation with a 27-fold increased infection prevalence. Here, we need 475 tests per 100’000 per day to reach R=1 (green dots in panel B). (A,C): Dashed blue lines: recovered plus infected plus deceased; dashed-dotted green lines: infected; solid red lines: deceased people. Inserts show the same data, but with a linear scale for Y-axes. All presented cases assume false negative rate = 5%; test speed = 1 day. An additional mitigation measure is serological testing to identify the subpopulation that has already recovered from the disease and is now immune (*15*). This could be combined with mass testing to reduce the number of tests by removing the immune subpopulation from the pool of candidates to be tested. However, we are still at a stage of this pandemic where the fraction of immune people is likely quite small (see **Fig. 2A**). Therefore, serological testing could only minimally reduce the number of needed tests per 100’000 people per day. In no way are we implying, however, that serological testing is not useful: especially in the case of healthcare workers, and other essential personnel that is in close contact with the public, it is of great value to know who is immune to the disease. Also, serological testing could complement mitigation at later phases of such an epidemic when larger fractions of the population have become immune (**Fig. 2A**; >10% immune or deceased beyond day 140). Finally, we assessed how virus testing could be combined with Bluetooth-based tracing applications (e.g. PEPP-PT (*16*)), which are designed to identify individuals who had recent contacts with infected people and might therefore be infected, too. Contact tracing alerts the detected contacts to encourage isolation or virus testing. However, as this strategy misses most transmission events by asymptomatically infected cases and cases with mild symptoms that remained undetected (**Fig. 1**, dark blue), it is unlikely to suffice alone, given the epidemiological parameters in Europe (*3*) (**Fig. S5**, **supplementary text section 4**). Therefore, we explored an alternative strategy that relies on smart testing. In theory, smart testing would use additional information to focus virus tests on subpopulations having a higher prevalence than the overall population. We will first analyze the theoretical principle and then discuss a practical implementation. For smart testing to work, two crucial requirements have to be fulfilled. The prevalence in the tested subpopulation has to be high enough, so that a limited number of tests will suffice to detect the necessary number of infected people. Furthermore, the traced subpopulation featuring the desired prevalence has to be large enough. We describe the detailed derivation of the corresponding mathematical expressions in the **supplementary text** section 5. Smart testing alone can substantially reduce the number of virus tests needed to achieve R= 1 (**Fig. 3B**, dashed blue line). In fact, if we assume that smart testing identifies a subpopulation with 27-fold higher prevalence than the overall population, we can achieve R=1 with 475 tests per 100’000 per day (**Fig. 3B**, dashed blue line, green dot) and thus keep the number of infections constant as long as smart testing is maintained (**Fig. 3C**, days 50–300). Fewer tests would be needed, if we could further increase the virus prevalence in the tested subpopulation. Alternatively, one could combine smart testing with mild social distancing, e.g. to stop a surge of incidence when strong distancing measures are eased. For example, if social distancing would reduce the infection rate by merely 33%, we would only need 167 tests per 100’000 per day to reach R=1 (**Fig. 3B**, red line, green dot). Sufficient testing capacity for this combined strategy would already be available today in several European countries. We conclude that smart testing could work: it could achieve the same result as mass testing (**Fig. S4B**), but with much fewer tests (**Fig. 4B**). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F4) Figure 4. STeCC can achieve R=1 using a realistic number of tests and realistic rates of app-usage. **(A)** STeCC approach vs contact tracing. Contacts of a focal individual are logged using software. Using contact tracing, all contacts would be quarantined if the focal individual shows symptoms. Using STeCC, the focal individual is tested for infection if their contacts exceed a threshold number (e.g., if they visit a sporting event and a bar, in addition to being at home and at work), and quarantined if tested positive. (**B**) Progressive test cycling approach. Example for 80% app users among smart-phone owners and the availability of 165 tests per 100’000 people per day. (**C**) Covid-19 mitigation by STeCC alone with quarantine of contacts of virus-positive cases. Graph showing how the number of tests per 100’000 per day relates to the percentage of app users. We assume R=2.4 for the unmitigated case, test speed of 1 day, immediate notification and quarantine, false negative ratio 5%, asymptomatic cases (1/3 of all) are 50% as infective as symptomatic cases. Solid line: values to achieve R=1; dashed line: values to achieve R=0.8; dashed dotted line: values to achieve R=1.2. (**D**) Covid-19 mitigation by combining STeCC with app-based contact tracing. Graph showing how the number of tests per 100’000 per day relates to the percentage of app users. We assume R=2.4 for the unmitigated case, test speed of 1 day, immediate notification and quarantine, false negative ratio 5%, asymptomatic (1/3 of all) are 50% as infective as symptomatic cases. Solid line: values to achieve R=1; dashed line: values to achieve R=0.8; dashed dotted line: values to achieve R=0.6. The coloring indicates the achievable R with the respective mitigation strategy in place (see **supplementary text**, sections 6 and 7). Can we use Bluetooth-based tracing apps for identifying high prevalence subpopulations? A first attempt at achieving this would be an adaptation of app-based contact tracing to identify individuals with recent infection-relevant contact to detected Covid-19 cases (*4*). This contact subpopulation should have a higher prevalence than the overall population. However, such a "trace+test" approach may not suffice to achieve R=1, at least when applied without any other mitigation strategies (*3*). This is explained by its focus on contacts of symptomatic individuals (while contacts to asymptomatic cases are missed), as well as the virus’ infection and transmission dynamics (for details, see the **supplementary text** section 6). Contact tracing is therefore not sufficient to reach the requirements for smart testing to work as a standalone mitigation strategy. As many asymptomatic SARS-CoV2 carriers contribute to the transmission, we reasoned that smart testing should cover symptomatic and asymptomatic cases, alike. To solve this problem, we suggest a different strategy for identifying high-prevalence subpopulations. Importantly, this can be achieved using the same Bluetooth based technology as contact tracing. Our strategy relies on the fact that some individuals will have many more infection-relevant contacts than most others (*17*). For example, with realistic scale-free network model assumptions, ≈1% of the population with the most infection-relevant contacts (at least 10 times more than the average) make the difference between R=1 and R=2.4 (see **supplementary text** section 6). These individuals can be identified by counting of infection-relevant contacts (i.e. contacts which could facilitate transmission), regardless, if contacts were infected or not. Their sheer number of contacts makes high-contact individuals much more likely than others to become infected and to transmit the virus, and they are known to be highly important for epidemiological dynamics ((*18*), referred to as super-spreaders). To test if smart testing with contact counting (STeCC; **Fig. 4A**) works in realistic scenarios, we used a scale free network model, which can account for the heterogeneity of the number of contacts in a population (**supplementary text** section 6). The fraction of app-users in the population will also affect the success of STeCC. We assume that children under age 10 contribute little to transmission in the general population (*19*, *20*) and would not use smartphones, as suggested by others (*3*). High-risk persons (i.e. >70 years old; 20% of the overall population) would remain shielded in isolation for safety reasons and would therefore not contribute to the infection dynamics. People aged 10–70 (in particular the high-contact subgroup) would use smartphones, and we assessed how their app-usage would affect mitigation. Then, we tested three different scenarios how STeCC could be implemented to reach R=1. We used a progressive testing cycle, as testing the same high-contact individuals every day would seem unrealistic. On the first day, based on the number of available tests, we would identify the optimal subgroup with the highest frequency of infection-relevant contacts, and invite them to be tested (**Fig. 4B**; dark blue; **supplementary text** section 6). During the next days, we would repeat the same procedure, but exclude the groups previously identified. After 7 days, we would begin the cycle again and identify the optimal high-contact group from the entire susceptible population. Using this testing cycle, we observed that STeCC alone could reduce R, but cannot reach R=1 with realistic numbers of app-users and tests per day (**Fig. S6**). Next, we assessed a STeCC strategy, which would refer two groups of people to quarantine, i.e. the high-contact individuals that tested positive and their recent contacts. This strategy would reach R=1, if 90% of smartphone users would use the app (63% of the overall population) and ≈397 tests per 100’000 people per day were carried out (**Fig. 4C**, red lines). Finally, we tested a combination of the latter STeCC approach with Bluetooth app-based contact tracing (similar to (*4*)), which may be possible using the same data. This would refer three groups of people to quarantine, i.e. the high-contact individuals testing positive, their recent contacts and the contacts of symptomatic individuals. If 72% of smartphone users would use the app (≈50% of the overall population) we would need ≈166 tests per 100’000 people per day to reach R=1 (**Fig. 4D**, red lines). This testing capacity is currently available in several European countries, like Switzerland. Please note that our example in **Fig. 4D** assumes an R=2.4, if no mitigation was applied. An identical 2.4-fold reduction is achievable for any other basic virus reproduction rate, as long as 72% of smartphone users would use the app and ≈166 tests per 100’000 people per were carried out. We have conducted sensitivity studies to study effects of the false negative rate, network parameters, test processing times and R for the unmitigated case. In all considered variations, STeCC with contact tracing could achieve R=1 with reasonable fractions of app-users and test numbers (see **supplementary text** section 6). Thus, STeCC with contact tracing could add to any mitigation policy. What are the advantages of a STeCC-based mitigation strategy? First, it includes detection and removal of asymptomatic cases and their contacts. In contrast to contact tracing alone (**Fig. S6**), the combination with STeCC would work when >72% of the smartphone users (≈50% of the overall population) would use the app, or if unfavorable assumptions about the pandemic’s parameters would turn out to be true (**Fig. S8**). This synergy is attributable to the different selection processes of both approaches. Furthermore, STeCC could be implemented within the same apps as developed for contact tracing. This could offer additional choices for app-based mitigation, i.e. the parallel use of STeCC and contact tracing, which is particularly powerful (**Fig. 4D**). As STeCC focusses on a small high-contact group, it works with relatively small numbers of tests and would only quarantine those that are virus positive (and their contacts, e.g. **Fig. 4C,D**). It is realistic given the available test capacities in several countries and would affect smaller fractions of the population than other mitigation strategies (i.e. contact tracing or social distancing). Thereby, one could enhance the impact of Bluetooth-based tracing applications. Due to its key function in virus transmission, the identified high-contact subgroup would also be a promising priority for vaccination, once limited supplies of a vaccine become available. What are the limitations of STeCC? First, we cannot exclude that some assumptions used in our model may be too optimistic or that more precise information might be derived later from alternative, more detailed simulation approaches. However, modeling the impact with less favorable parameters verified that STeCC would still provide substantial benefits (**Fig. S7**, **Fig. S8**). Second, contact counting has not been a focus during the development of Bluetooth-based proximity testing applications. Thus, small adaptations might be needed, in order to enable efficient detection and notification of high-contact individuals. Strict social distancing has been successful in achieving R<1 in many countries, but at a high economical and societal cost. Easing of these measures is presently being discussed or implemented. However, if a large fraction of the population has remained susceptible, a second wave of disease is bound to occur in the absence of effective alternative mitigation strategies. We suggest using a combination of contact tracing and STeCC, as a simple mitigation approach which relies on identifying high-contact individuals, testing them for infection, and quarantining positive cases. STeCC requires the same information as contact tracing (a list of unique contacts in a given period of time), and can be implemented using the same information as the contact tracing apps that are currently being developed or already in use. STeCC would be achievable using a number of tests that is realistic today in several European countries, and can be achieved in many others with appropriate efforts. STeCC adds to the portfolio of mitigation strategies for the Covid-19 pandemic. It could be deployed quickly in countries with sufficient testing capacities like Switzerland (capacity ≈230 tests per 100’000 per day). A combination with contact tracing might be particularly powerful. Our study provides quantitative estimates for the number of tests needed by starting out with realistic assumptions about the relevant parameters, like number of app users. Once STeCC is applied, one can adjust the strategy flexibly in order to ensure the desired performance. STeCC offers a realistic approach to help relaxing broad social distancing policies in the near future without compromising health, while at the same time providing public health officials with much needed actionable information on the success of their interventions. This will be an important prerequisite for reclaiming our normal public life and initiating economic recovery. ## Data Availability All data used for preparation of this manuscript is publicly available. ## Methods The dynamic model was implemented with Maple 2018. The calculations for mass testing, contact tracing and smart testing were implemented with MATLAB and the Statistics Toolbox Release 2018b. They will be made available at GitHub. ## Supplementary Text for: STeCC: Smart Testing with Contact Counting Enhances Covid-19 Mitigation by Bluetooth App Based Contact Tracing ### 1. Model A model is proposed to compute the numbers of infected people whose infection has not (yet) been detected and the numbers of infected persons with a detected infection (![Graphic][1] and ![Graphic][2]), respectively (Fig. 1). Note that detected here refers to persons being isolated, which comprises not only those being tested positive, but also those who feel strong symptoms and thus stay in self-quarantine. It is further important to notice that in the case of SARS-CoV2, the undetected infected people are main contributors to the spread of the pandemic [13]. The exact definitions of detected and undetected, as well as those of all other variables and model parameters are found in Table S1. Furthermore, we compute the number of fatalities (*nd*) and the number of people who recovered after a detected or an undetected infection. Importantly, we assume that these people will have developed protective immunity and we assume that they cannot be infected again in the considered time frame. The initial susceptible population ![Graphic][3] is naive (i.e. it lacks immunity against the infection) and *ns* is the number of persons who are susceptible at a given time *t*. In our model we assume that the virus is mainly passed on by undetected asymptomatic and mild symptomatic persons; the detected population with mild symptoms transmits at a much lower rate (because of self-isolation, hygiene precautions in hospitals and/or quarantine). While the model does not consider age dependency, it accounts for higher mortality rates due to temporary shortage of intensive care units. The graph in Fig. 1 (main part) shows the dynamic dependencies. The Covid-19 specific parameters have to be estimated from the available data; their values are listed below. It should be noted, that the implementation of our model allows updating our current estimates with more precise values, as new data come in. Initially, the entire population is susceptible and can get infected. Infected persons first get exposed and are not infectious before the latency time has passed. Then, they either become asymptomatic or mild symptomatic. Asymptomatic persons eventually recover without symptoms, while the others develop symptoms approximately half a day after the end of the latency time. We assume that persons with mild symptoms isolate themselves approximately one day after onset of symptoms and then either recover or become strong symptomatic, which requires hospitalization. Hospitalized individuals either recover or die. Once *ns* becomes smaller, which happens quickly without any measures, the infection rate slows down by a factor of ![Graphic][4]. This mechanism of slowing down spread of the epidemic due to a shrinking susceptible population is equivalent to herd immunity. It is crucial for the system dynamics that detected persons are isolated (either by self-isolation at home, by hygienic isolation in hospital setting or in other care facilities, or by organized isolation programs for detected infected people, e.g. in hotel rooms) and thus participate at a much smaller rate or not at all in spreading the disease. We assume that these detected infected people have a 10-fold lower likelihood of infecting others than undetected infected people. All this leads to a dynamic system, which is governed by the following ordinary differential equations: View this table: [Table S 1.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/T1) Table S 1. Terminology and nomenclature of model parameters and variables ![Formula][5] ![Formula][6] ![Formula][7] ![Formula][8] ![Formula][9] ![Formula][10] ![Formula][11] ![Formula][12] ![Formula][13] Note that ∊ ∈ [0,1] is the transmission reduction factor of the self-isolated individuals. It is assumed that those infected persons who were detected by testing or hospitalized infect much less due to strong isolation and other precautions. Therefore, their effect on the infection rate is neglected here. In order to parametrize the model, besides the rates *α*, *βa, βs, βa, γss, γss*, *ξms, ξss* and *θ*, also the relative rate *Q* of infections from outside, i.e., by travel or from the animal world, has to be determined. Further, the initial values of *ns*(*t*)*, ne*(*t*), *nia*(*t*)*, nim*(*t*), *nra*(*t*)*, nms*(*t*)*, nss*(*t*),*nrs*(*t*) and *nd*(*t*) have to be chosen. The variables ![Graphic][14], ![Graphic][15], ![Graphic][16], ![Graphic][17], ![Graphic][18], ![Graphic][19], ![Graphic][20] and ![Graphic][21] denote the respective numbers of persons who were tested positive and thus are removed from transmission. The detection rates of exposed (*ne*), asymptomatic *(nia)* and mild symptomatic *(nim)* persons due to testing are proportional to *ke, ka* and *ks,* respectively. These individuals are then accounted for by the respective numbers ![Graphic][22], ![Graphic][23] and ![Graphic][24]; see Fig. S1. Note that the graph in Fig. S1 is very similar as the one in Fig. 1, except that there is no node for susceptible persons (since by definition a susceptible person can not be detected infected) and that there exist sources due to testing (dotted arrows) instead of sinks. To account for the dynamics with testing, the system (1)-(9) has to be augmented by the ordinary differential equations ![Formula][25] ![Formula][26] ![Formula][27] ![Formula][28] ![Formula][29] ![Formula][30] ![Formula][31] ![Formula][32] The effect of testing is further discussed in Section 3. Next it is described how the parameters can be estimated based on literature data. ### 2. Parameter Estimation Our devised generalized SEIR model becomes closed once we tune the rate coefficients. These coefficients were computed mainly based on data provided in recently published reports [6, 27]. Before giving the values for transfer rates between different compartments, let us analyze the basic reproduction number ![Graphic][33] of this virus infection with *Q* = 0 and ![Graphic][34]. Note that ![Graphic][35] represents ”the expected number of secondary cases produced, in a completely susceptible population, by a typical infective individual” [4]. If ![Graphic][36] becomes < 1, virus spread will decline, and if ![Graphic][37], virus spread will increase. To compute ![Graphic][38], we split the dynamics of the infected population into the infection driven propagation *f* and the remainder *V*, i.e., ![Formula][39] ![Fig. S 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F5.medium.gif) [Fig. S 1.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F5) Fig. S 1. Graph showing the dependencies of the compartments describing the dynamics of the positively tested people. Note that testing is not considered here, i.e., *ke*, *ka* and *ks* are zero. The ![Graphic][40] of this system is the spectral radius of *fV*−1, that is, ![Formula][41] By inspecting Eq. (19), we observe that we can move towards a stable state (corresponding to ![Graphic][42]) by reducing the infection rate *α* via mitigation policies such as social distancing. Importantly, as shown later, ![Graphic][43] can be reduced as well by introducing mass testing, contact tracing, smart testing and subsequent isolation of detected infected individuals. Next, to clarify our choice of model coefficients, we discuss the rates which appear in transmissive and non-transmissive compartments separately. Finally, the increase of mortality due to lack of intensive care units is modeled. 1. *Transmissive:* We model the incubation time to be log-normally distributed with mean 5.84 (day) and standard deviation 2.98 (day) [25, 18]. In accordance with [6] we take the latency time *xl* such that in average it becomes half a day shorter than the incubation time. Similar to the incubation time, we adopt a log-normal distribution for *xl* but with mean 5.34 (day) and standard deviation 2.7249 (day). We assume that 1/3 of the cases won’t have noticeable symptoms and 2/3 become symptomatic half a day after latency [6]. This leads to ![Graphic][44] (1/day), where ![Graphic][45] denotes the expectation which gives us the average latency rate. Due to the ratio of 1/3 to 2/3 between asymptomatic and symptomatic cases we get the transfer rate from being exposed to infectious symptomatic as ![Graphic][46] (1/day). We suppose that it takes around 1 day from onset of symptoms to self-isolation [6]. Since it takes half a day time delay from becoming infectious to symptomatic, we get ξ*ms* = 1/1.5 = 0.6667 (1/day). The average onset to discharge time of clinical cases is around 22 days [27]. We assume that for mild-symptomatic cases the onset to recovery time would be half of this amount, i.e., around 11 days. Therefore the average recovery time from end of the latency period becomes 11.5 days for mild-symptomatic cases. We set the same recovery time for asymptomatic cases which leads to *γa* = 1/11.5 = 0.087 (1/day). A range of values have been suggested for infectiousness of asymptomatic cases; one finds 0.1 in [7], 2/3 in [6] and 1 in [12]. We assume that the asymptomatic cases are 50% less infectious. Furthermore we consider the self quarantined patients to be 90% less infectious, i.e., *∊* = 0.1 is adopted. To compute the infection rate *α*, we assume ![Graphic][47] [6]. Following Eq. (19), the infection rate becomes *α* = 0.6711 (1/day). Since the basic reproduction number is the most important single parameter of the system, we performed sensitivity studies by changing ![Graphic][48] [14]. 2. *Non-transmissive:* The mean delay time from appearance of symptoms to hospitalization has been reported to be around 11 days [27]. However, 80% of symptomatic cases would not require hospitalization [15]. For those who develop strong symptoms, the delay from self-isolation to hospitalization then becomes 11 − 1 = 10 days. Hence we get *ξss* = 0.2 × 1/10 = 0.02 (1/day) and *γms* = 0.8 × 1/10 = 0.08 (1/day). Note that the latter gives onset to recovery time of 11 days for mild cases consistent with our earlier assumption. The average hospital treatment time is 11 days [27, 6]. In case of availability of intensive care units we assume that 20% of hospitalized cases die [15]. Accordingly, we get ![Graphic][49] (1/day) and *θ*(0) = 0.2 × 1/11 = 0.0182 (1/day). 3. *Fatality increase:* We assume that the case fatality ratio increases by two-fold in saturation of the health system. This is justified by noting that the case fatality ratio has increased from approximately 5% in China [2] to roughly 10% in Wuhan while it was the epicentre of the outbreak [20]. By taking this factor into account, and assuming that the average time of hospital treatment remains 11 days, we can compute the death rate of hospitalized cases once saturation of intensive care units is reached as *θ*(*sat*) = 0.4 × 1/11 = 0.0364 (1/day). Note that consistently one obtains ![Graphic][50] (1/day). It is assumed that there exist eight intensive care beds per 100’000 persons [6] and that 40% of the hospitalized cases need such treatment. Accordingly, saturation is reached once the number of hospitalized cases, i.e., *nss,* exceeds *C*(*sat*) = 0.02% of the total population. The adjusted rate ![Formula][51] then quantifies the death rate of hospitalized cases as the weighted average of *θ*(0) and *θ*(sat); the consistently adjusted recovery rate becomes ![Formula][52] All estimates here are summarized in Table S2; note that these values can easily be adapted, if more reliable data becomes available. The resulting parameter values for our base case are provided in Table S3. Figure 2A shows the model results without mitigation for a period of 600 days with ![Graphic][53]. Results with ![Graphic][54] are shown in Fig. S2A,B. Dashed lines represent the immune plus deceased plus infected ![Graphic][55], dash-dotted lines the infected ![Graphic][56] and solid lines the deceased ![Graphic][57] population. For bigger values of ![Graphic][58], a larger immune population is needed to achieve herd immunity (right half of the graphs), and the peak in the number of infections is higher and sharper. The plot in Fig. S2C shows the case with ![Graphic][59] without intensive care unit limitation; compare with Fig. 2A, which shows the same case with intensive care unit limitation (our base case). In both cases 87% of the population will become immune, which compares well with 81% infected people predicted by [6] for the UK and US populations in the absence of mitigation plans. Without intensive care the chance of dying is roughly twice as high for strong symptomatic people than with proper treatment (4.6% vs. 2.3%). While these numbers are subject to errors (mainly due to uncertainties in the parameter values and efforts to increase intensive care and respirator availability), it can be expected that the relevant dynamics is captured to a high degree. If the results are regarded with respect to the base case, much insight can be gained, e.g. how social distancing, mass testing and smart testing can be combined most effectively. View this table: [Table S 2.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/T2) Table S 2. Estimations made for the model closure, social distancing, contact tracing and testing. View this table: [Table S 3.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/T3) Table S 3. List of estimated parameters and initial values. Note that our model allows to easily replace any of these parameters by more precise estimates, as more data become available. The initial values of all numbers except *ne* are set to zero. ![Fig. S 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F6.medium.gif) [Fig. S 2.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F6) Fig. S 2. Alternative model outcomes when changing ![Graphic][60] of the pandemic, or relaxing the assumption that ICU beds are limiting. Dashed blue lines represent recovered and deceased, dashed-dotted green lines infected, and solid red lines deceased people. Changing ![Graphic][61] to (A) 1.9 and (B) 2.9 changes the outcomes quantitatively, but does not change the overall picture. (C) Model outcome if the assumption that full ICUs increase the death rate is dropped. ### 3. Mass Testing - How the Number of Tests Relates to ![Graphic][62] Here we analyze how many tests are needed to mitigate the Covid-19 pandemic if no other mitigation strategies were applied. If we can use more tests than that, ![Graphic][63] will drop below one, the incidence would decline and the pandemic would eventually end, even if no vaccines or infection therapies become available. Therefore we studied how ![Graphic][64] varies once confirmed cases get isolated (in addition to self-quarantined and hospitalized individuals). Of particular value is the relationship between ![Graphic][65] and the interval of testing the susceptible population (i.e., the frequency of testing needed for reducing ![Graphic][66] below one). Thereby we can determine the key technical parameter of interest, i.e., the number of tests per 100’000 people that must be tested per day in order to achieve the desired ![Graphic][67] value; we chose ![Graphic][68] as the target value for our analyses (if not indicated otherwise), which would suffice to keep the number of infected people constant. To be realistic, we suppose that the processing time *τproc* of mass testing would be somewhere between half a day and two days. Furthermore, a fraction *η* = 0.05 of false negative test results is taken into account [17]. It should be noted, however, that this is a rough estimate. The use of standards allows for very high reproducibly of the virus RNA detection results even between different laboratories [3]. The true rate of false negatives and false positives is currently not known. We assume that a false negative rate of 5% is a conservative estimate. The current virus RNA testing capacities in continental Europe reach up to 230 tests per 100’000 people per day (e.g. in Switzerland). If equipment and supplies are not limiting for testing, e.g. by using a quantitative polymerase chain reaction (qPCR) method1, we estimate that up to around 1’000 samples within a time frame of eight hours can be analyzed per machine. Mass testing (i.e., if > 500 − 1ʹ000 tests per 100’000 perople per day would be required) could be realizable by taking advantage of next-generation RNA extraction, reverse transcription and sequencing (combined with reverse transcription and PCR) to detect the virus RNA of infected people. For example, in [11] a massively parallel diagnostic assay is described for testing up to 19’200 patient samples per work flow. In principle, such very high throughput approaches can be parallelized (and potentially optimized) to provide millions of tests per day. In reality, the logistics of collecting these millions of samples would however be a major hurdle. First, we want to assess how many tests were indeed necessary to stop the virus spread (i.e., to reach ![Graphic][69] if no other mitigation strategies were applied). Obviously, the scenario that the whole susceptible population is tested perfectly at once would lead to a trivial disease free state. However, this is an unrealistic scenario, not only because of a lack of test capacity, but also due to logistic and compliance concerns. Therefore, it is only realistic to assume that individuals would be tested at different schedules. Let us consider a situation where each person is tested once every *N* days. Note that this is equivalent to testing a random fraction of 1/*N* of the susceptible population every day. Therefore we focus on the set {1,..., *N*} of days. An individual is infected at some random time *x*. To characterize *x*, we assume that the likelihood of getting infected does not vary much during these days, which is justified if testing is applied at an intensity such that ![Graphic][70]. Hence *x* becomes uniformly distributed in the interval [1, *N*]. The time delay between infection and detection would then be *τdet* = *N* − *x* + *τproc*. Finally, we sample the latency time *xl* from a log-normal distribution with 5.34 ± 2.7249 (see Section 2). Intuitively, by conducting mass testing on individuals who are neither self-quarantined nor hospitalized, positive cases will be detected from exposed, asymptomatic and mild-symptomatic compartments. To quantify each detection rate, it is essential to compare detection time versus the latency period; therefore we consider the effect of testing on these three compartments individually: 1. *Exposed:* Once testing occurs during the latency period of an infected individual, they would be detected from the exposed compartment. This translates into an event set ![Graphic][71]. The rate of detecting individuals by testing from the exposed population then reads ![Formula][72] where ![Graphic][73] and ![Graphic][74] denote frequency of such events and conditional expectation, respectively. 2. *Mild-symptomatic:* Once testing occurs after the latency period, one has to distinguish between two types of infection developments. According to our setting, two thirds of the infected individuals would develop symptoms that will lead them to self-isolate. Please note, that the fraction of infected individuals that remain asymptomatic may range between 33%-50% [6, 7, 24]; see Table S2. Therefore, we have also tested scenarios where 60% or 50% of infected will will progress to develop symptoms; see Fig. S5. These individuals may in fact turn to testing centers in order to detect the virus and incentivize the decision for self-quarantine. They can be detected by testing and therefore sent into quarantine in the span of one and a half days after becoming infectious. The relevant event set is ![Graphic][75], from which one obtains ![Formula][76] for the test detection rate of mild-symptomatic persons. 3. *Asymptomatic:* This is arguably the most important group to consider, because they will not know that they are infected and they can make up as much as 50% of the entire group of the infected people. Importantly, they are very hard to identify using the current mitigation strategies, including app based contact tracing [7]. We will discuss this in more detail, below. Testing may catch asymptomatic cases. These individuals won’t have symptoms and would then recover after 11.5 days. In their first one and a half days they share same the time-line as mild symptomatic ones. Therefore, in this scenario two event sets ![Graphic][77] and ![Graphic][78] become relevant. We obtain the detection rate from the asymptomatic compartment as ![Formula][79] Before finding the map from testing frequency *N*−1 to ![Graphic][80], we need to find out how ![Graphic][81] varies with respect to the test detection rates *ke, ka* and *ks.* Hence, let us define ![Graphic][82] as the reproduction number subject to testing. To compute R“, the main dynamics of the infected population, which can be described by ![Graphic][83], is split into the rate of appearance *f* of new infected individuals and transfer ![Graphic][84] of infected ones across different compartments, i.e., ![Formula][85] and ![Formula][86] which leads to ![Formula][87] It is important to emphasize that in contrast to social distancing the effect of testing on the reproduction number is not through reducing the infection rate (i.e., *f* remains the same), but rather through transfer of infected individuals to quarantine, either at home or in hotel rooms, or in care units with strict hygiene barriers, all of which would reduce the likelihood of transmission (i.e., changing *V* to ![Graphic][88]). Using Monte-Carlo to estimate the test detection rates given by Eqs. (22) and (23) one can compute the ratio ![Graphic][89] with respect to the testing frequency *N*−1. The plots in Figs. S3A,B show the number of tests required to reach ![Graphic][90]. It depends on the time from sampling to result (note: we assume immediate notification and immediate implementation of quarantine measures upon notification) and on the false negative rate (5% in Fig. S3A and 15% in Fig. S3B). The horizontal green dashed lines indicate ![Graphic][91], if the virus reproduction rate in the case without any mitigation is ![Graphic][92]. ![Fig. S 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F7.medium.gif) [Fig. S 3.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F7) Fig. S 3. Mass testing: The number of tests required to reach ![Graphic][93] depends on the time from sampling to result. A mitigation strategy relying on mass testing alone is assumed and we computed the number of tests performed per day, which are needed to achieve a particular test-speed dependent ![Graphic][94] ratio; for (A) 5% and (B) 15% false negative test results are assumed. The green lines indicate ![Graphic][95]. Test speeds were: 0.5 days (black line), 1 day (orange line), 1.5 days (blue line) and 2 days (purple line). One can see for example that reducing the initial ![Graphic][96] to one would require to test the entire sceptible population roughly once every eight days, if the false negative rate is 5%, if testing results are available after one day and quarantine of the detected individuals commences immediately. Corresponding results are depicted in Fig. S4B, which shows model outcomes if ![Graphic][97] is reduced to one by testing 12’600 per 100’000 people per day. Note that this corresponds to a testing interval of 7.92 days or equivalently to a fraction of 1/7.92 which has to be tested every day. Figure S4A shows the outcome when social distancing is applied to reduce ![Graphic][98] to one, which is equivalent to a 58% reduction of the infection rate. We define this as moderate social distancing; see Table S1. By comparing Figs. S4A and S4B one observes that mass testing and social distancing yield qualitatively equivalent results. For half a day delay time the testing interval can be increased to roughly twelve days and for a delay time of one and a half days it would be five days. This information is important for optimizing technical development decisions. The overall testing capacity needs to be larger, if the testing method requires more time (Fig. 3) or if delays in case-notification or implementation of quarantine measures occur. This could be justified, if a low-cost technique (such as next generation sequencing) could be devised [11]. Alternatively, one could aim for fewer tests, if they are completed in less than half a day. If combined with social distancing, the basic reproduction number can be reduced by the factor ![Graphic][99], where *λ* ∊ [0,1] is the intensity of social distancing *(λ* = 0 means no social distancing and *λ* = 1 means complete isolation of everybody). ![Fig. S 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F8.medium.gif) [Fig. S 4.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F8) Fig. S 4. Moderate social distancing vs. mass testing: Model outcomes if ![Graphic][100] is reduced to 1 by mass testing or social distancing. (A) Moderate social distancing (![Graphic][101], equivalent to a 58% reduction in infection rate) and (B) testing 12’600 per 100’000 people per day yield qualitatively equivalent results. ### 4. Contact Tracing - How Size and Infectiousness of the Asymptomatic Population Relates to ![Graphic][102] Contact tracing has been proposed to slow down or even stabilize the pandemic [7, 10, 8]. The strategy is that symptomatic individuals who go into self quarantine will use an App to alert all proximity contacts of the past two weeks. Subsequently, these identified individuals self-isolate themselves. For the following analysis we assume that the notice-to-quarantine time delay is negligible [7, 10]. To analyze the effect of contact tracing on the basic reproduction number, the expected infectiousness *A*(*τ*) at time *τ* after infection plays a central role. By knowing *A*(*τ*) we can extract the basic reproduction number as ![Formula][103] Note that the above equation is consistent with our previous computation of ![Graphic][104] in the linear regime [4]. In our compartmental setting shown in Fig. 1 one obtains the expression ![Formula][105] where ![Graphic][106] is the probability that an infected individual is in compartment *ci,* ![Graphic][107] is the index set of all compartments and *Ai*(*τ*) denotes the infectiousness at time *τ* after infection, if the individual is in *ci.* Correspondingly, we denote the time spent in each compartment as ![Graphic][108]. From Eqs. (28) and (29), if we assume that infectiousness is constant inside each compartment, we obtain ![Formula][109] where *r*1 = 1 − *S*(*m*) is the fraction of exposed people who develop no symptoms (thus remain asymptomatic), and *r*2 is the factor by which asymptomatic people are less infectious than symptomatic ones. For both *r*1 and *r*2 different values are suggested in the literature. For *r*1 one finds 1/3 in [6], 0.4 in [7] and 0.5 in [24], and for *r*2 one finds 0.1 in [7], 2/3 in [6] and 1 in [12]. Based on these published numbers we consider *r*1 ∊ [0.3, 0.5] and *r*2 ∊ [0.1, 0.5]. From Eq. (30) with our base case values *r*1 = 1/3 and *r*2 = 0.5 one obtains *α =* 0.67 (1/day), which is consistent with our previous parameter estimation. Based on our previous assumptions and on the values in Table S2 we obtain ![Graphic][110] days, ![Graphic][111] days and ![Graphic][112]. Interesting here is the fraction ![Formula][113] of infected people who got infected by symptomatic cases, since this is the maximum relative reduction of ![Graphic][114] which can be achieved by tracing contacts of symptomatic individuals (by classical contact tracing or by using an app) with a success rate of *ζ* = 1. Ignoring secondary infections (their probability becomes around 0.8%) we obtain the approximation ![Formula][115] for the basic reproduction number, if contact tracing is employed. Figures S5A and S5B show the performance of contact tracing for different *r*1-, *r*2- and *ζ*-values. The numbers attached to the isolines refer to the ratio ![Graphic][116], and the bold contours depict combinations of *r*2 and *ζ* for which ![Graphic][117] is reduced from 2.4 to ![Graphic][118]. One can conclude that the effect of contact tracing, if only applied to identify contacts with symptomatic persons, is limited to optimistic assumptions concerning the parameters dictating the Covid-19 pandemic and strongly depends on size and infectiousness of the asymptomatic population relative to size and infectiousness of the symptomatic one. Even if the most optimistic assumptions would hold, more substantial reductions of ![Graphic][119] would be desirable in order to accelerate the end of the pandemic (e.g. even in the absence of effective therapies or vaccines). For our base parameters of *r*1 = 1/3 and *r*2 = 1/2, the contact tracing alone would not lead to basic reproduction number of one (see Fig. S 5A at 50% infectiousness). ![Fig. S 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F9.medium.gif) [Fig. S 5.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F9) Fig. S 5. Contact tracing: Effectiveness of contact tracing as function of relative infectiousness of asymptomatic cases (*r*2) and the success rate of contact tracing (*ζ*). (A): 33% of infected ones are asymptomatic (*r*1 = 0.33); (B): 50% of infected ones are asymptomatic (*r*1 = 0.5). The numbers attached to the isolines refer to the ratio ![Graphic][120], and the bold contours depict combinations of *r*2 and *ζ* for which ![Graphic][121] is reduced from 2.4 to ![Graphic][122]. In all computations, notice-to-quarantine time was neglected and it was assumed that quarantined contacts are not infectious. ### 5. Smart Testing - How Selectivity Relates to ![Graphic][123] Here it is studied how the number of required tests can be reduced, if one is able to identify (and propose testing to) a subpopulation with a prevalence higher than the overall population. Without any additional knowledge, to achieve the discussed detection rates *ke*, *ka* and *ks* and corresponding reductions in ![Graphic][124], one has to test the entire undetected population once every *N* days (or equivalently every day a random fraction of 1*/N*). To improve the efficiency of testing, i.e., the probability per test of getting a positive result by avoiding unnecessary testing of people who have a low likelihood of being infected, one can reduce the sample population by the following approaches: 1. *Serological testing:* With serological testing one can remove the recovered population from the pool of undetected individuals, and thus the same number of positive test results can be achieved with fewer tests. However, during an early stage of the pandemic the relative size of the undetected recovered compared to the whole undetected population is very small, and therefore the gain would be negligible. Nevertheless, this approach can easily be integrated into any mass testing strategy once reliable serological tests become available. Here, we can assume that immune individuals will remain immune for an extended period of time (e.g. up to 1–2 years; however, this is still subject to verification). If immune people cannot be infected for a second time (or are infected at a much lower rate than susceptible individuals), one can collect the information on positive test results to exclude the immune individuals from the Covid-19 RNA testing. 2. *Inference from contact tracing:* A more effective approach would be based on an inference model (e.g. by using contact tracing of infected individuals), which allows to divide the sample population ![Graphic][125] into one subpopulation ![Graphic][126] with a higher and the remainder with a lower percentage of infected individuals. For the following analysis we denote the size of ![Graphic][127] with *n* and that of ![Graphic][128] with *ń*. Further, ![Graphic][129] is the fraction of infected persons in ![Graphic][130] and *p* that in ![Graphic][131]. Testing every person in ![Graphic][132] at a frequency of 1/*N* would require *ń/N* tests per day; opposed to *n/N* tests, if the whole sample population was tested. The respective numbers of positive test results per day, on the other hand, would be ![Graphic][133] opposed to *P*(*N*) = *pn/N*. In order to obtain the same number of positive results from the subpopulation ![Graphic][134] as one would get from ![Graphic][135], one has to reduce the test interval *N* to *Ñ*, such that ![Graphic][136]. From this one obtains ![Graphic][137] and one can conclude that the number of tests required to achieve the same overall quota reduces by the factor ![Formula][138] In order for this result to be practically meaningful, *Ñ* has to be at least one, which translates into the requirement that ![Formula][139] In short, if one can identify a subpopulation ![Graphic][140] for which the percentage of infections is higher by a factor of *r*−1 than in ![Graphic][141], and which is larger than *nr/N,* then the number of tests needed to obtain the same reproduction number reduces by the factor *r*. The curves in Fig. 3B show the relationship between number of tests per 100’000 people per day needed to achieve ![Graphic][142] and the prevalence ratio between sub- and overall population; in combination with mild social distancing (solid line) and without social distancing (dashed line). ### 6. Contact Counting - How to Screen Large Enough Subpopulations with High Prevalence Here we devise a way to screen large enough subpopulations with a much higher prevalence and infectiousness than the overall population, which is a prerequisite for smart testing. We study three subpopulations as potential candidates for our smart testing mitigation approach; of interest are their prevalence, their infectiousness and their size. Next we describe them and provide quantitative estimates for the most relevant subpopulation. 1. *Contacts of symptomatic cases:* A straight-forward approach would be to choose the contacts of symptomatic cases as our subpopulation. While this group is highly likely to be infected, this approach has one major drawback. In fact, the outcome of testing contacts of symptomatic cases would not be much different than that of contact tracing mentioned before (and discussed in detail by Ferretti et al.[7]). Therefore, it suffers from the same limitation of not catching sufficient numbers of asymptomatic infections. Besides tracing contacts which potentially got infected by a symptomatic individual, one may also find the contact by whom it got infected. That person has most likely recovered, since he/she got infected roughly 10–14 days ago. Therefore this contact would not be tested positive (as virus titers may already be low and as we still lack reliable serological tests) and hence testing contacts of symptomatic ones would not lead us to a larger group with a sufficient number of asymptomatic cases. 2. *Direct and indirect contacts of symptomatic cases:* One way to cope with the issue arising from lack of enough asymptomatic cases in the contacts of symptomatic ones is to enlarge our sample population and include also indirect contacts of symptomatic individuals in the past two weeks. This strategy, while most probably catching enough asymptomatic cases, may not reduce the burden of mass testing, since now the size of the subpopulation becomes simply too large. This problem has also been noted by others [10]. For example at a prevalence of 1% in the total population, and assuming 10 contacts per person, the size of this subpopulation becomes almost as large as the whole population. 3. *High-contact individuals* (Fig. 4A): In this scenario we only test those with significantly more contacts than the average. In the following we show that indeed this strategy allows to screen a high prevalence subpopulation which is also large enough to stop the pandemic. It is also important to emphasize that the prerequisite of this strategy is to utilize a contact counter, which may be integrated into an existing contact tracing app that uses bluetooth technology. While the improvements resulting from contact counting can be estimated based on the prevalence ratio, the contact counting scheme has a more fundamental feature that exhibits itself directly in the basic reproduction number. In fact by cutting out the highly transmissive parts of the population network, we reduce the basic reproduction number significantly. This reduction in ![Graphic][143] can be evaluated by considering the transmissibility of the disease *T*. In short, considering a normalized recovery rate, *T* is the probability that an infected person infects one of their contacts per unit of time. Consider ![Graphic][144] to be the degree of connectivity of a person; therefore we can compute *T* for a heterogeneous network as [19, 22] ![Formula][145] Now imagine a scenario where we halt the virus-spread among all individuals with a degree of connectivity above ![Graphic][146]; for example via vaccinating every person who has contact numbers above ![Graphic][147]. Therefore, we get a reduction in the virus reproduction number ![Formula][148] where ![Graphic][149] denotes the basic reproduction number of the new network. For a specified size of the new network this would be the maximum reduction one could achieve by reshaping the network, as the tail of the contact distribution is removed (Fig. 4B). However, in practice there are failures in containing the virus-spread through highly connected people. The efficacy ![Graphic][150] of a smart testing based on contact counting thus depends on the number of highly connected people who would employ the contact counting app, as well as the accuracy of the tests. In the following we compute how these boundary conditions affect STeCC. However, before proceeding, notice that we suppose that the isolated individuals have negligible contributions to the virus-spread. Furthermore, we assume that elderly people (above 70 years of age) are shielded by isolation and that children below 10 years would not contribute to the infection dynamics [26, 9]. Therefore our target subpopulation is considered to be in possession of smart-phones. Let us define a testing regime, where we screen through app users with number of connections ![Graphic][151] at day 1, ![Graphic][152] at day 2 and so-forth until testing ![Graphic][153] at day *Tt.* Consequently, based on the test results, we ask the positively tested individuals to quarantine themselves. We fix the testing cycle to *Tt* = 7 days (see Fig. 4). Consider *ζ* to be the fraction of smart-phone owners who utilize the app. The portion of the network besides the fraction *ζ* that we can disconnect from the population depends on the probability of the event that a positively tested highly connected individual could pass on the virus to at least one person in the past *Tt + τproc* days. Let us denote such an event by ![Graphic][154] and consider the test processing time *τproc.* Therefore the reduction in ![Graphic][155] during the whole cycle of STeCC becomes ![Formula][156] where ![Graphic][157] is the efficacy of STeCC and can be computed via ![Formula][158] Note that *ŋ* is the fraction of false negative test results, ![Graphic][159] is the probability of the event ![Graphic][160], ![Graphic][161] is the probability that an individual has contacts below and ![Graphic][162]. It is evident that ![Graphic][163] and ![Graphic][164] would now depend on the network topology and latency time, respectively. #### 6.1. Network Topology In order to model the heterogeneity of a population relevant for disease modeling, scale-free networks offer appropriate features [5, 19, 16, 22]. Scale-free networks are characterized by a power-law distribution which determines the probability density function *Pk* of the degree of connectivity *k* per node. We employ a continuous approximation [23], that is, ![Formula][165] where ![Formula][166] Notice that by contacts we mean disease relevant contacts. We adopt *γ* = 0.3 and fix the upper and lower cut-offs by *kc =* 552 and *m =* 4, respectively. These choices were made in order to have a realistic range for number of contacts and obtain an average number of contacts of *µ* = 13.4, which is consistent with the data provided in [21]. Two extreme choices of *γ* include *γ* = 1 and *γ* = 0. While the former leads to the celebrated Barabási and Albert (BA) model [1]; the latter has been used in disease spread models, e.g. see [19]. Note that we conducted a sensitivity analysis of our results with respect to *γ*, see Figs. S7 and S8. With the choice of the degree of connectivity distribution we can compute the relative size of the population with a degree of connectivity above a certain ![Graphic][167]: ![Formula][168] Before proceeding further, let us mention that for the perfect efficacy ![Graphic][169], we can reach ![Graphic][170] with ![Graphic][171]. This accounts for 0.83% of the population with highest degree of connectivity. Now, to translate ![Graphic][172] of the population into the equivalent number of tests, suppose we conduct *Ntest* per day per 100’000 people. Thus we obtain ![Formula][173] #### 6.2. Uncontained Virus Spread Since we test these high-contact individuals once every *Tt* = 7 days (test cycle depicted in Fig. 4B), there is a chance that a virus-positive person has already transmitted the virus between two successive tests. The probability of such events depends on the latency time *τl* and average number of contacts of the individual. The infection time *τ* is uniformly distributed between 0 and *Tt*. Now, since these individuals have contacts way above average, we take a conservative estimation that ![Formula][174] which becomes 0.4296 for the case of *Tt* = 7 (day) and *τproc* = 1 (day). #### 6.3. Scenario Analysis Now we are ready to consider the following scenarios and compute the reduction of the basic reproduction number as a function of number of app users. 1. *Scneario A:* In this STeCC alone scenario we identify high-contact individuals in cycles of 7 days (Fig. 4B), ask them to be tested and ask the positively tested individuals to go into self-quarantine. We compute the reduced basic reproduction number ![Graphic][175] using Eq. (37) with the efficacy equation (38). For our base parameters together with the network topology with *γ* = 0.3, *m* = 4 and *kc* = 552, Fig. S 6 shows the performance of STeCC-A. While up to 30% reduction of ![Graphic][176] can be achieved with a combination of 90% app users among smart phone users (which would only be achieved in very optimistic scenarios) and 200 tests per 100’000 per day, it is clearly insufficient to halt the pandemic. The reason lies in the fact that we would not be able to contain the virus spread from almost 40% of highly connected people resulting from the transmission events that occur between their infection and the date of the virus test. Next, we introduce a combination that can significantly reduce this virus spread. 2. *Scenario B:* Here we consider a variant of STeCC in which besides the positively tested ones also their contacts are asked to quarantine. This results in a much lower basic reproduction number of ![Graphic][177], which is obtained from Eq. (37) with the improved efficacy ![Formula][178] The results shown in Fig. 4C predict a much stronger effect than STeCC alone (Fig. S 6). The combination of 90% user percentage and almost 400 tests per 100’000 per day now lead to ![Graphic][179]. Again, this approach alone will likely not suffice to halt the pandemic under realistic conditions. 3. *Scenario C:* Since the app already provides the contacts, we can further improve the STeCC based mitigation by combining it with conventional contact tracing (see §4). This is especially interesting since neither contact tracing nor STeCC alone are effective enough to halt the pandemic under realistic conditions. Therefore, on top of scenario B, we also ask the contacts of symptomatic cases to quarantine, which leads to a basic reproduction number of ![Formula][180] where ![Graphic][181] and ![Graphic][182] (see §4). The result corresponding to this scenario is shown in Fig. 4D. Accordingly, we predict that ![Graphic][183] can be achieved with 72% app users among smart phone users and 166 tests per 100’000 per day. This is very encouraging, since 72% app users among smart phone users corresponds to only about 50% app users of the whole population. Furthermore, a testing capacity of 166 per 100’000 per day already is available in several developed countries, including Switzerland. ![Fig. S 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F10.medium.gif) [Fig. S 6.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F10) Fig. S 6. STeCC alone scenario (STeCC-A): ![Graphic][184] as function of number of smart tests per 100’000 per day and the percentage of smartphone owners who participate in STeCC (Scenario A). For ![Graphic][185], the dotted line indicates ![Graphic][186], the dashed line ![Graphic][187] and the dashed-dotted line ![Graphic][188]. #### 6.4. Sensitivity Study In order to gain further confidence in our STeCC related mitigation scenarios, we conducted studies to investigate the sensitivity of the basic reproduction number ![Graphic][189]. The parameters which we varied are *y* ∊ {2, 2.5, 3.5} in the exponent of the power-law distribution of the degree of connectivity, the ratio *η* ∊ {0.1, 0.15} of false negatives and the test processing time *τproc* ∊ {0.5,1.5} (day). Note that our base setting is the combination of ![Graphic][190], *γ =* 0.3, *η =* 0.05 and *τproc =* 1 (day). Figures S 7 and 8 show the sensitivity of STeCC-B and-C scenarios, respectively, for varying ![Graphic][191] and *γ* values. Figure S9 depicts the sensitivity of STeCC-C with respect to the fraction of false negatives and test processing time. We observe that for a large range of parameters considered the combination of STeCC and conventional contact tracing leads to stopping the pandemic with realistic app user percentage (i.e., 60% to 85%) and number of tests per day (i.e., 50 to 350 tests per 100’000). ### 7. Model Implementation The dynamic model was implemented with Maple 2018. The calculations for mass testing, contact tracing and smart testing were implemented with MATLAB and the Statistics Toolbox Release 2018b. ![Fig. S 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F11.medium.gif) [Fig. S 7.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F11) Fig. S 7. STeCC plus isolation of contacts with positively tested individuals (STeCC-B): ![Graphic][192] as function of number of smart tests per 100ʹ000 per day and the percentage of smart-phone owners who participate in STeCC (scenario B). For (*γ* = 0.3, *kc* = 552) and (A) ![Graphic][193], (C) ![Graphic][194] and (E) ![Graphic][195]; for ![Graphic][196] and (B) (*γ* = 0.2, *kc* = 227), (D) (*γ* = 0.25, *kc* = 327) and (F) (*γ* = 0.35, *kc* = 1325). A corresponding map with ![Graphic][197] and *γ* = 0.3 is shown in Fig. 4C. The bold lines indicates the combinations for which ![Graphic][198]. ![Fig. S 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F12.medium.gif) [Fig. S 8.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F12) Fig. S 8. STeCC plus isolation of contacts with positively tested individuals plus classical contact tracing (STeCC-C): ![Graphic][199] as function of number of smart tests per 100ʹ000 per day and the percentage of smart-phone owners who participate in STeCC (scenario C). For (*γ* = 0.3, *kc =* 552) and (A) ![Graphic][200], (C) ![Graphic][201] and (E) ![Graphic][202]; for ![Graphic][203] and (B) (*γ* = 0.2, *kc =* 227), (D) (*γ* = 0.25, *kc =* 327) and (F) (*γ* = 0.35, *kc =* 1325). A corresponding map with ![Graphic][204] and *γ* = 0.3 is shown in Fig. 4D. The bold lines indicates the combinations for which ![Graphic][205]. ![Fig. S 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/07/2020.03.27.20045237/F13.medium.gif) [Fig. S 9.](http://medrxiv.org/content/early/2020/05/07/2020.03.27.20045237/F13) Fig. S 9. STeCC plus isolation of contacts with positively tested individuals plus classical contact tracing (STeCC-C): ![Graphic][206] as function of number of smart tests per 100’000 per day and the percentage of smart-phone owners who participate in STeCC (scenario C). For (*η* = 0.05) and (A) *τproc* = 0.5 (day) and (C) *τproc* = 1.5 (day); for *τproc* = 1 (day) and (B) *η* = 0.1, (D) *η* = 0.15. A corresponding map with *η* = 0.05 and *τproc* = 1 (day) is shown in Fig. 4D. The bold lines indicates the combinations for which ![Graphic][207]. ## Acknowledgements The authors are very thankful to Dario Ackermann, who created the website with the simulation tool based on the model presented in this paper. The authors would like to thank Emma Slack, Erik Bakkeren and Noemi Santamaria for helpful comments on the manuscript. HG acknowledges funding from the Swiss National Science foundation (grant number 174060). ## Footnotes * 1 [https://www.roche.com/media/releases/med-cor-2020-03-13.htm](https://www.roche.com/media/releases/med-cor-2020-03-13.htm) * Received March 27, 2020. * Revision received April 28, 2020. * Accepted May 7, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Johns Hopkins University, Worldometer (2020) (available at [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/)). 2. 2. R. Li et al., Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science. 6, eabb3221–9 (2020). 3. 3. R. Hinch et al., Effective Configurations of a Digital Contact Tracing App: A report to NHSX. A report to NHSX, 1–29 (2020). 4. 4. L. Ferretti et al., Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science, eabb6936–13 (2020). 5. 5. N. Ferguson et al., Impact of non-pharmaceutical interventions (NPIs) to reduce COVID19 mortality and healthcare demand. Imperial College COVID-19 Response Team (2020), doi:10.25561/77482. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.25561/77482&link_type=DOI) 6. 6. E. Lavezzo et al., Suppression of COVID-19 outbreak in the municipality of Vo, Italy. medRxiv, 1–23 (2020). 7. 7.Robert Koch Institut, Täglicher Lagebericht des RKI zur Coronavirus-Krankheit-2019 (COVID-19), 1–10 (2020). 8. 8. J. Scire et al., “Reproductive number of the COVID-19 epidemic in Switzerland with a focus on the Cantons of Basel-Stadt and Basel-Landschaft” (2020), pp. 1–13. 9. 9. D. Baud et al., Real estimates of mortality following COVID-19 infection. Lancet Infect Dis, 1–1 (2020). 10. 10. Y. Liu, A. A. Gayle, A. Wilder-Smith, J. Rocklöv, The reproductive number of COVID-19 is higher compared to SARS coronavirus. Journal of Travel Medicine. 27, 1–4 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taaa020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) 11. 11. Z. Liu, X. Bing, Z. Zhi, The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (covid-19) in China. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. 41 (2020). 12. 12. K. Men et al., Estimate the incubation period of coronavirus 2019 (COVID-19). medRxiv, 1–13 (2020). 13. 13. V. Virlogeux et al., Incubation Period Duration and Severity of Clinical Disease Following Severe Acute Respiratory Syndrome Coronavirus Infection. Epidemiology. 26, 666–669 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/EDE.0000000000000339&link_type=DOI) 14. 14. F. Zhou et al., Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 395, 1054–1062 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30566-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) 15. 15. C. Sheridan, Fast, portable tests come online to curb coronavirus pandemic. Nat Biotechnol (2020), doi:10.1038/d41587-020-00010-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/d41587-020-00010-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32203294&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) 16. 16.Pan-European Privacy-Preserving Proximity Tracing, HOME | Pepp-Pt (2020) (available at [www.pepp-pt.org](http://www.pepp-pt.org)). 17. 17. L. A. Meyers, B. Pourbohloul, M. E. J. Newman, D. M. Skowronski, R. C. Brunham, Network theory and SARS: predicting outbreak diversity. Journal of Theoretical Biology. 232, 71–81 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2004.07.026&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15498594&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225304900008&link_type=ISI) 18. 18. M. J. Keeling, K. T. D. Eames, Networks and epidemic models. J. R. Soc. Interface. 2, 295–307 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsif.2005.0051&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16849187&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234342000003&link_type=ISI) 19. 19. D. F. Gudbjartsson et al., Spread of SARS-CoV-2 in the Icelandic Population. N Engl J Med, NEJMoa2006100–14 (2020). 20. 20. Y. Xu et al., Characteristics of pediatric SARS-CoV-2 infection and potential evidence for persistent fecal viral shedding. Nature Medicine, 1–9 (2020). ## References 1. [1]. Albert-László Barabási and Réka Albert. Emergence of scaling in random networks. science, 286(5439):509–512, 1999. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIyODYvNTQzOS81MDkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8wNy8yMDIwLjAzLjI3LjIwMDQ1MjM3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 2. [2]. David Baud, Xiaolong Qi, Karin Nielsen-Saines, Didier Musso, Léo Pomar, and Guillaume Favre. Real estimates of mortality following covid-19 infection. The Lancet Infectious Diseases, 2020. 3. [3]. Victor M Corman, Olfert Landt, Marco Kaiser, Richard Molenkamp, Adam Meijer, Daniel KW Chu, Tobias Bleicker, Sebastian Brünink, Julia Schneider, Marie Luisa Schmidt, Daphne GJC Mulders, Bart L Haagmans, Bas van der Veer, Sharon van den Brink, Lisa Wijsman, Gabriel Goderski, Jean-Louis Romette, Joanna Ellis, Maria Zambon, Malik Peiris, Herman Goossens, Chantal Reusken, Marion PG Koopmans, and Christian Drosten. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Eurosurveillance, 25(3), January 2020. doi: 10.2807/1560-7917.es.2020.25.3.2000045. URL [https://doi.org/10.2807/1560-7917.es.2020.25.3.2000045](https://doi.org/10.2807/1560-7917.es.2020.25.3.2000045). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.es.2020.25.3.2000045&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31992387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) 4. [4]. Odo Diekmann, Johan Andre Peter Heesterbeek, and Johan AJ Metz. On the definition and the computation of the basic reproduction ratio *r* in models for infectious diseases in heterogeneous populations. Journal of mathematical biology, 28(4):365–382, 1990. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00178324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2117040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DG35000001&link_type=ISI) 5. [5]. Victor M Eguiluz and Konstantin Klemm. Epidemic threshold in structured scale-free networks. Physical Review Letters, 89(10):108701, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevLett.89.108701&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12225235&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) 6. [6]. Neil M Ferguson, Daniel Laydon, Gemma Nedjati-Gilani, Natsuko Imai, Kylie Ainslie, Marc Baguelin, Sangeeta Bhatia, Adhiratha Boonyasiri, Zulma Cucunubá, Gina Cuomo-Dannenburg, et al. Impact of non-pharmaceutical interventions (npis) to reduce covid-19 mortality and healthcare demand. London: Imperial College COVID-19 Response Team, March, 16, 2020. 7. [7]. Luca Ferretti, Chris Wymant, Michelle Kendall, Lele Zhao, Anel Nurtay, Lucie Abeler-Dörner, Michael Parker, David Bonsall, and Christophe Fraser. Quantifying sarscov-2 transmission suggests epidemic control with digital contact tracing. Science, 2020. 8. [8]. Christophe Fraser, Steven Riley, Roy M Anderson, and Neil M Ferguson. Factors that make an infectious disease outbreak controllable. Proceedings of the National Academy of Sciences, 101(16):6146–6151, 2004. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAxLzE2LzYxNDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8wNy8yMDIwLjAzLjI3LjIwMDQ1MjM3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 9. [9]. Daniel F Gudbjartsson, Agnar Helgason, Hakon Jonsson, Olafur T Magnusson, Pall Melsted, Gudmundur L Norddahl, Jona Saemundsdottir, Asgeir Sigurdsson, Patrick Sulem, Arna B Agustsdottir, et al. Spread of sars-cov-2 in the icelandic population. New England Journal of Medicine, 2020. 10. [10]. Robert Hinch, Will Probert, Anel Nurtay, Michelle Kendall, Chris Wymant, Matthew Hall, Katrina Lythgoe, Ana Bulas Cruz, Lele Zhao, Andrea Stewart, Luca Ferretti, Michael Parker, Ares Meroueh, Bryn Mathias, Scott Stevenson, Daniel Montero, James Warren, Nicole K Mather, Anthony Finkelstein, Lucie, Abeler-Dörner, David Bonsall, and Christophe Fraser. Effective configurations of a digital contact tracing app: A report to NHSX, April 2020. URL [https://045.medsci.ox.ac.uk/files/files/report-effective-app-configurations.pdf](https://045.medsci.ox.ac.uk/files/files/report-effective-app-configurations.pdf). 11. [11]. Ayaan Hossain, C. Alexander Reis, Sarthok Rahman, and M. Howard Salis. A massively parallel covid-19 diagnostic assay for simultaneous testing of 19200 patient samples, 2020. 12. [12]. Enrico Lavezzo, Elisa Franchin, Constanze Ciavarella, Gina Cuomo-Dannenburg, Luisa Barzon, Claudia Del Vecchio, Lucia Rossi, Riccardo Manganelli, Arianna Loregian, Nicolò Navarin, et al. Suppression of covid-19 outbreak in the municipality of vo, italy. medRxiv, 2020. 13. [13]. Ruiyun Li, Sen Pei, Bin Chen, Yimeng Song, Tao Zhang, Wan Yang, and Jeffrey Shaman. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (sars-cov2). Science, 2020. 14. [14]. Ying Liu, Albert A Gayle, Annelies Wilder-Smith, and Joacim Rocklöv. The reproductive number of covid-19 is higher compared to sars coronavirus. Journal of travel medicine, 2020. 15. [15]. Zhonghua Liu, Xing Bing, and Xue Za Zhi. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (covid-19) in china. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team, 41(2):145, 2020. 16. [16]. Robert M May and Alun L Lloyd. Infection dynamics on scale-free networks. Physical Review E, 64(6):066112, 2001. 17. [17]. Kenneth McIntosh, Martin S Hirsch, and Allyson Bloom. Coronavirus disease 2019 (covid-19). UpToDate, 5, 2020. 18. [18]. Ke Men, Xia Wang, Yihao Li, Guangwei Zhang, Jingjing Hu, Yanyan Gao, and Henry Han. Estimate the incubation period of coronavirus 2019 (covid-19). medRxiv, 2020. 19. [19]. Lauren Ancel Meyers, Babak Pourbohloul, Mark EJ Newman, Danuta M Skowronski, and Robert C Brunham. Network theory and sars: predicting outbreak diversity. Journal of theoretical biology, 232(1):71–81, 2005. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2004.07.026&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15498594&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F07%2F2020.03.27.20045237.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225304900008&link_type=ISI) 20. [20]. Kenji Mizumoto and Gerardo Chowell. Estimating risk for death from 2019 novel coronavirus disease, china, january-february 2020. Emerging infectious diseases, 26(6), 2020. 21. [21]. Joël Mossong, Niel Hens, Mark Jit, Philippe Beutels, Kari Auranen, Rafael Mikolajczyk, Marco Massari, Stefania Salmaso, Gianpaolo Scalia Tomba, Jacco Wallinga, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS medicine, 5(3), 2008. 22. [22]. Mark EJ Newman. Spread of epidemic disease on networks. Physical review E, 66(1):016128, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevE.66.016128&link_type=DOI) 23. [23].Romualdo Pastor-Satorras and Alessandro Vespignani. Epidemic dynamics in finite size scale-free networks. Physical Review E, 65(3):035108, 2002. 24. [24]. Robert Verity, Lucy C Okell, Ilaria Dorigatti, Peter Winskill, Charles Whittaker, Natsuko Imai, Gina Cuomo-Dannenburg, Hayley Thompson, Patrick GT Walker, Han Fu, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet Infectious Diseases, 2020. 25. [25]. Victor Virlogeux, Vicky J Fang, Joseph T Wu, Lai-Ming Ho, JS Malik Peiris, Gabriel M Leung, and Benjamin J Cowling. Incubation period duration and severity of clinical disease following severe acute respiratory syndrome coronavirus infection. Epidemiology (Cambridge, Mass.), 26(5):666, 2015. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/EDE.0000000000000339&link_type=DOI) 26. [26]. Yi Xu, Xufang Li, Bing Zhu, Huiying Liang, Chunxiao Fang, Yu Gong, Qiaozhi Guo, Xin Sun, Danyang Zhao, Jun Shen, et al. Characteristics of pediatric sars-cov-2 infection and potential evidence for persistent fecal viral shedding. Nature medicine, pages 1–4, 2020. 27. [27]. Fei Zhou, Ting Yu, Ronghui Du, Guohui Fan, Ying Liu, Zhibo Liu, Jie Xiang, Yeming Wang, Bin Song, Xiaoying Gu, et al. Clinical course and risk factors for mortality of adult inpatients with covid-19 in wuhan, china: a retrospective cohort study. The Lancet, 2020. [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/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/graphic-11.gif [11]: /embed/graphic-12.gif [12]: /embed/graphic-13.gif [13]: /embed/graphic-14.gif [14]: /embed/inline-graphic-5.gif [15]: /embed/inline-graphic-6.gif [16]: /embed/inline-graphic-7.gif [17]: /embed/inline-graphic-8.gif [18]: /embed/inline-graphic-9.gif [19]: /embed/inline-graphic-10.gif [20]: /embed/inline-graphic-11.gif [21]: /embed/inline-graphic-12.gif [22]: /embed/inline-graphic-13.gif [23]: /embed/inline-graphic-14.gif [24]: /embed/inline-graphic-15.gif [25]: /embed/graphic-15.gif [26]: /embed/graphic-16.gif [27]: /embed/graphic-17.gif [28]: /embed/graphic-18.gif [29]: /embed/graphic-19.gif [30]: /embed/graphic-20.gif [31]: /embed/graphic-21.gif [32]: /embed/graphic-22.gif [33]: /embed/inline-graphic-16.gif [34]: /embed/inline-graphic-17.gif [35]: /embed/inline-graphic-18.gif [36]: /embed/inline-graphic-19.gif [37]: /embed/inline-graphic-20.gif [38]: /embed/inline-graphic-21.gif [39]: /embed/graphic-23.gif [40]: /embed/inline-graphic-22.gif [41]: /embed/graphic-25.gif [42]: /embed/inline-graphic-23.gif [43]: /embed/inline-graphic-24.gif [44]: /embed/inline-graphic-25.gif [45]: /embed/inline-graphic-26.gif [46]: /embed/inline-graphic-27.gif [47]: /embed/inline-graphic-28.gif [48]: /embed/inline-graphic-29.gif [49]: /embed/inline-graphic-30.gif [50]: /embed/inline-graphic-31.gif [51]: /embed/graphic-26.gif [52]: /embed/graphic-27.gif [53]: /embed/inline-graphic-32.gif [54]: /embed/inline-graphic-33.gif [55]: /embed/inline-graphic-34.gif [56]: /embed/inline-graphic-35.gif [57]: /embed/inline-graphic-36.gif [58]: /embed/inline-graphic-37.gif [59]: /embed/inline-graphic-38.gif [60]: F6/embed/inline-graphic-39.gif [61]: F6/embed/inline-graphic-40.gif [62]: /embed/inline-graphic-41.gif [63]: /embed/inline-graphic-42.gif [64]: /embed/inline-graphic-43.gif [65]: /embed/inline-graphic-44.gif [66]: /embed/inline-graphic-45.gif [67]: /embed/inline-graphic-46.gif [68]: /embed/inline-graphic-47.gif [69]: /embed/inline-graphic-48.gif [70]: /embed/inline-graphic-49.gif [71]: /embed/inline-graphic-50.gif [72]: /embed/graphic-31.gif [73]: /embed/inline-graphic-51.gif [74]: /embed/inline-graphic-52.gif [75]: /embed/inline-graphic-53.gif [76]: /embed/graphic-32.gif [77]: /embed/inline-graphic-54.gif [78]: /embed/inline-graphic-55.gif [79]: /embed/graphic-33.gif [80]: /embed/inline-graphic-56.gif [81]: /embed/inline-graphic-57.gif [82]: /embed/inline-graphic-58.gif [83]: /embed/inline-graphic-59.gif [84]: /embed/inline-graphic-60.gif [85]: /embed/graphic-34.gif [86]: /embed/graphic-35.gif [87]: /embed/graphic-36.gif [88]: /embed/inline-graphic-61.gif [89]: /embed/inline-graphic-62.gif [90]: /embed/inline-graphic-63.gif [91]: /embed/inline-graphic-64.gif [92]: /embed/inline-graphic-65.gif [93]: F7/embed/inline-graphic-66.gif [94]: F7/embed/inline-graphic-67.gif [95]: F7/embed/inline-graphic-68.gif [96]: /embed/inline-graphic-69.gif [97]: /embed/inline-graphic-70.gif [98]: /embed/inline-graphic-71.gif [99]: /embed/inline-graphic-72.gif [100]: F8/embed/inline-graphic-73.gif [101]: F8/embed/inline-graphic-74.gif [102]: /embed/inline-graphic-75.gif [103]: /embed/graphic-39.gif [104]: /embed/inline-graphic-76.gif [105]: /embed/graphic-40.gif [106]: /embed/inline-graphic-77.gif [107]: /embed/inline-graphic-78.gif [108]: /embed/inline-graphic-79.gif [109]: /embed/graphic-41.gif [110]: /embed/inline-graphic-80.gif [111]: /embed/inline-graphic-81.gif [112]: /embed/inline-graphic-82.gif [113]: /embed/graphic-42.gif [114]: /embed/inline-graphic-83.gif [115]: /embed/graphic-43.gif [116]: /embed/inline-graphic-84.gif [117]: /embed/inline-graphic-85.gif [118]: /embed/inline-graphic-86.gif [119]: /embed/inline-graphic-87.gif [120]: F9/embed/inline-graphic-88.gif [121]: F9/embed/inline-graphic-89.gif [122]: F9/embed/inline-graphic-90.gif [123]: /embed/inline-graphic-91.gif [124]: /embed/inline-graphic-92.gif [125]: /embed/inline-graphic-93.gif [126]: /embed/inline-graphic-94.gif [127]: /embed/inline-graphic-95.gif [128]: /embed/inline-graphic-96.gif [129]: /embed/inline-graphic-97.gif [130]: /embed/inline-graphic-98.gif [131]: /embed/inline-graphic-99.gif [132]: /embed/inline-graphic-100.gif [133]: /embed/inline-graphic-101.gif [134]: /embed/inline-graphic-102.gif [135]: /embed/inline-graphic-103.gif [136]: /embed/inline-graphic-104.gif [137]: /embed/inline-graphic-105.gif [138]: /embed/graphic-45.gif [139]: /embed/graphic-46.gif [140]: /embed/inline-graphic-106.gif [141]: /embed/inline-graphic-107.gif [142]: /embed/inline-graphic-108.gif [143]: /embed/inline-graphic-109.gif [144]: /embed/inline-graphic-110.gif [145]: /embed/graphic-47.gif [146]: /embed/inline-graphic-111.gif [147]: /embed/inline-graphic-112.gif [148]: /embed/graphic-48.gif [149]: /embed/inline-graphic-113.gif [150]: /embed/inline-graphic-114.gif [151]: /embed/inline-graphic-115.gif [152]: /embed/inline-graphic-116.gif [153]: /embed/inline-graphic-117.gif [154]: /embed/inline-graphic-118.gif [155]: /embed/inline-graphic-119.gif [156]: /embed/graphic-49.gif [157]: /embed/inline-graphic-120.gif [158]: /embed/graphic-50.gif [159]: /embed/inline-graphic-121.gif [160]: /embed/inline-graphic-122.gif [161]: /embed/inline-graphic-123.gif [162]: /embed/inline-graphic-124.gif [163]: /embed/inline-graphic-125.gif [164]: /embed/inline-graphic-126.gif [165]: /embed/graphic-51.gif [166]: /embed/graphic-52.gif [167]: /embed/inline-graphic-127.gif [168]: /embed/graphic-53.gif [169]: /embed/inline-graphic-128.gif [170]: /embed/inline-graphic-129.gif [171]: /embed/inline-graphic-130.gif [172]: /embed/inline-graphic-131.gif [173]: /embed/graphic-54.gif [174]: /embed/graphic-55.gif [175]: /embed/inline-graphic-132.gif [176]: /embed/inline-graphic-133.gif [177]: /embed/inline-graphic-134.gif [178]: /embed/graphic-56.gif [179]: /embed/inline-graphic-135.gif [180]: /embed/graphic-57.gif [181]: /embed/inline-graphic-136.gif [182]: /embed/inline-graphic-137.gif [183]: /embed/inline-graphic-138.gif [184]: F10/embed/inline-graphic-139.gif [185]: F10/embed/inline-graphic-140.gif [186]: F10/embed/inline-graphic-141.gif [187]: F10/embed/inline-graphic-142.gif [188]: F10/embed/inline-graphic-143.gif [189]: /embed/inline-graphic-144.gif [190]: /embed/inline-graphic-145.gif [191]: /embed/inline-graphic-146.gif [192]: F11/embed/inline-graphic-147.gif [193]: F11/embed/inline-graphic-148.gif [194]: F11/embed/inline-graphic-149.gif [195]: F11/embed/inline-graphic-150.gif [196]: F11/embed/inline-graphic-151.gif [197]: F11/embed/inline-graphic-152.gif [198]: F11/embed/inline-graphic-153.gif [199]: F12/embed/inline-graphic-154.gif [200]: F12/embed/inline-graphic-155.gif [201]: F12/embed/inline-graphic-156.gif [202]: F12/embed/inline-graphic-157.gif [203]: F12/embed/inline-graphic-158.gif [204]: F12/embed/inline-graphic-159.gif [205]: F12/embed/inline-graphic-160.gif [206]: F13/embed/inline-graphic-161.gif [207]: F13/embed/inline-graphic-162.gif