Skip to main content
  • Research article
  • Open access
  • Published:

Socioeconomic status and survival after stroke – using mediation and sensitivity analyses to assess the effect of stroke severity and unmeasured confounding

Abstract

Background

Although it has been established that low socioeconomic status is linked to increased risk of death after stroke, the mechanisms behind this link are still unclear. In this study we aim to shed light on the relationship between income level and survival after stroke by investigating the extent to which differences in stroke severity account for differences in survival.

Methods

The study was based on patients registered in Riksstroke (the Swedish stroke register) with first time ischemic stroke (n = 51,159) or intracerebral hemorrhage (n = 6777) in 2009–2012. We used causal mediation analysis to decompose the effect of low income on 3-month case fatality into a direct effect and an indirect effect due to stroke severity. Since causal mediation analysis relies on strong assumptions regarding residual confounding of the relationships involved, recently developed methods for sensitivity analysis were used to assess the robustness of the results to unobserved confounding.

Results

After adjustment for observed confounders, patients in the lowest income tertile had a 3.2% (95% CI: 0.9–5.4%) increased absolute risk of 3-month case fatality after intracerebral hemorrhage compared to patients in the two highest tertiles. The corresponding increase for case fatality after ischemic stroke was 1% (0.4–1.5%). The indirect effect of low income, mediated by stroke severity, was 1.8% (0.7–2.9%) for intracerebral hemorrhage and 0.4% (0.2–0.6%) for ischemic stroke. Unobserved confounders affecting the risk of low income, more severe stroke and case fatality in the same directions could explain the indirect effect, but additional adjustment to observed confounders did not alter the conclusions.

Conclusions

This study provides evidence that as much as half of income-related inequalities in stroke case fatality is mediated through differences in stroke severity. Targeting stroke severity could therefore lead to a substantial reduction in inequalities and should be prioritized. Sensitivity analysis suggests that additional adjustment for a confounder of greater impact than age would be required to considerably alter our conclusions.

Peer Review reports

Background

There is a well-documented social gradient in health, where socially underprivileged groups have an increased risk of disease and adverse health outcomes [1,2,3]. This is also true for stroke, the second most common cause of death worldwide according to the World Health Organization (WHO) and a leading cause of disability among adults. Socioeconomic inequalities in stroke incidence, quality of and access to care, and outcome after stroke have been well established globally where low- and middle-income countries, as well as socially disadvantaged groups within high-income countries have been found to be disproportionately affected by stroke [4].

Patients with low socioeconomic status (SES) are at an increased risk of death after stroke, even in countries with universal access to healthcare [5,6,7,8,9,10]. The causes behind these differences have not been established but are likely multifaceted, with aspects related both to patient characteristics (risk factors, comorbidities, etc.) and structural problems within the health care system. Patients with low SES have been found to have more severe strokes and higher incidence of cardiovascular risk factors than high SES patients [4, 6, 11]. Examining how much of the effect of SES on survival that is attributable to differences in disease severity may shed light on the extent to which differences are related to pre-stroke (e.g. comorbidities, lifestyle, primary prevention) vs. post stroke factors such as structural differences in acute stroke management and secondary prevention.

By using causal mediation analysis we can separate the effect of SES on post-stroke survival into a direct effect and an indirect effect operating through stroke severity to gain insight into the relative importance of each pathway [12, 13]. The estimation of these effects from observational data requires strong assumptions regarding residual confounding of the relations involved. These assumptions are often difficult or impossible to verify from the observed data and it is therefore crucial to perform sensitivity analyses to assess the robustness of the results to violations [14,15,16,17,18]. The mediation analysis literature has focused almost exclusively on methods to gauge sensitivity to unobserved confounding of the relation between the intermediate variable (mediator) and the outcome, assuming that randomization of the exposure or adjustment for baseline confounders adequately addresses confounding of the exposure-mediator and exposure-outcome relations [14,15,16, 18, 19]. However, in observational studies no such randomization takes place and it is difficult to guarantee that a sufficiently rich set of baseline confounders has been adjusted for.

The aim of this study is to shed light on the relationship between socioeconomy and survival after stroke by investigating to which extent it can be explained by differences in stroke severity. We focused on short term survival (up to 3 months) after stroke and used individual patient data from Riksstroke, the national Swedish stroke register, in combination with causal mediation analysis. Recent techniques for sensitivity analysis were used to assess the robustness of the results to different types of residual confounding [17].

