Health related quality of life and satisfaction with care of stroke patients in Budapest: A substudy of the EuroHOPE project

Background Disadvantaged socioeconomic status is associated with higher stroke incidence and mortality, and higher readmission rate. We aimed to assess the effect of socioeconomic factors on case fatality, health related quality of life (HRQoL), and satisfaction with care of stroke survivors in the framework of the European Health Care Outcomes, Performance and Efficiency (EuroHOPE) study in Hungary, one of the leading countries regarding stroke mortality. Methods We evaluated 200 consecutive patients admitted for first-ever ischemic stroke in a single center and performed a follow-up at 3 months after stroke. We recorded pre- and post-stroke socioeconomic factors, and assessed case fatality, HRQoL and patient satisfaction with the care received. Stroke severity at onset was scored by the National Institutes of Health Stroke scale (NIHSS), disability at discharge from acute care was evaluated by the modified Rankin Score (mRS). To evaluate HRQoL and patient satisfaction with care we used the EQ-5D-5L, 15D and EORTC IN PATSAT 32 questionnaires. Results At 3 months after stroke the odds of death was significantly increased by stroke severity (NIHSS, OR = 1.209, 95%CI: 1.125–1.299, p<0.001) and age (OR = 1.045, 95%CI: 1.003–1.089, p = 0.038). In a multiple linear regression model, independent predictors of HRQoL were age, disability at discharge, satisfaction with care, type of social dwelling after stroke, length of acute hospital stay and rehospitalization. Satisfaction with care was influenced negatively by stroke severity (Coef. = -1.111, 95%C.I.: -2.159- -0.062, p = 0.040), and positively by having had thrombolysis (Coef. = 25.635, 95%C.I.: 5.212–46.058, p = 0.016) and better HRQoL (Coef. = 22.858, 95%C.I.: 6.007–39.708, p = 0.009). Conclusion In addition to age, disability, and satisfaction with care, length of hospital stay and type of social dwelling after stroke also predicted HRQoL. Long-term outcome after stroke could be improved by reducing time spent in hospital, i.e. by developing home care rehabilitation facilities thus reducing the need for readmission to inpatient care.


Introduction
Stroke has a major impact on the quality of life of patients [1,2], with considerable differences between populations with various socioeconomic background [3]. People of lower socioeconomic status show higher stroke incidence and mortality, higher severity in the acute phase, higher readmission rate to hospital and lose more disability-adjusted life-years [4][5][6][7]. Disparities in clinical or behavioral risk factors or differential access to healthcare facilities do not completely explain regional differences in stroke incidence and outcome [8][9][10][11][12].
The impact of differences in social background on post-stroke health-related quality of life (HRQoL) has not been sufficiently evaluated in Central-Eastern European countries. Although there are some data on the consequences of social inequalities on stroke features among neighborhoods within the capital city of Hungary [13], reports evaluating HRQoL data after stroke are scarce from Central-Eastern European countries [14].
The European Health Care Outcomes, Performance and Efficiency (EuroHOPE) project compared health system performance in 6 European countries (Italy, Finland, Sweden, Scotland, the Netherlands and Hungary) [15], evaluating stroke, acute myocardial infarction, hip fracture, breast cancer, and low-birth-weight. In the EuroHOPE study Hungary had the highest stroke incidence, the largest all-cause case fatality 1 year after stroke and the largest regional differences [16,17].
To improve long-term outcome after stroke, it would be of major importance to identify social groups and features of post-stroke care that are associated with fatality or reduced HRQoL. We previously found worse stroke outcome in patients residing in the poorest district of Budapest compared to those of the wealthiest neighborhood [18,19], and also found that the lower the annual taxable income, the lower the age at stroke onset in the 23 districts of Budapest [13].
The aim of the present study was to assess the impact of stroke-related, demographic and socioeconomic factors on acute and 3-months case fatality, on HRQoL and on satisfaction with care of patients after stroke.

