Evidence of COVID-19 fatalities in Swedish neighborhoods from a full population study

The COVID-19 pandemic has highlighted a debate about whether marginalized communities suffered the disproportionate brunt of the pandemic’s mortality. Empirical studies addressing this question typically suffer from statistical uncertainties and potential biases associated with uneven and incomplete reporting. We use geo-coded micro-level data for the entire population of Sweden to analyze how local neighborhood characteristics affect the likelihood of dying with COVID-19 at individual level, given the individual’s overall risk of death. We control for several individual and regional characteristics to compare the results in specific communities to overall death patterns in Sweden during 2020. When accounting for the probability to die of any cause, we find that individuals residing in socioeconomically disadvantaged neighborhoods were not more likely to die with COVID-19 than individuals residing elsewhere. Importantly, we do find that individuals show a generally higher probability of death in these neighborhoods. Nevertheless, ethnicity is an important explanatory factor for COVID-19 deaths for foreign-born individuals, especially from East Africa, who are more likely to pass away regardless of residential neighborhood.

The ongoing COVID-19 pandemic has reminded us that the features that endow cities with distinctive socioeconomic advantages, such as their vast socioeconomic networks and their density, also make them vulnerable to the spread of pathogens, especially emerging infectious diseases [1][2][3] .It is by now well known in the social sciences that the effects of urban environments on individual-level outcomes are mediated, to a large extent, by local residential contexts, i.e., neighborhoods [4][5][6][7][8][9][10][11] .There were concerns from the very start of the COVID-19 pandemic that the manner and pace of contagion and the impacts of the disease would be greatly magnified by neighborhood characteristics, with neighborhood-level inequalities expressed as concentrated poverty, segregation, and deprivation amplifying the virus's infectivity and deadly consequences.These expectations support a set of hypotheses that we test here in the context of Sweden, using the country's full population records of mortality.The aim of the paper is thus to analyze whether individuals residing in ethnically segregated and socioeconomically weak neighborhoods were more likely to die with COVID-19.We specifically contribute to existing literature by conditioning COVID-19 deaths on overall rates of mortality.Given our access to geocoded individual micro-level data, we are thus able to trace out whether individuals in disadvantaged neighborhoods had a higher-than-expected risk of dying due to the pandemic.
Expectations for neighborhood effects resulting in disproportionate impacts of the COVID-19 pandemic on disadvantaged local communities follow from decades of observational and empirical work in sociology, economics, and public health sciences 12,13 .This literature gained a renewed focus since the 1980s with the work of Wilson, who characterized processes of local and cumulative disadvantage in deindustrializing US cities, which created mass unemployment among working class Black communities and facilitated the formation of inner-city marginalized neighborhoods 14 .Analyses of neighborhood effects have grown substantially in recent years supporting a multidisciplinary account of the origins and consequences of spatially concentrated local disadvantage [15][16][17][18] .More recent empirical work, using detailed data sources, such as tax records, has better established the critical importance of temporal exposure to disadvantaged neighborhoods, especially during childhood, with life course consequences affecting future income, family structure, health, and life expectancy [19][20][21] .
The empirical richness and methodological diversity which characterize the investigation of neighborhood effects studies suggest complicated chains of causality underlying local disadvantage involving interlocked social, Insurance Agency.The database is updated on a yearly basis and applies to a whole year as the temporal measurement period, using the population on the 31st of December.The variables in the database are grouped by demographic characteristics (age, gender), education and training attainment, employment status, income and to what extent an individual has been financially compensated via the social insurance system, family characteristics (e.g., civil status, number of children in the household), and workplace establishment and firm information (e.g., industry, number of employees, occupational structures).Information is also available about individuals' residential housing arrangement such as the number of rooms in their dwelling.A detailed description of the variables utilized in this study is available in the supplementary material (Tables S1-S3).
To the already existing full-population data set, we have linked additional data from the "Cause of Death Register" constructed by the Swedish National Board of Health and Welfare, covering all deaths in the country during 2020.This dataset includes information on causes of death (including multiple causes), date of death as well as information on diabetes diagnosis.Diabetes is important to control for as it is highlighted as a risk factor for COVID-19 deaths.Since individuals with underlying diseases face a higher risk of severe infection and death in COVID-19 and are more likely to reside in disadvantaged areas 51,52 , we control for multiple causes of death to not overestimate the impact of neighborhood characteristics.In total, 97,374 individuals passed away in Sweden during 2020, of which 97,197 can be matched to the 2019 micro-data.We matched against the year 2019 since the individual characteristics register data from Statistics Sweden is based on the last day of the year and thus the individuals who died would not be registered during the year 2020.Of those, 10,168 individuals had COVID-19 as a cause of death, 4.8 percent had COVID-19 as the sole cause of death, 90 percent of the dead with COVID-19 during 2020 were aged 70 and above, and 53 percent were male.To get a first glimpse of the spatial variation of individuals dying with COVID-19 in Sweden, Fig. 1 shows the distribution of COVID-19 deaths across municipalities, relative to the overall number of deaths in each municipality.COVID-19 cases are mainly overrepresented in the Stockholm region, as shown by the darker red colors in this area.
The European Union statistical agency reported that Sweden had an excess mortality rate of 7.7 percent in 2020 in a comparison with the previous 4 years 53 .Of the 30 countries for which data are available, 21 had a higher excess mortality rate than Sweden.However, Sweden had significantly higher excess mortality than its Nordic neighbors Denmark (1.5 per cent) and Finland (1.0 per cent).
In the estimations of COVID-19 deaths, we control for potential virus exposure in each municipality by counting the number of COVID-19-infected in the municipality and region (weighted by distance to the municipality in question) in the week preceding the individual's death.As COVID-19 hit different municipalities at different times during the year, this is a more relevant measure than simply controlling for time (for example, weeks or months) in general.
Data processing and all statistical analysis, including the generation of Figs. 4, 5 and 6 are performed using STATA, version 18 54 .Maps are created using ArcGIS software by Esri, version 10.8.1 55 .). Vol:.(1234567890)