Methods

Data and variable definitions

Patients included in the study were registered in Riksstroke with intracerebral hemorrhage, ICH (ICD-10 code I61) or ischemic stroke, IS (ICD-10 code I63) between January 1, 2009 and October 1, 2012. Patients with previous stroke, living in institution, or dependent in ADL (activities of daily living, i.e. unable to manage clothing, toileting or walking unassisted) were excluded from the study. A total of 6777 patients with ICH and 51,159 patients with IS were included in the analyses (Fig. 1). We performed separate analyses for the two main types of stroke, IS and ICH, as these patient groups tend to differ on background variables and present with different severity levels and short term mortality rates [20, 21].

Fig. 1
figure 1

Flow-chart for inclusion of patients in the study

The purpose of Riksstroke is to monitor and support improvement of quality of stroke care in Sweden [22]. Patient-level information is collected during the acute phase at all Swedish hospitals that admit acute stroke patients (currently 72 hospitals) as well as at follow-up 3 months and 1 year after stroke. For this study Riksstroke data were linked with data from the Swedish Cause of Death Register (managed by the Swedish National Board of Health and Welfare) and the LISA database (managed by Statistics Sweden) using the personal identity numbers (Swedish national identification numbers) of the patients. Information on survival was retrieved from the Swedish Cause of Death Register and information on socioeconomic status was retrieved from LISA.

The outcome was death 0–3 months after stroke (defined as death within 90 days after stroke onset).

The choice of proxy for low SES was motivated by previous studies. In a study on a similar cohort we found income to be more closely associated with post-stroke survival than education [5]. This has also been found in register based studies in Denmark and Norway [6, 9]. Based on this we chose to use low income as a proxy for low SES. As in our previous study, we define low income as being in the lowest tertile of the individual’s part of the family disposable income the year before the stroke.

The mediator of interest, stroke severity, was measured using level of consciousness on admission to hospital, based on the Re-action Level Scale (RLS) [23]. An RLS of 1 corresponded to fully conscious and an RLS 2–8 to lowered consciousness.

The covariates included for confounding adjustment were selected based on medical knowledge and availability in Riksstroke. In addition, the definition of a confounder in this scenario is a variable that affects at least two of the exposure, mediator and outcome. Based on these criteria our analyses were adjusted for a set of baseline confounders consisting of patient age (years), sex (male/female), whether or not the patient was living alone at the time of stroke, and the cardiovascular risk factors atrial fibrillation, diabetes and smoking history (smoker/non-smoker/unknown).

We opted not to adjust our main analyses for other SES variables (e.g. education level) as we were interested in a general effect of SES, not specifically a pure income effect. Supplementary analyses with additional adjustment for highest attained education level were performed.

This study was approved by the Ethical Review Board in Umeå.

Statistical methods

The distributions of characteristics among low and mid/high income patients within each stroke type were compared using Pearson χ2 tests (categorical variables) and independent samples t tests (continuous variables).

The directed acyclic graph (DAG) in Fig. 2 illustrates the causal relations assumed between the variables in the study. Our exposure of interest was low income, which we assumed had a causal effect on lowered consciousness upon hospital arrival and death 0–3 months after stroke for both stroke types. The baseline confounders were age, sex, living alone, atrial fibrillation, diabetes and smoking history. We decomposed the total effect of low income on death 0–3 months into two parts, an indirect effect of low income on death 0–3 months that acts through lowered consciousness upon arrival and a direct effect that does not act through lowered consciousness. Here we focus on a decomposition of the total effect into the natural indirect effect and the natural direct effect [12, 13]. In our study the indirect effect contrasts the risk of death for patients with low income with their observed level of consciousness to a counterfactual risk where level of consciousness is set to levels that would have been observed had the low income patient been mid/high income. The direct effect contrasts the risk of death for patients with low income and patients with mid/high income if, for all patients, level of consciousness were set to levels that would have been observed had the patient been mid/high income. A more formal definition of these effects is given in Section 1 of Additional file 1.

Fig. 2
figure 2

Directed acyclic graph showing the causal relations assumed between the variables in the study

