key: cord-308069-iydjrmhh
authors: Contreras, Sebastián; Biron-Lattes, Juan Pablo; Villavicencio, H. Andrés; Medina-Ortiz, David; Llanovarced-Kawles, Nyna; Olivera-Nappa, Álvaro
title: Statistically-based methodology for revealing real contagion trends and correcting delay-induced errors in the assessment of COVID-19 pandemic
date: 2020-07-03
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.110087
sha: 
doc_id: 308069
cord_uid: iydjrmhh

COVID-19 pandemic has reshaped our world in a timescale much shorter than what we can understand. Particularities of SARS-CoV-2, such as its persistence in surfaces and the lack of a curative treatment or vaccine against COVID-19, have pushed authorities to apply restrictive policies to control its spreading. As data drove most of the decisions made in this global contingency, their quality is a critical variable for decision-making actors, and therefore should be carefully curated. In this work, we analyze the sources of error in typically reported epidemiological variables and usual tests used for diagnosis, and their impact on our understanding of COVID-19 spreading dynamics. We address the existence of different delays in the report of new cases, induced by the incubation time of the virus and testing-diagnosis time gaps, and other error sources related to the sensitivity/specificity of the tests used to diagnose COVID-19. Using a statistically-based algorithm, we perform a temporal reclassification of cases to avoid delay-induced errors, building up new epidemiologic curves centered in the day where the contagion effectively occurred. We also statistically enhance the robustness behind the discharge/recovery clinical criteria in the absence of a direct test, which is typically the case of non-first world countries, where the limited testing capabilities are fully dedicated to the evaluation of new cases. Finally, we applied our methodology to assess the evolution of the pandemic in Chile through the Effective Reproduction Number R(t), identifying different moments in which data was misleading governmental actions. In doing so, we aim to raise public awareness of the need for proper data reporting and processing protocols for epidemiological modelling and predictions.

Since the outbreak of novel SARS-CoV-2 in late 2019, the spread of COVID-19 has changed nearly every aspect of our daily life, challenging modern society to find a way to function under conditions never seen before. Governmental plans on public health have played a crucial role in its control in the absence of an effective treatment to cure COVID-19 or a vaccine to prevent it. Typically reported variables in the COVID-19 pandemic, available in public repositories (as Worldometers.info, 02 Jun, 2020; Dong et al., 2020, among others) , are the incidence of new cases ∆T , total cases T , discharged/recovered cases R, and deaths D. These variables serve as input for the development and evaluation of governmental plans and to fit the vast variety of SIR-like mathematical models recently proposed (see, e.g., Contreras et al., 2020a; Yang et al., 2020; Pang, 2020; Aguilar et al., 2020 , and references therein for a brief review of them). Among the several factors conditioning the quality/reliability of the variables mentioned above are those related to the sensitivity and specificity of diagnostic tests, delays between sampling and diagnosing, and the delay between contagion, development of symptoms and testing. The latter varies from country to country, depending heavily on the local government's testing strategy and resources. Even though SIR-like models have been reported to be not-suitable to predict turning points or any other quantitative insight on this pandemic (Castro et al., 2020) , to the best knowledge of the authors, delays have not been pointed out as the cause, as raw data has always been used to fit them. Different parameters can be used to evaluate the evolution of the SARS-CoV-2 outbreak. Among them we may find the documentation rate (Wilder et al., 2020) , secondary infection rate, serological response to infection (WHO, 2020b), number of vacancies at ICU (Liew et al., 2020) , and the Basic Reproduction Number R 0 (Allieta et al., 2020; Gevertz et al., 2020) , which is one of the most widely used. This parameter (R 0 ) represents the number of persons a single infected individual might infect before either recovering or dying, when the population is entirely susceptible (Perasso, 2018) . Traditional forms to estimate the R 0 are rather complex, heavily depending on the fitting of SIR models to local data (Heesterbeek, 2002; Delamater et al., 2019; Wang et al., 2020; Ma, 2020) . In a previous work , we proposed a methodology to obtain real-time estimations of the Effective Reproduction Number R t directly from raw data, which was satisfactorily applied to evaluate the panorama of the COVID-19 spread in different countries and to forecast its evolution (Medina-Ortiz et al., 2020a) . Nevertheless, its heavy dependence on reported data required the study of common error sources affecting this parameter, and the development of methodologies to control, correct, and quantify their impact .