Spatial units of analysis
To analyze how the neighborhood of residence affected the likelihood of dying with COVID-19, we take advantage of the geographic information in the data from Statistics Sweden to construct neighborhood characteristics based on ethnic background, education level, and incomes.There is no unique neighborhood definition in Sweden, but several definitions have been developed during different time periods.The SAMS definition (Small Areas for Market Statistics) was introduced in 1994 and has often been used in Swedish studies of neighborhood effects 56 .However, this definition has been criticized because these neighborhood areas are not as homogeneous as it is sometimes assumed.In 2018, SAMS was replaced by a new neighborhood definition known as DeSo (Demographic Statistical Areas) that took into account changes in buildings and infrastructure, since SAMS was developed, which means that the SAMS definition was no longer considered to meet the existing needs.We define neighborhoods according to Statistics Sweden's recent (2020) classification of Regional Statistical Areas (RegSO), which builds on the DeSo definition and delineates 3,363 units of analysis covering the whole nation.This geographic specification is especially advantageous for studying neighborhood effects since the purpose of RegSO is to allow for statistical studies that assess socioeconomic segregation 67 .The final RegSo classification was introduced in 2020 and is intended to remain unchanged over time.A benefit of this classification of neighborhoods is that location-specific demographic conditions are defined with spatial (neighborhood) borders following, e.g., streets, waterways, and railways.RegSO areas thus represent actual neighborhoods rather than arbitrary administrative constructs.
Based on the residential geo-coded information, we can construct the socio-economic characteristics of the neighborhood in which an individual resides.We measure neighborhood characteristics in terms of the share of residents in the neighborhood who (1) were born in a non-Nordic/EU15 country (population aged 15 +), (2)  had only an elementary school degree as highest educational attainment (population aged 20-64, thus excluding individuals in school age and the older generation who were overall less likely to attend high school and higher education levels), and (3) are at the risk of poverty, defined as having a disposable income below 60 percent of the corresponding national median (population aged 20 + , thus excluding individuals in school age) 57 .Figure 2 illustrates the neighborhood characteristics focusing on the areas in and around Stockholm. Figure 3 shows the overall deaths and COVID-19 deaths for the same geographic area.Further descriptive statistics for the neighborhood variables are available in the supplementary material (Fig. S1-S7, Table S4).