Study setting and recruitment of participants
Data presented here were collected within the framework of the EuroHOPE project [15], a retrospective study based on national hospital healthcare administrative records in 6 European countries, with the index ischemic stroke between 2006-2008 [16,20]. The present study was conducted according to Work Package 2 of the EuroHOPE project ("HRQol protocol and patient satisfaction questionnaire-Stroke"). We prospectively included patients with ischemic stroke (International Classification Diseases, ICD-10, I63). According to the study protocol, inclusion of consecutive patients with first-ever ischemic stroke was planned up till including 200 patients. That was intended to be an exploratory study with no formal sample size estimation. Exclusion criteria were history of previous ischemic stroke and refusal by the patient to participate in the survey.
Patients were recruited between September 2012 and May 2013, and evaluated at the Department of Neurology of Semmelweis University, a stroke center with a catchment area covering Districts 8 and 9 of Budapest, and the Ráckeve region. The catchment area of this stroke center was defined by geographic availability, i.e. to be the nearest stroke center for the regions of the catchment area. From the socioeconomic point of view, District 8 ranks the last of the 23 districts of Budapest on the list of income; the other regions of the catchment area are in the medium or high income categories. Income constitutes the major socioeconomic difference between the regions in the catchment area. Due to consecutive patient inclusion, our set of patients may be considered representative of ischemic stroke patients of this stroke center's catchment area.

Study design
Ethical approval was obtained from the Regional and Institutional Committee of Science and Research Ethics of Semmelweis University (no. 98/2012) and written consent was obtained from each patient or from the carer in case of patients who were unable to consent. Initial evaluation of the patients was performed at admission by the physician on duty. The research investigators re-evaluated each case to ascertain that it is indeed an acute first-ever ischemic stroke. As the next step, clinical features, pre-stroke demographic and socioeconomic factors ( Table 1) were recorded by a standard questionnaire. Follow-up was performed three months after the onset of stroke by personal visit or by phone interview to assess post-stroke socioeconomic factors, HRQoL and patient satisfaction with care. Length of post-acute hospital stay (i.e. institutional rehabilitation or rehospitalization for any other reason after discharge from

PLOS ONE
Stroke outcome, quality of life and patient satisfaction in residents of Budapest the acute setting) was recorded. Post-acute hospital stay was defined as the number of all inpatient days after discharge from the index hospital admission within the three months of follow-up. Outcome measures were case fatality during the initial hospitalization and during followup, HRQoL and satisfaction with care of the survivors at follow-up.

Data collection instruments
During admission and at follow-up we have used detailed, standardized questionnaires to assess sociodemographic and socioeconomic factors listed in Table 1.
Initial stroke severity was assessed by the National Institutes of Health Stroke Scale (NIHSS), [21] on admission, whereas disability at discharge was scored by the modified Rankin Scale (mRS), [22]. The NIHSS is a 15-item neurological examination scale used to evaluate the effect of stroke on the level of consciousness, language, neglect, visual-field loss, extraocular movement, motor strength, coordination, dysarthria, and sensory loss. The NIHSS score ranges from 0 to 42, higher scores representing more severe states. The mRS addresses disability ranging from 0 (asymptomatic) to 6 (death). The classification of ischemic stroke was performed using Trial of Org 10172 in Acute Stroke Treatment (TOAST) criteria [23]. According to TOAST, ischemic stroke can be classified into large-artery atherosclerosis, cardioembolism, small-vessel disease, stroke of other determined etiology, and stroke of undetermined etiology.
Stroke and its consequences have significant effects on health related quality of life (HRQoL). We applied indirect methods to evaluate HRQoL in our study with the use of multiattribute utility (MAU) instruments: by the EuroQoL-5 Dimension-5 Levels (hereafter EQ-5D), developed by the EuroQoL group [24], and by the HRQoL 15D instrument (hereafter 15D) [25,26]. Both questionnaires are indirect generic instruments to evaluate HRQoL. EQ-5D addresses 5 dimensions: mobility, self-care, usual activities, pain/discomfort and anxiety/ depression, grading these on 5 levels and a visual analog scale. The overall HRQoL described by the EQ-5D can be used to derive health state utility values, and has been found to be a valid descriptive system as a generic health outcome measure in patients with acute stroke [27]. The 15D is a comprehensive, self-administered health-related QoL instrument, which consists of 15 dimensions: breathing, mental function, speech, vision, mobility, usual activity, vitality, hearing, eating, elimination, sleeping, distress, discomfort and symptoms, depression and sexual activity. While the EQ-5D is widely used [28], the results might vary depending on the regression method used [29]. EQ-5D seems somewhat less sensitive when it comes to evaluating several chronic diseases [30][31][32][33]. The 15D performs better in terms of reliability, discrimination and responsiveness [34], and performs well after critical illness [28]. The agreement between the two utility measures was only moderate [28], and both the mean utilities and the standard deviation differs between EQ-5D and 15D [35]. For these reasons-as recommended by Feeny et al. [36]-we decided to use multiple questionnaires in this study.
We have also assessed the experience of patients with the care received during the hospital stay by EORTC IN-PATSAT32 (hereafter PATSAT), a tool that was developed by the European Organization for Research and Treatment of Cancer [37]. It evaluates 32 items regarding the satisfaction with technical competence, information provision, interpersonal skills, availability, waiting time, access, comfort and overall care perception. The following terms can be used: global patient satisfaction expressed by the PATSAT score of the questionnaire; three subcategories of patient satisfaction (with doctors, with nurses and with services and care organization) and the item number 32 of the questionnaire which is the general satisfaction with care received during the hospital stay of the patient. A higher scale score represents a higher level of satisfaction with care.

