Identifying Risk of Depressive Symptoms in Adults With Physical Disabilities Receiving Rehabilitation Services: Propensity Score Approaches
Article information
Abstract
Objective
To compare the risk of depressive symptoms in adults with physical disabilities between rehabilitation-receivers and non-receivers.
Methods
A total of 3,568 adults with physical disabilities were retrieved from the 2014 Korean National Survey on People with Disabilities database. Four covariate adjustment methods (a multivariable regression model, inverse probability of treatment weighting [IPTW] adjusted for normalized weight, IPTW with stabilized weight, and greedy algorithm with 1:1 propensity score matching) were used to estimate the odds of having depressive symptoms. The dependent variable was depressive symptoms and the independent variable was the use of rehabilitation services. Baseline covariates were 19 demographic variables and 10 chronic condition variables.
Results
The four covariate adjustment methods revealed that adults with physical disabilities receiving rehabilitation services had a higher risk of depressive symptoms than those who did not receive these services (adjusted odds ratio, 1.191–1.294).
Conclusion
Our findings suggest that adults with physical disabilities receiving rehabilitation services have higher risk of developing depressive symptoms. Therefore, rehabilitation professionals need to pay attention to depressive symptoms and establish therapeutic strategies that can reduce such risk in rehabilitation settings.
INTRODUCTION
Physical disability defined as a condition in which an individual experiences deviation or loss in their body function or structure [1] can result in difficulty, restriction, or dependence in performing basic and instrumental activities of daily living (ADLs and IADLs) [2,3]. In South Korea, with increasing number of people with physical disabilities [4], the proportion of population with physical disabilities receiving rehabilitation services is also increasing. Therefore, developing a deeper understanding of people with physical disabilities is indispensable for rehabilitation professionals to provide more successful interventions for this population.
The International Classification of Functioning, Disability, and Health (ICF) based on a bio-psycho-social model recommends that rehabilitation professionals need to consider not only physical functions, but also psychological symptoms such as depressive symptoms or anxiety to better understand people with physical disabilities and deliver better interventions [1]. To be specific, depressive symptoms can exacerbate individuals’ physical disability symptoms [5] which may result in low health-related quality of life [6] and higher healthcare expenditures [7]. Considering adverse effects of depressive symptoms, rehabilitation professionals should monitor and track the risk of depressive symptoms and understand impacts of rehabilitation services on depressive symptoms to successfully manage health of people with physical disabilities.
Previous studies have acknowledged the importance of identifying the risk of depressive symptoms and reported a higher prevalence rate of depressive symptoms or depression in people with disabilities compared to the general population [5,8,9]. These previous studies have revealed that disabilities can increase the risk of depressive symptoms itself. To investigate whether receiving rehabilitation services can affect the risk of depressive symptoms among adults with physical disabilities, it is imperative to compare the risk of depressive symptoms between rehabilitation receivers and non-receivers within the same patient group. However, to the best of our knowledge, no research has compared the risk of depressive symptoms between these two groups.
To accurately compare the risk of depressive symptoms, risk factors for depressive symptoms should be controlled to minimize effects of confounding factors. For instance, there are well-known risk factors for depressive symptoms or depression, including low income [10], social support [11], osteoarthritis [12], low back pain [13], age [14], and sex [15]. Among patients receiving rehabilitation, depressive symptoms are also associated with higher comorbidity and impairment, lower social support, lower self-efficacy, and specific work-related problems [16]. Since previous studies have demonstrated various factors associated with depressive symptoms, it is important to effectively control for those factors related to depressive symptoms when investigating the risk of depressive symptoms.
The ICF model is an effective model for understanding and controlling for related risk factors of depressive symptoms. According to the ICF model, depressive symptoms are defined as impairments in an individual’s emotional function which belongs to the category of body function. Since impairments in body function can be affected by various factors such as health condition, personal and environmental factors, activities, and participation, related factors of depressive symptoms can be categorized according to the ICF model. Thus, information about relationships between depressive symptoms and its risk factors identified from the ICF model will be useful when estimating the risk of depressive symptoms by controlling for these risk factors.
Methodologically, propensity score (PS) matching methods are robust statistical methods effectively accounting for covariates. Although traditional regression adjustment methods are effective for controlling for confounders, these methods have limitations in addressing selection bias which is a typical limitation of observational studies [17-20]. In contrast, PS matching methods can mimic a randomized clinical trial by matching two comparison groups using conditional probability of receiving exposure given a set of covariates. For instance, various PS matching methods have been demonstrated to be robust after controlling for selection bias and confounders, such as inverse probability of treatment weighting adjustment (IPTW), IPTW with stabilized weight, and 1:1 matching [18-22].
Therefore, the objective of this study was to compare the risk of depressive symptoms between adults with physical disabilities who received rehabilitation services and those who did not receive these services by utilizing four different covariate adjustment methods.
MATERIALS AND METHODS
Study sample
We retrieved records of 104,703 adults with physical disabilities from the 2014 Korea National Survey on People with Disabilities (NSPD) database [23]. The inclusion criteria for our study were: (1) the presence of physical impairments in the upper and/or lower extremities and/or spinal cord injury and (2) age older than 20 years. Physical impairments were defined as the state of being physically limited due to amputation, a medical disorder, or impairments of joints, extremities, or the spinal column. The study data set is de-identified and publicly available.
Definition of depressive symptoms and rehabilitation services
The presence of depressive symptoms was defined as a depressive mood or hopelessness that continued for 2 or more weeks and affected the participant’s daily activities during the past 12 months (September 2013 to August 2014). The survey participants were also asked if they had received rehabilitation services (treatment group vs. control group), including physical, occupational, speech, music, play, art, behavioral, and/or other therapies (e.g., exercise, horticultural therapy, therapeutic riding, and pottery psychotherapy). In other words, if participants received any of these rehabilitation services, we considered them to be receiving rehabilitation services.
Covariates based on the ICF model
Included covariates were categorized based on the ICF model.
Health condition
Health condition includes issues associated with chronic or acute illnesses [1]. In this study, history of falls, self-rated health (1 ‘very poor’ to 5 ‘very good’), and chronic conditions present for 3 months or more were included as factors for health condition.
Body function
Body function was classified as physiological functions of body systems according to the ICF model [1]. We selected physical function score as a factor for body function. It was estimated using an ADL instrument validated by an item response theory model [24].
Activity
The ICF model defines activity as the execution of a task or action by an individual [1]. Participation in physical activity (at least 3 days/week) was defined as a factor of activity.
Participation
Participation refers to involvement in life situations as noted in the ICF [1]. Disability grades (1 ‘severe’ to 6 ‘mild’) assigned by the Korean government include assessment of participation levels. Thus, disability grade was categorized as a factor for participation.
Personal factors
Personal factors are defined as characteristics of an individual, including race, age, and health conditions [1]. Thus, demographic variables such as age, years since primary impairment onset, sex, marital status (single or married/living with a partner), employment (currently working), low income (less than $1,388/month), education, smoking, and drinking were included as personal factors.
Environmental factors
The ICF model defines environmental factors as physical, social, and attitudinal environments in which people live and interact, including personal care providers and personal assistants, health professionals, and education and training services [1]. Thus, in this study, outpatient visits (past 2 weeks including local hospitals, health centers, or dental clinics), preventative health services (education programs related to exercise, nutrition, smoking, and drinking in the past 12 months), mental health services (in the past 12 months), currently receiving depressive symptoms treatment, and social support (living with family member(s) or partner) were defined as environmental factors.
Statistical analysis
Descriptive analyses were conducted to describe demographic characteristics and covariates. To compare the two groups (rehabilitation receivers and non-receivers), chi-square test was used for categorical variables and Wilcoxon rank sum non-parametric test was used for continuous variables. Odds of depressive symptoms were estimated by using the four covariate adjustment methods (a traditional multivariable logistic regression and three PS matching methods) to examine the robustness of study results. The three PS matching methods were: (1) IPTW using propensity score with normalized weight, (2) IPTW with stabilized weight, and (3) greedy algorithms with 1:1 matching of the PS [18,21,25].
First, we used a complete pooling approach in PS analyses [26]. PS was estimated based on the likelihood of having depressive symptoms from a logistic regression model that accounted for participant’s characteristics (Table 1). Normalized weights for the IPTW matching method were applied to avoid extreme values of weight by dividing each individual propensity score by the mean of all propensity scores. This allows a comparison to samples representative of specific populations [18]. Therefore, normalized weight in the IPTW matching method does not result in a loss of observations.
While the robustness of IPTW with normalized weight has been demonstrated in observational studies, the null hypothesis is frequently rejected due to inflated sample size [22]. In other words, this method often produces large variance estimates that result in high type I error rates. Thus, we also used IPTW with stabilized weight method to account for high propensity scores in treatment and control groups and validated point estimation from normalized weight for the IPTW method.
The most common PS matching method is the 1:1 matching method [18,19]. In this study, we used a greedy 1:1 matching algorithm without replacement using a caliper of width equal to 0.10 of the standard deviation (SD) of the logit of complete pooling PS [25]. Consequently, this method allows matched sets of survey participants between treatment and control groups to share similar PS values estimated using covariates listed in Table 1.
After performing the three PS matching methods, we re-calculated descriptive statistics between the two comparison groups to examine if covariates listed in Table 1 were balanced. We also performed two PS diagnostic tests for the 1:1 matching method, including absolute standardized differences in means of all covariates [19,27] and percentage bias reduction of covariates [26]. We considered an absolute standardized difference of less than 0.2 in means of all covariates as indicative of a good match between the two groups [27]. If there were unbalanced covariates from the three PS matching methods, we additionally controlled for these covariates using logistic regression models. Lastly, we conducted sensitivity analysis to examine the magnitude of hidden bias in the risk of depression symptoms caused by unobserved covariates [18,28,29]. Sensitivity analysis result was reported with a value of gamma at two-sided α=0.05 level. All analyses were performed using SAS version 9.4 software (SAS Institute Inc., Cary, NC, USA). Statistical significance was determined at an alpha level of 0.05 in two-sided tests.
RESULTS
We analyzed records of a total of 3,568 adults with physical disabilities. The average age was 62.7 years (SD=13.6). Among retained study sample, 775 participants were identified as having depressive symptoms, including 312 (27.3%) in the treatment group and 463 (18.9%) in the control group. The majority (n=1,126, 99.8%) of patients in the treatment group received physical therapy. Baseline covariate differences between treatment and control groups before matching are presented in Table 1. Before matching, 23 out of 29 covariates were significantly different between the treatment group (n=1,126, 31.5%) and the control group (n=2,442, 68.5%).
PS ranged from 0.049 to 0.841 in the treatment group and from 0.031 to 0.834 in the control group (Appendix 1). Baseline covariate differences between the two groups after performing the three matching methods are shown in Table 2. After applying the IPTW with normalized weight and stabilized weight methods, all baseline covariates were balanced between the two comparison groups (all p>0.05).
After applying the 1:1 PS matching method, seven covariates were significantly different between treatment and control groups, including age, marital status, low income, current smoker, social support, osteoarthritis, and back pain (all p<0.05). However, these variables satisfied the criterion of minimum reduction in absolute standardized differences (Fig. 1).
Total percentage bias reduction of covariates by the 1:1 PS matching method was 75.8%. We controlled for these seven significantly different variables when estimating odds of depressive symptoms. Table 3 shows risks of depressive symptoms in adults with disabilities who have received rehabilitation services after applying covariate adjustment and the PS method plus further covariate adjustments.
Unadjusted risk of depressive symptoms was 63.9% (odds ratio=1.639; 95% confidence interval [CI], 1.389–1.933). However, this risk decreased to 19.1%–29.4% according to the four covariate adjustment methods, including a traditional multivariable logistic regression (adjusted odds ratio [AOR]=1.233; 95% CI, 1.013–1.501), IPTW normalized with logistic regression adjustment (AOR=1.191, 95% CI, 1.014–1.400), IPTW stabilized weight with logistic regression adjustment (AOR=1.191; 95% CI, 1.002–1.417), and 1:1 PS matching method with logistic regression adjustment (AOR=1.294; 95% CI, 1.065–1.574). Sensitivity analysis for the matched sample resulted in a gamma value of 1.11 at two-sided α=0.05 level, indicating that unobserved covariates could change our study results.
DISCUSSION
In this study, we compared the risk of depressive symptoms in adults with physical disabilities who received rehabilitation services to those who did not receive such services by utilizing four different covariate adjustment methods (a traditional regression model and three different propensity score matching methods). Our study findings revealed that adults with physical disabilities receiving rehabilitation services had higher risk of developing depressive symptoms (AOR, 1.191–1.294). To the best of our knowledge, this is the first study to compare the risk of depressive symptoms between rehabilitation receivers and non-receivers among the same adult population with physical disabilities. For instance, a recent study in Germany has investigated the prevalence of depression in rehabilitation settings (i.e., neurological, orthopedic, cancer rehabilitation, and others) and reported a higher prevalence of depression in those who have received rehabilitation than healthy individuals [16]. However, this previous study only compared people with disabilities to healthy individuals.
In a rehabilitation setting, patients encounter unavoidable physical and emotional challenges such as pain, setbacks, and frustration [30]. In particular, if people experience permanent physical disabilities, they could be exposed to these types of personal challenges while receiving rehabilitation. For example, patients with physical disabilities are likely to experience challenges such as frustration and depression when their recovery does not progress as rapidly as they expect [30]. These challenges might increase the risk and likelihood of depressive symptoms in rehabilitation receivers. Our findings suggest that rehabilitation receivers require special care for depressive symptoms.
The ICF model proposes that rehabilitation professionals should follow a multidisciplinary approach to determining the best options for patients [1]. However, in South Korea, rehabilitation services are primarily designed to provide services targeting physical functions or functional performance, not psychological symptoms such as depressive symptoms. This current rehabilitation practice is partially due to the current reimbursement policy and the lack of attention to or care for depressive symptoms in rehabilitation settings. Nevertheless, rehabilitation professionals can indirectly provide therapeutic interventions such as recreational activities to ameliorate clients’ depressive symptoms. Therefore, the current reimbursement policy needs to be revised for practitioners to administer systematic professional therapeutic strategies for people with psychological symptoms.
It is important to note that the identified risk of depressive symptoms in adults with physical disabilities who received rehabilitation services might stem from particular depressive symptoms covariates. After conducting the 1:1 PS matching method, the seven covariates (age, sex, low income, current smoker, social support, osteoarthritis, and back pain) showed significant differences between rehabilitation receivers and non-receivers within adults having physical disabilities (all p<0.05). All of these seven uncontrolled covariates are associated with depressive symptoms in adults with physical disabilities. Numerous studies examining risk factors for depressive symptoms support these findings [10-15]. While some of these risk factors are non-modifiable personal factors (age, sex, low income), other relevant variables are potentially modifiable factors, such as environmental factors (social support) and health condition (osteoarthritis and back pain). Therefore, rehabilitation professionals might need to address modifiable risk factors to reduce the risk of depressive symptoms in rehabilitation settings.
Of note, gamma level of 1.11 was relatively small. In other words, there are high probability that hidden bias can affect our results. Thus, a prospective study that controls for various depressive symptom covariates that are currently unavailable in this study is needed in the future to validate our current study findings.
In conclusion, Korean adults with physical disabilities who receive rehabilitation services show risk of developing depressive symptoms. Because depressive symptoms can cause a myriad of challenges for both the society and individuals, depressive symptoms in rehabilitation receivers should be carefully monitored to provide appropriate interventions.
Notes
No potential conflict of interest relevant to this article was reported.
Acknowledgements
This research was supported in part by grant (No. K12 HD055929) from the National Institutes of Health (NIH). Its contents are solely the responsibility of the authors. They do not necessarily represent official views of the NIH.
Notes
Conceptualization: all authors. Methodology: all authors. Formal analysis: Hong I. Funding acquisition: Hong I, Park HY. Project administration: Park HY. Visualization: Hong I, Lee MJ. Writing – original draft: Lee Y, Hong I. Writing – review and editing: all authors. Approval of final manuscript: all authors.