Estimation of the probability of dying with COVID-19
Differences in individuals' propensity to die with COVID-19 across neighborhoods are hypothesized to be caused by different socio-economic conditions.At the same time, recent studies point to that also general death rates vary across neighborhoods, resulting in lower life expectancies in disadvantaged neighborhoods 58,59 .There may thus be an issue of selection bias when estimating neighborhood effects on the likelihood of dying with COVID-19 coming from the fact that individuals residing in certain types of neighborhoods are more likely to die overall.Hence, we estimate neighborhood effects on the probability of dying with COVID-19 conditional on the probability of dying in general, using the Heckman selection model [60][61][62] .Even though this may not be the typical selection bias problem identified in previous literature, the Heckman two-step approach is a feasible and general-enough method to deal with selection bias.Originally proposed to correct for non-randomly selected data samples, the technique explicitly models the selection probability of each observation (the selection equation) together with the conditional expectation of the dependent variable (the outcome equation).
That selection bias does indeed exist in our case is supported by the estimations showing a rho, ρ, significantly different from zero and a rejection of the hypothesis of independent equations (using the Wald test).Additionally, running probit estimations on COVID-19 deaths without selection bias correction produce vastly different estimates of several variables (Table S6).Of particular importance for this paper is that ordinary probit estimations tend to overestimate the importance of neighborhood effects on COVID-19 deaths, confirming the need for the two-step procedure.
In the first step of the Heckman selection model, the probability of dying for a population on a given year is estimated.The canonical specification for this relationship is a probit regression of the form: where s * is an unobserved latent variable, Z is a vector of explanatory variables displayed in Table S1, γ is a vector of unknown parameters (regressors), and u 1 is standard normally distributed error term.The binary outcome variable, s = 0, 1(meaning death or no death over the time period), are observed if The probability of being selected (that is, dying) is specified by where is the cumulative distribution function of the standard normal distribution.The model estimation, in terms of the coefficients γ , yields results that can be used to predict the death probability for each individual.
In the second step of estimation, we correct for "self-selection" (the fact that COVID-19 deaths are only observed in "the population of dead") by incorporating a transformation of these predicted individual probabilities as an additional explanatory 61 .Analogously to the first step, the probability of dying with COVID-19 may be specified using another latent variable (1) where X is a vector of explanatory variables (Table S1) and where we can observe the binary outcome y = 1, 0 (specific death with COVID-19), given that an individual has died such that The assumptions of the error terms are u 1 ∼ N(0, 1) , u 2 ∼ N(0, 1) and corr(u 1 , u 2 ) = ρ � = 0. Thus, we get the conditional expectation of dying with COVID-19 according to Under the error terms assumptions, we have that where ρ is the correlation between unobserved determinants of the probability of dying u 1 , and unobserved determinants of the probability of dying with COVID-19, u 2 , and is the inverse Mills ratio, evaluated at Zγ .In other words, sample selection is viewed as a form of omitted-variable bias, as conditional on both X and on λ as if the sample is randomly selected.(As usual, φ denotes the standard normal density function.)(4) We estimate the Heckman selection model on the total population of Sweden aged 15 or above as well as the population aged 70 and above since the latter comprise the bulk of COVID-19 deaths.Additionally, we perform separate estimations for different regional types, by partly identifying municipalities belonging to the Stockholm region and partly distinguishing between metropolitan, urban, and rural municipalities following the classification by the Swedish Agency for Economic and Regional Growth 63 .

