Abstract
Objective
Identify subgroups of postoperative older adults using a latent class analysis (LCA) of routinely collected patient data, including preoperative frailty and postoperative mobility.
Methods
Retrospective analysis of routinely collected data for adults 65 years and older (N=2,036) admitted to the hospital following surgery at our institution. We utilized LCA to identify patient subgroups based on routinely collected measures of frailty, mobility, activities of daily living, and general health status. We compared hospital outcomes between the identified groups, including extended length of stay (LOS) (>0.5 SD beyond mean LOS by surgical category), discharge disposition (i.e., home versus non-home discharge), and utilization (weekly visit frequency) of physical therapy (PT) and occupational therapy (OT).
Results
We identified 3 subgroups of patients that we labeled Low Frailty-High Mobility (LF-HM), High Frailty-Low Mobility (HF-LM), and Low Frailty-Low Mobility (LF-LM), representing 15.3%, 27.6%, and 57.1% of the cohort, respectively. Discharge to home was highest among the LF-HM group (99%), followed by LF-LM (96%), and HF-LM (77%). Extended LOS was most common among the HF-LM group (27%), followed by LF-LM (18%), and LF-HM (6%). PT and OT visit frequencies were highest in the HF-LM group followed by the LF-LM and LF-HM groups.
Conclusions
This study identified 3 subgroups of postoperative older adults using routinely collected patient data. These groups may help to explain variations in hospital outcomes that have not been adequately explained using measures of preoperative or postoperative measures alone. These results may be useful in developing tailored postoperative care for older adults.
Introduction
More than half of the 20 million surgeries performed each year in the US are for adults who will require an inpatient hospital stay. 1 Among these, older adults are at increased risk for extended hospital lengths of stay (LOS) and discharge to post-acute care facilities,2–7 which conflict with patient preferences and contribute to increased cost of care.8–10 Previous studies have identified individual patient-level factors associated with postoperative outcomes. For example, studies have found that preoperative frailty is associated with increased risk for loss of independence following surgery, as well as longer LOS, among older adults.7,11 Another study found that patients with lower levels of mobility one day after gastrointestinal surgery were more likely to experience an extended LOS. 12 Additional studies have observed associations between age, gender, body mass index, preoperative function and American Society of Anesthesiologists (ASA) rating with postoperative outcomes following surgery.2,6,13,14 However, few studies have examined preoperative and postoperative patient data in combination, which is likely to account for variability in patient outcomes better than individual preoperative or postoperative measures alone.
Latent class analysis (LCA) is a model-based statistical approach used to identify hidden subgroups (i.e., “latent classes”) within a population based on patterns across multiple observed variables. 15 As a type of finite mixture modeling, LCA identifies groups of individuals who share similar characteristics or response patterns, allowing researchers to characterize heterogeneity that may not be directly observable. Identifying subgroups in this manner is commonly used to inform the development of tailored or “precision” treatment approaches.16–20 The purpose of this study was to identify subgroups of postoperative older adults using a latent class analysis of routinely collected patient-level data from a large academic health care system. We also examined differences in key postoperative outcomes among those in each of the identified latent subgroups.
Methods
We utilized a retrospective cohort design to identify latent subgroups of older adults admitted to the hospital following surgery and to examine differences in hospital outcomes among those in each latent subgroup. Data for this study was collected during routine care encounters and extracted from the electronic health record (EHR) for analysis.
Patient cohort
We included older adults (≥65) hospitalized at Johns Hopkins Hospital or Johns Hopkins Bayview Medical Center following inpatient surgery between September 2016 and March 2020. Eligible patients included those with data available for each of our class indicators (described below) and having any of the following surgical procedures: abdominal surgery, gynecologic surgery, pancreatic surgery, head & neck surgery, neurosurgery, orthopaedic surgery, plastic surgery, spine surgery, thoracic surgery, urologic surgery, or vascular surgery. Individuals undergoing outpatient surgery were not included in our cohort.
Latent class indicators
Based on previous evidence and standard data collection practices at our institution, we utilized patient age, general health status, frailty, mobility, and activities of daily living to identify latent subgroups of older adults after surgery.2,6,7,11–14 General health status was measured prior to surgery by the anesthesiology team using the American Society of Anesthesiologists Physical Status Classification System (ASA PS Classification), a method used to succinctly assess and communicate patients’ pre-anesthesia medical comorbidities. 21 ASA PS Classification scores range from 1 (normal health patient) to 6 (brain-dead). Frailty was measured at preoperative appointments by nursing using the Edmonton Frail Scale (EFS), a valid and reliable measure of frailty with scores ranging from 0-17 and higher scores indicative of greater frailty.22,23 Patient mobility function was assessed using the Activity Measure for Post-Acute Care (AM-PAC) Inpatient Basic Mobility “6-Clicks” Short Form.24,25 Patient function in activities of daily living was measured using the AM-PAC Inpatient Daily Activity “6-Clicks” Short Form. Both AM-PAC short forms include 6 items scored by clinicians on a scale of 1-4 with lower scores indicating greater impairment. Raw scores were converted to t-scores and used in all analyses; these t-scores range from 16.6 to 57.7. 26 In our hospitals, AM-PAC short forms are scored daily by nursing and at all physical therapy (PT) visits (mobility short form) and occupational therapy (OT) visits (daily activity short form). Previous studies have shown high interrater reliability across and within disciplines for these measures.27,28 We utilized the first observed AM-PAC mobility and daily activity scores recorded by either discipline following surgery (maximum 48 hours after surgery) for analysis.
Hospital outcomes
Hospital outcomes for this study included LOS, discharge disposition, and utilization of PT and OT services. LOS was categorized as an extended LOS, defined as a LOS (days) that was ≥0.5 standard deviations above the mean LOS for each surgical category, or a non-extended LOS, defined as any LOS less than the extended LOS definition. This approach was used to define extended LOS as LOS varies substantially depending on the type of surgery a patient receives. Defining extended LOS in this manner accounts for differences between surgical categories while still allowing our team to examine the number of patients experiencing an extended LOS across the full cohort. Discharge disposition was categorized as home discharge or non-home discharge (e.g., skilled nursing facility, inpatient rehabilitation facility). PT and OT utilization were examined separately based on weekly visit frequency, which was identified based on the number of PT or OT visits received, divided by patients’ LOS, and multiplied by 7.
Statistical analysis
Descriptive statistics were used to describe the sample. Latent class analysis was used to identify unobserved subgroups within the study population based on patterns across multiple class indicator variables. This model-based clustering approach estimates the probability that individuals belong to distinct latent classes defined by similar combinations of indicator values. Class indicators used in this study included a combination of continuous and categorical variables. Continuous class indicators included: age, AM-PAC mobility and AM-PAC daily activity scores, and EFS total score. Categorical class indicators included ASA ratings and individual EFS items (i.e., general health status, functional dependence, functional performance). ASA ratings were dichotomized to those patients who were healthy or with mild systemic disease (i.e., scores of 1-2) or patients with severe disease (i.e., scores of 3 and above). The general health status item on the EFS was dichotomized to those patients in good to excellent health or patients in fair to poor health. The functional dependence criterion of the EFS was dichotomized to those who required help in 0-1 activities or those who required help in 2 or more activities. The functional performance criterion of the EFS was dichotomized to those who completed the timed up and go test in 0-10 seconds or for those who completed it in 11 or more seconds. All LCA models were estimated using Mplus, Version 8.6. 29 Fit statistics, classification accuracy and relative class size were considered in selecting the best fitting latent class model. Guidelines suggest a minimum of 500 observations for LCA as used in this study. 15
Model fit was measured using the Bayesian information criterion (BIC) and the sample size-adjusted BIC (aBIC). 30 For the BIC and aBIC, smaller values indicate better model fit. The Vuong-Lo-Mendell-Rubin Likelihood Ratio Test (VLMR-LRT) was used to compare the relative fit of a model with k classes to a model with one fewer class. Significant p values for the VLMR-LRT indicate that the latent class solution with fewer classes should be rejected in support of the solution with more classes. 30 Entropy, ranging from 0 to 1, was also examined with larger values indicating better classification accuracy. Finally, class sizes were examined for each solution as simulation studies indicate that small or uncommon classes can be difficult to identify reliably and it is important to avoid over extracting classes. 30
Once the unconditional latent class model was established (step 1), we used the auxiliary function in Mplus with Bolck, Croon and Hagenaars (BCH) weights to examine differences across classes in the prevalence of two dichotomous distal hospital outcomes (i.e., discharge disposition, extended length of stay) and in the means of two continuous distal hospital outcomes (i.e., PT and OT weekly frequencies). In this automatic 2-step BCH approach, weights were applied to individual participants based on their posterior probabilities of class membership to account for classification uncertainty and an overall test of equality as well as pairwise tests of differences in each outcome across classes using one degree of freedom were conducted.29,31 P-values of less than 0.01 were considered statistically significant.
The surgical heterogeneity present in our sample has the potential to obscure meaningful patterns among more homogenous surgical populations within our cohort. To examine the influence of including patients with multiple types of surgeries in our analysis, we conducted two sensitivity analyses that limited our analysis to patients with spine surgery or abdominal surgery. These populations were selected as they represented homogenous patient groups and were among the largest subgroups.
Results
Cohort description.
aASA Ratings for those in this study included 1 (normal healthy patient), 2 (patient with mild systemic disease without significant functional limitation or end-organ involvement), 3 (patient with one or more severe systemic diseases causing substantive functional limitation) and 4 (patient with severe systemic disease that poses a constant threat to life).
bBased on number of hospitalizations over past year and self-reported general health (higher scores = worse health, instrument score range 0-2).
cBased on the number of daily activities that an individual requires assistance with (higher scores= more dependence, instrument score range 0-2).
dBased on impact of baseline health on ability to complete physical activity (higher scores = more impaired function, instrument score range 0-2).
EFS: Edmonton Frail Scale; AM-PAC: Activity Measure for Post-Acute Care; ASA: American Society of Anesthesiologists.
Fit indices for latent class models.
Note. LL-value = Log-likelihood value. BIC = Bayesian information criterion. aBIC = sample size adjusted BIC. LMR = Lo-Mendell-Rubin likelihood ratio test. LMR-LRT and entropy are not applicable for the 1-class model.
Description of identified classes and outcomes of class membership.
N=2,036.
Note. For categorical hospital outcomes, values by class represent proportions. For continuous outcomes values represent means.
†Significantly different compared to other classes in row (p<0.001).
αDo not sum to 100% due to rounding.
EFS: Edmonton Frail Scale; AM-PAC: Activity Measure for Post-Acute Care; ASA: American Society of Anesthesiologists.
*Represents a length of stay ≥0.5 SD above mean length of stay (days), by surgical category.
There were notable differences in class indicators between the three classes (Table 3). Age was similar between the LF-HM and LF-LM classes (71.6 and 72.5 years, respectively), but was higher in the HF-LM class (75.0 years). Frailty scores were similar in the LF-HM class (2.6) and LF-LM class (2.2) but were notably higher among those in the HF-LM class. Mean AM-PAC mobility scores were highest in the LF-HM class (54.2), followed by the LF-LM class (41.6), and HF-LM class (37.7). AM-PAC daily activity scores followed a similar pattern with mean scores of 54.8, 40.7, and 37.5 in the LF-HM, LF-LM, and HF-LM classes, respectively. ASA also differed by class with the HF-LM class having the highest probability (0.87) of receiving an ASA rating of 3 or higher and the LF-HM and LF-LM having similar but lower probabilities of receiving an ASA rating of 3 or higher (0.64 and 0.67, respectively). Individual EFS items also differed by class. Specifically, individuals in the HF-LM had a much higher probability of having fair to poor health status, being functionally dependent, and requiring more than 10 seconds to complete the up and go test as compared to the LF-HM and LF-LM classes.
Class membership varied by surgical category (Figure 1). For example, the LF-HM class was most represented among those undergoing head & neck surgery (32%) and was least represented among those undergoing hepatobiliary/pancreatic surgery (2%). The HF-LM class was most represented among those undergoing vascular surgery (50%) and least represented by those undergoing urologic surgery (12%). Finally, the LF-LM class was most represented by those undergoing hepatobiliary/pancreatic surgery (69%) and least among those undergoing plastic surgery (38%). Class Membership by Surgical Category HB: hepatobiliary.
All hospital outcomes of interest were significantly different by class (p<0.001) (Table 3). The proportion of those discharged home was 99% in the LF-HM class, 96% in the LF-LM class, and 77% in the HF-LM class. We observed that the percentage of patients experiencing an extended LOS was highest in the HF-LM class (27%), followed by the LF-LM class (18%), and the LF-HM class (6%). PT and OT visit frequencies were highest in HF-LM class (PT: 2.1 visits/week; OT: 1.1 visits/week), followed by the LF-LM class (PT: 1.2 visits/week; OT: 0.5 visit/week), and the LF-HM class (PT: 0.3 visits/week; OT: 0.1 visits/week).
We conducted two sensitivity analyses where we repeated each step of our LCA among 1 patients who had abdominal surgery (n=437) and 2 patients who had spine surgery (n=287). For the abdominal surgery group, we identified a 3-class model with results that were very similar to the full cohort with EFS scores within 1 point and AM-PAC scores within 2 points of those reported for the full cohort. For the spine surgery group, we also selected a 3-class model. While this model was not statistically superior to the 2-class model based on the VLMR p-value, we viewed this as a product of low sample size versus poor model fit, especially because other fit statistics, including smallest class size (12.1%), supported this model. This model was similar to the 3-class model identified using the full cohort with some notable differences. The spine HF-LM class included higher frailty scores than the equivalent class in the full cohort model. We also observed lower mobility scores across all 3 classes compared to the full cohort. Taken together, these sensitivity analyses suggest that the 3-class model identified using the full cohort is valid across different surgical categories but that the levels of frailty, mobility, and other class indicators may vary slightly by surgical category. Complete results of the sensitivity analyses are available as Supplemental Material.
Discussion
The results of this study indicate that LCA may be a useful approach for identifying latent subgroups of postoperative older adults, a large patient population vulnerable to poor outcomes. We specifically identified three latent subgroups of postoperative older adults, which we have labeled as Low Frailty-High Mobility (LF-HM) (15.3%), High Frailty-Low Mobility (HF-LM) (27.6%), and Low-Frailty-Low Mobility (LF-LM) (57.1%). These subgroups were present within each of the included surgical categories, although the proportion of patients belonging to each subgroup varied by type of surgery. Furthermore, patients in each subgroup differed significantly from one another in discharge disposition, LOS, and utilization of PT and OT services.
The classes identified in this study align with clinical reasoning and highlight the importance of physical function in the International Classification of Functioning, Disability and Health (ICF) domains of mobility and self-care (activities of daily living). 32 It was not surprising to us that our analysis identified a subgroup of patients (HF-LM) with poorer general health, higher frailty scores, and lower physical function and that these patients were more commonly discharged to post-acute care and experienced higher rates of extended LOS. Likewise, it stands to reason that our analysis identified another class of patients (LF-HM) that were generally healthy, had low frailty scores, had high mobility and daily activity function scores, were nearly always discharged home and rarely experienced an extended LOS. However, less expected was the identification of the LF-LM class, which is similar to the LF-HM group for all class identifiers, with the exception of AM-PAC scores which are both lower in the LF-LM class compared to the LF-HM class. Despite only this difference between groups, the LF-LM class experienced extended LOS rates 3 times higher and utilized 4-5 times as many PT and OT visits per week compared to the LF-HM class. These findings align with previous studies that have highlighted the influence of patient mobility and daily activity function on hospital outcomes and the importance of measuring these items as part of routine care.12,33–35
The findings of this study may inform strategies for developing tailored postoperative care for postoperative older adults by providing insights into how preoperative and postoperative factors may influence patient outcomes in combination. For example, previous studies suggest immediate postoperative mobility levels are closely tied with LOS and discharge disposition.12,33 In this study, we identified two groups of patients with relatively low mobility score (i.e., HF-LM, LF-LM). However, the group with higher preoperative frailty experienced extended LOS and non-home discharges significantly more often than those with lower frailty – suggesting that both of these measures should be considered when determining what resources (e.g., PT, OT) patients may require during their postoperative recovery. While the results of this study do immediately provide a deployable classification tool or new treatment strategy, they do suggest that tailored or “precision” postoperative care strategies for older adults should account for both preoperative frailty and postoperative mobility as these factors provide a more comprehensive view of the patient than either measure in isolation.
The classes identified in this study, while corresponding to one institution, are likely generalizable to other institutions performing surgery among older adults seeking to identify clinical subgroups. Importantly, this study was conducted exclusively using routinely collected data extracted from the EHR. While this increases the overall feasibility of the study performed, it does place additional emphasis on the identification and selection of clinical measures, as well as systematic collection and documentation of these measures. Previous publications can provide guidance on the process of selecting outcome measures for healthcare systems seeking to implement or refine data collection strategies.36,37
Strengths of this study include a large patient sample, well beyond the sample size of ≥500 observations recommended for LCA studies. 15 This study also utilized data elements that are routinely collected at many hospitals, increasing the generalizability of our results. There are also limitations to this study that should be considered when interpreting results. Our study used a retrospective design and was limited to data available in our EHR. It is likely that there are factors not measured routinely in our healthcare system that would have altered the subgroups identified in this study had they been included in our analysis. For example, social determinants of health (e.g., caregiver availability, home accessibility, family support structures) and mental health (i.e., depression) have been shown to influence LOS and discharge disposition among hospitalized patients and may have altered the composition of the classes we identified.38–42 Our analysis was also performed specifically among older adults after surgery and it is possible that this affects the generalizability of our results to other populations, such as postoperative adults under 65 years of age. It is also possible that missing data that affected eligibility for our analysis (e.g., AM-PAC) was not missing at random and could have systematically biased our results in ways unknown to the authors, affecting generalizability.
Conclusion
Using LCA to explore subgroups of postoperative adults, we identified 3 subgroups of patients based on routinely collected data extracted from the EHR. We observed different proportions of patients with extended LOS and discharged home across the identified subgroups. We also observed different rates of PT and OT utilization among these groups. Our results suggest that combining preoperative and postoperative measures may better account for heterogeneity of outcomes among postoperative older adults than preoperative or postoperative measures along. Using these routinely collected measures may help to inform care decisions following surgery for older adults. However, further research is needed to validate the findings of our analysis at other institutions, as well as to examine the influence of social factors on the classes identified in this study and the hospital outcomes observed among these groups.
Supplemental material
Supplemental material - Using latent class analysis to identify subgroups of postoperative older adults
Supplemental material for Using latent class analysis to identify subgroups of postoperative older adults by Kevin McLaughlin, Amie Bettencourt, Daniel Young, Erik Hoyer, Michael Friedman, Elizabeth Colantuoni, Lee A. Goeddel and Pedro Gozalo by Sage Open Medicine.
Footnotes
Ethical considerations
This study was acknowledged as exempt from review by the Johns Hopkins Medicine Institutional Review Board (IRB00337243).
Consent to participate
A waiver of written informed consent was formally granted by the Johns Hopkins Medicine Institutional Review Board.
Author contributions
Conceptualization: KM, DY, EH, MF, LG. Methodology: EC, AB. Data curation: KM, EH, LG. Formal Analysis: AB, KM. Writing (original draft): KM, AB, PG. Writing (reviewing and editing): KM, AB, DY, EH, MF, EC, LG, PG.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Funding for this study was provided by the Learning Health Systems Rehabilitation Research Network (LeaRRn) (NIH 5P2CHD101895-03)
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