Statistical modeling
The relationship between variables was first assessed independently of any other effects (i.e. without control variables), followed by multivariate analysis.
We have built two types of statistical models, the first using variables related to pre-stroke features, characteristics of the acute stroke, stroke management and post-stroke variables as well. Including stroke related and post-stroke variables in the analysis was especially important when inspecting patient satisfaction with the care received, as the care received is necessarily applied after the onset of stroke. In a supplementary analysis we excluded variables that describe the patients after the acute stroke.
Logistic regression was used to evaluate acute and follow-up fatality. Linear regression was used during the analysis of the aggregated, continuous scores of HRQoL and patient satisfaction. The ordinal elements and scores of HRQoL and patient satisfaction were analyzed with ordered logistic regression. While HRQoL can be defined for the deceased, we chose to create HRQoL-related models only for the survivors, as the imputation of values for deceased patients could have introduced a bias. We also conducted analysis using post-stroke information, which were only available among the survivors.
Linear and logistic regression tables include the intercept, which is just the mean of the dependent variable or odds of the event when the values of the explanatory variables are zero or at the reference value. The ordered logistic regression models include multiple intercepts, but these are omitted from the tables since their value is not important for model interpretation.
In multivariate model building, first all potential covariates were included in the regression model alone to check their effect on the dependent variable. The multiple regression models then were built in the following way: we defined a set of control variables which were included in all models of a given type of dependent variable. Age, sex and NIHSS were control variables in all models. Discharge mRS was also included as control variable in all models that allowed post-stroke information. Additional control variables in the PATSAT models were education and marital status, in the HRQoL models education and in the death models dwelling district (8 th or not). We then assessed the effect of all other possible covariates one-by-one, controlled by the previously set variables. In the next step, we placed the control variables and all the other variables which in the previous step had a coefficient significant at the 0.1 level into the model. We then checked which variables were significant at the 0.05 level and eliminated those which were not. In this last step we also considered the viability of the model and checked the effect of pre-stroke and post-stroke version and the original and pooled version of the same variables and also interaction terms. We also aimed to synchronize models with similar dependent variables. The reason behind this algorithm is that it allowed us to check the effect of the variables together with the control variables and then all together. This way we were able to discover relationships which might have remained hidden if the variable selection was based only on univariate analysis. Automatized tools such as stepwise regression were considered but discarded due to low patient number.