Neighborhoods and COVID-19 mortality
We estimated the probability of dying with COVID-19 conditional on the probability of dying in general using a two-step Heckman Probit regression model [60][61][62] .In the estimation, we examine to what extent neighborhood characteristics, in terms of ethnicity, education level, and income level, significantly explain the probability of dying with COVID-19 (step 2) conditional the probability of dying in general (step 1), while still controlling for other individual and regional characteristics.The following results Tables 1, 2, 3, 4, 5 and 6 are organized based on neighborhood characteristics, in terms of the share of residents in the neighborhood who are (1) foreign born, (2) low educated, and (3) at the risk of poverty.All coefficients in the tables are expressed as average marginal effects (AMEs), in other words the average change in probability of dying (step 1) and dying with COVID-19 (step 2) when the explanatory variable in question changes by one unit.For categorical variables, the AMEs show the change in probability in relation to the base category.
Table 1 presents the estimation results of various socioeconomic characteristics on the general likelihood of death (all causes), controlling for individual and regional characteristics.
The results on the neighborhood variables are in line with the expectations, showing that individuals residing in more ethnically segregated and socioeconomically weaker neighborhoods have a greater probability of dying.Starting with ethnicity, neighborhoods with higher shares of foreign born (> 75%) have a marginal effect of 0.0052, which is about double the marginal effect of all other ethnic share categories.The marginal effect on general death rates is insignificant for neighborhoods with the lowest shares of foreign born but increases to 0.0004 and 0.0009 for neighborhoods with 15-35 percent foreign born, increasing further with higher shares of foreign born.We find similar patterns for the influence of education, where neighborhoods with higher shares of low-educated individuals display a higher probability of death from all causes.The marginal effect on mortality for the group with the highest shares of low educated is 0.0043, compared to 0.0003 and 0.0007 for the groups with the lowest (better educated).For poverty, the patterns are similar: The higher the share of individuals in poverty, the higher the marginal effects on the general probability of dying.However, for the category with the highest shares of the poor (> 55%), we observe a decline in the death probability because this category features low-income student areas in college towns (this is discussed further below in regards to COVID-19 fatalities).
The marginal effects for disadvantaged neighborhoods reported in Table 1 may seem relatively small.For example, individuals residing in the neighborhoods with the highest shares of low educated (> 35%) are 0.43  www.nature.com/scientificreports/percentage points more likely to die than individuals residing in the neighborhoods with the lowest shares of low educated (0-7%).To illustrate that this difference is not only statistically significant but also non-negligible from a socioeconomic perspective, Figs. 4, 5 and 6 show the predicted probability of dying for each socioeconomic category and neighborhood group.These figures illustrate how the probability to die increases as a function of the share of vulnerable individuals, for example, Fig. 4 shows that the probability to die is almost 50 percent higher in the most ethnically segregated neighborhoods (> 35% non-Nordic/EU15 immigrants), as compared to the least ethnically segregated neighborhoods (0-7% non-Nordic/EU15 immigrants).The predictive probabilities shown in Figs. 4, 5 and 6 are generated from the average marginal effects estimated for Table 1.
Table 2 shows the marginal effects of each of these local environment characteristics on the likelihood of COVID-19 death for an individual in each type of neighborhood.Note that the size of the estimates in Table 2 should not be compared to the estimates in Table 1 since the baseline probability of dying with COVID-19, given death, is greater than the overall probability of dying.
Surprisingly perhaps, we see that the general effects of distinct neighborhood characteristics -whether ethnicity, education attainment or incomes levels-on individual risk of dying with COVID-19 are statistically insignificant in 17 out of 22 cases.In other words, once we control for individual and locational characteristics associated with higher likelihood of death from all causes, we find little or no support that any neighborhood characteristics moderate the proportion of dead dying with COVID-19 specifically.The richness and completeness of the data for Sweden allows us to discuss these effects in great detail.
We find a positive and significant result for COVID-19 mortality for neighborhoods with lower shares of foreign born (> 7-15% and > 15-25% in Table 1), and also some support for the educational level effect as individuals in more educated neighborhoods (with up to 10-25% share of less educated) were less likely to die.Most importantly, we find no significant results for an increased risk of COVID-19 deaths in the most extreme communities in terms of education or ethnicity (> 7-10% nor > 35%), so all effects are found at more intermediate levels of mixing.These more extreme communities comprise of the most marginalized neighborhoods or the neighborhoods that are most affluent.In terms of income, we find a higher marginal effect on the COVID-19 share of mortality for the highest poverty neighborhoods (> 55%).Further analysis shows that several of these neighborhoods are student areas in college towns (that is, neighborhoods where many university students reside) with a low share of non-Nordic/EU immigrants and a low share of the less educated.These neighborhoods are not considered marginalized, which explains their lower and insignificant marginal effect on the probability for overall death (selection equation).On the other hand, there have been reports of COVID-19 outbreaks among specific student groups (for example, in the university town of Uppsala), which implies that students may have contributed to a greater spread of infection into the local general population.Indeed, the results of the outcome equation show that individuals residing in student areas have a 3% higher probability of dying with COVID-19 than from other causes.These deaths are not infected students, but rather the more vulnerable and/or older individuals living in the same areas, as we discuss below in terms of age profiles.
At the individual level (see "Individual characteristics"), it is interesting to note the negative significant results for older age groups relative to dying with COVID-19 (Table 2).While this may sound counterintuitive, it is simply the result of the selection bias that we correct for in Step 1 of the analysis (Table 1).Older persons are more likely to die to start with (general causes) but, given death at an advanced age, are actually less likely to have COVID-19 as their cause of death.For individuals who are foreign-born, we find an opposite pattern.Individuals with a non-Swedish ethnic background are in most cases less likely to die from general causes (Table 1), but among those who died, the cause of death is more likely to be COVID-19 (Table 2).The ethnicity most at risk in our analysis is East African, showing the highest marginal effect: approximately 7 times higher than for individuals from EU15.These results are in line with previous studies on COVID-19 deaths across ethnic groups.
As expected, we also find a higher probability of dying with COVID-19 the more deleterious health conditions an individual has.This confirms clinical observations that people with multiple serious conditions are at a greater risk of becoming fatally ill from a COVID-19 infection.There is a significant marginal effect of up to 13 percentage points for individuals with five or more causes of death.On the other hand, after controlling for all other factors, individuals diagnosed with diabetes are less likely to die with COVID-19, but the effect is small since the decrease in the probability of dying with COVID-19 is less than one percentage point.Since diabetes has been brought up as a risk factor in COVID-19 deaths, our results may appear counterintuitive.They may however be explained by the fact that diabetes patients often suffer from several potentially life-threatening conditions 64  www.nature.com/scientificreports/and the risk of diabetes may thus be captured by the effect of multiple causes of death.Indeed, when excluding multiple causes of death diabetes turns positive although insignificant.Regarding regional characteristics, individuals residing in Stockholm (the Swedish capital) and those exposed to a greater number of COVID-19 infected in their municipality have a higher probability of having COVID-19 as a cause of death.On the other hand, the results on population density are significantly negative, which goes against the expectations and is likely due to collinearity with the variable on virus exposure.Indeed, when excluding virus exposure population density turns significantly positive, despite a less than 0.5 bivariate correlation between the variables.The remaining variables are robust to the exclusion of virus exposure.
Taken together, we find little support for neighborhood effects mediated by socioeconomic segregation on the likelihood of dying with COVID-19, but strong support that individual-level characteristics, especially related to specific ethnic backgrounds, may have played a role.