In this work, we analyze the sources of error in the typically reported epidemiological variables and their impact on our understanding of COVID-19 spreading dynamics. We address the existence of different delays in the report of new cases, induced by the incubation time of the virus and testing-diagnosis time gaps, and provide a straightforward methodology to avoid the propagation of delay-induced errors to model-derived parameters. Using our statistically-based algorithm, we perform a temporal reclassification of individuals to the day where they were -statistically-most likely to have acquired the virus, building a new smooth curve with corrected variables. We postulate this new temporally modified curve as the proper one for fitting purposes in SIR-like models. We present an analogous methodology to estimate the number of discharged/recovered individuals, based on the reported evolution of the viral infection, the performance of the different tests for its diagnosis, and the case fatality, which can be easily adapted for a particular country.

We used our methodology to assess the evolution of the pandemic in Chile, identifying different moments in which the use of raw data was misleading governmental actions.

Different methods for diagnosing COVID-19 have been developed and reported in the literature, with real-time RT-PCR being the standard applied globally (Lippi et al., 2020) . Nevertheless, techniques such as the IgG and IgM rapid tests, the Chest Computed Tomography (Chest CT), and CRISPR-Cas systems are also being used. In this section we provide a brief analysis of them, highlighting the different characteristics of both the techniques and their basic approach to detect viral infection.

Real-time reverse transcription polymerase chain reaction (RT-PCR) is a mechanism for amplification and detection of RNA in real time (see Gibson et al., 1996 , for an exhaustive description of the technique).

Initially, RNA obtained from samples is retrotranscribed to DNA using a reverse transcriptase enzyme. By applying temperature cycles, the conditions are created for new copies of the DNA to be synthesized from the initial one. The lower the initial DNA concentration, the lower the probability of a synthesis reaction in a given cycle. It is assumed that the minimum time from contagion until testing positive in the RT-PCR test is 2.3 days (95 % CI, 0.8-3.0 days) before the onset of symptoms (He et al., 2020) , which typically appear 5.2 days (95 % CI, 4.1-7.0 days) after contagion . He et al. (2020) ; Ganyani et al. (2020) showed that 44 % to 62 % of the total contagions occur in the pre-symptomatic period. It has been inferred that the viral load reaches a peak value before 0.7 days from the onset of symptoms (95 % CI, 0.2-2.0 days), when it starts falling monotonically together with the infectivity rate (He et al., 2020) . Finally, the virus has been detected for a median of 20 days after the onset of symptoms (Zhou et al., 2020) , but infectivity may decrease significantly eight days after symptom appearance.

A high false-negative rate (Li et al., 2020c; Xiao et al., 2020) and a sensitivity of 71 to 83 % (Fang et al., 2020; Long et al., 2020) have been reported for the real-time RT-PCR technique, and several vulnerabilities of it have been identified and quantified (Lippi et al., 2020) . Considering sampling, handling, testing, and reporting, the total time necessary to get RT-PCR results for an individual may range between 2 to 3 days 4 (Nguyen et al., 2020) . However, the time it takes to perform the RT-PCR experiment takes about 2 to 3 hours (Li et al., 2020d) .

Part of the immune system response to the SARS-CoV-2 infection is the production of specific antibodies against it, including IgG and IgM (Li et al., 2020b) . Serological tests detect the presence of those antibodies and, unlike the other detection methods, take only 15 minutes to produce results (Li et al., 2020d) . These tests have a sensitivity of 88.6 % and a specificity of 90.6 % (Li et al., 2020d) . This technology was developed for the SARS-CoV epidemic, which was caused by a virus belonging to the same family of coronaviruses as SARS-CoV-2, providing satisfactory results after 2-3 days from the onset of symptoms (for IgG), and after eight days (for IgM) (Woo et al., 2004) .

The principle behind the Chest Computed Tomography (Chest CT) is the analysis of cross-sectional lung images to identify viral pneumonia features, like ground-glass opacity, consolidation, reticulation/thickened interlobular septa or nodules (Ai et al., 2020) . Chest CT has shown a sensitivity between 97 % and 98 % (Long et al., 2020; Fang et al., 2020) . However, due to the similarities between CT images accounting for COVID-19 and CT images for other viral types of pneumonia, false-positives are likely to occur. Compared to RT-PCR, Chest CT tends to be more reliable, practical, and quick to diagnose COVID-19 (Ai et al., 2020) . Nevertheless, requiring the presence of the potentially infected patient in a health center lacks the flexibility that rapid tests provide, and can backfire on movement restriction measures. COVID-19 pneumonia manifests with abnormalities on computed tomography images of the chest, even in asymptomatic patients (Shi et al., 2020) .

