記住我
Musculoskeletal (MSK) conditions are highly prevalent worldwide, resulting in significant disability and suffering [], and were associated with up to US $380.9 billion of total medical expenditures in 2016 in the United States alone []. Exercise-based physical therapy is the mainstay treatment for such conditions [,-]. Recently, telerehabilitation and digital physical therapy have emerged as alternative care delivery systems for a wide range of MSK conditions [-]. These alternative care delivery systems have shown to be effective and feasible compared to traditional physical therapy [-] while increasing access and affordability to patients and easing the burdens of conventional programs []. Accessibility is increased by reducing travel limitations and time barriers while eliminating any geographic restrictions. Additionally, digital therapy can increase compliance by allowing patients to undergo treatment at their convenience and at their own pace, increasing patient empowerment and self-management [,,].
Despite the many benefits of telehealth, inequities in health other than underlying health status still exist based on age, geography, respective availability of health care facilities, and socioeconomic factors [-]. In fact, compared to those in urban areas, patients from rural areas tend to be older, are more likely to be obese, have higher rates of disability, have more chronic health conditions, and have higher fall rates [,]. Rural areas are known to have higher proportions of uninsured and underinsured individuals and higher costs of health care services when compared to urban areas []. Overall, 65% of rural US counties are designated as health professional shortage areas [], and rural areas have lower patient-to-primary care physician ratios []. Patients in rural areas of the United States have fewer opportunities for in-person physical activity programs due to limited access to indoor facilities and limited transportation when compared to urban patients []. These inequities are further compounded by lower educational levels, higher rates of poverty, and lower rates of internet access in rural areas [,]. Also, individuals with limited or no digital literacy or with limited access to digital technology may not have the means to pursue and maintain a telehealth intervention []. It is therefore crucial to identify strategies for improved access and quality of physical therapy in these historically disinvested areas.
To our knowledge, no study has been conducted on the impact of urban or rural location on engagement and clinical outcomes following a telerehabilitation program for MSK conditions. We have previously reported a multimodal digital care program (DCP) combining exercise-based physical therapy with psychoeducational components, which provided a comprehensive approach to pain management. This program encourages patients to develop strategies and self-management skills to manage their pain and has been validated in several acute and chronic MSK conditions [-,-]. Additionally, the impact of race and ethnicity [], as well as baseline mental health [] and fear-avoidance beliefs [], on final clinical outcomes have also been explored. The purpose of this study was to assess the impact of geographical location on engagement and clinical outcomes after a multimodal DCP, with the hypothesis being that patients from both rural and urban areas would have similar engagement and significant improvement in outcomes after program completion.
This study is an ad hoc analysis of a decentralized, single-arm investigation into clinical and engagement-related outcomes following a multimodal DCP in patients with musculoskeletal (MSK) pain conditions. The DCP was administered at the patients’ homes and delivered between March 1, 2021, and March 10, 2022.
Ethics ApprovalThis study is part of a trial that was prospectively registered on ClinicalTrials.gov (NCT04092946) on September 17, 2019, and approved by the New England Institutional Review Board (120190313) on June 18, 2020.
PopulationThe study population included adults (≥18 years of age) who were beneficiaries of employer health plans from 50 US states and the District of Columbia. Employees and their dependents who reported either acute or chronic MSK pain in the spine, upper limbs, or lower limbs were eligible and were invited to apply to the DCP of Sword Health (located in Draper, Utah) through a dedicated website. Throughout enrollment, participants were asked to provide demographic data, including zip codes and baseline clinical information (eg, initial pain levels). Participants were informed about the study and invited to provide consent. The exclusion criteria were as follows: (1) a health condition (eg, cardiac or respiratory) not allowing a participant to engage in at least 20 minutes of light to moderate exercise, (2) being under treatment for active cancer, and (3) rapid loss of strength or numbness in the arms or legs or change in bowel or urinary function in the previous 2 weeks.
InterventionThe DCP has been described previously [-,-]. In brief, this multimodal program consisted of 4-, 8-, or 12-week telerehabilitation interventions comprising exercise, education, and cognitive behavioral therapy (CBT). This program digitally interfaced between the patient and an assigned physical therapist (PT), who monitored the patient for the study duration. Participants who lacked internet access at home were given a Wi-Fi hotspot. A US Food and Drug Administration–listed class II medical device that consisted of inertial motion trackers, a mobile app in a dedicated tablet, and a cloud-based portal was made available to all patients. Briefly, the personalized exercises were displayed on the tablet, with trackers allowing real-time video and audio biofeedback on performance. At session end, the data related to the exercise sessions, such as compliance, presence or absence of movement errors, and level of pain and fatigue during the exercise, were registered and stored in a cloud-based portal. This portal enabled remote and asynchronous monitoring by the assigned PT, who revised the prescribed exercises if needed. Patients were recommended a frequency of 3 exercise sessions per week. The education and CBT components of the program were developed by a multidisciplinary team following current clinical guidelines and state-of-the-art research [-]. The education component delved into topics focused on anatomy, physiology, symptoms, evidence-based treatments, fear avoidance, and active coping skills (including managing feelings of anxiety and depression). The CBT program was based on third-generation techniques—mindfulness, acceptance, and commitment therapy; empathy-focused therapy; fear-avoidance behavior; and constructive coping. The education and CBT materials were delivered to the patients through written articles, audio content, and interactive modules. Bidirectional communication with the assigned PT was ensured through built-in secure chat within a smartphone app and through video calls. Participants who did not perform any exercise session for 28 consecutive days were considered dropouts.
Demographic DataDemographic data included age, BMI, patient gender, educational level, and employment status. The gender category included “man,” “woman,” “nonbinary,” “other,” and “prefer not to specify.” The employment status categories were defined as the following: full-time employed, part-time employed, or not employed. The educational levels were (1) high school or less (including technical or vocational training), (2) some college, including a bachelor’s degree, community college, or an associate degree, (3) some graduate school, including a master’s or doctoral degree, and (4) “not available” or “prefer not to answer.”
Patients were coded according to their zip codes to a specific rural-urban commuting area (RUCA) code []. RUCA codes characterize all census areas regarding their rural and urban status and relationships. This classification system uses the standard Bureau of Census urbanized area and urban cluster definitions in combination with work-commuting information [,]. Rural areas have been defined as having an urban core of 50,000 people or less []. Therefore, using primary RUCA codes, we defined urban areas by scores from 1 to 3, and rural areas by aggregating codes 4 to 10 (, Table S1 []).
OutcomesOutcomes were collected at baseline and 4, 8, and 12 weeks, and mean changes were calculated between baseline and program end. Engagement and clinical outcomes are described in .
Table 1. Engagement and clinical outcomes in this study.OutcomeDescriptionEngagementMeasured through the following:aWPAI: Work Productivity and Activity Impairment.
Statistical AnalysisAnalyses of baseline characteristics (demographics and clinical data), as well as engagement metrics, were performed using a 2-tailed, 2-sample t test, a 2-way analysis of variance (ANOVA) with Bonferroni post hoc test, a chi-square test, or a 2-proportion z test. Patients who completed the program were defined as “completers” and those that did not were defined as “noncompleters.”
Latent growth curve analysis (LGCA) was used to estimate trajectories of each outcome over time [,]. The LGCA has been recognized as one of the most powerful methods to analyze longitudinal data, since it provides a measure of fitness and addresses missing data through full information maximum likelihood (FIML) [-]. FIML uses all available data at each time point from all participants to calculate maximum likelihood estimates, outperforming multiple imputation by chained equations or listwise deletion [,]. In addition, the LGCA uses a structural equation model to define trajectories through intercept, slope, and curvature for each variable, allowing analysis of the recovery pace and leveling of the effect for each outcome. In order to account for unbalanced group sizes, a multiple-group LGCA was conducted. This allows for creating separate models for rural and urban groups while simultaneously performing intergroup comparisons (eg, mean change). A conditional analysis was conducted to assess the influence of age, gender, BMI, employment status, and education level and was fitted as a random effect. Additionally, analysis of subpopulations was performed by focusing on participants who met the following criteria at baseline: Generalized Anxiety Disorder 7-item (GAD-7) and Patient Health Questionnaire 9-item (PHQ-9) scores equal or greater than 5 points [,] and a Work Productivity and Activity Impairment (WPAI; comprising overall, work, time, and activity) score greater than 0 points. A robust sandwich estimator was used in all models for standard errors.
All statistical analyses were conducted using commercially available software (SPSS version 22; IBM Corp), and the level of significance was set at P<.05 for all tests. The LGCA was coded using R (version 4.2.2; R Foundation for Statistical Computing).
A total of 14,754 participants were screened for eligibility (). Of these, 2151 were excluded, for a total of 12,603 (85.4%) eligible patients, of whom 2611 were excluded due to unavailable RUCA data or not starting the program, resulting in a total of 9992 patients at program start. A total of 7378 of 9992 (73.8%) patients completed the program. The study flow diagram is presented in .
Figure 1. Study flow chart showing the number of participants who were excluded, included, and dropped out. RUCA: rural-urban commuting area. Baseline CharacteristicsPatients’ baseline demographics grouped by urban and rural areas are presented in , while baseline characteristics stratified by completers and noncompleters can be found in , Table S2.
Table 2. Baseline characteristics for urban and rural groups following an intention-to-treat analysis. Filtered cases correspond to participants who reported relevant impairment at baseline (>0 or ≥5 points). Statistically significant P values are italicized.aA total of 1.1% (114/9992) of patients were postsurgical.
bGAD-7: Generalized Anxiety Disorder 7-item scale.
cPHQ-9: Patient Health Questionnaire 9-item scale.
dWPAI: Work Productivity and Activity Impairment scale.
Patients from rural areas were significantly older than patients from urban areas (51.85, SD 12.57 years vs 48.11, SD 12.37 years, respectively; P<.001), had higher BMI (30.83, SD 7.48 kg/m2 vs 28.96, SD 6.60 kg/m2, respectively), had a higher percentage of unemployed individuals (214/1183, 18.1% vs 853/8809, 9.7%, respectively; P<.001), and had a higher percentage of individuals with a lower educational level (166/1183, 14% vs 700/8809, 7.9%, respectively; P<.001). No statistically significant differences in gender distribution were found between groups.
Patients in rural areas also presented a higher prevalence of chronic pain and low back pain conditions than patients from urban areas (535/1183, 45.2% vs 3441/8809, 39.1%, respectively). In opposition, a higher prevalence of knee-related conditions was observed in patients from urban areas than from rural areas (1292/8809, 14.7% vs 146/1183, 12.3%, respectively; P<.001).
Overall, similar clinical metrics were observed between patients from rural and urban areas at baseline. The statistically significant differences found were in baseline depression (rural PHQ-9 score 2.70, SD 4.41 vs urban PHQ-9 score 2.33, SD 4.22; P=.004) and percentage of patients with at least mild or moderate depression at baseline (281/1183, 23.8% vs 1790/8809, 20.3%, respectively; P=.006), where more patients from rural areas were depressed.
When comparing completers with noncompleters, the latter were younger (46.23, SD 12.61 years vs 49.38, SD 12.29 years, respectively; P<.001) and reported higher BMI (30.22, SD 7.36 kg/m2 vs 28.81, SD 6.47 kg/m2, respectively; P<.001). A larger proportion of noncompleters were employed full-time (1732/2614, 66.3% vs 4539/7378, 61.5%, respectively; P<.001) and reported a lower educational level. Noncompleters also reported higher levels of impairment in productivity (WPAI overall, P=.004 and WPAI work, P=.001) and non–work-related activities (WPAI activity, P=.01) at baseline.
EngagementIndividuals from rural areas were more likely to complete the program than patients from urban areas (906/1183, 76.6% vs 6472/8809, 73.5%, respectively; P=.02). However, independently of dropout rates, both groups had a similar pattern of engagement. Engagement data stratified by patients from urban and rural areas is presented in . The 2 groups had similar time dedicated to exercise (P=.48), number of sessions (P=.77), sessions per week (P=.11), and interactions with the PT (P=.14). Average satisfaction scores were similarly high in both groups (rural score 8.6, SD 1.7 and urban score 8.6, SD 1.8; P=.95). The single significant difference in engagement between groups was the number of educational articles consulted, with patients from rural areas reading more articles than patients from urban areas (P<.001; ).
Table 3. Engagement data across the groups. Statistically significant P values are italicized.Engagement outcomesUrban, mean (SD)Rural, mean (SD)P valueSessions, n33.88 (32.23)34.18 (30.94).77Sessions per week, n2.76 (1.14)2.78 (1.12).11Training time, minutes472.02 (485.56)482.54 (485.96).48Articles read, n2.72 (5.27)3.44 (6.18)<.001Interactions with physical therapist, n11.79 (12.54)12.37 (13.90).14Average satisfaction score8.6 (1.7)8.6 (1.8).95Clinical OutcomesMean changes in clinical outcomes for both urban and rural groups following an intent-to-treat analysis are presented in , while the corresponding model estimates and model fitness are presented in , Tables S3 and S4, respectively [,]. Change trajectories of each outcome are depicted in . The impact of the covariates in clinical outcomes is presented in , Table S5. The same analysis following a per-protocol approach is presented in , Tables S6-S8. Similar results were observed from both intention-to-treat and per-protocol approaches. Since intention-to-treat analysis offers an overview change of the entire cohort, the following section is focused on these results.
Table 4. Mean changes between baseline and program end and mean differences between groups for the studied clinical outcomes following an intent-to-treat analysis. Statistically significant P values are italicized.ScoresUrbanRuralMean differenceaGAD-7: Generalized Anxiety Disorder 7-item scale.
bPHQ-9: Patient Health Questionnaire 9-item scale.
cWPAI: Work Productivity and Activity Impairment scale.
Figure 2. Longitudinal changes across time using intent-to-treat analysis. A: pain level; B and C: mental health (GAD-7-and PHQ-9 scores, respectively) for cases with at least mild or moderate anxiety or depression at baseline; D-F: work productivity (WPAI overall, WPAI work, and WPAI activity scores, respectively) for cases reporting impairment at baseline. The shaded areas are the 95% CI. GAD-7: Generalized Anxiety Disorder 7-item scale; PHQ-9: Patient Health Questionnaire 9-item scale; WPAI: Work Productivity and Activity Impairment scale. PainPain levels decreased similarly in both groups from baseline to program end (P=.62), corresponding to a significant mean change of 2.2 points (95% CI 2.2-2.3) in the urban cohort and 2.3 points (95% CI 2.1-2.5) in the rural cohort ( and A). The rate of patients meeting pain MCID at program end was not statistically different between urban and rural cohorts (67.1% versus 68.3%, P=.30).
Mental HealthAmong those who reported at least mild or moderate anxiety (GAD-7 score ≥5) at baseline, we observed a significant decrease in anxiety in both groups (urban score 4.5, 95% CI 3.7-5.4; P<.001 and rural score 4.6, 95% CI 4.3-4.9; P<.001); this decrease was similar in the groups (P=.85). Likewise, among those reporting at least mild or moderate depression at baseline (ie, PHQ-9 score ≥5), we observed a significant decrease in depression scores in both groups (urban score 4.5, 95% CI 3.39-5.53; P<.001 and rural score 4.9, 95% CI 4.5-5.2; P<.001); this decrease was also similar in the groups (P=.48). In patients from urban areas, chronic MSK conditions were associated with steeper initial recovery of depression, followed by a stronger leveling effect. In contrast, in patients from rural areas, chronic pain was associated with slower improvement that was more sustained over time (C, difference in slope between groups –0.07, P=.004; difference in curve between groups 1.06, P<.001; , Table S5).
ProductivityProductivity improvements were observed in both groups with no differences between them (D-F). For WPAI overall, mean changes of 14.6 points (95% CI 11.7-17.4) and 15.6 points (95% CI 14.5-16.7) were reported for urban and rural groups, respectively (P=.50). Presenteeism, measured through WPAI work, improved by 14.3 points (95% CI 11.6-17.0) in the urban group and by 15.4 points (95% CI 14.3-16.4) in the rural group (P=.46). Absenteeism, measured through WPAI time missed, was reduced by 11.4 points (95% CI 6.9-16.0) and 11.8 points (95% CI 9.9-13.7) in the urban and rural groups, respectively. Similarly, impairment in non–work-related activities (ie, WPAI activity) had improvements of 19.4 points (95% CI 18.6-20.3) and 18.4 points (95% CI 16.1-20.7) in the urban and rural groups, respectively.
In urban areas, individuals with higher BMI reported a greater leveling effect on absenteeism improvement than those from rural areas (the difference in curve between groups was 0.14, P=.04; , Table S5). In patients from rural areas, the presence of chronic pain was associated with a faster recovery in absenteeism in comparison with that in patients from urban areas (the difference in slope between groups was 0.16, P=.03; , Table S5).
The multimodal DCP herein reported was able to reach all US states in both urban and rural locations and had a completion rate of 73.8% (7378/9992), which is similar to previous studies reporting the use of digital interventions for MSK pain management [,]. The percentage of participants in urban areas was significantly higher than in rural locations, which was expected given that only approximately 19.6% of the US population is located in rural areas according to the US Census Bureau [].
Health inequities between urban and rural populations are prevalent in the United States. []. Rural populations have been reported to have demographics associated with a poorer prognosis for MSK pain [-]. In accord with this, in this study, patients from rural areas were older [], had higher BMI [], had lower educational levels [], and had a higher prevalence of depression []. Some of these factors have previously been associated with lower chances of receiving physical therapy []. Lower education (and consequently lower digital literacy) have been associated with lower adherence, for example [-].
Studies have shown that patients in rural areas of the United States may face additional difficulties in recovery due to fewer opportunities for in-person physical activity programs as a consequence of limited access to indoor facilities, limited transportation, and a lower overall health status when compared to urban patients [,]. Additionally, rural residents are less likely to report having home broadband than those living in urban or suburban areas [], which seriously impacts their access to digital health care tools and electronic communication with health providers [].
In this study, engagement was similar between both rural and urban areas (eg, the number of sessions and interactions with a PT), and completion rates were higher in the rural cohort. The reasons behind these observations may be multifactorial, but one can speculate that the lack of access to alternative health care resources, as well as the provision of a Wi-Fi hotspot to those without internet, might have prompted patients from rural areas to not only engage with the exercise sessions but also to achieve higher completion rates [,]. Also, despite lower educational levels, patients from rural areas engaged more with curated health educational articles advocating for telerehabilitation programs as enablers of health literacy.
Despite the worse clinical outcomes reported at baseline by those in rural communities, in line with what has been described before [,-], similar improvements in pain, mental health, and productivity impairment were observed in both groups, again reinforcing the notion that higher MSK pain burden in rural areas may be associated with lack of access to care. Pain improvements were above a 2-point change independently of the studied group, with 67.1% to 68.3% of participants meeting the MCID for pain [,]. The percentages of patients meeting the MCID were within the ranges previously reported for digital interventions (49%-75.6%) [-] and in-person physical therapy [].
The prevalence of depression and anxiety has been reported to be higher in residents of rural areas compared to urban areas [], with those from rural areas facing a shortage of mental health services []. Since mental health and MSK pain are tightly associated [,], the scarcity of psychological support can seriously impact the recovery rates of rural populations. This study confirms a higher depression burden in patients from rural areas but found similar improvement in mental health scores in both rural and urban patients, reinforcing the notion that lack of access to mental health resources may be the main driver for the higher burden of disease in rural areas. Additionally, MSK pain has been reported to be a main driver for loss of work productivity [,]. The factors weighing on absenteeism recovery were BMI and the presence of chronic pain, both previously reported to negatively impact MSK pain recovery [,]. Nevertheless, we did not observe significant differences in improvement in any productivity domains between urban and rural groups; all of these domains showed significant improvement following the DCP.
Despite the wide reach of telerehabilitation, many areas across the United States are still facing unmet needs. The results observed herein support the need for further research and investment in digital rehabilitation to mitigate inequities in health care access and care delivery optimization.
Strengths and LimitationsThere are many strengths to this study, namely the novelty of investigating the urban-rural dichotomy within a digital therapy program in a large sample size from a real-world context, including patients from 50 US states and the District of Columbia, which allows for a diverse population and thus better generalizability. Another strength is the DCP itself, which uses a multimodal approach that includes exercises with real-time biofeedback, mental support, regular communication with the PT, and a digital format. All these components favor accessibility and maximize engagement and clinical outcomes, allowing us to study different aspects of the problem, from pain to mental health to productivity.
The classification of rural and urban areas is a challenging topic considering the multitude of factors that can highly influence the obtained readings. Despite the application of a recognized classification system [-], we cannot rule out the existence of other confounding factors with contributions that were not taken into account during this exploratory analysis, including desirability bias. Other limitations include the lack of control groups (to account for nonspecific treatment effects) and the lack of long-term outcome and objective outcome measures (ie, through activity trackers). Nevertheless, this exploratory study may lay a foundation for future work in this field, identifying areas in need of improvement for future telerehabilitation programs. Further prospective controlled studies are warranted to better characterize the effect of rural and urban inequities on digital therapy outcomes.
ConclusionThis study provides important insights regarding the impact of a multimodal digital program for MSK pain management in rural and urban settings. The DCP was able to reach all areas across the United States with high completion rates in both settings. Despite the inherent health inequities between patients from rural and urban areas, similarly high satisfaction and engagement, alongside significant improvements in pain, mental health, and productivity, were observed in both groups. This showcases the potential of the DCP to mitigate inequities by improving the accessibility of MSK care independently of geographic location.
The authors acknowledge the team of physical therapists responsible for managing the participants. The authors also acknowledge the contributions of João Tiago Silva and Guilherme Freches in data validation (both employees of Sword Health). Critical revision of the manuscript for important intellectual content was done by all authors. All authors were involved with the final approval of the manuscript. The study sponsor, Sword Health, was involved in the study design, data collection, and interpretation and writing of the manuscript.
The data sets generated during or analyzed during this study are available from the corresponding author upon reasonable request.
All authors made a significant contribution to the work. FC, FDC, and JL were responsible for the study concept and design. MM acquired the data. RM performed the statistical analysis. JS, FC, ACA, DJ, MM, and FC interpreted the data. JS was responsible for drafting the work. VB was responsible for funding.
Sword Health owns the Digital Care Program and as such this program is only available to those covered by a commercial agreement with Sword Health. FC, DJ, AA, MM, FDC, VB, and VY are employees of Sword Health, the sponsor of this study. DJ, FC, FDC, VY, and VB also hold equity in Sword Health. RM is an independent scientific consultant responsible for the statistical analysis. JS and JL are independent scientific and clinical consultants who received adviser honorariums from Sword Health.
Edited by G Eysenbach, L Buis; submitted 17.11.22; peer-reviewed by R Jamison; comments to author 05.01.23; revised version received 12.01.23; accepted 02.02.23; published 16.03.23
©Justin Scheer, Anabela C Areias, Maria Molinos, Dora Janela, Robert Moulder, Jorge Lains, Virgílio Bento, Vijay Yanamadala, Fernando Dias Correia, Fabíola Costa. Originally published in JMIR mHealth and uHealth (https://mhealth.jmir.org), 16.03.2023.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR mHealth and uHealth, is properly cited. The complete bibliographic information, a link to the original publication on https://mhealth.jmir.org/, as well as this copyright and license information must be included.
留言 (0)