Statistical analysis by age groups and urban areas
To further investigate other important effects, we focus our estimation procedures on more specific subgroups since (1) the majority of deaths are in the age group 70 and older, and (2) the role of neighborhoods may matter differently in urban versus rural settings.Our objective with this more disaggregated analysis is to investigate possible confounding results due to averaging.For example, it may be the case that the selection bias estimation (step 1) generates a different result if we target specific groups instead of comparing the probability of dying for the whole population.Starting with an analysis where we only include individuals aged 70 and over, Tables 3  and 4 present the average marginal effects for neighborhood characteristics (full tables are available from the authors upon request).As noted, 90% of deaths with COVID-19 in Sweden belong to this age group.The number of observations decreases relatively more in the selection equations (compare Table 3 with Table 1) than in the outcome equations (compare Table 4 with Table 2), reflecting the higher probability of dying in the older age group (70 +) than in the overall population (15 +).
We find that neighborhood characteristics did not add any significant risk to the likelihood of dying with COVID-19 for the oldest age groups and, specifically, that neither ethnically nor educationally segregated neighborhoods add to this risk.However, we do find a weakly significant (10 percent level) increased risk for elderly individuals who live in neighborhoods where more than 55% of the population is categorized as poor.Recall that these are primarily student residential areas in college towns.The elderly may have been adversely affected by what was reported in the media as an increased spread of COVID-19 resulting from student gatherings and lack of compliance with social distancing recommendations in several of these towns during the pandemic.While the students themselves may not have been at any severe risk of death, we find that the elderly residing in these neighborhoods are more likely to die specifically with COVID-19.
We also performed the estimation analysis separately for urban and rural regions since neighborhood effects may play out differently depending on this type of context.For example, recent work indicates that bigger, more urban regions also tend to be more segregated 65,66 .We also ran a separate analysis for Stockholm, Sweden's capital city, as well as for the metropolitan regions of Stockholm, Gothenburg, and Malmö combined (metropolitan regions being the equivalent of spatially extended labor markets).Table 5 and 6 present the average marginal effects for the neighborhood segregation variables for metropolitan municipalities, urban municipalities, and rural municipalities, respectively.The separate estimation for individuals residing in the Stockholm region includes 19 municipalities classified as metropolitan plus six urban municipalities and one rural municipality.
Table 5 summarizes these results showing that, regardless of regional type, we observe that more segregated neighborhoods in terms of ethnicity, education levels, and poverty increase the marginal effects of the probability of dying from all causes.The exception is rural regions, with fewer significant but still positive coefficients.As per Table 6, living in marginalized neighborhoods, in this sense, does not add significantly to the risk of dying specifically with COVID-19.There are some significant estimates however, in particular for Stockholm

