**OSF_Neurological_Psychiatric_Sequelae_COVID.Rdata** contains the data for the 6-month neurological and psychiatric sequelae of COVID-19. Loading the Rdata object in R results in one variable called `res` that contains the results for the comparison of neurological and psychiatric sequelae between cohorts. `res` is an array of 7 lists and each list corresponds to the results for the comparison between two specific cohorts (after propensity score matching). Specifically: - `res[[1]]`: results for the comparison between COVID-19 and influenza (as presented in Table 3 of the manuscript) - `res[[2]]`: results for the comparison between COVID-19 and other respiratory tract infections (as presented in Table 3 of the manuscript) - `res[[3]]`: results for the comparison between COVID-19 and influenza in patients without hospitalisation (as presented in Table 4 of the manuscript) - `res[[4]]`: results for the comparison between COVID-19 and other respitatory tract infections in patients without hospitalisation (as presented in Table 4 of the manuscript) - `res[[5]]`: results for the comparison between COVID-19 with vs. without hospitalisation (as presented in Table 5 of the manuscript) - `res[[6]]`: results for the comparison between COVID-19 with vs. without ITU admission (as presented in Table 5 of the manuscript) - `res[[7]]`: results for the comparison between COVID-19 with vs. without encephalopathy (as presented in Table 5 of the manuscript) Each element of the array contains information about cohorts being compared, the outcomes being compared, and the results of the comparison. Specifically they contain: - `cohort1`: the name of the first cohort (usually *COVID-19*) - `cohort2`: the name of the second cohort (usually the control health event) - `outcomeNames`: the name of the outcomes - `outcomes`: an array with the results for each outcome The results of the Kaplan-Meier analysis for the *i*th outcome are stored in `outcomes[[i]]$KM` which contains the following elements: - `endSurvivals`: the survival (1-event probability) probability at the end of the follow-up period for the two cohorts - `Chi2`: the Chi squared statistics for the statistical test - `p`: the corresponding p-value - `survival$time`: a vector of time in days corresponding to the x-axis of all KM figures - `survival$values1`: the survival (1-event probability) for the first cohort - `survival$values2`: the survival (1-event probability) for the second cohort - `survival$CI1`: the 95% confidence interval for the survival for the first cohort - `survival$CI2`: the 95% confidence interval for the survival for the second cohort - `HR`: the hazard ratio for the comparison - `HR_ci`: the 95% CI for the hazard ratio As an illustration the following few lines of code would plot, for the comparison between COVID-19 and other respiratory tract infections (`res[[2]]`), the outcome probability for "Intracranial haemorrhage (any)" (`outcomes[[8]]`) for the COVID-19 cohort (i.e. first cohort of the comparison so using `values1`): `load('OSF_Neurological_Psychiatric_Sequelae_COVID.Rdata')` `kmResults=res[[2]]$outcomes[[8]]$KM$survival` `plot(kmResults$time,1-kmResults$values1,type='l')`