To estimate the direct and indirect effects we used a regression model based approach [17, 24]. For each stroke type two regression models were fitted, one for the mediator (lowered consciousness upon arrival) given the exposure (low income) and the observed confounders, and one for the outcome (death 0–3 months) given exposure, mediator and observed confounders. As both the mediator and outcome are binary we fit probit regression models. The estimated regression parameters from these models were then used to estimate the direct and indirect effects (see Additional file 1, Section 1 for details).

The estimates of the direct and indirect effects are valid under the assumption that there is no residual confounding of the exposure-mediator, mediator-outcome or exposure-outcome relations and provided that the regression models are correctly specified [24]. See Fig. 1 in Additional file 1 for a DAG illustrating the three types of residual confounding. We also assume that there is no confounder of the mediator-outcome relation (observed or unobserved) that is caused by the exposure.

To assess the robustness of our results to the presence of residual confounding we performed sensitivity analyses using a method introduced in Lindmark et al. [17] We use a sensitivity parameter to quantify the effect of a possible unobserved confounder. The goal is then to capture the effect that would have been obtained had we been able to adjust for the unobserved confounder.

For our sensitivity parameter we use the fact that omitting a confounder from our regression models means that the error terms in these models will depend on the omitted confounder and therefore be correlated. Negative correlations correspond to an unobserved confounder that affects both variables of interest in opposite directions (e.g. increasing the risk of the mediator but decreasing the risk of the outcome). Positive correlations correspond to the unobserved confounder affecting both variables of interest in the same direction, either increasing the risk of both or decreasing the risk of both. The larger the absolute value of the correlation (closer to − 1 or + 1), the stronger the effect of the unobserved confounder, with a correlation of 0 indicating no unobserved confounding.

We performed three separate sensitivity analyses, one for each type of confounding. For mediator-outcome confounding the sensitivity parameter is the correlation between the error terms of the mediator and outcome models. For unobserved exposure-mediator and exposure-outcome confounding the sensitivity parameters are the correlations between a model for the exposure given the observed confounders and the mediator and outcome models, respectively.

In the sensitivity analysis we obtained estimates of the direct and indirect effects under a given level of unobserved confounding (a given value of the sensitivity parameter, for details see Section 2 of Additional file 1). Sensitivity analyses are sometimes criticized for being subject to the investigators choice of sensitivity parameter values. To increase objectivity we report the results of the sensitivity analyses based on recommendations in VanderWeele [25]:

  1. i.

    Estimates are reported across a wide range of the sensitivity parameter (correlations ranging from − 0.9 to 0.9).

  2. ii.

    The value of the sensitivity parameter required to reduce the estimated effect to 0 is reported.

  3. iii.

    The most important observed confounder is identified and the results are compared to the impact of omitting this confounder from the estimation.

All analyses were performed using R (version 3.6) [26]. Both the estimation and the sensitivity analyses were performed using the sensmediation package [27]. See Additional file 2 for R code for the mediation and sensitivity analyses.

Results

Descriptive statistics by stroke type and income level are displayed in Table 1. Patients with ICH had higher short term case fatality and more severe stroke than patients with IS, but were on average younger and had a lower proportion of risk factors.

Table 1 Descriptive statistics by stroke type. Number (%)

Comparing low income and mid/high income patients within each stroke type we see that the proportion of patients who died 0–3 months and the proportion of patients with lowered level of consciousness was higher among low income patients for both ICH and IS (Table 1). The low income patient group was also to a larger extent female, living alone and on average older than the mid/high income group for both IS and ICH. For IS the low income group also exhibited higher proportions of atrial fibrillation and diabetes compared to the mid/high income group, while the proportion of smokers was lower among low income patients.

The estimated probit regression models for the mediator lowered level of consciousness and the outcome death 0–3 months used for estimation of the indirect and direct effects are shown in Section 1 of Additional file 3. Following recommendations by VanderWeele [25] we initially included an interaction term between income and level of consciousness in the outcome models. As these interaction terms were small in magnitude and the estimated direct and indirect effects were virtually identical with and without interactions we did not include them in further analyses.

The estimated total effect of low income on death within 3 months after stroke was larger for ICH than for IS. For ICH, low income patients had a 3.2% (95% CI: 0.9–5.4%) increased absolute risk of death 0–3 months compared to mid/high income patients even after adjustment for observed confounders. The corresponding increase for IS was 1% (0.4–1.5%).