Discussion and conclusions
The COVID-19 pandemic brought to the forefront of both public discussions and scientific research inequalities in the burden of disease and mortality, especially affecting marginalized local communities.Much research reported so far, found excessive deaths in ethnically segregated and socioeconomically disadvantaged local communities 67,68 .There are also several previous empirical studies connecting COVID-19 deaths to individual and neighborhood characteristics.However, these studies have not adjusted for the underlying selection bias resulting from spatial variations in general mortality.
The most ambitious studies have looked at the likelihood of dying given infection 69 , but not if this is a group that is likely to die to start with.A typical approach has been to examine "excess mortality" but at the national level, not at the individual level 70 .Here, we used a full-population, geocoded micro-dataset for Sweden to analyze in detail if and how the characteristics of residential neighborhoods affected the likelihood of dying with COVID-19 after controlling for a rich set of individual characteristics and comparing with background death rates for the same groups in Sweden during 2020.The scope of our datasets and the design of our statistical approach allow us therefore to account explicitly for selection bias reflecting the probability of overall death by all causes, which proved critically important in our findings.At the same time, our two-step modelling approach, which is unique in studies on COVID-19 deaths, also comes with the limitation that our estimates are not directly comparable to previous studies on this topic.Nevertheless, we argue that to understand and to not overestimate the role played by residential contexts during the pandemic, it is necessary to account for individual and geographic heterogeneity in the general risk of dying.
Our results show that neighborhood characteristics played at best an extremely marginal role in explaining the probability of dying with COVID-19 in Sweden.Neither ethnicity nor lower levels of education at the neighborhood scale add to the likelihood of dying with COVID-19.On the other hand, we do observe an increased risk of COVID-19 deaths in neighborhoods with a high proportion of residents categorized as "poor."On closer inspection, however, we found that these neighborhoods are mainly student residential areas in college towns, not classically considered as marginalized communities.In these local contexts, we show that high COVID-19 death rates are the result of specific population mixing patterns, where a less compliant young population in terms There may be several reasons why ethnicity influenced the probability of dying during the pandemic in the context of how Swedish authorities responded to the pandemic.Instead of relying on a mandatory lockdown, the Swedish COVID-19 strategy was based on the assumption that individuals would follow the recommendations of public health authorities 71 .International surveys, such as Inglehart's World Value Survey, show that residents of Sweden generally have high levels of trust in their authorities and in researchers.This trust in government experts in turn was taken to imply trusting public health advice [72][73][74] .Studies of migrant groups in Sweden, however, indicate that specific groups can have significantly lower levels of trust in public authorities 75 .As a result, information about what to do in response to the pandemic might have been obtained from sources in the migrants' countries of origin or through informal channels within their own networks.Swedish medical doctors raised concerns that this strategy would have been better suited for a more homogeneous society 76 .Contrary to prevalent views about Sweden, the country is significantly heterogeneous, with 20% of the resident population being foreign-born.Language may also have been a barrier as the initial contact with Swedish healthcare often takes place by telephone and some individuals may have been limited in their ability to describe their medical situation.Further, there could be cultural differences that affect the extent to which one seeks care for different types of symptoms 77 .Studies have also shown how immigrants in Sweden are three to four times more likely than native Swedes to say they suffer from poor or very poor health, often due to factors related to place of origin, the refugee process, as well as the challenge of being new to a country and in a socially vulnerable situation 78 .Foreign born individuals are also less likely to use health screening programs 79,80 .Female foreign-born individuals from some countries are also more likely to suffer from obesity than native Swedish women 81 .
We did find, however, that neighborhood characteristics play a clear and important role in the general likelihood of death from all causes both pre-and during the pandemic.Individuals living in marginalized areas display a significantly higher risk of dying, whatever their individual characteristics.Reconciling these two results means that this elevated likelihood of death is shared among all causes and is not specific to COVID-19.A major conclusion from these results is that the role of neighborhood and individual vulnerability-captured here by the risk of death-is in part "ecological", in that it cannot be fully explained by individual characteristics.
It is important to note that Sweden has an extensive public (i.e., free) healthcare system.It is organized at the provincial level ("regions" in Swedish) and divided into two levels: primary and secondary care.Primary care consists of slightly over 1000 health centers distributed throughout the country.This is normally the first contact made with healthcare.If necessary, the patient is referred to secondary or more specialized tertiary care, which consists mainly of hospitals that cover larger geographical areas.It could be the case that health centers in marginalized neighborhoods were unable to provide care of the same quality as the ones located in better off neighborhoods.It has been documented that health disparities are greater within regions than between them, and that this has a negative effect on the disadvantaged groups in society 82 .Furthermore, the work situation at health centers in marginalized neighborhoods is reported to be more challenging as the need is often greater, but resources are limited 83 .Even though Swedish healthcare is free and of a high quality by international standards, there are still reports of discrimination against certain groups when it comes to treatment and a lower probability of gaining access to certain types of care 84 .
Although the datasets available for Sweden are special in terms of their details and completeness, our statistical procedures and findings should be replicable in other nations with population registries, such as other Scandinavian nations, the Netherlands or Austria, as well as in other contexts where detailed death records are available.The hypothesis advanced by our study for these contexts is that apparent elevated burdens of mortality with COVID-19 fall along the same statistical lines as burdens of general mortality at the individual and community levels, once these specific contexts can be accounted for.We encourage further studies to test this hypothesis beyond the case of Sweden to allow for cross-country comparisons.
A limitation of our study following the official definition of COVID-19 deaths in Sweden is however that an individual passing away within 30 days after a confirmed positive COVID-19 test is classified as a COVID-19 death even though the death may have been primarily due to other causes.Even though we do control for multiple causes of death, this definition may bias our estimates of neighborhood effects on COVID-19 deaths somewhat, most likely downwards.For future studies a more detailed coding of deaths from COVID-19 (and other infectious virus diseases) by the Swedish National Board of Health and Welfare would be beneficial.Analyses of determinants of mortality once vaccination became available in 2021 are also an important theme for future studies.
The results presented here might have generalizable and practical policy implications, specifically in contributing methods to better assess the social determinants of health and emphasizing the importance of holistic approaches to human development.They show that statistical analyses of health inequalities in marginalized neighborhoods conducted in connection with the COVID-19 pandemic must be examined in broad terms, beyond the analysis of a singular emergency.In our study, COVID-19 appears as a new stress creating additional mortality along statistically engrained population risk patterns, expressed as systemic community level dynamics affecting all mortality.These findings suggest therefore that while emergencies such as the COVID-19 pandemic necessarily drive attention and policy over the short term, there might be systemic factors that need be better understood and mitigated over the longer term associated with higher risk of death in marginalized neighborhoods.These factors cannot be explained by an individual's age, income, or education.While the COVID-19 pandemic has highlighted the vulnerability of marginalized neighborhoods at a time of crisis, we are learning that the best defense against the next pandemic may well be the systemic improvement of health conditions in disadvantaged communities at all times.