In CRISPR-Cas systems, a guide RNA (gRNA) is designed to recognize a specific RNA sequence, like any particular gene or partial ARN sequence of SARS-CoV-2 coronavirus. Endonuclease enzymes of the Cas family and the specific gRNA will search for the sequence match. This match will deliver a signal that confirms the presence of SARS-CoV-2 RNA in the sample (LeMieux, 2020). DETECTR and SHERLOCK

are two examples of CRISPR-Cas technologies for the detection of SARS-CoV-2, being able to obtain results in less than 1 hour at a significantly lower cost compared to the RT-PCR technique. DETECTR showed a 95% positive predictive agreement and 100% negative predictive agreement (Broughton et al., 2020) , while SHERLOCK has not been validated using real patient samples and is not suitable for clinical use at this time (Jin et al., 2020) .

Our work aims to expose and quantify, both theoretically and in a case study, the impact of different sources of error in commonly reported data of the COVID-19 spread, such as newly reported cases ∆T , total cases T , infected I, and recovered R fractions of the population.

First, we define random variables associated with the delay in both sampling and diagnosing new cases, and by modeling their probability distribution functions, we derive a method to re-classify the newly reported cases accordingly. As reclassification occurs backwardly, for re-evaluating the current scenario through our methodology, we cast predictions on the reported new cases using an ARIMA (Autoregressive Integrated Moving Average) model. Having the corrected variables, we evaluate differences on the effective reproduction number R t , following the methodology presented by Contreras et al. (2020b) . Assuming the outbreak is well represented with a SIR model in a given timeframe i, R t can be estimated as βi γi , which can be directly estimated from the discrete version of the SIR differential equations. Decoupling the "recovered" fraction in the SIR model as clinically recovered individuals R and deaths D, R t would be given by equation 1:

Given the existence of a latent and incubation period for COVID-19, which will be carefully described in the next section, SIR models might not capture all the features of the outbreak and SEIR models would be more suitable. In such models, only infected individuals can propagate the infection and, equation-wise, estimations of R t remain more or less the same (Chowell et al., 2007) . As data of exposed individuals is not widely available or entirely reliable, its contribution was not included in the estimations of R t to keep its straightforward use.

Autoregression models for the forecast of ∆T were implemented using the statsmodels Python library (Seabold and Perktold, 2010) . All other calculations and visualizations were made using DMAKit-lib Python library (Medina-Ortiz et al., 2020b) and MATLAB R2018a.

There exists an incubation period t c for the development of SARS-CoV-2-related symptoms, which average has been reported to range between 5.1 days (Li et al., 2020a; He et al., 2020 ) and 6.4 days (Backer et al., 2020) . He et al. (2020) also reported the existence of a latency period, which is the time required before an infected individual could spread the infection, which was extended until 2.3 days (95% CI, 0.8-3.0 days) before the onset of symptoms, which occurred t c days after contagion.

The incubation time t c is especially relevant in the case of a symptoms-based testing strategy or when the spread has reached the non-traceability stage. To model it, we use τ 1 with log-normal distribution, as 6 reported in and suggested in Nishiura (2007):

According to the values reported in , µ = 1.63 = ln(t c ) and σ = 0.41 retrieve the expected values of 5.1 for the median, and 11.5 for the 98'th percentile. Even though the required time for performing the test is short (Li et al., 2020d) , delays between testing and diagnosis have been reported (Nguyen et al., 2020) . We will sum up secondary delays, such as the symptom-to-testing and testing-to-diagnosis time gaps, into a random variable τ 2 , which, for the sake of simplicity, will be assumed to follow a uniform distribution between t min and t max :

Further knowledge on the nature of factors affecting τ 2 could lead to a different distribution. We may postulate a reclassification for obtaining the real number of new contagions occurred in a day t as the sum of contributions of cases reported with a delay of k days:

where w i represents the fraction of patients that were notified at time t = t + k but had acquired the virus at t. Note that the different delays are referred to the random variable Z = τ 1 + τ 2 . The probability distribution function for Z is obtained by the convolution method, assuming τ 1 and τ 2 are independent and combining equations 2 and 3:

Assuming that data is reported on a daily basis, we can calculate the probability associated to having a delay of k days:

For practical reasons, we can define a threshold U < 1 for truncating the probability mass distribution (equation 6), which otherwise would assign a probability to every k ∈ N. Let n 1 be the first positive integer for which equation 7 holds,