Statistical tests
Distribution of continuous variables was checked by the Shapiro-Wilk test. In case of non-normal distribution, non-parametric tests (Mann-Whitney test and Kruskal-Wallis ANOVA) were used for comparisons. Univariate relationship between two continuous variables were analyzed by correlation analysis (Pearson, Kendall, Spearman), Kruskal-Wallis test, and the Mann-Whitney test. The relationship between two discrete variables was assessed using the chi-squared test of independence. The overall significance of a categorical variable with more than 2 levels in a multiple regression model was checked using likelihood-ratio test.
Goodness of fit of regression models was tested using statistical and visual tools. The Hosmer-Lemeshow test and separation plot [38] were used for logistic regression. Shapiro-Wilk test and quantile-quantile plot were used for linear regression. The ordinal models were also inspected by the Hosmer-Lemeshow test. Multicollinearity was checked using the variable inflation factor. R version 3.5.2 was used for data analysis with packages ggplot2, MASS, separationplot, generalhoslem, car, rms.

Participant characteristics
In total 200 patients participated in the study. The general description of the study population can be seen in Table 2. Patients from District 8 had more severe strokes (higher NIHSS on admission) and more severe disability at discharge (higher mRS scores). Despite an older age at stroke, fatal cases of dwellers of District 8 were younger both at the initial hospitalization and at 3 months after stroke, although none of these differences were statistically significant.

Case fatality
Overall hospital case fatality was 7.5% (12% among District 8 residents and 6% among those living outside District 8). At 3 months after stroke overall case fatality was 16.5% (22% vs 15% among District 8 residents and those living outside District 8). In a multiple logistic regression model admission NIHSS and living alone prior to stroke (vs. living with partner) increased the odds of acute death. No other variables-including the sociodemographic variables-had a statistically significant effect. At follow-up at 3 months, age at stroke onset, and stroke severity by the NIHSS were significant predictors of case fatality, while living alone vs. living with partner lost its significance. For further information, see Table 3. Our study was underpowered to make reliable comparisons between features of inhabitants of the poorer District 8 and other regions.

Questionnaire completion
Assessment of HRQoL and patient satisfaction at follow-up could have been possible in 81% of the survivors (136 cases alive, able and willing to express themselves, 7 unable-aphasic, demented or comatose, 24 patients declined to answer). The response rates among the 167 patients surviving 3 months were 122 (73%) for PATSAT score, 126 (75%) for EQ-5D utility index and 124 (74%) for 15D utility index. The difference from 136 is due to incomplete answers at item level. The response rate for HRQoL and patient satisfaction questionnaires were almost the same (73%-75%), while initial stroke severity and disability data (i.e. NIHSS and mRS) could be gathered from all 167 survivors. The 15D missing data were imputed at the item level using an SPSS algorithm for those with less than 4 missing items, based on the age and gender of the patients [25,26].

PLOS ONE
Stroke outcome, quality of life and patient satisfaction in residents of Budapest

Health-related quality of life
A. Main analysis. During the statistical model building we have initially included at least in the univariate analysis all the factors assessed in the study (see Table 1). The control variables were sex, age, education, NIHSS and discharge mRS. The effect of patient satisfaction was also checked in every HRQoL model. In the multiple linear regression model, we have found that out of all the factors assessed (Table 1) the following had a significant independent impact on the HRQoL measured by EQ-5D index at 90-days after stroke: age, mRS at discharge from the acute hospitalization, stroke type, patient satisfaction, acute and post-acute length of hospital stay, and the type of dwelling after stroke (Table 4). These last two features were related but caused no multicollinearity issues during statistical analysis. The effect of employment status prior to stroke and admission source were only marginally significant predictors of the EQ-5D index, while admission from hospital entailed significantly lower 15D utility scores on average.
The association of age and disability (i.e. mRS at discharge) with HRQoL is already known: the older and the more disabled report poorer quality of life. Some effects are worthy of closer inspection, though. TOAST classification influenced HRQoL in the study population. Cardioembolism was associated with better HRQoL compared to large-artery atherosclerosis

