Six-month Neurological and Psychiatric Outcomes in 236,379 Survivors of COVID-19 ================================================================================ * M. Taquet * J.R. Geddes * M. Husain * S. Luciano * P.J. Harrison ## Abstract **Background** Neurological and psychiatric sequelae of COVID-19 have been reported, but there are limited data on incidence rates and relative risks. **Methods** Using retrospective cohort studies and time-to-event analysis, we estimated the incidence of ICD-10 diagnoses in the 6 months after a confirmed diagnosis of COVID-19: intracranial haemorrhage; ischaemic stroke; Parkinsonism; Guillain-Barré syndrome; nerve/nerve root/plexus disorders; myoneural/muscle disease; encephalitis; dementia; mood, anxiety, and psychotic disorders; substance misuse; and insomnia. Data were obtained from the TriNetX electronic health records network (over 81 million patients). We compared incidences with those in propensity score-matched cohorts of patients with influenza or other respiratory infections using a Cox model. We investigated the effect on incidence estimates of COVID-19 severity, as proxied by hospitalization and encephalopathy (including delirium and related disorders). **Findings** 236,379 patients survived a confirmed diagnosis of COVID-19. Among them, the estimated incidence of neurological or psychiatric sequelae at 6 months was 33.6%, with 12.8% receiving their first such diagnosis. Most diagnostic categories were commoner after COVID-19 than after influenza or other respiratory infections (hazard ratios from 1.21 to 5.28), including stroke, intracranial haemorrhage, dementia, and psychotic disorders. Findings were equivocal for Parkinsonism and Guillain-Barré syndrome. Amongst COVID-19 cases, incidences and hazard ratios for most disorders were higher in patients who had been hospitalized, and markedly so in those who had experienced encephalopathy. Results were robust to sensitivity analyses, including comparisons against an additional four index health events. **Interpretation** The study provides evidence for substantial neurological and psychiatric morbidity following COVID-19 infection. Risks were greatest in, but not limited to, those who had severe COVID-19. The information can help in service planning and identification of research priorities. **Funding** National Institute for Health Research (NIHR) Oxford Health Biomedical Research Centre. Since the COVID-19 pandemic began, there has been concern that survivors might be at an increased risk of neurological disorders. This concern, initially based on findings from other coronaviruses,1 was followed rapidly by case series,2-4 emerging evidence of COVID-19 CNS involvement,5-7 and the identification of mechanisms by which this could occur.8-11 There have been similar concerns regarding psychiatric sequelae of COVID-19,12,13 with evidence showing that survivors are indeed at increased risk of mood and anxiety disorders, as well as dementia, in the three months after infection.14 However, large-scale, robust, and longer-term data are needed if the consequences of the COVID-19 pandemic on brain health are to be identified and quantified. Such information is required both to plan services and to identify research priorities. We used an electronic health records network to investigate the incidence of neurological and psychiatric diagnoses in survivors in the six months after documented clinical COVID-19 infection, and the hazard ratios compared to other health conditions. We explored whether hospitalization and encephalopathy during acute COVID-19 illness affect these risks. The trajectory of hazard ratios across the six-month period was also investigated. ## Methods ### Data and study design The TriNetX Analytics Network ([www.trinetx.com](http://www.trinetx.com)) is a federated network recording anonymized data from electronic health records in 62 healthcare organizations, totalling 81 million patients. The network and its functionalities have been described elsewhere.12 Available data include demographics, diagnoses (using ICD-10 codes), medications, procedures, and measurements (e.g. blood pressure). The healthcare organizations are a mixture of hospitals, primary care, and specialist providers and contribute data from uninsured as well as insured patients. Using the TriNetX user interface, cohorts can be created based on inclusion and exclusion criteria, matched for confounding variables using a built-in propensity score matching algorithm, and compared for outcomes of interest over specified time periods. TriNetX has a waiver from the Western Institutional Review Board since only aggregated counts and summaries of de-identified information are used. Further details are provided in the Appendix (p. 1). ### Cohorts The primary cohort was defined as all patients who had a confirmed diagnosis of COVID-19 (ICD-10 code U07.1). We also produced two matched control cohorts: patients diagnosed with influenza (ICD-10 codes J09-J11) and patients diagnosed with any respiratory tract infection (ICD-10 codes J00–J06, J09–J18, or J20–J22). We excluded patients with a diagnosis of COVID-19 or a positive test for COVID-19 from the control cohorts. We refer to the diagnosis of COVID-19 (in the primary cohort), influenza, or other respiratory tract infections (in the control cohorts) as index events. All cohorts included all patients over the age of ten who had the index event on or after January 20, 2020 (the date of the first recorded COVID-19 case in the USA) and who were still alive at the time of the analysis (December 13, 2020). Further details on cohorts are provided in the Appendix (p. 1). ### Covariates A set of established and suspected risk factors for COVID-19 and for more severe COVID-19 illness was used:15,16 age, sex, race, ethnicity, obesity, hypertension, diabetes, chronic kidney disease, asthma, chronic lower respiratory diseases, nicotine dependence, substance misuse, ischaemic heart disease and other forms of heart disease, socioeconomic deprivation, cancer (and haematological cancer in particular), chronic liver disease, stroke, dementia, organ transplant, rheumatoid arthritis, lupus, psoriasis, and other immunosuppression. To capture these risk factors in patients’ health records, 55 variables were used. More details are provided in the Appendix (pp. 2-3). Cohorts were matched for all these variables, as described below. ### Outcomes We investigated neurological and psychiatric sequelae of COVID-19 in terms of 14 outcomes occurring in the period from 1 to 180 days after the index event: intracranial haemorrhage (ICD-10 codes I60–I62); ischaemic stroke (I63); Parkinson’s disease and Parkinsonism (G20– G21); Guillain-Barré syndrome (G61.0); nerve, nerve root, and plexus disorders (G50–G59); myoneural junction and muscle disease (G70–G73); encephalitis (G04, G05, A86, or A85.8); dementia (F01–F03, G30, G31.0, or G31.83); psychotic, mood, or anxiety disorder (F20–F48) as well as each category separately; substance use disorder (F10–F19), and insomnia (F51.0 or G47.0). For outcomes that are chronic illnesses (e.g. dementia, Parkinson’s disease), patients who had the diagnosis before the index health event were excluded. For outcomes that tend to recur or relapse (e.g. ischaemic strokes, psychiatric diagnoses), we estimated separately the incidence of *first* diagnoses (i.e. excluding those who had a diagnosis before the index event) and the incidence of *any* diagnosis (i.e. including patients who had a diagnosis at some point before the index event). For other outcomes (e.g. Guillain-Barré syndrome), the incidence of *any* diagnosis was estimated. More details are provided in the Appendix (p. 3). Finally, to assess the overall risk of neurological and psychiatric outcomes after COVID-19, we estimated the incidence of any of the 14 outcomes, and the incidence of a first diagnosis of any of the outcomes. Note that this is less than the sum of the incidences of each outcome since some patients have more than one diagnosis. ### Secondary analyses We investigated whether the neurological and psychiatric sequelae of COVID-19 is affected by the severity of the illness in three ways. The incidence of outcomes was estimated separately among patients who: (a) had required hospitalization, within a time window from 4 days before their COVID-19 diagnosis (taken to be the time it might take between clinical presentation and confirmation) to 2 weeks afterwards; (b) had not required hospitalization during that window, and (c) patients who were diagnosed with delirium or other forms of altered mental status during that window; here, we use the term ‘encephalopathy’ to describe this group of patients (see Appendix p. 3 for list of ICD-10 codes).17,18 Differences in outcome incidence between these subgroups might reflect differences in their baseline characteristics. For each outcome, we therefore estimated the HR between patients requiring hospitalization and a matched cohort of patients not requiring hospitalization, and between patients with encephalopathy and a matched cohort of patients without encephalopathy. Finally, HRs were calculated for patients who had not required hospitalization for COVID-19, influenza, or other respiratory infections. To provide benchmarks for the incidence of neurological and psychiatric sequelae, HRs for COVID-19 were compared to four additional matched cohorts of patients diagnosed with health events selected to represent a range of acute presentations during the same time period. The health events were: (i) skin infection, (ii) urolithiasis, (iii) fracture of a large bone, and (iv) pulmonary embolism. For details, see Appendix (p. 4). To test whether differences in sequelae between cohorts could be accounted for by differences in follow up, we counted the average number of health visits that each cohort had during the follow-up period. ### Statistical analyses Propensity score matching19 was used to create cohorts with matched baseline characteristics, and carried out within the TriNetX network. Propensity score 1:1 matching used a greedy nearest neighbour matching approach with a caliper distance of 0.1 pooled standard deviations of the logit of the propensity score. Any characteristic with a standardized mean difference (SMD) between cohorts lower than 0.1 is considered well matched.20 The incidence of each outcome was estimated using the Kaplan-Meier estimator. Comparisons between cohorts were made using a log-rank test. Hazard ratios (HR) were calculated using a proportional hazard model wherein the cohort to which the patient belonged was used as the independent variable. The proportional hazard assumption was tested using the generalized Schoenfeld approach. When the assumption was violated, the time-varying HR was assessed using natural cubic splines fitted to the log cumulative hazard.21 See Appendix (p. 4) for further details. Statistical analyses were conducted in R version 3.4.3 except for the log-rank tests which were performed within TriNetX. Statistical significance was set at two-sided p-values < 0⋅05. ## Results Using electronic health records, the incidence of neurological and psychiatric diagnoses was measured in the six months after COVID-19 infection, and compared to propensity-score-matched cohorts of patients with influenza or other respiratory tract infections. We explored how incidences and HRs differed according to hospitalization and encephalopathy during the acute illness, and how HRs changed over the six months. Key findings are summarised here, with additional details, results, and sensitivity analyses, shown in the Appendix. Adequate propensity-score matching (SMD≤0.1) was achieved for all comparisons and baseline characteristics. Table 1 (and Table S1 in the Appendix) summarises the main demographic features and comorbidities of the COVID-19 cohort (n=236,379), as well as the subgroups who were not-hospitalized (n=190,077), hospitalized (n=46,302), or who had a diagnosis of encephalopathy (n=6,229). Table 1 also presents the estimated diagnostic incidence of the major neurological and psychiatric outcomes over the following six months. All diagnoses are commoner in those who had been hospitalized, and markedly so in those who were encephalopathic during the illness. Results according to gender, race and age are shown in the Appendix (Table S2). View this table: [Table 1:](http://medrxiv.org/content/early/2021/01/24/2021.01.16.21249950/T1) Table 1: Baseline characteristics and major outcomes for the whole COVID-19 cohort, the non-hospitalized and hospitalized sub-groups, and those with encephalopathy during the illness. Only characteristics with a prevalence higher than 5% in the whole population are displayed. For additional baseline characteristics and outcomes see Table S1. Table 2 summarises the HRs for COVID-19 compared to the matched cohorts with influenza (n=105,579; matched baseline characteristics in Table S3) and with other respiratory infections (n=236,038; matched baseline characteristics in Table S4); HRs for diagnostic subcategories are given in Table S5). Kaplan-Meier curves for COVID-19 compared with the other respiratory tract infections cohort are illustrated in Figure 1 (and Figures S1–S3). HRs are statistically significantly higher than 1 for all diagnoses for COVID-19 compared to influenza, except for Parkinsonism and Guillain-Barré syndrome, and are statistically significantly higher than 1 for all diagnoses relative to other respiratory tract infections. Similar results are observed when COVID-19 is compared to the four other index events (see Appendix, Table S6 and Figures S4–S7), except where an outcome has a predicted relationship to the comparator condition (e.g. stroke incidence is higher in the pulmonary embolism cohort; intracranial haemorrhage is commoner in association with fracture of a large bone). View this table: [Table 2:](http://medrxiv.org/content/early/2021/01/24/2021.01.16.21249950/T2) Table 2: Hazard ratios for the major outcomes after COVID-19 compared to influenza and other respiratory infections (RTI). For additional details on cohort characteristics and diagnostic subcategories, see Appendix Tables S3–S5. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/24/2021.01.16.21249950/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/01/24/2021.01.16.21249950/F1) Figure 1: Kaplan-Meier estimates for the incidence of major outcomes after COVID-19 (pink) compared with other respiratory tract infections (blue). 95% confidence intervals are shaded. For diagnostic subcategories and additional details, see Appendix. There were no violations of the proportional hazards assumption for most of the neurological outcomes relative to influenza or other respiratory infections over the six months (Table S7 and Figures S8–S9). The only exception is for intracranial haemorrhage or stroke when compared to other respiratory tract infections (p=0.012 and p=0.032 respectively). For the overall psychiatric disorder category (F20-F48), the HR did vary with time, declining but remaining significantly above 1, indicating that the risk is attenuated but maintained at six months. HRs for COVID-19 compared to the other four index events showed more variation with time, in part reflecting the natural history of the comparator condition (Table S8 and Figures S10–S13). The effect of COVID-19 severity on outcomes was explored in three ways. First, by restricting analyses to matched cohorts of patients who had not required hospitalization (matched baseline characteristics in Tables S9 and S10). As shown in Table 3 (and Figure S14–S15), HRs remain significant in this subgroup. Second, by calculating the HRs for the matched hospitalized vs non-hospitalized COVID-19 cohorts (n=44,927; matched baseline characteristics in Table S11). This shows greater HRs for all outcomes in the hospitalized group, except for nerve/nerve root/plexus disorders (Table 3 and Figure 2). Third, the HRs for all diagnoses were greater for the group who had encephalopathy diagnosed during acute illness compared to a matched cohort who did not (n=6,221; matched baseline characteristics in Table S12) (Table 3 and Figure 2). HRs for the psychiatric diagnoses were less elevated in association with COVID-19 severity than were the neurological diagnoses. View this table: [Table 3:](http://medrxiv.org/content/early/2021/01/24/2021.01.16.21249950/T3) Table 3: Hazard ratios for the major outcomes in non-hospitalized patients after COVID-19 compared to influenza or other respiratory tract infections, and for hospitalized versus non-hospitalized COVID-19 patients, and for encephalopathic versus non-encephalopathic COVID-19 patients. For details on cohort characteristics, see Appendix Tables S9–S12. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/01/24/2021.01.16.21249950/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/01/24/2021.01.16.21249950/F2) Figure 2: Kaplan-Meier estimates for the incidence of major outcomes after COVID-19 comparing patients requiring hospitalization (solid blue line) with matched patients not requiring hospitalization (dashed blue line), and comparing those who had encephalopathy (solid red line) with matched patients who did not have encephalopathy (dashed red line). For clarity, confidence intervals are omitted, but are shown in the Appendix, Figure S16. We inspected for other factors that might influence the findings. (1) The results regarding hospitalization (which we had defined as occurring up to 14 days after diagnosis) could be confounded by admissions due to an early complication of COVID-19 rather than to COVID-19 itself. This was explored by excluding outcomes during this period, and findings remained, albeit with many HRs being reduced (Table S13). (2) COVID-19 survivors had fewer healthcare visits during the six-month period compared to the other cohorts (Table S14). Hence their higher incidence of many diagnoses is not simply due to having had more diagnostic opportunities. ## Discussion A range of adverse neurological and psychiatric outcomes after COVID-19 have been predicted and reported. Data from a large electronic health records network confirm this prediction, and provide estimates of their incidence (Table 1), and the HRs compared to other matched health conditions occurring contemporaneously with the COVID-19 pandemic (Table 2 and Figure 1). The severity of COVID-19 infection had a clear effect on subsequent neurological diagnoses (Table 3 and Figure 2). Overall, COVID-19 was associated with increased risk, but the incidences and the HRs were greater in patients who required hospitalization, and markedly so in those who developed encephalopathy, even after extensive propensity score matching for other factors (e.g. age, prior cerebrovascular disease). Potential mechanisms for this association include viral invasion of the central nervous system,10,11 hypercoagulable states,22 and neural effects of the immune response.9 On the other hand, it is also notable that incidence of these diagnoses was increased even in the COVID-19 cases who had not required hospitalization. Some specific neurological diagnoses merit mention. Consistent with several other reports,23,24 the risk of cerebrovascular events (stroke and intracranial haemorrhage) is elevated after COVID-19, with the incidence of stroke rising to almost 1 in 10 (or 3 in 100 for a first stroke) in those who had been encephalopathic (i.e. had experienced delirium or other altered mental status; see Appendix p. 3). A similarly increased risk for stroke relative to influenza was recently reported.25 Whether COVID-19 is associated with Guillain-Barre syndrome has been unclear;26 our data indicate that the incidence is increased in hospitalized and encephalopathic patients, but not significantly in the whole cohort. There have been concerns about post-COVID-19 Parkinsonian syndromes, driven by the encephalitis lethargica epidemic that followed the 1918 influenza pandemic.27 Our data provide some support for this possibility although the incidence was low and not all HRs significant. It is possible that Parkinsonism is a delayed outcome, in which case a clearer signal may emerge with longer follow-up. Our previous study reported preliminary evidence for an association between COVID-19 and dementia.14 The present data confirm this association, with a HR of 2.33 compared to influenza, and 1.71 compared to other respiratory infections. Although the estimated incidence is modest in the whole COVID-19 cohort (0.67%), 1.46% of hospitalized cases, and 4.72% who were encephalopathic receive a first diagnosis of dementia within six months. Our findings regarding psychiatric disorders are consistent with the 3-month outcome data we reported in a smaller number of cases using the same network,14 and show the HR remains elevated, although decreasing, at the six-month period. Unlike the earlier study, there is also an increased risk of psychotic disorders (likely reflecting the larger sample size and longer duration of follow-up), and substance use disorders are also more common (with the exception of the incidence of a first diagnosis of substance misuse compared to other respiratory infections). Hence, like the neurological outcomes, the psychiatric sequelae of COVID-19 appear widespread, and persist to and probably beyond six months. Compared to the neurological disorders, HRs for the common psychiatric disorders (mood and anxiety disorders) show a weaker relationship with hospitalization and encephalopathy (Table 1). This may indicate that their occurrence reflects, at least partly, the psychological and other implications of a COVID-19 diagnosis rather than being a direct manifestation of the illness. HRs for most neurological outcomes were constant, and hence the risks associated with COVID-19 persist until the six-month time point. Longer-term studies will be needed to ascertain the duration of risk, and the trajectory for individual diagnoses. Our findings are relatively robust, given the sample size, the propensity score matching, and the results of the sensitivity and secondary analyses. Nevertheless, they have weaknesses inherent to an electronic health records study,28 such as the unknown completeness of records, lack of validation of diagnoses, and sparse information on socioeconomic and lifestyle factors. These issues primarily affect the incidence estimates rather than hazard ratios, but the choice of cohorts against which to compare COVID-19 outcomes influences the HRs (see Appendix Table S6). The analyses regarding encephalopathy deserve a note of caution. Even amongst hospitalized patients, only about 11% received this diagnosis, even though much higher rates would be expected.18,29 Under-recording of delirium and other altered mental states during acute illness is well known, and likely means that the diagnosed cases had prominent and/or sustained features; as such, results for this group should not be generalised to all patients with delirium. Finally, a study of this kind can only demonstrate associations; efforts to identify mechanisms and assess causality will require prospective cohort studies or more elaborate study designs. ## Supporting information Appendix [[supplements/249950_file03.pdf]](pending:yes) ## Data Availability The TriNetX system returned the results of these analyses as csv files which were downloaded and archived. Data presented in this paper and the Appendix can be freely accessed at [URL to be added on publication]. Additionally, TriNetX will grant access to researchers if they have a specific concern (via the third-party agreement option). ## Declarations of interest SL is an employee of TriNetX Inc. ## Data sharing The TriNetX system returned the results of these analyses as csv files which were downloaded and archived. Data presented in this paper and the Appendix can be freely accessed at [URL to be added on publication]. Additionally, TriNetX will grant access to researchers if they have a specific concern (via the third-party agreement option). ## Acknowledgments PJH and MT were granted unrestricted access to the TriNetX Analytics network for the purposes of research, and with no constraints on the analyses performed nor the decision to publish. Work supported by the National Institute for Health Research (NIHR) Oxford Health Biomedical Research Centre (grant BRC-1215-20005). The views expressed are those of the authors and not necessarily those of the UK National Health Service, NIHR, or the UK Department of Health. MT is an NIHR Academic Clinical Fellow. * Received January 16, 2021. * Revision received January 16, 2021. * Accepted January 24, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Rogers JP, Chesney E, Oliver D, et al. Psychiatric and neuropsychiatric presentations associated with severe coronavirus infections: a systematic review and meta-analysis with comparison to the COVID-19 pandemic. Lancet Psychiatry 2020; 7: 611–27. 2. 2.Ellul MA, Benjamin L, Singh B, et al. Neurological associations of COVID-19. Lancet Neurol 2020; 19: 767–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(20)30221-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32622375&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 3. 3.Varatharaj A, Thomas N, Ellul MA, et al. Neurological and neuropsychiatric complications of COVID-19 in 153 patients: a UK-wide surveillance study. Lancet Psychiatry 2020; 7: 875–82. 4. 4.Paterson RW, Brown RL, Benjamin L, et al. The emerging spectrum of COVID-19 neurology: clinical, radiological and laboratory findings. Brain 2020; 143: 3104–20. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 5. 5.Kremer S, Lersy F, Anheim M, et al. Neurologic and neuroimaging findings in patients with COVID-19. A retrospective multicenter study. Neurology 2020; 95: e1868–82. 6. 6.Pezzini A, Padovani A. Lifting the mask on neurological manifestations of COVID-19. Nat Rev Neurol 2020; 16 : 636–44. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 7. 7.Raman B, Cassar MP, Tunnicliffe EM, et al. Medium-term effects of SARC-CoV-2 infection on multiple vital organs, exercise capacity, cognition, quality of lige and mental health, post- hospital discharge. EclinicalMedicine (in press) 8. 8.Iadecola C, Anrather J, Kamel H. Effects of COVID-19 on the nervous system. Cell 2020; 183: 16–27. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.08.028&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32882182&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 9. 9.Kreye J, Reincke SM, Pruss H. Do cross-reactive antibodies cause neuropathology in COVID- 19? Nat Rev Immunol 2020; 20: 645–6. 10. 10.Meinhardt J, Radke J, Dittmayer C, et al. Olfactory transmucosal SARS-CoV-2 invasion as a port of central nervous system entry in individuals with COVID-19. Nat Neurosci 2020. 11. 11.Rhea EM, Logsdon AF, Hansen KM, et al. The S1 protein of SARS-CoV-2 crosses the blood- brain barrier in mice. Nat Neurosci 2020. 12. 12.Holmes EA, O’Connor RC, Perry VH, et al. Multidisciplinary research priorities for the COVID-19 pandemic: a call for action for mental health science. Lancet Psychiatry 2020;7:547–60. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 13. 13.Vindegaard N, Benros ME. COVID-19 pandemic and mental health consequences : Systematic review of the current evidence. Brain Behav Immun 2020; 89: 531–42. 14. 14.Taquet M, Luciano S, Geddes JR, Harrison PJ. Bidirectional associations between COVID-19 and psychiatric disorder: retrospective cohort studies of 62 354 COVID-19 cases in the USA. Lancet Psychiatry 2020 (AOL 9 Nov 2020). 15. 15.de Lusignan S, Dorward J, Correa A. Risk factors for SARS-CoV-2 among patients in the Oxford Royal College of General Practitioners Research and Surveillance Centre primary care network : a cross-sectional study. Lancet Infect Dis 2020; 20: 1034–42. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 16. 16.Williamson EJ, Walker AJ, Bhaskaran K, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature 2020; 584: 430–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2521-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 17. 17.Slooter AJL, Otte WM, Devlin JW, et al. Updated nomenclature of delirium and acute encephalopathy: statement of ten Societies. Intensive Care Med 2020: 46: 1020–2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00134-019-05907-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 18. 18.Oldham MA, Slooter AJC, Cunningham C, et al. Characterising neuropsychiatric disorders in patients with COVID-19. Lancet Psychiatry 2020; 7: 932–3. 19. 19.Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res 2011; 46: 399–424. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00273171.2011.568786&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21818162&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291533400002&link_type=ISI) 20. 20.Haukoos JS, Lewis RJ. The propensity score. JAMA 2015; 314: 1637–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2015.13480&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26501539&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 21. 21.Royston P, Parmar MKB. Flexible parametric proportional-hazards and proportional-odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects. Stat Med 2002; 21: 2175–97. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.1203&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12210632&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177051300004&link_type=ISI) 22. 22.Panigada M, Bottino N, Tagliabue P, et al. Hypercoagulability of COVID-19 patients in intensive care unit: A report of thromboelastography findings and other parameters of hemostasis. J Thromb Haemost 2020; 18: 1738–42. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 23. 23.Siow I, Lee, KS, Zhang JJY, et al. Stroke as a neurological complication of COVID-19: a systematic review and meta-analysis of incidence, outcomes and predictors. J Stroke Cerebrovasc Dis 2020; 30: 105549. 24. 24.Hernandez-Fernandez F, Valencia HS, Barbella-Aponte R, et al. Cerebrovascular disease in patients with COVID-19: neuroimaging, histological and clinical description. Brain 2020 (AOL 9 July 2020). 25. 25.Xie Y, Bowe B, Maddukuri G, Al-Aly Z. Comparative evaluation of clinical manifestations and risk of death in patients admitted to hospital with covid-19 and seasonal influenza: cohort study. BMJ 2020; 371: m4677. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE4OiIzNzEvZGVjMTRfMTcvbTQ2NzciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMS8yNC8yMDIxLjAxLjE2LjIxMjQ5OTUwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 26. 26.Keddie S, Pakpoor J, Mousele C, et al. Epidemiological and cohort study finds no association between COVID-19 and Guillain-Barre syndrome. Brain (AOL 14 December 2020). 27. 27.Hoffman LA, Vilensky JA. Encephalitis lethargica: 100 years after the epidemic. Brain 2017; 140: 2246–51. 28. 28.Casey JA, Schwartz BS, Walter F, et al. Using electronic health records for population health research: a review of methods and applications. Annu Rev Public Health 2016; 37: 61–81. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-publhealth-032315-021353&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26667605&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F01%2F24%2F2021.01.16.21249950.atom) 29. 29.Docherty AB, Harrison EM, Green CA, et al. Features of 20 133 UK patients in hospital with COVID-19 using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. BMJ 2020; 369: m1985. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNjkvbWF5MjJfMS9tMTk4NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAxLzI0LzIwMjEuMDEuMTYuMjEyNDk5NTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)