we may rewrite equation 4 as:

The magnitude of the total delay between infection and diagnosis can be estimated through the expected value of equation 5 (or equivalently, equation 6). A schematic representation of the proposed methodology is presented in Figure 1 . Assuming the lowest reported value for the average incubation time, t c = 5.2, and a conservative timeframe for the delay between the appearance of symptoms, testing and diagnosing, (2-5 days), the expected delay E (Z) is about 7.5 days.

Algorithm 1: Statistically-based temporal reclassification algorithm for correcting delay-induced errors in the report of new COVID-19 infections Result: Corrected differential and cumulative total cases T .

dT : Officially reported new cases per day. dT ∈ R n ; λ, a, b : Parameters of the τ 1 and τ 2 distributions;

U : Probability threshold for the maximum possible delay; Calculate [P k, F ], respectively the probability and cumulative distribution functions for Z = τ 1 + τ 2 ; idx = F < U ; logical indexes of the values of vector F smaller than U ; P k = P k(idx); P k ∈ R m ; Define K, delay indexes of P k in idx. Renormalize P k;

Forecast the next m values for dT , using an autoregression (ARIMA) model;

Define an empty vector dT corr ∈ R 2m+n−1 to record the corrections; Round dT corr values to the nearest integer; dT corr = round(dT corr);

As discussed previously, errors in the amount of discharged/recovered patients R are likely to be greater only when no quantitative criteria are applied. In such cases, some countries (like Chile) have adopted the following criteria (officially reported in the MINSAL (14 Apr, 2020) report), possibly based on the recommendations published by the WHO (2020a). • If there were no previously existing pathologies, a patient should be discharged 14 days after testing positive for COVID-19.

• If there were previously existing pathologies, the patient should be discharged 28 days after testing positive for COVID-19.

This criterion turns out to be quite simplistic, especially considering the existence of uncertainties regarding the diagnosis and contagion days. If we try to model the probability of recovering from COVID-19, some assumptions are necessary. Let τ r be the random variable for the time of discharge/recovery such that:

• P (τ r ≤ 14) = 0

• P (τ r ≤ 28) ≈ 1

Further assumptions are necessary to estimate the probability distribution function f τr , as it depends on local diagnosis criteria, testing strategy, and the fraction of the population having preexisting pathologies. In particular, depending on the test applied for diagnosis and its sensitivity/specificity -which were carefully described in Section 2-the probability profile would change. The simplest form that can be assumed, and which, for clarity reasons, is adopted herein, is a triangular distribution:

Further knowledge on the nature of factors affecting τ r could lead to a different distribution.

Algorithm 2: Statistical estimation of dR, daily patients that have been discharged/recovered from

Result: Estimated discharged/recovered cases per day R.

dT corr : corrected new cases per day. dT corr ∈ R n ; dD : Officially reported daily deaths due to COVID-19. dD ∈ R n ; θ : Parameters of the probability distribution function for the time of discharge/recovery τ r .;

t min , t max : time frame for discharging an infected patient.; w = t max −t min +1, width of the time-frame of discharge; K = linspace(t min , t max , w); Calculate the probability distribution function for τ r ; P kR = f ( θ, w);

Define an empty vector dRcorr ∈ R w+n to record the corrections;

Subtract the daily deaths due to COVID-19 to obtain the final estimator, round to the nearest integer.

Confirm the result is positive or zero.; dRcorr = max(round(dRcorr − dD), 0) ;

The spread of COVID-19 in Chile is far from being controlled, as shown by the exponential growth that new cases have had in recent weeks (Worldometers.info, 02 Jun, 2020) . In order to apply our methodology, we need to cast predictions on the trends of ∆T . Figure 2 presents the current and forecast trends, using an autoregression ARIMA model. These models were chosen over others because of the great success they have had on representing COVID-19 dynamics in other relevant studies (Chakraborty and Ghosh, 2020; Benvenuto et al., 2020) . (02 Jun, 2020), while the red curve was generated using an autoregression ARIMA model. As these reported data are likely to be affected by several exogenous factors (Ribeiro et al., 2020) , the best performance metric for the generated forecast is trend consistency.

First, we performed a temporal reclassification of new cases to obtain ∆T corr , presented in Figure 3 . It can be seen that our methodology, besides exhibiting an horizontal semi-displacement, generates a smooth curve, which is due to the inclusion of a probability function that spreads the influence of daily informed cases in previous days (where they were most likely to have occurred). The last part of the red curve is dashed because it partially contains contributions of the forecast of ∆T , and therefore might change in the upcoming days, when the required data for completing the reclassification of cases would be available.