PLOS ONE
Stroke outcome, quality of life and patient satisfaction in residents of Budapest expressed in EQ-5D, with marginal effect on 15D (see Tables 4 and 5). Out of the socioeconomic factors assessed in this study, the type of social dwelling after stroke proved to be an independent predictor of EQ-5D index, meaning that the HRQoL at follow-up was highest in patients living alone independently in their home, somewhat lower in those at home with company, even lower in case of patients receiving family help; the poorest results were reported in patients dwelling with a professional carer (Fig 1). The length of post-acute hospital stay also proved to be an independent predictor of both the EQ-5D index and the 15D utility score: the longer the duration of post-acute hospitalization, the lower the HRQoL (Fig 2). We have also inspected the correlation between the scores of different instruments (EQ-5D utility index, 15D utility index, PATSAT score and its subcategories) and we have found that there is a significant positive correlation between each pair (see S4 Table). When examining the sub-item of general patient satisfaction (item 32 of the PATSAT questionnaire), we have also found a significant relationship with both HRQoL indices (S5 Table, Kruskal-Wallis test p-value = 0.019 for EQ-5D and 0.040 for 15D).

Patient satisfaction with care
Although 81% of the survivors were able and willing to respond to the questionnaires at follow-up, the response rate of PATSAT was 73%. Even less patients chose the categories of poor,

PLOS ONE
Stroke outcome, quality of life and patient satisfaction in residents of Budapest fair or good satisfaction. Thus, when performing ordinal logistic regression analysis, these categories had to be merged due to low patient numbers.
In the main analyses of patient satisfaction (see Table 6), the control variables were sex, age, education, NIHSS, discharge mRS and marital status. The effect of HRQoL was also checked in every patient satisfaction model, even if it was not included by the algorithm. Severity of stroke at onset measured by the NIHSS significantly influenced the global patient satisfaction. The EQ-5D index was found in the multiple regression model to be an independent predictor of the global PATSAT score and all its sub-categories at least at a marginally significant level (p<0.1). Thrombolysis also predicted better satisfaction globally, with the care provided by nurses and with organization. Altogether, the multiple linear regression models showed that the patient satisfaction is reduced by stroke severity at admission and is improved by better EQ-5D scores and by having had thrombolysis. The rest of the parameters described in Table 1 showed no effect on global patient satisfaction score or its subcategories. When we performed the multiple ordinal logistic regression model analysis on the question of general patient satisfaction (item 32 of the PATSAT), we have found that both mRS at discharge and EQ-5D index had a marginally significant effect on this item (see Table 6). All patient satisfaction linear regression models showed some deviation from normality (Shapiro-Wilk p<0.05), this somewhat weakens the reliability of the results. There were no significant problems with the fit of the ordinal logistic regression models.

PLOS ONE
Stroke outcome, quality of life and patient satisfaction in residents of Budapest B. Analysis excluding post-stroke variables. We excluded post-stroke variables during this type of model building: we included only age, sex, education, social and anamnestic factors, admission NIHSS and TOAST category. The control variables were sex, age, education and NIHSS. We present the results of these analyses for EQ-5D, 15D and PATSAT in S1-S3 Tables. The multiple linear regression models showed that more advanced age and higher stroke severity were associated with poorer utility indices of both the EQ-5D and 15D (see S1 and S2 Tables). It also showed that ischemic stroke subtype, determined by TOAST criteria might influence 3-month-post-stroke HRQoL. We have found that in this restricted model not only patients who have suffered cardioembolism but also those with small vessel occlusion reported better HRQoL at 3 months after stroke (expressed as higher utility indices of both EQ-5D and 15D) compared to patients who have suffered large vessel atherosclerosis. The stroke subtypes of other determined or undetermined etiology also showed a marginally significant effect on EQ-5D utility index, while this relationship was not detectable in terms of 15D utility index (see S1 and S2 Tables).
We have also found two associations detectable on only one of the two HRQoL instruments used. By the EQ-5D utility index, we have detected that patients who have not been employed prior to their stroke reported higher HRQoL compared to those patients who were employed at time of the stroke. By the 15D utility index, we have seen that patients admitted from another hospital at the onset of the index stroke reported poorer HRQoL compared to their