In the decomposition of the total effect for ICH, the indirect effect of low income, mediated by lowered level of consciousness, was 1.8% (0.7–2.9%) and the direct effect 1.4% (− 0.6–3.3%). For IS the indirect effect was 0.4% (0.2–0.6%) and the direct effect 0.6% (0.1–1.1%). Looking at effect sizes, a large part of the total effect of low income on death 0–3 months after stroke was an indirect effect working through lowered level of consciousness. For ICH over half the effect (57.0%) of income on death within 3 months was explained by lowered consciousness upon arrival. This proportion was smaller for IS but still close to two fifths (38.5%) of the effect of low income was mediated by lowered level of consciousness.

Effect of possible unmeasured confounding

The results of sensitivity analyses of the estimated indirect effect to residual confounding are displayed in Figs. 3, 4 and 5, respectively. The sensitivity parameter (denoted ρ) ranged from − 0.9 to 0.9. To assess the sensitivity to residual confounding involving the exposure, an additional probit regression model for low income given the observed confounders was estimated (see Additional file 3, Section 1).

Fig. 3
figure 3

Results of sensitivity analyses to residual confounding of the income-stroke severity relation on the estimated natural indirect effect for ICH and IS

Fig. 4
figure 4

Results of sensitivity analyses to residual confounding of the stroke severity-death 0–3 months relation on the estimated natural indirect effect for ICH and IS

Fig. 5
figure 5

Results of sensitivity analyses to residual confounding of the income-death 0–3 months relation on the estimated natural indirect effect for ICH and IS

In Fig. 3 we see similar overall patterns across the range of the sensitivity parameter for both ICH and IS. Adjusting for a confounder which affects the risk of low income and death 0–3 months in opposite directions (negative ρ) would increase the estimated indirect effect, while adjusting for a confounder which affects the risk of low income and death 0–3 months in the same direction (positive ρ) would decrease the estimated indirect effect. Moderate positive values of ρ (between 0 and 0.1) would yield a zero effect.

We see a similar but weaker pattern in the sensitivity analysis to residual confounding between the level of consciousness-death 0–3 months relation (Fig. 4). Here, an unobserved confounder would need to have a strong impact (ρ around 0.7 for ICH and ρ around 0.6 for IS) in order for additional adjustment to reduce the estimated indirect effect to 0.

In Fig. 5 we see that the estimated indirect effect does not appear to be sensitive to residual confounding of the income-death 0–3 months relation as it remained relatively constant over the range of the sensitivity parameters for both stroke types.

In the sensitivity analyses for the estimated direct effect we see that the direct effect is not sensitive to residual confounding of the income-level of consciousness relation (Fig. 6). The estimated direct effect increased over the range of the sensitivity parameter in the sensitivity analysis to residual confounding of the level of consciousness-death 0–3 month relation (Fig. 7) but decreased over the range of the sensitivity parameter in the sensitivity analysis to residual confounding of the income-death 0–3 month relation (Fig. 8). For the latter, moderate positive values of ρ (close to 0) would yield a zero effect.

Fig. 6
figure 6

Results of sensitivity analyses to residual confounding of the income-stroke severity relation on the estimated natural direct effect for ICH and IS

Fig. 7
figure 7

Results of sensitivity analyses to residual confounding of the stroke severity-death 0–3 months relation on the estimated natural direct effect for ICH and IS

Fig. 8
figure 8

Results of sensitivity analyses to residual confounding of the income-death 0–3 months relation on the estimated natural direct effect for ICH and IS

Finally, we looked at the effect of omitting age, the most important observed predictor of income, level of consciousness and death 0–3 months, from our models. When not adjusting for age the estimated indirect effect was 2.0% (95% CI: 0.9–3.1%) and the direct effect 2.1% (0.1–4.1%) for ICH. For IS the indirect effect was 0.6% (0.4–0.8%) and the direct effect was 1.4% (0.9–1.9%) without adjustment for age. That is, additional adjustment for age yielded a modest decrease of the estimated indirect effect and a more substantial decrease of the estimated direct effect and therefore an increase in the proportion mediated for both stroke types. These values cannot be directly translated to correlations in Figs. 3, 4, 5, 6, 7 and 8, since correlation values correspond to looking at each type of confounding in isolation, assuming that the other two types are not present. This assumption does not hold for age, which affects all three relations, but we still get an indication that it would take an unobserved confounder of greater impact than age to considerably alter our conclusions.

