Using Generalized Structured Additive Regression Models to determine factors associated with and clusters for COVID-19 hospital deaths in South Africa ====================================================================================================================================================== * Innocent Maposa * Richard Welch * Lovelyn Ozougwu * Tracy Arendse * Caroline Mudara * Lucille Blumberg * Waasila Jassat ## Abstract **Background** The first case of COVID-19 in South Africa was reported in March 2020 and the country has since recorded over 3.6 million laboratory-confirmed cases and 100 000 deaths as of March 2022. Transmission and infection of SARS-CoV-2 virus and deaths in general due to COVID-19 have been shown to be spatially associated but spatial patterns in in-hospital deaths have not fully been investigated in South Africa. This study uses national COVID-19 hospitalization data to investigate the spatial effects on hospital deaths after adjusting for known mortality risk factors. **Methods** COVID-19 hospitalization data and deaths were obtained from the National Institute for Communicable Diseases (NICD). Generalized structured additive logistic regression model was used to assess spatial effects on COVID-19 in-hospital deaths adjusting for demographic and clinical covariates. Continuous covariates were modelled by assuming second-order random walk priors, while spatial autocorrelation was specified with Markov random field prior and fixed effects with vague priors respectively. The inference was fully Bayesian. **Results** The risk of COVID-19 in-hospital mortality increased with patient age, with admission to intensive care unit (ICU) (aOR=4.16; 95% Credible Interval: 4.05-4.27), being on oxygen (aOR=1.49; 95% Credible Interval: 1.46-1.51) and on invasive mechanical ventilation (aOR=3.74; 95% Credible Interval: 3.61-3.87). Being admitted in a public hospital (aOR= 3.16; 95% Credible Interval: 3.10-3.21) was also significantly associated with mortality. Risk of in-hospital deaths increased in months following a surge in infections and dropped after months of successive low infections highlighting crest and troughs lagging the epidemic curve. After controlling for these factors, districts such as Vhembe, Capricorn and Mopani in Limpopo province, and Buffalo City, O.R. Tambo, Joe Gqabi and Chris Hani in Eastern Cape province remained with significantly higher odds of COVID-19 hospital deaths suggesting possible health systems challenges in those districts. **Conclusion** The results show substantial COVID-19 in-hospital mortality variation across the 52 districts. Our analysis provides information that can be important for strengthening health policies and the public health system for the benefit of the whole South African population. Understanding differences in in-hospital COVID-19 mortality across space could guide interventions to achieve better health outcomes in affected districts. Key words * COVID-19 * spatial effects * health systems * hospitalizations * nonlinear effects * clusters * deaths ## BACKGROUND Over the past two years, the world has grappled with COVID-19, a disease caused by a novel coronavirus, SARS-CoV-2. This was first noted in December 2019 in Wuhan, China, where a cluster of patients with pneumonia of unknown cause was identified (1–3) South Africa recorded its first confirmed COVID-19 case in March 2020 and the government took decisive action to mitigate the spread of the disease through implementing a state of disaster and then adjusted mitigation levels(4) The first wave (peak) of COVID-19 infections occurred in July 2020 in South Africa and four more waves followed (5–8). The most recent waves dominated by Omicron sub-variants have been less severe due to high levels of immunity from vaccination and prior infection in South Africa (5,9,10). As of March 2022, there has been more than 3.6 million confirmed COVID-19 cases and 100 000 confirmed deaths in South Africa, making it the most affected country in Sub-Saharan Africa. Low- and middle-income countries reported large numbers of COVID-19 related hospitalizations and deaths. New Variants of Concern (VOC) emerged including Alpha, Beta and Delta, that were more transmissible and associated with more severe disease (9). The highly transmissible Omicron BA.1 (B.1.1.529; hereafter BA.1) VOC heralded a new surge of infections in South Africa from November 2021 (10). Two Omicron sub-lineages (BA.4 and BA.5) dominated a surge of new infections from April 2022 (10). The Omicron variant, though highly transmissible, was not as virulent (5)hence had lower hospitalizations and deaths comparably (9,10). COVID-19 mortality rates were highly disproportionate in South Africa, with the elderly, males, people of colour, those with comorbidities and those admitted at public health facilities and in certain provinces, being at higher risk (6). Several studies have shown that the risk for severe COVID-19 disease were disproportionally born among minority communities (11–13). A South African study highlighted that older age, male sex, minority race groups and lower socioeconomic status (SES) are associated with severe COVID-19 disease and deaths (6). In addition to these factors, patient outcomes generally differ by health sector (public or private hospital) and facility type (district, regional or tertiary hospital) (14). Jassat et al described the demographic and clinical characteristics of individuals admitted to hospital with laboratory confirmed COVID-19 throughout South Africa in first and second waves using the DATCOV surveillance data and also assessed risk factors for in-hospital mortality(5,6,8). In addition, Jassat et al 2022 described and compared admissions and deaths by age, sex, race and health sector as proxy for SES using same database(6). Waasila et al allude to the role of socio-economic status, race and health care facility type in COVID19 mortality. These can all latently be reflected spatially with similar race groups more likely to be clustered in same neighbourhoods and sharing the same socio-economic status. South Africa being a highly unequal society has a Gini coefficient of 63% (15). Access to efficient health care was reported to be inequitable during the pandemic, and poor communities relied mainly on public health facilities as a consequence of structural inequalities still prevalent in the country (6). The authors also observe a relationship between structural inequality and COVID-19 susceptibility and severity. The COVID-19 epidemic revealed strong spatial heterogeneity in the spread of infections in countries such as France and Italy (1). A study by Sannigrahi et al 2020 assessed spatial association between socio-demographic variables and COVID-19 cases and deaths and showed distribution of cases and deaths were spatially heterogenous across Europe(16). In England, Sartorius et al 2021 showed that a number of administrative areas (small areas or contiguous small areas) appeared to be at a significantly elevated risk of high COVID-19 transmission and also at increased risk for higher mortality rates(17). Jassat et al 2022, in South Africa, showed that hospitalized COVID-19 patients in the Eastern Cape and Limpopo provinces were 60% and 50% more at risk of in-hospital deaths respectively compared to those hospitalized in the Western Cape province highlighting spatial heterogeneity in hospital deaths(6). It is evident from these studies that variations in COVID-19 hospital mortality can be found across the socio-economic spectrum with vulnerable communities being at highest risk. This latter study, although it assessed the regional effect, did not account for spatial autocorrelation and also did not model the spatial effect at a finer resolution which can be helpful in defining interventions and deciding effective public health policy. In recent years, there has been a growing interest in the application of spatial analysis and modelling techniques as a tool for in-depth understanding of public health problems including identifying hotspots, spatial distribution, patterns and effects (1,16–20). SARS-CoV-2 transmission and infections, and COVID-19 deaths, have been shown to be spatially distributed in Europe(16), England(17), and France(1). In Afghanistan, COVID-19 cases were shown to be spatially distributed (21). However, in South Africa, though several studies have been done to determine factors associated with hospital COVID-19 mortality (6), none have considered simultaneously modelling spatial effects, fixed effects and nonlinear effects in investigating factors associated with COVID-19 hospital deaths. In addition, to the best of our knowledge, none of the studies have used flexible structured additive logistic regression models within the Bayesian framework to estimate district spatial effects adjusting for other known COVID-19 mortality risk factors. This study seeks to determine 1) the district level clusters and spatial effects or variability of COVID-19 hospital deaths and 2) identify other factors associated with hospital deaths adjusting for spatial correlation. ## METHODS ### Source of data and sample The data used in this study was obtained from the NICD. The data was collected using DATCOV an active national COVID-19 hospital surveillance system. The database contains 486 344 hospitalizations with corresponding minimal key individual level data points including age, month and year of hospitalization, patient gender and clinical markers like being on ventilation, on oxygen and admission in ICU. The facility, subdistrict, district and province of the hospital where the patient was admitted are also recorded. We assessed district level spatial effects together with other important covariates that were chosen on the basis of biological plausibility and available evidence. These spatial effects need to be modelled and estimated simultaneously with linear and possibly nonlinear effects. We used district level effects to allow for spatial correlation and any other unknown regional heterogeneity of COVID-19 hospital deaths. ### Bayesian Structured Additive Logistic Regression Model Let *y**ij* be the hospital death status for a hospitalized patient *i* in desctrict. *y**ij* = 1 if the patient *i* in district *j* died in hospital and *y**ij* = 0 otherwise. A vector *X**ij* *=* (*x**ij*1, *x**ij*2,*…, x**ijp*)′ contains *p* continuous covariate random variables and *Z**ij* *=* (*z**ij*1, *z**ij*2,*…, z**ijr*)′ contains some r categorical variables. In our study, *p*= 2 and *r*= 4. This study assumes that the dependent variable, *y**ij* is a Bernoulli distributed random variable with *y**ij* |*p**ij* *∼ Bernoulli* (*p**ij*) with an unknown *E* (*y**ij*) *= p**ij*, being related to the covariates through the link function ![Formula][1] The link function in this equation is known as the logit link, *β* is the *p* dimensional vector of coefficients for the continuous random variables, and *θ* is an *r* dimensional vector of coefficients for categorical random variables. In order to assess for both non-linear effects of continuous random variables and spatial autocorrelation in our data we employed a semi-parametric model which utilizes a penalized regression approach(22,23). The penalized regression approach is a non-parametric method of ordinary least squares (OLS) which relaxes the highly restrictive linear predictor for a versatile semi-parametric predictor. The flexible semi-parametric predictor is defined by: ![Formula][2] where *f**v* (·) represents the non-linear twice differentiable smooth function for the continuous covariates and *f**spat* (*s**j*) is the variable that denotes the spatial effects for each district. In our study, as in Ngesa et al(24), we consider a convolution approach to the spatial effects. The assumption is that the spatial effects can be decomposed into two pure components, that is, spatially structured and spatially unstructured effects given as ![Graphic][3]. The final model for our study then becomes: ![Formula][4] We fit a generalized structured additive logistic regression model for COVID-19 hospital deaths using Markov chain Monte Carlo (MCMC) simulations. ### Data analysis Statistical analysis was performed in R/RStudio v 4.1.0. Spatial distribution and patterns including Global Moran I and local Moran I for assessing autocorrelation and local clustering were assessed for district level aggregated rates. A conditional autoregressive (CAR) generalized structured additive logistic regression model with binomial link was fit using BayesX R package accounting for spatial effects (25). In addition, non-linear effects of hospital COVID-19 deaths for some continuous covariates were also assessed. The structured additive logistic regression fits a multivariable model with patient sex, facility type, having been on oxygen, ventilator or in ICU as fixed effects, patient age and month (proxy for temporal COVID-19 evolution) as non-linear effects and district as spatial effects. The analysis included data from all waves up to March 2022. All tests were two-sided and a p-value of less or equal to 0.05 was considered to indicate statistical significance. The 95% credible intervals were reported with adjusted Odds Ratios for the full Bayesian inference. Ethical approval was obtained from the Human Research Ethics Committee (Medical) of University of the Witwatersrand for the DATCOV surveillance programme (M2010108). ## RESULTS Our analysis included data from all provinces in South Africa over the period March 2020 to March 2022. The database contained 484 699 COVID-19 patients who were hospitalized in either private or public health care facilities. The patient data including hospitalization facility was reported at the individual and subdistrict level. The median age for patients who died was 63 years (IQR: 53-73) while for those who survived or were discharged alive was 48 years (IQR: 34-62). Over half of the patients, 60.9% (n=62928), who died in hospital were admitted and received care from public health care facilities (Table 1). View this table: [Table 1:](http://medrxiv.org/content/early/2022/09/17/2022.09.16.22280020/T1) Table 1: COVID19 hospitalized patient characteristics and outcomes, South Africa, March 2020 – March 2022 District level spatial distribution of hospital mortality rates are described in Figure 1 below. Figure 1A shows the distribution of observed crude COVID-19 hospital deaths rates. Most districts in the Eastern Cape including Amathole, Buffalo City, O.R. Tambo, Joe Gqabi, Chris Hani, Cacadu, Nelson Mandela, Alfred Nzo and those in the Limpopo province such as Vhembe, Capricorn, Mopani and Greater Sekhukhune showed elevated observed rates of COVID-19 hospital deaths. Districts such as Xhariep, Mangaung, John Taolo Gaetsewe, Bojanala, Dr Kenneth Kaunda and Sisonke had low observed rates of COVID-19 hospital deaths. Figure 1B shows the significant hot and cold spots (clusters) for hospital COVID-19 deaths using local indicators of spatial autocorrelation (LISA). Most districts in the Eastern Cape show above expected high rates of hospital deaths in neighbouring areas and similar pattern also observed in some districts (eg Mopane) in Limpopo province. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/17/2022.09.16.22280020/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/09/17/2022.09.16.22280020/F1) Figure 1: Hospital COVID-19 deaths rates distribution, South Africa, March 2020 – March 2022 Table 2 shows patient and clinical care factors associated with COVID-19 deaths adjusting for period of hospitalization (month) and patient age as nonlinear effects and district level spatial autocorrelation. Adjusting for spatial correlation and nonlinear effects, male gender (adjusted odds ratio [aOR]=1.22; 95% Credible Interval [CI]: 1.20-1.24), admission in a public hospital (aOR=3.17; 95% CI: 3.11-3.23), being on oxygen (aOR=1.37; 95% CI: 1.35-1.39), admitted in ICU (aOR=4.29; 95% CI: 4.18-4.40) and invasive mechanical ventilation (aOR=3.51; 95% CI: 3.40-3.64) were all significantly associated with hospital deaths among admitted COVID-19 patients. View this table: [Table 2:](http://medrxiv.org/content/early/2022/09/17/2022.09.16.22280020/T2) Table 2: Patient and clinical care factors adjusted effects associated with mortality, South Africa, March 2020 – March 2022 ### Nonlinear effects of patient age and admission month Figure 2 shows the nonlinear effect of time in months on the likelihood of hospital COVID-19 deaths in South Africa over the years 2020 to 2022. The figure provides posterior mean of the smooth time function and its corresponding 80% as well as 95% Credible Intervals. It is clear from the figure that the risk of hospital deaths fluctuated significantly over the peak of the COVID-19 waves. The association between time in months and COVID-19 hospital deaths was nonlinear and assuming linear effect would bias the result and lead to incorrect interpretation of time effect. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/17/2022.09.16.22280020/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/09/17/2022.09.16.22280020/F2) Figure 2: Nonlinear effect of time in months on COVID-19 in-hospital deaths, South Africa, March 2020 – March 2022 Figure 3 shows the nonlinear effect of age on hospital COVID-19 deaths. The likelihood of hospital deaths increased with increasing age in South Africa over the peak periods of the epidemic. The figure provides posterior mean of the smooth age function and its corresponding 80% as well as 95% Credible Intervals. The association between patient age and COVID-19 hospital deaths was nonlinear and assuming linear or peace-wise effect, as with month of hospitalization, would bias the result and lead to incorrect interpretation of age effects. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/17/2022.09.16.22280020/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/09/17/2022.09.16.22280020/F3) Figure 3: Nonlinear effect of patient age on COVID-19 in-hospital deaths, South Africa, March 2020 – March 2022 We investigated the spatial effects using the structured additive logistic regression model and figure 4 shows the residual spatially associated hospital deaths after adjusting for fixed and nonlinear effects. Districts with blue colour show lower odds of hospital deaths while those with red colour indicate significant higher odds of hospital deaths. From Figure 4, there is clear evidence of spatial variation of hospital COVID-19 deaths after controlling for some known risk factors. The spatial random effect, that is, the heterogeneity between districts in COVID-19 deaths dominates the residual spatial variation explaining about 83% of the variance showing differences in in-hospital COVID-19 deaths by district. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/17/2022.09.16.22280020/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/09/17/2022.09.16.22280020/F4) Figure 4: Residual spatial odds of in-hospital deaths, South Africa, March 2020 – March 2022 ## DISCUSSION We examined factors associated with in-hospital COVID-19 mortality in South Africa between March 2020 and March 2022, adjusting for spatial autocorrelation at district level, fixed effects and some non-linear temporal and age effects. Besides confirming that older age, male gender, admission in ICU, being treated with supplemental oxygen and invasive mechanical ventilation, predisposed hospitalized COVID-19 patients to a higher risk of in-hospital mortality, our study highlighted substantial heterogeneity in mortality across districts. Our study describes the spatial distribution of COVID-19 in-hospital deaths in the whole population in South Africa adjusting for several factors including month of admission as a proxy for epidemic waves. The impact of spatial neighbourhood dependents and heterogeneity have not previously been explored at higher resolution in South Africa regarding COVID-19 in-hospital mortality. To respond effectively to new epidemics, it is essential to unpack the spatial dynamics of in-hospital deaths at higher spatial resolution such as at subdistrict or district level in order to identify priority areas for intervention to improve health outcomes. Our spatial model adjusted for temporal effects, suggested that age had a nonlinear effect on in-hospital mortality with increasing risk for increasing age as was determined in other previous studies (6,12,26–29). The study also provides evidence that in-hospital mortality was clustered geographically with several districts in Eastern Cape and Limpopo showing elevated risk. Previous studies showed that hospitalized COVID-19 patients in the Eastern Cape and Limpopo provinces were at higher odds of death (6), however, our study goes further in identifying the districts in those provinces which predisposed patients to higher odds of in-hospital mortality. A study in Brazil similarly explored the spatial distribution of COVID-19 cases and deaths in paediatric population and showed that forty municipalities had higher mortality, and these were observed in regions with poor socio-economic indicators and greatest health disparities (19). Another study in South Africa assessed spatial distribution of deaths at provincial level and other COVID-19 outcomes showed that Eastern Cape had higher risk for deaths (30,31). Understanding spatial heterogeneity in relation to the socio-environmental determinants and COVID-19 related outcomes is central to targeting interventions for vulnerable populations (1). Our results indicate substantial geographical variations in the distribution of COVID-19 in-hospital mortality across South Africa districts. In South Africa, 21% of the population have access to private health care while the remainder rely on public health care support (6). There are also inequities in health care access in South Africa, with poorer access and availability of health facilities in rural provinces such as Eastern Cape, KwaZulu-Natal and Limpopo. Quality of care has also historically been poor in many public health facilities (14). A study that assessed inequalities in access to health care in South Africa showed large differentials in private health care admissions between people living in the poorer, more rural provinces like Limpopo, Eastern Cape and Mpumalanga compared to more urban provinces like Gauteng (32). We provide evidence of spatial clustering and spatial heterogeneity in in-hospital mortality in South Africa. Highlighting districts at excess risk of COVID-19 in-hospital mortality can help guide local health care system policies to better protect vulnerable population subgroups. Health care access inequality is not a new aspect in South Africa, however, when it contributes to some groups of the population left at higher risk of in-hospital mortality, this calls for immediate action for implementation of policies that enhance equal access to and improved delivery of health care(6,14,32). Besides showing spatial structure in the COVID-19 in-hospital mortality in South Africa, our study using the non-linear effects of months as a proxy for temporal evolution of the epidemic highlight higher risk of deaths during the second and third waves. Jassat et al similarly reveal in their study that in-hospital case-fatality ratio during the Omicron wave was 10·7%, compared with 21·5% during the first wave, 28·8% during the second wave, and 26·4% during the third wave (7,8,33). The Beta and Delta VOCs drove the second and third COVID-19 waves in South Africa. Our study reveals an ebbing off in the in-hospital deaths during the last months which were driven by the latter variants of Omicron. Omicron marked a change in the SARS-CoV-2 epidemic curve, clinical profile, and deaths in South Africa (33) Davies et al showed similar patterns in the Western Cape province (9). Major strengths of this study include the application of flexible structured additive logistic regression model within the Bayesian framework which allows for exploration of spatial association with health outcomes at higher spatial resolutions as well as allow for inclusion of non-linear effects and fixed effects. In addition, this study utilizes a COVID-19 hospital admissions high-quality national data repository with their associated individual level outcomes ensuring ample admission sample data as well as outcomes for analysis at district level. This study is not without some limitations. For this analysis, we adjusted for some of the known risk factors for COVID-19 in-hospital deaths, however, there may be other covariates which we may not have included. A further limitation may be that the DATCOV database does not distinguish between patients who were admitted with coincidental positive SARS-CoV-2 test and those with COVID-19. As highlighted in the preceding limitations, there was limited information on other important variables and ancillary predictors that can be used to model in-hospital mortality. Therefore, there is a need to utilize the same modelling framework to account for the impact of these and other factors especially aggregated at subdistrict or district levels. As in one Brazilian study (19), our findings confirmed the higher burden of COVID-19 in-hospital mortality in some districts in Limpopo and Eastern Cape provinces. This may be linked to poor health care access including limited access to private health care (32) due to unaffordability of medical insurance. To reduce the burden of in-hospital deaths due to potentially new epidemics and COVID-19 in particular, the South African government and private stakeholders need to build strong and resilient health care infrastructure which can be accessed equally by both the poor and rich independent of class or race especially in vulnerable districts. Modelling the effects of underlying factors and in-hospital disease mortality at small spatial scale is essential to plan effective control strategies for disease outcome risks (20). The relationship between COVID-19 in-hospital mortality rates and sociodemographic, clinical and districts spatial dependents can consequently guide the development of specific intervention actions for these places. ## Conclusions The results from this study reveal substantial COVID-19 in-hospital mortality variation across the districts in South Africa. This highlights the importance of modelling spatial patterns simultaneously with fixed and nonlinear effects of continuous covariates to identify clusters at high risk of health outcome. The flexible approach to modelling data that has spatial patterns helps to account for possible loss of efficiency due to spatial correlation that spatial patterns can induce in data. Our analysis suggests notable COVID-19 hospital deaths clustering in some districts in Limpopo and Eastern Cape provinces and this information can be important in strengthening the public health and health systems policies for the benefit of the whole South African population. Understanding differences in in-hospital COVID-19 mortality across space could guide interventions to achieve better health outcomes. The findings of this analysis reveal that strengthening health systems, specifically critical intensive care, oxygen support, and availability of health care workers, could lead to substantial reductions in mortality rates in South Africa in the case of a resurgence of COVID-19 epidemic or new epidemics. It is also important for public health authorities to manage medical care effectively and in a timely manner, as well as to direct the intensity and type of interventions needed to overcome the pandemic. ## Data Availability The dataset analysed for the manuscript is available upon reasonable request. The data dictionary is available at request to the co-author: waasilaj{at}nicd.ac.za ## FUNDING DATCOV as a national surveillance system, was initially funded by the NICD and the South African National Government, and subsequently by the support of the American people through the United States Agency for International Development (USAID) via the mechanism awarded to Right to Care. The contents of this study are the sole responsibility of the authors and do not necessarily reflect the views of USAID, PEPFAR, or the United States Government. The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication. ## DATA SHARING AGREEMENT The dataset analysed for the manuscript is available upon reasonable request. The data dictionary is available at request to the co-author: waasilaj{at}nicd.ac.za ## DECLARATION OF INTEREST The authors declare that there are no conflicts of interest. ## ACKNOWLEDGEMENTS We acknowledge the NICD team responsible for reporting hospitalisation data. Thanks to the National Department of Health and the NICD for support and oversight. Our gratitude to the laboratories, clinicians and data teams in all public and private sector hospitals throughout the country reporting cases and hospitalisation data, who are acknowledged and listed as the DATCOV author group: [https://www.nicd.ac.za/diseases-a-z-index/covid-19/surveillance-reports/daily-hospital-surveillance-datcov-report/](https://www.nicd.ac.za/diseases-a-z-index/covid-19/surveillance-reports/daily-hospital-surveillance-datcov-report/) * Received September 16, 2022. * Revision received September 16, 2022. * Accepted September 17, 2022. * © 2022, 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.Deguen S, Kihal-Talantikite W. Geographical pattern of COVID-19-related outcomes over the pandemic period in France: A nationwide socio-environmental study. Int J Environ Res Public Health. 2021 Feb 2;18(4):1–16. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph18168831&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F17%2F2022.09.16.22280020.atom) 2. 2.Zhu Z, Xu S, Wang H, Liu Z, Wu J, Li G, et al. COVID-19 in Wuhan: Immediate Psychological Impact on 5062 Health Workers. medRxiv [Internet]. 2020;2020.02.20.20025338. Available from: [https://doi.org/10.1101/2020.02.20.20025338](https://doi.org/10.1101/2020.02.20.20025338) 3. 3.Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. New England Journal of Medicine. 2020 Feb 20;382(8):727–33. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F17%2F2022.09.16.22280020.atom) 4. 4.Giandhari J, Pillay S, Wilkinson E, Tegally H, Sinayskiy I, Schuld M, et al. Early transmission of SARS-CoV-2 in South Africa: An epidemiological and phylogenetic report. International Journal of Infectious Diseases. 2021 Feb 1;103:234–41. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.11.128&link_type=DOI) 5. 5.Jassat W, Abdool Karim SS, Mudara C, Welch R, Ozougwu L, Groome MJ, et al. Clinical severity of COVID-19 in patients admitted to hospital during the omicron wave in South Africa: a retrospective observational study. Lancet Glob Health. 2022 Jul 1;10(7):e961–9. 6. 6.Jassat W, Ozougwu L, Munshi S, Mudara C, Vika C, Arendse T, et al. The intersection of age, sex, race and socio-economic status in COVID-19 hospital admissions and deaths in South Africa. S Afr J Sci. 2022;18(5–6). 7. 7.Jassat W, Abdool Karim SS, Ozougwu L, Welch R, Mudara C, Masha M, et al. TRENDS IN CASES, HOSPITALISATION AND MORTALITY RELATED TO THE OMICRON BA.4/BA.5 SUB-VARIANTS IN SOUTH AFRICA. Available from: [https://doi.org/10.1101/2022.08.24.22279197](https://doi.org/10.1101/2022.08.24.22279197) 8. 8.Jassat W, Mudara C, Ozougwu L, Tempia S, Blumberg L, Davies MA, et al. Difference in mortality among individuals admitted to hospital with COVID-19 during the first and second waves in South Africa: a cohort study. Lancet Glob Health. 2021 Sep 1;9(9):e1216–25. 9. 9.Davies MA, Morden E, Rosseau P, Arendse J, Bam JL, Boloko L, et al. Outcomes of laboratory-confirmed SARS-CoV-2 infection during resurgence driven by Omicron lineages BA.4 and BA.5 compared with previous waves in the Western Cape Province, South Africa. medRxiv [Internet]. 2022 Jul 1; Available from: [http://www.ncbi.nlm.nih.gov/pubmed/35794899](http://www.ncbi.nlm.nih.gov/pubmed/35794899) 10. 10.Madhi SA, Kwatra G, Myers JE, Jassat W, Dhar N, Mukendi CK, et al. South African Population Immunity and Severe Covid-19 with Omicron Variant. medRxiv [Internet]. 2022 Jan 1;2021.12.20.21268096. Available from: [http://medrxiv.org/content/early/2022/01/27/2021.12.20.21268096.abstract](http://medrxiv.org/content/early/2022/01/27/2021.12.20.21268096.abstract) 11. 11.Nguyen LH, Drew DA, Graham MS, Joshi AD, Guo CG, Ma W, et al. Risk of COVID-19 among front-line health-care workers and the general community: a prospective cohort study. Lancet Public Health. 2020 Sep 1;5(9):e475–83. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-2667(20)30164-X&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F17%2F2022.09.16.22280020.atom) 12. 12.Sasson I. Age and COVID-19 mortality. Source: Demographic Research [Internet]. [cited 2022 Sep 13];44:379–96. Available from: [https://www.jstor.org/stable/10.2307/27032918](https://www.jstor.org/stable/10.2307/27032918) 13. 13.Sze S, Pan D, Nevill CR, Gray LJ, Martin CA, Nazareth J, et al. Ethnicity and clinical outcomes in COVID-19: A systematic review and meta-analysis. EClinicalMedicine. 2020 Dec 1;29–30. 14. 14.Bor J, Gage A, Onoya D, Maskew M, Tripodis Y, Fox MP, et al. Variation in HIV care and treatment outcomes by facility in South Africa, 2011–2015: A cohort study. PLoS Med. 2021 Mar 31;18(3). 15. 15.The World Bank. Inequality In Southern Africa: An Assessment of the Southern African Customs Union [Internet]. 2022 [cited 2022 Sep 15]. Available from: [https://documents1.worldbank.org/curated/en/099125303072236903/pdf/P1649270c02a1f06b0a3ae02e57eadd7a82.pdf](https://documents1.worldbank.org/curated/en/099125303072236903/pdf/P1649270c02a1f06b0a3ae02e57eadd7a82.pdf) 16. 16.Sannigrahi S, Pilla F, Basu B, Basu AS, Molter A. Examining the association between socio-demographic composition and COVID-19 fatalities in the European region using spatial regression approach. Sustain Cities Soc. 2020 Nov 1;62. 17. 17.Sartorius B, Lawson AB, Pullan RL. Modelling and predicting the spatio-temporal spread of COVID-19, associated deaths and impact of key risk factors in England. Sci Rep. 2021 Dec 1;11(1). 18. 18.Quiliche R, Rentería R, de Brito Junior I, Luna A, Chong M. Using Spatial Patterns of COVID-19 to Build a Framework for Economic Reactivation. 2021; Available from: [www.preprints.org](http://www.preprints.org) 19. 19.Santana Santos V, Santos Siqueira T, Cubas Atienzar AI, Augusta Ricardo da Rocha Santos M, Cristina Fontes Vieira S, de Siqueira Alves Lopes A, et al. Spatial clusters, social determinants of health and risk of COVID-19 mortality in Brazilian children and adolescents: A nationwide population-based ecological study. 2022; Available from: [https://doi.org/10.1016/j](https://doi.org/10.1016/j). 20. 20.de Souza APG, de Miranda Mota CM, Rosa AGF, de Figueiredo CJJ, Candeias ALB. A spatial-temporal analysis at the early stages of the COVID-19 pandemic and its determinants: The case of Recife neighborhoods, Brazil. PLoS One. 2022 May 1;17(5 May). 21. 21.Haider MS, Salih SK, Hassan S, Taniwall NJ, Moazzam MFU, Lee BG. Spatial distribution and mapping of COVID-19 pandemic in Afghanistan using GIS technique. SN Social Sciences. 2022 May;2(5). 22. 22.Brezger A, Kneib T, Lang S. BayesX: Analyzing Bayesian Structured Additive Regression Models [Internet]. Vol. 14, JSS Journal of Statistical Software. 2005. Available from: [http://www.jstatsoft.org/](http://www.jstatsoft.org/) 23. 23.Kazembe LN, Chirwa TF, Simbeye JS, Namangale JJ. Applications of Bayesian approach in modelling risk of malaria-related hospital mortality. BMC Med Res Methodol. 2008;8. 24. 24.1. Paraskevis D, editor Ngesa O, Mwambi H, Achia T. Bayesian Spatial Semi-Parametric Modeling of HIV Variation in Kenya. Paraskevis D, editor. PLoS One [Internet]. 2014 Jul 25 [cited 2022 Aug 16];9(7):e103299. Available from: [https://dx.plos.org/10.1371/journal.pone.0103299](https://dx.plos.org/10.1371/journal.pone.0103299) 25. 25.Umlauf N, Adler D, Kneib T, Lang S, Zeileis A. Structured additive regression models: An R interface to Bayes X. J Stat Softw. 2015;63(21):1–46. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v063.i18&link_type=DOI) 26. 26.Goldstein JR, Lee RD. Demographic perspectives on the mortality of COVID-19 and other epidemics. Proc Natl Acad Sci U S A. 2020 Sep 8;117(36):22035–41. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzM2LzIyMDM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDkvMTcvMjAyMi4wOS4xNi4yMjI4MDAyMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 27. 27.Yanez ND, Weiss NS, Romand JA, Treggiari MM. COVID-19 mortality risk for older men and women. BMC Public Health. 2020 Dec 1;20(1). 28. 28.Ho FK, Petermann-Rocha F, Gray SR, Jani BD, Vittal Katikireddi S, Niedzwiedz CL, et al. Is older age associated with COVID-19 mortality in the absence of other risk factors? General population cohort study of 470,034 participants. PLoS One. 2020 Nov 1;15(11 November). 29. 29.1. S. 2. Karger AG Biswas M, Rahaman S, Biswas TK, Haque Z, Ibrahim B. Association of Sex, Age, and Comorbidities with Mortality in COVID-19 Patients: A Systematic Review and Meta-Analysis. Vol. 64, Intervirology. S. Karger AG; 2021. p. 36–47. 30. 30.Arashi M, Bekker A, Salehi M, Millard S, Botha T, Golpaygani M. Evaluating prediction of COVID-19 at provincial level of South Africa: a statistical perspective. Environmental Science and Pollution Research. 2022 Mar 1;29(15):21289–302. 31. 31.Arashi M, Bekker A, Salehi M, Millard S, Erasmus B, Cronje T, et al. Spatial analysis and prediction of COVID-19 spread in South Africa after lockdown. 32. 32.Harris B, Goudge J, Ataguba JE, McIntyre D, Nxumalo N, Jikwana S, et al. Inequities in access to health care in South Africa. J Public Health Policy. 2011;32(SUPPL. 1). 33. 33.Jassat W, Abdool Karim SS, Mudara C, Welch R, Ozougwu L, Groome MJ, et al. Clinical severity of COVID-19 in patients admitted to hospital during the omicron wave in South Africa: a retrospective observational study. Lancet Glob Health. 2022 Jul 1;10(7):e961–9. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/graphic-3.gif