PLOS ONE
Stroke outcome, quality of life and patient satisfaction in residents of Budapest fellows admitted from their home. The only model that showed significant fitting problems (Shapiro-Wilk p = 0.009) among the HRQoL models was the EQ-5D multiple linear regression model excluding post-stroke variables (S1 Table). The result of the quantile-quantile plot was more dubious. As this was just a sub-analysis and the variable effects are similar in the other models this does not undermine our overall conclusions.
In the supplementary set of analyses (excluding post-stroke variables) we have found that the PATSAT score, satisfaction with services and care organization and the sub-item 32 are significantly influenced by admission NIHSS (see S3 Table). Namely, more severe stroke at admission results in poorer patient satisfaction. We found no significant predictors in these models for satisfaction with physicians and satisfaction with nurses. The control variables were sex, age, education, NIHSS and marital status.

Discussion
In a consecutive sample of 200 patients with acute ischemic stroke we evaluated predictors of case fatality, health-related quality of life and patient satisfaction. Predictors of acute in-hospital case fatality were stroke severity and living alone prior to stroke. At 3 months after stroke, similarly to previous reports, we also found that case fatality related to initial stroke severity and age. Health utility index was diminished after stroke compared to the population norm for EQ-5D for Hungary, and utility indices tended to be lower in the disadvantaged District 8 than in wealthier regions both by EQ-5D and 15D, suggesting that in addition to the consequences of stroke, the general living standard of the patient also affects post-stroke HRQoL. Independent predictors of HRQoL at 3 months after stroke, similarly to other reports, were age, disability at discharge from the acute hospitalization, stroke type, and patient satisfaction. As far as we know, it has not been reported previously that acute and post-acute length of hospital stay after stroke also affect HRQoL. We found evidence that post-acute inpatient management might adversely affect HRQoL: longer than 5 days post-acute LOS was associated with poorer HRQoL. In our study the single independent predictor of all subcategories of patient satisfaction was HRQoL. Global patient satisfaction was influenced negatively by initial stroke severity and positively by thrombolysis and HRQoL. According to our knowledge, this positive effect of the provision of thrombolysis on patient satisfaction has not been described before.

Case fatality
For fatality in the acute phase and at follow-up, the number of patients was low in the different subgroups. This restricted the detailed data analysis and conclusion regarding socioeconomic factors. Nevertheless, we have found that stroke severity and living alone prior to the stroke increase the odds of dying. These associations are well known. Advanced age and case severity are the most important predictors of early and late case fatality [39,40]. Living alone was also documented to increase long-term mortality after stroke [41].