Discussion

Our results showed an increased risk of death within 3 months after stroke for low vs. mid/high income patients with first-time stroke and who were independent in ADL and living at home at the time of stroke. The gradient was especially pronounced for patients with ICH, corresponding to 89 additional deaths per 10,000 patients among low income patients compared to mid/high income patients. For IS the corresponding number was 10 per 10,000 patients. We found that for ICH nearly three fifths and for IS nearly two fifths of the effect operated through stroke severity (lowered consciousness upon hospital arrival).

An effect of low income/SES on death after stroke has been shown in studies from Denmark, England, Canada, Norway, Wales and the USA [6,7,8,9,10,11, 28, 29]. Several studies on the SES-stroke survival relation have adjusted the analyses for stroke severity, either as confounding adjustment [5,6,7, 11] or specifically focusing on the attenuation of the estimated SES-survival effect [8]. Earlier studies have also shown that patients with low SES tend to have more severe strokes than high SES patients [4, 6, 11]. However, ours is the first study to estimate the proportional contribution of stroke severity of the effect of SES on survival; we found that the contribution of stroke severity was substantial even after adjustment for other established prognostic factors such as age and atrial fibrillation. Further studies are needed to determine how SES is pathophysiologically linked to stroke severity. Potential factors may e.g. be differences in the early recognition of stroke, time to hospital, and other pre-hospital factors. Further studies should also include further subtyping of ischemic and hemorrhagic stroke, and results from neuroimaging and other diagnostic procedures.

Our results also indicate that a not unsubstantial part of the effect of income on death within 3 months is not mediated by stroke severity. The remaining association could be due to other pathways that link income with survival, e.g. differences in acute stroke management and/or secondary prevention.

In all observational studies there is a risk of residual confounding. Here, we used a recent method for sensitivity analyses [17] to assess the effect of violations of the crucial assumptions about residual confounding of not only the relation between the mediator (lowered consciousness) and the outcome (death within 3 months) but also the relations involving the exposure (low income).

We found that the indirect effect was most sensitive to unobserved confounders that increase both the risk of having low income and having lowered consciousness upon arrival to hospital. The direct effect was most sensitive to unobserved confounders that increase both the risk of having low income and death within 3 months after stroke. Individual linkage of several national registers can provide information on a rich set of possible confounders. Still, genetic factors, life style (e.g. diet, exercise habits, alcohol consumption), and demography (e.g. distance from the patient’s residence to the hospital) may contribute to residual confounding in studies such as ours. Additional adjustment for such confounders could potentially explain the observed effects. We did observe that additional adjustment for the most important observed confounder, age, attenuated the estimated direct effect more than the indirect effect leading to an increased proportion of the total effect mediated by stroke severity. This was also observed with additional adjustment for education level (see below).

In order to maintain a high coverage and participation rate in Riksstroke, the amount of information collected in the register has been kept at a minimum, and does not include detailed information on e.g. comorbidity, stroke location or size. In addition to the potential risk of residual confounding, there are also other limitations to the study. SES is a multidimensional concept that can be measured in different ways. Several previous studies, including one on a similar cohort, found income to be more strongly associated with post-stroke survival than education [5, 6, 9]. As we were interested in a general effect of SES, not specifically a pure effect of low income we opted not to adjust for other SES variables, such as education level, in our analyses. Supplementary analyses with adjustment for education showed an attenuation of the total effect but an increased proportion mediated (see Additional file 3, Section 2). We measure income the year before stroke, focusing on adulthood SES, but there is evidence in favor of a link between childhood SES and adulthood stroke [30] as well as a cumulative effect of SES over the life course [31] that is probably not reflected in our measurement.

We used level of consciousness, based on RLS, upon arrival to hospital as a proxy for stroke severity, our mediator of interest. This measurement was dichotomized into fully conscious and lowered level of consciousness. It is possible that the use of a more comprehensive measurement of stroke severity, such as NIHSS, would give different results. However, level of consciousness has been proven a good proxy for the full NIHSS in predicting death after stroke [32].