After obtaining statistically-based corrections of both ∆T and ∆R by following algorithms 1 and 2, and knowing the daily deaths due to COVID-19 ∆D, we proceed to calculate the variation on the active cases.

Based on the closure of the population-balance of the different considered classes, the cumulative new cases T will be the sum of active I, recovered/discharged R, and dead D, therefore their differences would follow: ∆T = ∆I + ∆R + ∆D ⇐⇒ ∆I = ∆T − (∆R + ∆D).

(10)

In the context of SEIR models, it is not possible to explicitly include exposed individuals E, as they are not being recorded separately from COVID-19 positive new cases.

After obtaining ∆I using equation 10, we proceed to calculate R t , using equation 1 with raw data, mobile averages of raw data, and the methodology proposed herein. As shown in Figure 4 , an abrupt growth in R t was evidenced around April 22nd, consistently with the relaxation of restrictive measures that were applied in Santiago, the capital and most populated city in the country, and the apogee of the governmental plan for a "safe return to work". Even though different trends seem to stabilise again in the second week of May, a lowering of the trend is not totally clear.

The different iconic dates highlighted in Figure 4 were obtained from the chronology presented in Wikipedia.org (23 May, 2020) and references therein. From this plot, we can assess both successful and unfortunate effects that the different governmental actions have had on the spreading dynamics of COVID-19

in Chile, and specifically in Santiago, by analyzing corrected R t trends and comparing them with daily reported cases. The various public-health actions were strongly based on locally-observed values from (raw) daily reported cases, yet appear to have been taken too late, or too early, according to the statisticallycorrected trend. In particular, the apogee of the governmental plan for a safe return to work occurred in a temporal window where raw-data-driven R t values were at a minimum, but the corrected contagion rates anticipated a steep growing trend for the following days. can be smoothed through mobile averages (dark red curve), but the trends are the same. A significantly different scenario is shown by the statistically-corrected Rt trend (green curve). Highlighted dates associated with iconic governmental actions in Chile: March 31st (second week of sectorised quarantine for the high-income districts of Santiago, capital of Chile), April 7th (compulsory use of facemask in public transport), April 24th (governmental call for a "safe return to work"), April 30th

(sectorised quarantine -low-income districts of Santiago-), May 15th (total quarantine in Santiago).

We have presented an exhaustive assessment of error sources in reported data of the COVID-19 pandemic and provided a methodology to minimize -and correct-their effect on both reported variables and modelderived parameters by applying a statistically-driven reclassification of newly reported cases, and corrections to discharge/recovery criteria. By using the corrected variables, SIR-like models could be fitted directly, as every value would represent the real contagion dynamics. We present the methodology as a general framework, aiming to provide a useful tool for researchers and decision-making actors looking to adapt it for their particular interests.

In a case study on the spreading dynamics of COVID-19 in Chile, we observed the effects that different iconic public health actions taken by the government had on SARS-CoV-2 spread rate R t , readily calculated using our previously reported approach, and discussed the reasons behind these effects interpreted under the light of raw data and corrected daily contagion rates, calculated with the methodology presented in this manuscript. The delay-induced error in raw data slowed-down the reaction time, so the actions were taken too late for restrictive actions or too early for the comeback to normality. Our statistically-driven method corrected such reporting errors, exposing the real contagion dynamics at a given time. The proposed methodology also serves as a non-invasive smoothing process, as it only temporally re-sorts daily reported 13 cases according to their most likely report delay from the real contagion day.

We expect our methodology to serve as a valuable input for researchers and public health practitioners trying to add statistical value to their calculations and predictions in order to vanquish the current SARS-CoV-2 pandemic. Finally, we also expect this work to contribute in raising public awareness on the need for a proper (and standardized) strategy for report and curation of data in the COVID-19 pandemic and other highly-contagious diseases.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. 

Investigating the impact of asymptomatic carriers on covid-19 transmission

Correlation of chest ct and rt-pcr testing in coronavirus disease 2019 (covid-19) in china: a report of 1014 cases

Covid-19 outbreak in italy: estimation of reproduction numbers over two months toward the phase 2. medRxiv

The incubation period of 2019-ncov infections among travellers from wuhan. China. medRxiv

Application of the arima model on the covid-2019 epidemic dataset

Crispr-cas12-based detection of sars-cov-2