Figure 2 .
Figure 2. Share of foreign-born (panel a), low educated (panel b), and in risk of poverty (panel c) per neighborhood in Stockholm and surrounding areas (blue color indicates water).

Figure 3 .
Figure 3. Share of dead (dead relative to total population) (panel a) and share of COVID-19 deaths (COVID-19 deaths relative to total deaths) (panel b) per neighborhood in Stockholm and surrounding areas (blue color indicates water).Panel b does not include data for neighborhoods with fewer than five dead in total.

Table 2 .
, Marginal effects for the probability of dying with COVID-19 at individual level, given that the individual has died (outcome equations).Notes: Dependent variable = 1 if dead in COVID-19, zero if other cause/s of death.* denotes significance at 1% level.

Table 3 .
Marginal effects for the probability of dying at individual level, aged 70 and above (selection equations).Notes: Dependent variable = 1 if dead, zero otherwise.* denotes significance at 1% level.
Vol:.(1234567890) Scientific Reports | (2024) 14:2998 | https://doi.org/10.1038/s41598-024-52988-3www.nature.com/scientificreports/and metropolitan municipalities and individuals there residing in neighborhoods with an intermediate share of foreign-born, low educated, and poor.However, the signs are negative, indicating that individuals residing in these neighborhoods are less likely to die with COVID-19 than individuals residing in the least ethnically segregated and socioeconomically strongest neighborhoods.

Table 4 .
Marginal effects for the probability of dying with COVID-19 at individual level, given that the individual has died, aged 70 and above (outcome equations).Notes: Dependent variable = 1 if dead with COVID-19, zero if dead from other cause/s.* denotes significance at 1% level.

Table 5 .
Marginal effects of the neighborhood characteristics for the probability of dying at individual level, per region type (selection equations).Notes: As for Table4above, the estimations including i) Ethnic, ii) Low education, and iii) Poverty (to describe neighborhoods) are run separately.All control variables included.* denotes significance at 1% level.

Table 6 .
Marginal effects of the neighborhood characteristics for the probability of dying with COVID− 19 at individual level, given that the individual has died, per region type (outcome equations).Notes: As for Table3above, the estimations including i) Ethnic, ii) Low education, and iii) Poverty (to describe neighborhoods) are run separately.All control variables included.*denotes significance at 1% level.ofsocial distancing interacts with an older set of individuals, who in turn bear the brunt of mortality.In other words, there is little evidence from our full population study that there are neighborhood effects connected to segregation or disadvantage mediating the probability of dying with COVID-19 in Sweden.Observed variations can be explained by differences of individual characteristics, where, e.g., specific ethnicity plays a significant role.