Even if it is evident that stroke severity plays an important role in the mechanism explaining the increased risk of death after stroke in low SES groups, the exact proportion mediated may not be generalized to countries with e.g. different demography, and health care settings.

We made assumptions about the directions of associations between the variables in the study. This cannot be directly ascertained from cross sectional data, and was based on clinical reasoning and prior knowledge. A statistical model will always be a simplification of the real world. Here we used mediation analysis to add information on the complex relation between SES and survival. Still, we did not consider e.g. intermediate confounders (confounders of the stroke severity- death within 3 months relation that are affected by income), or several mediating factors. This would require the use of alternative effect definitions and methods and may be a subject for future studies [33,34,35,36].

Conclusions

This study suggests that as much as half of the effect of low SES on increased risk of death operates through stroke severity. Hence, targeting stroke severity could potentially lead to a large reduction in inequalities and should therefore be prioritized.

Availability of data and materials

The data set comprises linked information from Riksstroke, Statistics Sweden and the National Board of Health and Welfare and cannot be made publically available due to Swedish legislation. Requests to access the dataset may be sent to Riksstroke, Statistics Sweden and the National Board of Health and Welfare after obtaining the appropriate ethics approval. Code for the main analyses are available in Additional file 2 and any additional code related to the results of the study are available from the corresponding author upon reasonable request.

Abbreviations

WHO:

World Health Organization

SES:

Socioeconomic status

ICH:

Intracerebral hemorrhage

IS:

Ischemic stroke

ADL:

Activities of daily living

LISA:

Longitudinal Integration Database for Health Insurance and Labor Market Studies

RLS:

Reaction Level Scale

DAG:

Directed acyclic graph

st. dev:

Standard deviation

CI:

Confidence interval

USA:

United States of America

NIHSS:

National Institutes of Health Stroke scale

