@article {Doblhammer2020.12.22.20248386, author = {Gabriele Doblhammer and Constantin Reinke and Daniel Kreft}, title = {Social disparities in the first wave of COVID-19 infections in Germany: A county-scale explainable machine learning approach}, elocation-id = {2020.12.22.20248386}, year = {2020}, doi = {10.1101/2020.12.22.20248386}, publisher = {Cold Spring Harbor Laboratory Press}, abstract = {Background Little is known about factors correlated with this geographic spread of the first wave of COVID-19 infections in Germany. Given the lack of individual-level socioeconomic information on COVID-19 cases, we resorted to an ecological study design, exploring regional correlates of COVID-19 diagnoses.Data and Method We used data from the Robert-Koch-Institute on COVID-19 diagnoses by sex, age (age groups: 0-4, 5-14, 15-34, 35-59, 60-79, 80+), county (NUTS3 region) differentiating five periods (initial phase: through 15 March; 1st lockdown period: 16 March to 31 March; 2nd lockdown period: from 1 April to 15 April; easing period: 16 April to 30 April; post-lockdown period: 1 May through 23 July). For each period we calculated age-standardized incidence of COVID-19 diagnoses on the county level, using the German age distribution from the year 2018. We characterized the regions by macro variables in nine domains: {\textquotedblleft}Demography{\textquotedblright}, {\textquotedblleft}Employment{\textquotedblright}, {\textquotedblleft}Politics, religion, and education{\textquotedblright}, {\textquotedblleft}Income{\textquotedblright}, {\textquotedblleft}Settlement structure and environment{\textquotedblright}, {\textquotedblleft}Health care{\textquotedblright}, {\textquotedblleft}(structural) Poverty{\textquotedblright}, {\textquotedblleft}Interrelationship with other regions{\textquotedblright}, and {\textquotedblleft}Geography{\textquotedblright}. We trained gradient boosting models to predict the age-standardized incidence rates with the macro structures of the counties, and used SHAP values to characterize the 20 most prominent features in terms of negative/positive correlations with the outcome variable.Results The change in the age-standardized incidence rates over time is reflected in the changing importance of features as indicated by the mean SHAP values for the five periods. The first COVID-19 wave started as a disease in wealthy rural counties in southern Germany, and ventured into poorer urban and agricultural counties during the course of the first wave. The negative social gradient became more pronounced from the 2nd lockdown period onwards, when wealthy counties appeared to be better protected. Population density per se does not appear to be a risk factor, and only in the post-lockdown period did connectedness become an important regional characteristic correlated with higher infections. Features related to economic and educational characteristics of the young population in a county played an important role at the beginning of the pandemic up to the 2nd lockdown phase, as did features related to the population living in nursing homes; those related to international migration and a large proportion of foreigners living in a county became important in the post-lockdown period.Discussion In the absence of individual level data, explainable machine learning methods based on regional data may help to better understand the changing nature of the drivers of the pandemic. High mobility of high SES groups may drive the pandemic at the beginning of waves, while mitigation measures and beliefs about the seriousness of the pandemic as well as the compliance with mitigation measures put lower SES groups at higher risks later on.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThis research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:This study bases on publicly available anonymised aggregated population data; hence, this study complied with the tenets of the Declaration of Helsinki and required no ethical approval.All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesAll data referred to in the manuscript are publicly available. https://npgeo-corona-npgeo-de.hub.arcgis.com/datasets/dd4580c810204019a7b8eb3e0b329dd6_0 https://www.regionalstatistik.de/genesis/online?operation=previous\&levelindex=0\&step=0\&titel=Themen+\%2F+Statistiken\&levelid=1608606041212\&acceptscookies=false$\#$abreadcrumb https://www.inkar.de/ https://www.umweltbundesamt.de/en/data https://statistik.arbeitsagentur.de/DE/Navigation/Statistiken/Interaktive-Angebote/Pendleratlas/Pendleratlas-Nav.html https://www.openstreetmap.de/karte.html}, URL = {https://www.medrxiv.org/content/early/2020/12/22/2020.12.22.20248386}, eprint = {https://www.medrxiv.org/content/early/2020/12/22/2020.12.22.20248386.full.pdf}, journal = {medRxiv} }