Predictability: Can the turning point and end of an expanding epidemic be precisely forecast?

Real-time forecasts and risk assessment of novel coronavirus (covid-19) cases: A data-driven analysis

Comparative estimation of the reproduction number for pandemic influenza from daily case notification data

Álvaro Olivera-Nappa, 2020a. A multi-group seira model for the spread of covid-19 among heterogeneous populations

Real-time estimation of r 0 for supporting public-health policies against covid-19

Complexity of the basic reproduction number (r0)

An interactive web-based dashboard to track covid-19 in real time. The Lancet infectious diseases

Sensitivity of chest ct for covid-19: comparison to rt-pcr

Estimating the generation interval for covid-19 based on symptom onset data

A novel covid-19 epidemiological model with explicit susceptible and asymptomatic isolation compartments reveals unexpected consequences of timing social distancing

A novel method for real time quantitative rt-pcr

Temporal dynamics in viral shedding and transmissibility of covid-19

A brief history of r 0 and a recipe for its calculation

Virology, epidemiology, pathogenesis, and control of covid-19

Qifang bi, forrest k jones, qulu zheng, hannah r meredith, andrew s azman, nicholas g reich, and justin lessler. the incubation period of coronavirus disease 2019 (covid-19) from publicly reported confirmed cases: estimation and application

The incubation period of coronavirus disease 2019 (covid-19) from publicly reported confirmed cases: estimation and application

Covid-19 drives crispr diagnostics: Crispr's role in dna detection, not editing, may fill the gap in covid-19 testing

Early transmission dynamics in wuhan, china, of novel coronavirus-infected pneumonia

Molecular immune pathogenesis and diagnosis of covid-19

Stability issues of rt-pcr testing of sars-cov-2 for hospitalized patients clinically diagnosed with covid-19

Development and clinical application of a rapid igm-igg combined antibody test for sars-cov-2 infection diagnosis

Preparing for covid-19: early experience from an intensive care unit in singapore

Potential preanalytical and analytical vulnerabilities in the laboratory diagnosis of coronavirus disease 2019 (covid-19)

Diagnosis of the coronavirus disease (covid-19): rrt-pcr or ct?

Estimating epidemic exponential growth rate and basic reproduction number

Country-wise forecast model for the basic reproduction number r 0 in the covid-19 outbreak. under review in

Dmakit: A user-friendly web platform for bringing state-of-the-art data analysis techniques to non-specific users

Criteria for discharging a covid-19 infected individual (criterios que se consideran para un paciente covid-19 sin riesgo de contagio)

novel coronavirus disease (covid-19): paving the road for rapid detection and point-of-care diagnostics

Early efforts in modeling the incubation period of infectious diseases with an acute course of illness

Public health policy: Covid-19 epidemic and seir model with asymptomatic viral carriers

An introduction to the basic reproduction number in mathematical epidemiology

Short-term forecasting covid-19 cumulative confirmed cases: Perspectives for brazil

Statsmodels: Econometric and statistical modeling with python

Radiological findings from 81 patients with covid-19 pneumonia in wuhan, china: a descriptive study. The Lancet Infectious Diseases

Zhonghua liu xing bing xue za zhi = Zhonghua liuxingbingxue zazhi

Considerations for quarantine of individuals in the context of containment for coronavirus disease (COVID-19): interim guidance

Protocol for assessment of potential risk factors for coronavirus disease 2019 (COVID-19) among health workers in a health care setting

Cronología de la pandemia de enfermedad por coronavirus de 2020 en chile

The role of age distribution and family structure on covid-19 dynamics: A preliminary modeling assessment for hubei and lombardy

Longitudinal profile of immunoglobulin g (igg), igm, and iga antibodies against the severe acute respiratory syndrome (sars) coronavirus nucleocapsid protein in patients with pneumonia due to the sars coronavirus

Official numbers for the coronavirus outbreak in chile

False-negative of rt-pcr and prolonged nucleic acid conversion in covid-19: Rather than recurrence

Modified seir and ai prediction of the epidemics trend of covid-19 in china under public health interventions

Clinical course and risk factors for mortality of adult inpatients with covid-19 in wuhan, china: a retrospective cohort study

The authors gratefully acknowledge support from the Chilean National Agency for Research and development through ANID PIA Grant AFB180004, and the Centre for Biotechnology and Bioengineering -CeBiB (PIA project FB0001, Conicyt, Chile). DM-O gratefully acknowledges Conicyt, Chile, for PhD fellowship 21181435.