References

  1. Pickett KE, Wilkinson RG. Income inequality and health: a causal review. Soc Sci Med. 2015;128:316–26. https://doi.org/10.1016/j.socscimed.2014.12.031.

    Article  PubMed  Google Scholar 

  2. Beckfield J, Olafsdottir S. Health inequalities in global context. Am Behav Sci. 2013;57(8):1014–39.

    Article  Google Scholar 

  3. Mackenbach JP, Stirbu I, Roskam AJ, Schaap MM, Menvielle G, Leinsalu M, et al. Socioeconomic inequalities in health in 22 European countries. N Engl J Med. 2008;358(23):2468–81. https://doi.org/10.1056/NEJMsa0707519.

    Article  CAS  PubMed  Google Scholar 

  4. Marshall IJ, Wang Y, Crichton S, McKevitt C, Rudd AG, Wolfe CD. The effects of socioeconomic status on stroke risk and outcomes. Lancet Neurol. 2015;14(12):1206–18. https://doi.org/10.1016/S1474-4422(15)00200-8.

    Article  PubMed  Google Scholar 

  5. Lindmark A, Glader EL, Asplund K, Norrving B, Eriksson M. Socioeconomic disparities in stroke case fatality - observations from Riks-stroke, the Swedish stroke register. Int J Stroke. 2014;9(4):429–36. https://doi.org/10.1111/ijs.12133.

    Article  PubMed  Google Scholar 

  6. Andersen KK, Dalton SO, Steding-Jessen M, Olsen TS. Socioeconomic position and survival after stroke in Denmark 2003 to 2012: nationwide hospital-based study. Stroke. 2014;45(12):3556–60. https://doi.org/10.1161/STROKEAHA.114.007046.

    Article  PubMed  Google Scholar 

  7. Bray BD, Paley L, Hoffman A, James M, Gompertz P, Wolfe CDA, et al. Socioeconomic disparities in first stroke incidence, quality of care, and survival: a nationwide registry-based cohort study of 44 million adults in England. Lancet Public Health. 2018;3(4):e185–e93. https://doi.org/10.1016/S2468-2667(18)30030-6.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Kapral MK, Fang J, Chan C, Alter DA, Bronskill SE, Hill MD, et al. Neighborhood income and stroke care and outcomes. Neurology. 2012;79(12):1200–7. https://doi.org/10.1212/WNL.0b013e31826aac9b.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Skyrud KD, Vikum E, Hansen TM, Kristoffersen DT, Helgeland J. Hospital variation in 30-day mortality for patients with stroke; the impact of individual and municipal socio-demographic status. J Am Heart Assoc. 2019;8(14):e010148. https://doi.org/10.1161/JAHA.118.010148.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Thorne K, Williams JG, Akbari A, Roberts SE. The impact of social deprivation on mortality following acute myocardial infarction, stroke or subarachnoid haemorrhage: a record linkage study. BMC Cardiovasc Disord. 2015;15:71. https://doi.org/10.1186/s12872-015-0045-x.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Hanchate AD, Schwamm LH, Huang W, Hylek EM. Comparison of ischemic stroke outcomes and patient and hospital characteristics by race/ethnicity and socioeconomic status. Stroke. 2013;44(2):469–76. https://doi.org/10.1161/STROKEAHA.112.669341.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Robins JM, Greenland S. Identifiability and exchangeability for direct and indirect effects. Epidemiology. 1992;3(2):143–55. https://doi.org/10.1097/00001648-199203000-00013.

    Article  CAS  PubMed  Google Scholar 

  13. Pearl J. Direct and indirect effects. Proceedings of the Seventeenth conference on Uncertainty in artificial intelligence. Seattle: 2074073: Morgan Kaufmann Publishers Inc.; 2001. p. 411–20.

    Google Scholar 

  14. VanderWeele TJ. Bias formulas for sensitivity analysis for direct and indirect effects. Epidemiology. 2010;21(4):540–51. https://doi.org/10.1097/EDE.0b013e3181df191c.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Imai K, Keele L, Tingley D. A general approach to causal mediation analysis. Psychol Methods. 2010;15(4):309–34. https://doi.org/10.1037/a0020761.

    Article  PubMed  Google Scholar 

  16. Tchetgen Tchetgen EJ, Shpitser I. Semiparametric theory for causal mediation analysis: efficiency bounds, multiple robustness and sensitivity analysis. Ann Stat. 2012;40(3):1816–45. https://doi.org/10.1214/12-AOS990.

    Article  Google Scholar 

  17. Lindmark A, de Luna X, Eriksson M. Sensitivity analysis for unobserved confounding of direct and indirect effects using uncertainty intervals. Stat Med. 2018;37(10):1744–62. https://doi.org/10.1002/sim.7620.

    Article  PubMed  Google Scholar 

  18. Vanderweele TJ, Arah OA. Bias formulas for sensitivity analysis of unmeasured confounding for general outcomes, treatments, and confounders. Epidemiology. 2011;22(1):42–52. https://doi.org/10.1097/EDE.0b013e3181f74493.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Albert JM, Wang W. Sensitivity analyses for parametric causal mediation effect estimation. Biostatistics. 2015;16(2):339–51. https://doi.org/10.1093/biostatistics/kxu048.

    Article  PubMed  Google Scholar 

  20. Andersen KK, Olsen TS, Dehlendorff C, Kammersgaard LP. Hemorrhagic and ischemic strokes compared: stroke severity, mortality, and risk factors. Stroke. 2009;40(6):2068–72. https://doi.org/10.1161/STROKEAHA.108.540112.

    Article  PubMed  Google Scholar 

  21. Henriksson KM, Farahmand B, Asberg S, Edvardsson N, Terent A. Comparison of cardiovascular risk factors and survival in patients with ischemic or hemorrhagic stroke. Int J Stroke. 2012;7(4):276–81. https://doi.org/10.1111/j.1747-4949.2011.00706.x.

    Article  PubMed  Google Scholar 

  22. Asplund K, Hulter Asberg K, Appelros P, Bjarne D, Eriksson M, Johansson A, et al. The Riks-stroke story: building a sustainable national register for quality assessment of stroke care. Int J Stroke. 2011;6(2):99–108. https://doi.org/10.1111/j.1747-4949.2010.00557.x.

    Article  PubMed  Google Scholar 

  23. Starmark JE, Stalhammar D, Holmgren E. The reaction level scale (RLS85). Manual and guidelines. Acta Neurochir. 1988;91(1–2):12–20.

    Article  CAS  Google Scholar 

  24. VanderWeele TJ, Vansteelandt S. Conceptual issues concerning mediation, interventions and composition. Stat Interface. 2009;2(4):457–68.

    Article  Google Scholar 

  25. VanderWeele TJ. Explanation in causal inference: methods for mediation and interaction. 1st ed. New York: Oxford University Press; 2015.

    Google Scholar 

  26. R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017.

    Google Scholar 

  27. Lindmark A. sensmediation: Parametric estimation and sensitivity analysis of direct and indirect effects. R package version 0.3.0. 2018. http://cran.R-project.org/package=sensmediation.

    Google Scholar 

  28. Brown AF, Liang LJ, Vassar SD, Merkin SS, Longstreth WT Jr, Ovbiagele B, et al. Neighborhood socioeconomic disadvantage and mortality after stroke. Neurology. 2013;80(6):520–7. https://doi.org/10.1212/WNL.0b013e31828154ae.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Elfassy T, Grasset L, Glymour MM, Swift S, Zhang L, Howard G, et al. Sociodemographic disparities in long-term mortality among stroke survivors in the United States. Stroke. 2019;50(4):805–12. https://doi.org/10.1161/STROKEAHA.118.023782.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Galobardes B, Smith GD, Lynch JW. Systematic review of the influence of childhood socioeconomic circumstances on risk for cardiovascular disease in adulthood. Ann Epidemiol. 2006;16(2):91–104. https://doi.org/10.1016/j.annepidem.2005.06.053.

    Article  PubMed  Google Scholar 

  31. Pollitt RA, Rose KM, Kaufman JS. Evaluating the evidence for models of life course socioeconomic factors and cardiovascular outcomes: a systematic review. BMC Public Health. 2005;5:7. https://doi.org/10.1186/1471-2458-5-7.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bray BD, Campbell J, Cloud GC, Hoffman A, James M, Tyrrell PJ, et al. Derivation and external validation of a case mix model for the standardized reporting of 30-day stroke mortality rates. Stroke. 2014;45(11):3374–80. https://doi.org/10.1161/Strokeaha.114.006451.

    Article  PubMed  Google Scholar 

  33. VanderWeele TJ, Vansteelandt S. Mediation analysis with multiple mediators. Epidemiol Methods. 2014;2(1):95–115. https://doi.org/10.1515/em-2012-0010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Vanderweele TJ, Vansteelandt S, Robins JM. Effect decomposition in the presence of an exposure-induced mediator-outcome confounder. Epidemiology. 2014;25(2):300–6. https://doi.org/10.1097/EDE.0000000000000034.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Daniel RM, De Stavola BL, Cousens SN, Vansteelandt S. Causal mediation analysis with multiple mediators. Biometrics. 2015;71(1):1–14. https://doi.org/10.1111/biom.12248.

    Article  CAS  PubMed  Google Scholar 

  36. Vansteelandt S, Daniel RM. Interventional effects for mediation analysis with multiple mediators. Epidemiology. 2017;28(2):258–65. https://doi.org/10.1097/Ede.0000000000000596.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to Riksstroke and the participating hospitals.