Health-related quality of life
As to the effect of location of residence, we have found that, although dwellers of the socioeconomically disadvantaged District 8 scored less by both health utility measures compared to residents of wealthier regions (Table 1), similarly to case fatality, the number of patients was not sufficient for this difference to reach statistical significance. We have considered the rate of respondents acceptable as it was similar or higher compared to the proportions reported after stroke [42,43]. Nonetheless we found that the average EQ-5D utility index in both subgroups (0.65±0.34 for District 8 and 0.75±0.28 for other regions) are reduced compared to the population norm of 0.82 [44]. Using EQ-5D the health state utility values after stroke were relatively stable in time: in a systematic review of a total of 199 publications on stroke the median value was 0.63 in studies before 2013 versus 0.65 in studies after 2013 [45]. For 15D, the difference in health utility values between the socioeconomically disadvantaged district and the other areas also did not reach the level of statistical significance in our study (0.73±0.17 for District 8 and 0.78±0.16 for other regions). These values for 15D are within the range reported for stroke patients with or without vision problems (0.73-0.89); [46]. It has been reported that EQ-5D and 15D should not be used interchangeably in economic evaluations after stroke, as the utility scores generated from the two instruments, although correlated well, they differed significantly from each other [47]. Using the same health state utility measures 3 months after intracerebral hemorrhage, predictors for lower HRQoL by both scales were higher NIHSS and older age, with similar ORs for EQ-5D-5L, and 15D [48]. Minimally important differences (MIDs) for health state utilities vary by measure and methods, and are not well established [49]. It has been suggested that for EQ-5D differences among health state utilities of at least 0.036 can be considered clinically important [50]. For EQ-5D, the mean MID among non-stroke patients was 0.074 (range -0.011-0.140) [51]. In a Korean version of the EQ-5D-3L questionnaire, MID in stroke patients ranged from 0.08 to 0.12 [52], and a similar value of 0.10 was reported by Chen et al [53]. Considering the population norms of EQ-5D reported for Hungary (0.82) [44] the mean state utility value for EQ-5D in our patients (0.65 for District 8 and 0.75 for the wealthier regions) is considerably reduced, i.e. HRQoL is obviously adversely affected in stroke survivors 3 months after the acute event with a signal for reduced HRQoL in District 8. As far as 15D is concerned, there is a lack of information on MID values in stroke patients, and also no population norms are known for Hungary. In other patient groups MID for 15D was estimated between 0.01-0.03 [54]. Although by standard statistical methods-similarly to our findings for EQ-5D -the difference between health state utilities by the 15D was not statistically significant, the difference between the mean values of District 8 and other regions (i.e. 0.78-0.73 = 0.05) suggests that HRQoL after stroke is reduced more in socioeconomically deprived regions than in wealthier areas.
It is understandable that graver disability after stroke has a major impact on numerous aspects of the patient's physical, mental, professional and social life which might all be translatable to HRQoL. In the statistical models which excluded post-stroke variables, we have also found that admission stroke severity and advanced age both reduce post-stroke HRQoL [55,56]. What seems more interesting, though, is the effect of ischemic stroke subtype as to TOAST classification on HRQoL. TOAST classification as a possible predictor of outcome (death or dependency) was extensively studied: only lacunar stroke was found to influence outcome independently from stroke severity measured by the NIHSS [23]. There is scarcity of data regarding TOAST classification as a predictor of HRQoL. Our finding that cardioembolism was associated with better HRQoL than large-artery atherosclerosis is not in line with previous studies stating that cardioembolism predicts graver physical dependency [57].
Including all variables, we found several independent predictors of HRQoL: age, disability at discharge from acute management, stroke type, patient satisfaction, longer duration of acute and post-acute care and social dwelling type after stroke.
The condition of the patients at discharge from acute management expressed in mRS predicts HRQoL of patients expressed both in EQ-5D and 15D. The more disabled the patients at discharge, the lower the HRQoL of patients at follow-up. This association is quite comprehensible and has been already described [55,56]. Age is also a well-known predictor of post-stroke HRQoL [58].
The social type of dwelling of patients after stroke also seemed to have an impact on HRQoL at follow-up. Those patients able to live alone independently reported better HRQoL compared to those living with company, and especially compared to those requiring professional help. The finding, that the type of social dwelling itself is associated with a certain HRQoL outcome, independently of the known confounders like stroke severity or disability, is unexpected. The other relationship-onset stroke severity predicting both post-stroke HRQoL and dependence on others-is well-documented [55] and understandable. The social aspects of the dwelling circumstances after stroke might contribute to post-stroke HRQoL above the immediate effects of the stroke related disability on HRQoL.
A strong association between HRQoL and acute length of stay was already described [59,60]. In our stroke population, above length of acute stay, the duration of post-acute management (i.e. length of stay in rehabilitation or other hospitalization after discharge from the acute care setting) was also associated to poorer post-stroke HRQoL: longer than 5 days (especially more than 30 days) of post-acute hospital-stay had a negative impact on HRQoL. Thus, longer post-acute hospital stay represents not only a costly intervention, but could also contribute to reduced HRQoL of patients. The length of post-acute management can be an expression of graver disability, i.e. definite need for treatment prolongation. On the other hand, it can also be the expression of poorer social background, as patients who are not entirely independent after stroke and are lacking family help will have longer hospital-stays to enable arranging the dwelling conditions necessary after discharge. The association of stroke severity at onset with length of stay and with discharge disposition was already described [61,62]. Still, length of post-acute hospital-stay was not mentioned yet among the many predictors of post-stroke HRQoL, as far as we know [63,64]. There is a large body of evidence supporting the benefits of inpatient rehabilitation [65,66]. What is more important, Stroke Units' Trialists Collaboration [67] states that "people with acute stroke are more likely to survive, return home and regain independence if they receive stroke unit care . . . that can offer a substantial period of rehabilitation if required". In spite of the latter recommendation, a decline in the long-term HRQoL after inpatient rehabilitation of stroke patients was also reported [68]. Early Supported Discharge Trialists bring convincing arguments in favor of shortening the length of inpatient care in selected patients in order to reduce long-term dependency [69,70]. Further research is needed regarding post-acute management of stroke patients to assess the long-term outcome of inpatient and home-based rehabilitation. Our results suggest that shorter hospital stays followed by well-organized regular outpatient care could contribute to better post-stroke HRQoL. These outpatient interventions could be motor or non-motor rehabilitation sessions, regular home-based check-ups for improving the adherence to secondary prevention, nurses available for administering drugs and assisting the patients in their every-day grooming and medical attendance.
Although non-married status and poorer social support was linked to poorer rehabilitation potential after stroke [71,72], marital status was not among the predictors of HRQoL in our study.