Funding

The study was supported by FORTE (Swedish Research Council for Health, Working Life and Welfare, grant 2018–00852) and the Swedish Research Council (grant 2018–02670). The funding bodies had no role in study design, data collection, analysis, interpretation, or writing of this manuscript. Open access funding provided by Umeå University.

Author information

Authors and Affiliations

Authors

Contributions

AL and ME conceived the original idea. AL conducted the data preparation and analysis. AL drafted the manuscript and BN and ME provided substantive revisions of previous versions. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Anita Lindmark.

Ethics declarations

Ethics approval and consent to participate

All patients are informed about registration in Riksstroke, that the register may be used for research purposes after ethical approval, and are informed of their rights to decline participation (verbal opt-out consent) and to withdraw their data at any time (this information can be found at www.riksstroke.org/patient-information/). The regional Ethical Review Board in Umeå Sweden, approved the study (reference number 2017–184-31). The committee waived the need for additional patient consent for this research study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Formal definitions and technical details.

Additional file 2.

R code.

Additional file 3.

Tables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lindmark, A., Norrving, B. & Eriksson, M. Socioeconomic status and survival after stroke – using mediation and sensitivity analyses to assess the effect of stroke severity and unmeasured confounding. BMC Public Health 20, 554 (2020). https://doi.org/10.1186/s12889-020-08629-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-020-08629-1

Keywords