Patient satisfaction
In our study, patients reporting better HRQoL also report better satisfaction with the care received and vice versa. There is evidence for each subcategory of patient satisfaction being influenced by the HRQoL, although in certain subcategories the evidence is weak. Our data suggest that patients with better HRQoL are more satisfied with physicians, nurses, organization of care and globally. Patient satisfaction was already known to be one of the most important predictors of HRQoL [73].
In addition to HRQoL, stroke severity also has a major impact on patient satisfaction: more severe strokes are associated with lower global satisfaction. On the other hand, the provision of thrombolysis had a significant positive impact on global patient satisfaction and the subcategory of satisfaction with nurses and care organization. We have not found any description of this association yet.
Patient satisfaction after stroke was not influenced by the rest of the numerous factors assessed. This is consistent with the known difficulty to influence patient satisfaction with direct measures [74,75] and questions patients' awareness of stroke management possibilities. Supposedly, patients socialized to have an impact on the care received would have valued their own satisfaction in a more differentiated manner. Patient satisfaction improved when they could actively influence the received care [76] and with better patient to nurse ratios [77]. Better patient satisfaction was also linked to outpatient rehabilitation [76][77][78], but in our study the patients discharged to their home, receiving professional help or living independently, have not valued the stroke care received any more than those discharged to rehabilitation or chronic care facilities.

Strengths and limitations of the study
The strengths of our study are the consecutive inclusion of patients thus reducing selection bias; the detailed systematic data collection, and the use of dual measures of HRQoL. Our study also has several limitations. First, the sample was relatively small, therefore does not have enough power to evaluate several possible predictors of HRQoL or patient satisfaction. Second, we have taken into account only the presence of pre-stroke depression but did not assess the patients for post-stroke depression, though according to published reports, depression has a significant impact on HRQoL and it can also influence negatively patient satisfaction [79]. Third, we have not assessed HRQoL in the acute phase of stroke, nor in a serial manner, therefore we cannot assess the change in HRQoL. Fourth, among those who did not return for an in-person check-up for the follow-up the rate of less satisfied patients may have been higher, and as the survey was done at 3 months after stroke, recall bias may have had an effect on patient satisfaction. Fifth, although EQ-5D is widely used, this might be a limitation considering its lower sensitivity in evaluating chronic diseases [30,32,33,47]. For this reason, we used 15D as well, adding higher sensitivity and reliability for evaluating HRQoL [30,[32][33][34]47].

Conclusions
Similarly to other studies, we found that initial stroke severity and age are the most important predictors of case fatality at 3 months after stroke. In addition to previously already identified predictors, we found that longer acute and post-acute hospital stay and the type of social dwelling were also associated with HRQoL of stroke survivors. These associations raise the issue whether long-term outcome could be improved by more intensive inpatient rehabilitation or rather by developing home care facilities. More research is needed regarding the effect of early discharge to patients' home supported by outpatient rehabilitation on HRQoL of the patients.