Abstract
Background:
Falls represent a major clinical and financial challenge for healthcare systems. Accurately predicting first falls remains challenging, especially when using routinely collected data.
Objective:
Develop and evaluate a predictive model to identify elderly people at risk of first fallādefined as a fall after a 90-day fall-free periodāin the Basque Country using routinely collected health records.
Method:
A retrospective study included patients aged ā„65 with at least two chronic conditions among heart failure, chronic obstructive pulmonary disease (COPD), and diabetes. Data on demographics, diagnoses, prescriptions and healthcare utilisation were obtained from Osakidetza-Basque Health Service databases. Patients were labelled as āfallersā if they fell during 2022ā2023 after a 90-day fall-free period rather than a true first-ever fall. Predictive modelsālogistic regression (LR), random forest (RF), and extreme gradient boosting (XGB)āwere trained using recursive feature elimination with cross-validation (RFECV). Shapley additive explanations (SHAP) enhanced model interpretability and explainability.
Results:
35,197 patients were included, with 10.6% experiencing a fall. All models achieved similar results, with an AUCROC score of 0.71, while precision remained low. Emergency room visits in the prior 3āÆmonths, presence of caregiver, and age consistently ranked the top predictors. The number of prescriptions and antidepressant use also emerged as relevant.
Discussion:
Models showed moderate predictive performance. Relying solely on routinely collected health records limited clinical applicability due to low precision. Future work should integrate diverse data sourcesāhealth records, real-time gait and balance metrics, environmental factorsāto improve fall prediction and support clinical decision-making.
1 Introduction
In recent years, the population over 65āÆyears and its life expectancy have increased (OECD, 2009; OECD, 2011). The growing longevity has led to a progressive rise in the prevalence of chronic diseases and associated polypharmacy (Allepuz Palau et al., 2015; Goodwin et al., 2012; Chang et al., 2020). As treatments become more complex and long-term (Brown and Bussell, 2011; PĆ©rez et al., 2002; Núñez Montenegro et al., 2014), older adults face a higher risk of adverse events such as falls (Hanlon et al., 1997).
Falls represent a major public health issue and a significant clinical and financial challenge for healthcare systems (Chang et al., 2020). About 30% of people over 65āÆyears and 50% of those over 80 experience a fall annually, with 75% of them likely to fall again (National Institute of Health and Care Excellence, 2013; Ministerio de Sanidad, Servicios Sociales e Igualdad, 2014). Falls are a leading cause of functional dependence in the elderly and the fifth leading cause of death among people over 65 years (Ministerio de Sanidad, Servicios Sociales e Igualdad, 2014). They are also associated with increased healthcare resource use (Campbell et al., 2004; INFAC, 2019; Steinman and Hanlon, 2010), accounting for 10% of emergency room visits and 6% of hospital admissions (Ministerio de Sanidad, Servicios Sociales e Igualdad, 2014), along with an associated increase in costs (INFAC, 2019). In addition, fear of falling is linked to negative psychological and behavioural effects, as well as reduced activity levels and lower quality of life (Scheffer et al., 2008; Schoene et al., 2019).
Given the negative health and economic impacts of falls, fall risk assessment and prevention have become a global healthcare priority (INFAC, 2019; Villafaina and GavilĆ”n, 2011; Orueta Mendia et al., 2014). Some risk factors can be controlled (LĆ”zaro del Nogal, 2009; MartĆnez et al., 2017; BalaguĆ© et al., 2015), such as the use of certain medications (Milos et al., 2014; de Vries et al., 2018; Fastbom and Schmidt, 2010), where timely interventions may contribute to fall prevention (Cameron et al., 2018; de Jong et al., 2013; Beers et al., 1991). In the Basque Country, 92% of accidents in people over 74āÆyears in 2012 were due to falls, a third of which were preventable (INFAC, 2019).
Consequently, early detection of individuals at risk of falling is gaining interest among healthcare organisations (National Institute of Health and Care Excellence, 2013; Villafaina and GavilƔn, 2011; Orueta Mendia et al., 2014). There is a pressing need for cost-effective prediction systems that can be applied in real-world settings (Rajagopalan et al., 2017). The secondary use of routinely collected data in clinical practise for this purpose could provide an effective solution (Lee et al., 2020), requiring no additional effort from already overburdened healthcare professionals (Duong and Vogel, 2023). Moreover, the digitalization of health records enables the integration of clinical and administrative data (Garrison et al., 2007; Menvielle et al., 2017; Noffsinger and Chin, 2000), facilitating the use of machine learning on large samples for the early detection of high-risk populations (Huang et al., 2021; Deschepper et al., 2019; Alowais et al., 2023). Combining electronic health records (EHR) with predictive models allows for the identification of key risk factors (Lee et al., 2020), providing healthcare professionals with interpretable models to support clinical decision-making (Orueta Mendia et al., 2014; Alowais et al., 2023).
Several studies have developed models to predict fall risk in elderly populations (Marier et al., 2016; Kuspinar et al., 2019; Lo et al., 2019; Millet et al., 2023; Lage et al., 2023; Lockhart et al., 2021; Noh et al., 2021; Dormosh et al., 2022)āusing techniques such as logistic regression (Marier et al., 2016; Millet et al., 2023; Lage et al., 2023), decision trees (Kuspinar et al., 2019; Millet et al., 2023), random forest (Lo et al., 2019; Millet et al., 2023; Lockhart et al., 2021), and extreme gradient boosting (Noh et al., 2021)ā, although few were based on routinely collected data (Marier et al., 2016; Dormosh et al., 2022; Dormosh et al., 2024). The multiple perspectives, methods, and predictors indicates that the field is still evolving (Deandrea et al., 2013). Most models include fall history among the analysed risk factors (Marier et al., 2016; Lo et al., 2019; Millet et al., 2023; Lage et al., 2023; Lockhart et al., 2021; Noh et al., 2021; Dormosh et al., 2022), which consistently emerges as the strongest predictor of a new fall (Pluijm et al., 2006; Stalenhoef et al., 2002), but limits their preventive value for first-time falls, an outcome that has been operationally defined in only one study as a fall occurring after a 90-day fall-free period (Kuspinar et al., 2019). Other models overlook prescriptions (Marier et al., 2016; Lo et al., 2019; Lockhart et al., 2021; Noh et al., 2021), despite the well-established role of certain drugs as modifiable fall risk factors (de Vries et al., 2018). Identifying individuals at risk of fall due to prescription-related issues is relevant for healthcare organisations to prevent adverse events and reduce costs (INFAC, 2019; Villafaina and GavilĆ”n, 2011). Some studies also fail to use a sufficiently large sample size (Marier et al., 2016; Millet et al., 2023; Lage et al., 2023; Lockhart et al., 2021; Noh et al., 2021). Furthermore, limited reporting on model metrics (Marier et al., 2016; Kuspinar et al., 2019), especially in terms of precision (Marier et al., 2016; Kuspinar et al., 2019; Lage et al., 2023; Lockhart et al., 2021), highlights the need for further robust research in the field (Seaman et al., 2022).
In this context, the European GATEKEEPER project was launched to harness new technologies to improve care for an ageing population (de Batlle et al., 2023). In the Basque Country, an intervention for chronic disease and polypharmacy management was implemented, with concern for falls linked to medication interactions. This initiative offered a unique opportunity to develop a predictive model to shed light on fall risk factors and assess its performance using routinely collected health data.
The present study contributes to the field by developing and evaluating a supervised machine learning predictive model that uses routinely collected health records to predict first fallsāoperationally defined as a fall occurring after a 90-day fall-free periodāand assess the role of polypharmacy in older adults with chronic conditions in the Basque Country, supporting early intervention and healthcare system sustainability.
2 Method
A retrospective study was conducted on the target population using anonymised data provided by Osakidetza, the Basque Health Service.
2.1 Target population
The target population was people registered in Osakidetza-Basque Health Serviceās databases that on 1 January 2022 were aged 65 or older and had 2 or more chronic diseases between heart failure, chronic obstructive pulmonary disease (COPD) and diabetes. Comorbidity criteria were defined using international classification of diseases, ninth revision (ICD-9) and tenth revision (ICD-10) diagnostic codes. The ICD-9 codes were 398.91, 402.*1, 404.*1, 404.*3, 425.*, 428.* for heart failure, 490.*-496.* for COPD and 250.* for diabetes. The ICD-10 codes were I09.9, I11.0, I13.0, I13.2, I25.5, I42.*-I43.*, I50.* for heart failure, J40.*-J47.* for COPD and E10.*-E13.* for diabetes.
Patients who experienced a first fall between 1 January 2022 and 31 December 2023 were classified as positive cases or āfallersā, with the fall date set as the reference. Falls were defined as events occurring after a 90-day fall-free period, rather than as guaranteed first-ever falls. This 90-day window was chosen pragmatically to align the outcome definition with the 90-day predictor window available from routinely collected data, consistent with previous studies using administrative records (Kuspinar et al., 2019). For ānon-fallersā, a random reference date within the same period was assigned. This case-reference design characterises the risk profile preceding falls that occur immediately after the 90-day observation window, rather than aiming to predict events within a predefined prospective time horizon.
2.2 Data source
All data were sourced from Osakidetza-Basque Health Serviceās anonymised corporative databases. The information was collected at the patient level and structured in three datasets (Figure 1). Patient demographics and disease registry (PDDR) included variables such as age, sex, socioeconomic status, deprivation index, Charlson index, ICD-9/10 diagnosis codes, fall history, and mortality. Patient medication history (PMH) recorded all prescriptions with start/end dates and anatomical therapeutic chemical (ATC) codes. Healthcare utilisation records (HUR) included, for primary care (PC), all contacts with medical and nursing staff, and for hospital care (HC), all contacts with outpatient services, emergency room, and hospitalisation. These data were obtained exclusively from the 90āÆdays preceding the reference date. All records dated on or after the reference date, including the index fall encounter itself, were excluded to ensure that the models could not capture the fall event and to prevent any information leakage from the outcome into the predictive models.
Figure 1
2.3 Predictive model
Following the GATEKEEPER project methodological framework for developing artificial intelligence (AI) systems in healthcare, this study applied established supervised machine learning methods across four main phases depicted in Figure 1: data preprocessing, feature selection, model development, and model evaluation. It also adhered to the transparent reporting of a multivariable prediction model for individual prognosis or diagnosis plus artificial intelligence (TRIPOD+AI) statement (Collins et al., 2024), adopted here as a reporting standard to ensure comprehensive and transparent reporting of the machine learning methods used, with further details provided in the Supplementary Table S1. All analyses were performed using Python (version 3.9.7) and the scikit-learn library (Scikit-Learn, n.d.).
2.3.1 Data preprocessing
First, a descriptive analysis of the sample was conducted. For categorical variables of two categories Fisherās exact test was used, whereas the chi-square test was applied to the others; for continuous variables with a normal distribution, Studentās t-test was employed. Given the large sample size, the standardised mean difference (SMD) was also calculated using Cohenās d for continuous variables and Cohenās h for categorical variables to quantify between-group differences independently of statistical significance.
Subsequently, to integrate the three aforementioned datasets a comprehensive data engineering process was conducted under the supervision of healthcare professionals to ensure scientific and technical rigour. The PDDR dataset, which captures data at a single time point, served as the core, where the definitions of all diagnostic predictor variables and their corresponding ICD-10 codes are provided in the Supplementary Table S2. All modelled variables were complete by construction, and therefore no patient was excluded due to missing data; the full eligible cohort was analysed. The PMH dataset provided prescription data, aggregated at the 3rd ATC levelāfirst four charactersāas previously in the literature (Milos et al., 2014). For each patient, the number of active medications per subgroup on the reference date was recorded. The HUR dataset offered insights into healthcare resource usage over the 90āÆdays before the reference date, aligning with first-fall definition. Pairwise correlations among variables were assessed considering all correlation values in absolute terms. For pairs with an absolute Pearson correlation coefficient greater than or equal to 0.8, the variable with the lower absolute correlation to the target was removed, thereby reducing multicollinearity and enhancing model interpretability.
The resulting final dataset preserved both clinical and temporal patterns while minimising data loss. It was split into training and test sets in an 80/20 ratio and standardised. Standardisation tools were trained on the training set and applied to the test set to prevent data leakage. Categorical variables with an ordinal structure were encoded using label encoding (LabelEncoder, n.d.), while numerical variables were scaled between 0 and 1 applying min-max normalisation based on the datasetās minimum and maximum values (MinMaxScaler, n.d.).
2.3.2 Feature selection
To reduce the number of variables and enhance model explainability, feature selection process was performed using recursive feature elimination with cross-validation (RFECV) (cv.āÆ=āÆ5) on the training set (RFECV, n.d.; Misra and Yadav, 2020). This method automatically determines the optimal number of features by adjusting an estimator across different subgroups of variables, which are then evaluated using a scorer. Three different algorithms were employed in this process: logistic regression (LR) (DomĆnguez-Almendros et al., 2011), random forest (RF) (Breiman, 2001), and extreme gradient boosting (XGB) (Chen and Guestrin, 2016). These models were selected for their widespread use in clinical risk prediction and prior applications in fall risk estimation works (Millet et al., 2023; Sitompul et al., 2025). Deep learning models were excluded, as XGB is recognised as the standard choice for classification tasks involving tabular data (Shwartz-Ziv and Armon, 2022). Models were configured with default hyperparameters and class weighting to address dataset imbalance, leading to an optimised evaluation using the area under the curve precision-recall (AUCPR) metric (Davis and Goadrich, 2006; McDermott et al., 2025). This process provided three selected variable subsets, each used to train the corresponding fall prediction model, enabling the evaluation of how different feature selection strategies influence model performance and interpretability.
2.3.3 Model development
The three algorithmsāLR, RF, and XGBāwere then trained on each of the three variable subsets identified during the feature selection phase. Before training, the hyperparameters were optimised using grid search with cross validation (cv.āÆ=āÆ5), with AUCPR serving as the target metric (LaValle et al., 2004). The different hyperparameters and their explored values for each algorithm are detailed in the Supplementary Tables S3āS5. To mitigate the effects of class imbalance, higher weights were assigned to the minority class during model training. The weight for each class was computed based on its relative frequency in the training data, giving more importance to underrepresented cases. Weighting was preferred over resampling because it addresses class imbalance without generating synthetic cases, preserving the real data distribution and allowing the models to account for the minority (faller) class through the AUCPR metric. This choice is supported by evidence that resampling-based corrections can impair model calibration and over-estimate minority-class risk without improving discrimination (van den Goorbergh et al., 2022; Carriero et al., 2025).
2.3.4 Model evaluation
Following the fine-tuning of hyperparameters, algorithms performance was assessed across all classes using several metrics (Scikit-Learn, n.d.): confusion matrix, area under the curve receiver operating characteristic (AUCROC) (Melo, 2013), AUCPR, accuracy, recall, and precision. Additionally, shapley additive explanations (SHAP) values were employed to quantify the contribution of each variable to the modelās predictions (Roth, 1988). SHAP values provide an interpretable measure of how much each variable increases or decreases the prediction for a given instance. This approach helps to identify the most influential variables, offering deeper insights into the modelās decision-making process and enhancing its interpretability and explainability.
3 Results
A total of 35,197 patients were included, with 10.6% experiencing falls (Table 1). No patient was excluded owing to missing or incomplete data. Women comprised 47% of the total, with a higher proportion in the fall group (51% vs. 46%). The mean age was 79.9āÆyears, significantly higher among fallers (82.6 vs. 79.6). The Charlson comorbidity index was also higher in this group (6.6 vs. 5.4), as was chronic medication use (10.3 vs. 8.8 drugs). Caregivers were present for 36% of patients, with a higher presence in the fall group (56% vs. 34%). Given the large sample size, all comparisons reached statistical significance. In terms of effect size, caregiver presence showed the largest effect size (0.44), followed by chronic prescriptions (0.39), age (0.37), and the Charlson index (0.36), all representing small-to-moderate differences, whereas sex and socioeconomic status showed only small differences (⤠0.12).
Table 1
| Variable | Total (N =āÆ35,197) | No fall (N =āÆ31,482) | Fall (N =āÆ3,715) | p-valuea | SMDb | ||||
|---|---|---|---|---|---|---|---|---|---|
| n | % (SE) or mean (SD) | n | % (SE) or mean (SD) | n | % (SE) or mean (SD) | ||||
| Sex | Women | 16,473 | 46.8 (0.3) | 14,591 | 46.3 (0.3) | 1,882 | 50.7 (0.8) | <0.01 | 0.09 |
| Men | 18,724 | 53.2 (0.3) | 16,891 | 53.7 (0.3) | 1,833 | 49.3 (0.8) | |||
| Age | Mean | 79.9 (8.1) | 79.6 (8.1) | 82.6 (7.5) | <0.01 | 0.37 | |||
| <80 | 17,433 | 49.5 (0.3) | 16,155 | 51.3 (0.3) | 1,278 | 34.4 (0.8) | <0.01 | 0.34 | |
| ā„ 80 | 17,764 | 50.5 (0.3) | 15,327 | 48.7 (0.3) | 2,437 | 65.6 (0.8) | |||
| Socioeconomic status | Low | 12,948 | 36.8 (0.3) | 11,775 | 37.4 (0.3) | 1,173 | 31.6 (0.8) | <0.01 | 0.12 |
| Medium | 15,313 | 43.5 (0.3) | 13,559 | 43.1 (0.3) | 1,754 | 47.2 (0.8) | |||
| High | 6,936 | 19.7 (0.2) | 6,148 | 19.5 (0.2) | 788 | 21.2 (0.7) | |||
| Charlson index | Mean | 5.5 (3.3) | 5.4 (3.3) | 6.6 (3.7) | <0.01 | 0.36 | |||
| 0 | 11 | 0.0 (0.0) | 9 | 0.0 (0.0) | 2 | 0.1 (0.0) | <0.01 | 0.30 | |
| 1ā2 | 5,979 | 17.0 (0.2) | 5,623 | 17.9 (0.2) | 356 | 9.6 (0.5) | |||
| 3ā4 | 10,668 | 30.3 (0.2) | 9,749 | 31.0 (0.3) | 919 | 24.7 (0.7) | |||
| ā„ 5 | 18,539 | 52.7 (0.3) | 16,101 | 51.1 (0.3) | 2,438 | 65.6 (0.8) | |||
| Chronic prescription | Mean | 9.0 (4.0) | 8.8 (3.9) | 10.3 (3.7) | <0.01 | 0.39 | |||
| 0ā9 | 20,339 | 57.8 (0.3) | 18,685 | 59.4 (0.3) | 1,654 | 44.5 (0.8) | <0.01 | 0.30 | |
| ā„10 | 14,858 | 42.2 (0.3) | 12,797 | 40.6 (0.3) | 2,061 | 55.5 (0.8) | |||
| Caregiver | No | 22,492 | 63.9 (0.3) | 20,837 | 66.2 (0.3) | 1,655 | 44.5 (0.3) | <0.01 | 0.44 |
| Yes | 12,705 | 36.1 (0.3) | 10,645 | 33.8 (0.3) | 2,060 | 55.5 (0.4) | |||
Patient characteristics of the sample.
Calculated using Fisherās exact test or chi square test for categorical variables and Studentās t-test for continuous variables.
Calculated using Cohen's d for continuous variables and Cohen's h for categorical variables.
SMD, standardised mean differences.
The dataset included 231 variables, with 77% from PMH, 21% from PDDR, and 2% from HUR, with a detailed breakdown by source provided in the Supplementary Table S6. The application of RFECV with the three estimators produced subsets of 183, 209, and 164 variables from LR, RF, and XGB, respectively. RFECV retained the feature subset maximising cross-validated AUCPR. The comparatively large subset reflects many low-magnitude but non-redundant signals, particularly in the medication (PMH) block. Details of these subsets and estimators are in the Supplementary Table S7. The outcomes of hyperparameter optimisation for the three models across the different database subsets are also presented in the Supplementary Table S8.
The evaluation results of the algorithms are summarised in Figures 2, 3, with a more detailed version available in the Supplementary Table S9. All models achieved similar performance, reaching an AUCROC value of 0.71 and an AUCPR value of 0.22. For recall, the XGB model with LR-selected features reached the highest value of 0.62, whereas precision peaked at 0.41 for the RF model with RF-selected features. Given the pronounced class imbalance, accuracy is not emphasised in interpreting these results.
Figure 2
Figure 3
The SHAP values for the 10 most influential variables in each model are illustrated in Table 2, while a detailed analysis is provided in the Supplementary Figures S1āS3. The colour intensity in the table reflects the influence of the variable, with darker shades indicating stronger influence and lighter shades indicating weaker influence. Emergency room visits in the past 3āÆmonths, caregiver presence, and age were consistently the most impactful variables across all models, followed by hospitalisations and persistent atrial fibrillation. Among drugs, antidepressants (N06A) had the highest impact, followed by anxiolytics (N05B) and high-ceiling diuretics (C03C), with the number of prescriptions also standing out as noteworthy.
Table 2
| Feature | Model | ||||||||
|---|---|---|---|---|---|---|---|---|---|
| Selec.: LR | Selec.: LR | Selec.: LR | Selec.: RF | Selec.: RF | Selec.: RF | Selec.: XGB | Selec.: XGB | Selec.: XGB | |
| Class.: LR | Class.: RF | Class.: XGB | Class.: LR | Class.: RF | Class.: XGB | Class.: LR | Class.: RF | Class.: XGB | |
| Number of ED visits in past 3āÆmonths | 0.146 | 0.033 | 0.253 | 0.133 | 0.037 | 0.239 | 0.145 | 0.039 | 0.256 |
| Caregiver | 0.203 | 0.033 | 0.234 | 0.201 | 0.034 | 0.210 | 0.200 | 0.037 | 0.234 |
| Age | 0.241 | 0.029 | 0.232 | 0.222 | 0.031 | 0.196 | 0.238 | 0.033 | 0.224 |
| Number of hospitalisations in past 3āÆmonths | 0.070 | 0.022 | 0.121 | 0.067 | 0.024 | 0.123 | 0.070 | 0.023 | 0.122 |
| Persistent atrial fibrillation | 0.123 | 0.017 | 0.096 | 0.124 | 0.015 | 0.104 | 0.120 | 0.016 | 0.092 |
| Number of PC nurse visits in past 3āÆmonths | ā | 0.021 | 0.112 | 0.047 | 0.019 | 0.112 | 0.054 | 0.021 | 0.104 |
| Sex | 0.116 | ā | 0.100 | 0.105 | ā | 0.081 | 0.110 | ā | 0.085 |
| Charlson points | ā | 0.019 | 0.075 | ā | 0.014 | 0.068 | ā | 0.016 | 0.069 |
| Mild liver disease | 0.081 | ā | ā | 0.083 | ā | ā | 0.081 | ā | ā |
| Antidepressants (N06A) | 0.059 | ā | 0.081 | 0.051 | ā | 0.070 | 0.056 | ā | 0.082 |
| Number of prescriptions | ā | 0.017 | 0.072 | ā | 0.013 | 0.047 | ā | 0.014 | 0.075 |
| Diabetes mellitus with chronic complications | 0.062 | ā | ā | 0.076 | ā | ā | 0.065 | ā | ā |
| Number of outpatient visits in past 3āÆmonths | ā | 0.016 | ā | ā | 0.012 | ā | ā | 0.012 | ā |
| Number of PC doctor visits in past 3āÆmonths | ā | 0.015 | ā | ā | ā | ā | ā | 0.012 | ā |
| Anxiolytics (N05B) | 0.056 | ā | ā | ā | ā | ā | ā | ā | ā |
| High-ceiling diuretics (C03C) | ā | ā | ā | ā | 0.012 | ā | ā | ā | ā |
Shapley additive explanations (SHAP) values of the 10 most influential features of each model.
Selec., selection algorithm; Class., classification algorithm; LR, logistic regression; RF, random forest; XGB, extreme gradient boosting. Darker shades: stronger influence. Lighter shades: weaker influence.
4 Discussion
Our findings revealed a moderate performance for first-fall prediction modelsāthat is, models predicting falls after a 90-day fall-free periodātargeting elderly individuals with chronic diseases and polypharmacy. While the predictive variables differed slightly across the models developed (Table 2), all aligned in identifying the same three variables as the strongest contributors to the modelās fall predictions: age, presence of a caregiver, and emergency room visits in the previous 3āÆmonths. Throughout the manuscript, these SHAP-based contributions should be interpreted as model-specific explanatory measures and not as evidence of causal or mechanistic effects on fall risk. Age is a well-established risk factor, as older adults tend to experience a decline in balance, strength, and mobility (Millet et al., 2023; Noh et al., 2021). The presence of a caregiver should be interpreted with caution, as it most likely acts as a proxy for underlying frailty, disability, or dependence, as well as increased healthcare utilisation and documentation, reflecting a greater need for assistance with activities of daily living rather than representing a direct risk factor (Rubenstein, 2006). Similarly, the number of recent ER visits may signal underlying health issues. All three of these factors are associated with the modelās fall predictions, but when combined, they suggest a specific patient profile: a dependent older individual who, due to residual effects or ongoing health issues, faces an elevated risk of falling. As healthcare utilisation was measured in the 90āÆdays before the fall, the patient profile identified by the models may partly reflect clinical deterioration preceding the fall rather than stable risk factors alone. However, this does not diminish the clinical relevance of these variables, as such deterioration represents a period of increased vulnerability during which the risk of falling is likely to be highest. Identifying this high-risk state may still provide a valuable opportunity for timely preventive interventions before the fall occurs, thereby averting its serious clinical and economic consequences. This suggests that the models may be better suited to identifying periods of acute, short-term fall riskāthat is, transient states of increased vulnerability rather than directly modifiable causes of fallsārather than capturing long-term susceptibility to falling.
When analysing drug effect, the number of prescriptions emerges as a relevant factor, underscoring that the consequences of polypharmacy on patient safety are significant, with risks increasing for each additional medication (Gnjidic et al., 2012). In particular, antidepressants (N06A) ranked among the top ten factors with the greatest contribution to the modelās fall prediction (Milos et al., 2014; Fastbom and Schmidt, 2010). Antidepressants can cause dizziness, sedation, and orthostatic hypotension, or interact with other drugs, especially in older adults, leading to cognitive impairment and drowsiness. These effects may compromise balance and coordination, increasing fall risk. High-ceiling diuretics (C03C) and anxiolytics (N05B) were also identified among the most relevant factors, in line with previous studies on medication-related falls (Supplementary material) (Milos et al., 2014; Fastbom and Schmidt, 2010). Nevertheless, moderate performance of the models suggests that not all relevant predictor variables were considered in this study, and that additional factors influencing the likelihood of a fall may exist but were not captured by the models.
In terms of performance, all analysed models showed similar results, with an AUCROC score of 0.71 in line with previous studies (Lo et al., 2019; Millet et al., 2023; Lage et al., 2023; Lockhart et al., 2021; Noh et al., 2021; Dormosh et al., 2022). However, a key requirement for reliable fall prediction is balancing recall and precision. Recall indicates modelās ability to identify true falls, with the highest recall of 0.62 achieved by the XGB model using LR-selected features. Precision reflects the proportion of fall predictions that correspond to true falls. High precision is critical for maintaining a manageable workload when reviewing at-risk patients and for minimising the costs associated with false positives. In this case, the precision value dropped significantly to 0.18 as in other studies (Lo et al., 2019; Noh et al., 2021; Dormosh et al., 2022), where the highest observed value is 0.41. Consequently, the high number of false positives hinders the practical implementation of the model in clinical practise settings, requiring substantial effort to review and validate predicted cases. At the observed prevalence of 10.6%, of every 1,000 patients screened approximately 106 would experience a fall and 894 would not. Given the precision achieved across models (0.18 to 0.41), this would translate into between roughly 259 and 589 patients flagged as at risk, of whom between approximately 153 and 483 would be false positives. Rather than reducing workload, this would increase the burden on healthcare professionals and compromise the cost-effectiveness of the models (Duong and Vogel, 2023). This highlights that relying solely on routinely collected health records is lacking when it comes to developing a robust and practically implementable fall prediction model, a challenge that could apply to any healthcare system that systematically collects similar information. Furthermore, it should also be considered that these models were developed within Osakidetzaāa publicly funded, vertically integrated single-payer healthcare system with a mature, shared electronic health recordāand consequently, their transportability is likely to be greatest in settings with similarly integrated records and comparable coding practises and data structures, whereas predictor availability and model performance may differ in more fragmented or insurance-based healthcare systems.
At this stage, it is important to recognise that fall prediction and prevention are complex and multifactorial challenges (Bargiotas et al., 2023). Intrinsic and extrinsic factors, including physiological, psychological, behavioural, and environmental aspects, interact to influence fall risk (Rajagopalan et al., 2017). Consequently, comprehensively managing all these factors remains highly challenging to date, as does gathering related relevant information to improve predictive model performance (Bargiotas et al., 2023). The absence of frailty measures, cognitive screening, comprehensive assessments of mobility and physical function, and environmental hazard data likely limited the modelsā ability to capture key mechanisms underlying falls. Consequently, the models relied on indirect, encounter-based proxies, which may partly explain their moderate discrimination and limited precision.
On one hand, relevant intrinsic data on gait and balance were unavailable. Although diagnoses of gait disorders, visual impairment, and/or vestibular dysfunction were available and included as predictors, they were binary indicators of prior diagnostic history rather than direct functional measures, likely constraining the discriminative and predictive performance achievable using routinely collected data alone. Technologies such as wearables and monitoring devices within the context of mobile health (mHealth) offer promising solution for gathering this information (Rajagopalan et al., 2017; Konara Mudiyanselage et al., 2025; Casilari et al., 2015), providing real-time insights into patientsā mobility patterns, balance stability, and overall physical activity levels (GonzĆ”lez-Castro et al., 2025). However, to date, the studies in the literature have only reported moderate performance in fall prediction models using these technologies (Konara Mudiyanselage et al., 2025). This suggests that integrating mHealth data with routinely collected health records may be key to making gait and balance information readily accessible and to enhancing model performance. Additionally, user acceptance and involvement in the monitoring process are crucial, as older adults often exhibit low engagement with wearable-based systems (Lage et al., 2023; MuƱoz Esquivel et al., 2023). As a more immediately implementable intermediate step, functional assessments already collected in routine primary careāsuch as frailty measures, mobility assessments, and/or standardised physical performance testsācould, where available, enrich routinely collected records at lower technological and adoption barriers than wearable-based systems.
On the other hand, another major gap lies in the absence of environmental data, such as geographic location, home layout, and other extrinsic factors affecting patients. While many predictive models are developed in controlled environments, their performance often differs substantially when applied to real-world settings (Rajagopalan et al., 2017). In this study, the inability to capture contextual dataādue to technological constraints such as limited access to new technologies or lack of integration with EHRsārestricted the understanding of patient mobility and living conditions, crucial for assessing external factors that may influence their health and well-being (Rajagopalan et al., 2017; Sczuka et al., 2023). Therefore, promoting the use and interoperability of data collected through technologies such as context-aware systems, monitoring devices, and ambient sensors is essential for identifying external fall risk factors, enhancing predictive models, and supporting targeted interventions for at-risk patients (Bargiotas et al., 2023; Job et al., 2020). Progress in this regard may also help overcome barriers to the adoption and widespread implementation of fall prevention systems in healthcare (Lee et al., 2020; Bargiotas et al., 2023; Ramadan et al., 2024).
Lastly, it is important to acknowledge that some falls occur purely by chance and are essentially unavoidable despite precautions, due to unpredictable factors like lapses in attention or sudden environmental changes. While prevention can reduce fall risks, not all falls can be prevented. However, given the high frequency and costs associated, even minor improvements in identifying high-risk individuals could result in significant cost savings.
Future work should prioritise user-centric longitudinal studies conducted in real-life conditions. This will require adaptable data-collection and fusion methods that leverage the strengths of each data source while ensuring compliance with patient-confidentiality regulations. It will also demand close interdisciplinary collaboration between IT and healthcare professionals to develop comprehensive fall prediction and prevention systems that support healthcare sustainability.
4.1 Limitations
The main limitation of this study lies in the inherent constraints of routinely collected data, which do not ensure complete and up-to-date information for every patient. While resource utilisation data are typically reliable due to consistent recording of service use, test and measurement data are only gathered when required for managing the patientās conditionāoften during in-person consultationsāso such information may not be available or up to date for all patients. Moreover, some measures, like the Barthel index and certain laboratory tests, have limited validity period and require frequent updates to remain accurate.
There is also the possibility that falls may be miscoded or not recorded at all, especially in cases where the fall did not result in injury or require medical attention (Hoffman et al., 2018). Information loss may be more substantial when integrating data from nursing homes and the private sector, due to shortcomings in data sharing (Vest and Kash, 2016). Consequently, this may have contributed to an underestimation of the total number of falls, with some true fallers likely misclassified as non-fallers, leading to conservative performance estimates.
Another limitation concerns the definition of first falls, based on a 90-day recall period. Defining falls as events occurring after a 90-day fall-free period is a pragmatic operational approach, and the authors acknowledge that it does not guarantee a true first-ever fall. While the timing and criteria used to define a first fall may be open to discussion, focusing on patients without a recent fall record remains valuable, as it enables predictive models to identify factors associated with falls and supports targeted preventive interventions. The way the first fall was defined in this study was intended to align the outcome definition with the predictor window available from routinely collected data, providing a practical basis for analysing factors associated with falls. The sensitivity of the findings to alternative washout periods was not assessed and remains an important area for future research, while a predefined prospective follow-up period also warrants further investigation.
A further limitation is the absence of external and temporal validation, as all models were developed and evaluated within a single healthcare system using an internal 80/20 split. Moreover, relying on a single random trainātest split may produce performance estimates that are more optimistic and less robust than those obtained through repeated cross-validation or evaluation in independent datasets, further highlighting the need for external and temporal validation. However, given the moderate predictive performance, limited precision, and high false-positive burden, external validation is not warranted at this stage. The priority is first to determine whether routinely collected data, complemented by additional sources (e.g., balance and gait, environmental, and lifestyle information), can achieve clinically useful predictive performance. For the same reason, a formal calibration assessment and a decision-curve (net-benefit) analysis were not performed, as they would not alter this conclusion. As the modelsā discrimination and precision are insufficient to support clinical use, a formal calibration assessment becomes informative only once a model achieves adequate discriminative performance, at which point both analyses would constitute relevant steps in future work. Given the modelsā insufficient discrimination and precision for clinical use, evaluating calibration would have provided limited additional value, as a formal calibration assessment becomes informative only once a model achieves adequate discriminative performance, at which point both calibration assessment and decision-curve analysis would constitute relevant steps in future work.
5 Conclusion
Relying solely on routinely collected health records makes it challenging to develop a robust and implementable fall prediction model for clinical practise. Given the multifactorial nature of falls, a comprehensive interdisciplinary approach is essential. Future models should incorporate up-to-date patient data and integrate diverse sources beyond medical records, such as real-time balance and gait, environmental, and lifestyle information. On the basis of routinely collected records alone, the performance achieved is insufficient to support clinical use, and authors do not propose these models for implementation or external validation at this stage. Only if future models, enriched with such additional data sources, achieve adequate predictive performance would external validation and, potentially, clinical deployment become justified, with the aim of helping healthcare providers identify high-risk individuals, improving quality of life and reducing fall-related healthcare costs.
Statements
Data availability statement
The datasets presented in this article are not readily available because given the potentially sensitive nature of the data, the Ethics Committees of the participating healthcare provider organisations did not authorise public access. Dataset access requests should be made through the official channels of each organisation. Requests to access these datasets should be directed to igor.larranagauribeetxebarria@bio-sistemak.eus.
Ethics statement
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. All procedures involving human subjects/patients were approved by the Clinical Research Ethics Committee of the Basque Country (PS2023033).
Author contributions
IL: Conceptualization, Data curation, Investigation, Resources, Validation, Visualization, Writing ā original draft, Writing ā review & editing. MR: Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing ā original draft, Writing ā review & editing. IA: Formal analysis, Methodology, Validation, Visualization, Writing ā review & editing. IE: Conceptualization, Investigation, Writing ā review & editing. MC: Methodology, Validation, Visualization, Writing ā review & editing. LL-P: Validation, Writing ā review & editing. MMe: Conceptualization, Investigation, Writing ā review & editing. MMa: Conceptualization, Investigation, Writing ā review & editing. HC-G: Conceptualization, Writing ā review & editing. EG: Writing ā review & editing. MH: Writing ā review & editing. LP: Writing ā review & editing. DF: Writing ā review & editing. GF: Project administration, Resources, Supervision, Writing ā review & editing. AF: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing ā review & editing.
Group member of GATEKEEPER Consortium
Claudio Caimi, Christian Tamporale, and Chiara Bonferini from Hewlett-Packard Italiana, Milan, Italy; Paolo Zampognaro and Federica Saca from Engineering Ingegneria Informatica S.p.A, Roma, Italy; Ioanna Drympeta and Konstantinos Votis from the Information Technologies Institute, Centre for Research and Technology Hellas, Thessaloniki, Greece; Frans Folkvord from PredictBy, Barcelona, Spain, and from Tilburg School of Humanities and Digital Sciences, Tilburg, Netherlands; Sergio Guillen and Juan-Carlos Naranjo from Mysphera S.L., Paterna, Spain; and Paula Curras, Jorge Posada, and German Gutierrez from Innova & European Projects Office, Integrated Health Solutions, Medtronic IbƩrica S.A., Madrid, Spain.
Funding
The author(s) declared that financial support was received for this work and/or its publication. The GATEKEEPER project has received funding from the European Unionās Horizon 2020 research and innovation programme under grant agreement no. 857223.
Acknowledgments
The authors would like to express their sincerest gratitude to all the partners and participants of the GATEKEEPER Consortium.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author GF declared that they were an editorial board member of Frontiers at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisherās note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frai.2026.1852439/full#supplementary-material
- AI
artificial intelligence
- ATC
anatomical therapeutic chemical
- AUCPR
area under the curve precision-recall
- AUCROC
area under the curve receiver operating characteristic
- COPD
chronic obstructive pulmonary disease
- EHR
electronic health records
- HC
hospital care
- HUR
healthcare utilisation records
- ICD-10
international classification of diseases, tenth revision
- ICD-9
international classification of diseases, ninth revision
- LR
logistic regression
- mHealth
mobile health
- PC
primary care
- PCM
patientsā clinical measurements
- PDDR
patient demographics and disease registry
- PMH
patient medication history
- RF
random forest
- RFECV
recursive feature elimination with cross-validation
- SHAP
shapley additive explanations
- SMD
standardised mean difference
- TRIPOD+AI
transparent reporting of a multivariable prediction model for individual prognosis or diagnosis plus artificial intelligence
- XGB
extreme gradient boosting
Glossary
References
1
Allepuz PalauA.PiƱeiro MĆ©ndezP.Molina HinojosaJ. C.Jou FerreV.Gabarró JuliĆ L. (2015). Evaluación económica de un programa de coordinación entre niveles para el manejo de pacientes crónicos complejos. Aten. Primaria47, 134ā140. doi: 10.1016/j.aprim.2014.05.002,
2
AlowaisS. A.AlghamdiS. S.AlsuhebanyN.AlqahtaniT.AlshayaA. I.AlmoharebS. N.et al. (2023). Revolutionizing healthcare: the role of artificial intelligence in clinical practice. BMC Med. Educ.23:689. doi: 10.1186/s12909-023-04698-z,
3
BalaguĆ©L.GaritanoB.MartĆnezJ.MayordomoM.PeƱaM. (2015). Recomendaciones basadas en evidencia para la prevención de caĆdas. Vitoria-Gasteiz: Departamento de Salud.
4
BargiotasI.WangD.MantillaJ.QuijouxF.MoreauA.VidalC.et al. (2023). Preventing falls: the use of machine learning for the prediction of future falls in individuals without history of fall. J. Neurol.270, 618ā631. doi: 10.1007/s00415-022-11251-3,
5
BeersM. H.OuslanderJ. G.RollingherI.ReubenD. B.BrooksJ.BeckJ. C. (1991). Explicit criteria for determining inappropriate medication use in nursing home residents. UCLA Division of Geriatric Medicine. Arch. Intern. Med.151, 1825ā1832.
6
BreimanL. (2001). Random forests. Mach. Learn.45, 5ā32. doi: 10.1023/a:1010933404324
7
BrownM. T.BussellJ. K. (2011). Medication adherence: WHO cares?Mayo Clin. Proc.86, 304ā314. doi: 10.4065/mcp.2010.0575,
8
CameronI. D.DyerS. M.PanagodaC. E.MurrayG. R.HillK. D.CummingR. G.et al. (2018). Interventions for preventing falls in older people in care facilities and hospitals. Cochrane Database Syst. Rev.9:CD005465. doi: 10.1002/14651858.CD005465.pub4,
9
CampbellS. E.SeymourD. G.PrimroseW. R. (2004). A systematic literature review of factors affecting outcome in older medical patients admitted to hospital. Age Ageing33, 110ā115. doi: 10.1093/ageing/afh036,
10
CarrieroA.LuijkenK.de HondA.MoonsK. G. M.van CalsterB.van SmedenM. (2025). The harms of class imbalance corrections for machine learning based prediction models: a simulation study. Stat. Med.44:e10320. doi: 10.1002/sim.10320,
11
CasilariE.LuqueR.MorónM. J. (2015). Analysis of android device-based solutions for fall detection. Sensors15, 17827ā17894. doi: 10.3390/s150817827,
12
ChangT. I.ParkH.KimD. W.JeonE. K.RheeC. M.Kalantar-ZadehK.et al. (2020). Polypharmacy, hospitalization, and mortality risk: a nationwide cohort study. Sci. Rep.10:18964. doi: 10.1038/s41598-020-75888-8,
13
ChenT.GuestrinC. (2016). XGBoost: A scalable tree boosting system. Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. 785ā794. doi: 10.1145/2939672.2939785
14
CollinsG. S.MoonsK. G. M.DhimanP.RileyR. D.BeamA. L.Van CalsterB.et al. (2024). TRIPOD+AI statement: updated guidance for reporting clinical prediction models that use regression or machine learning methods. BMJ385:e078378. doi: 10.1136/bmj-2023-078378,
15
DavisJ.GoadrichM.The relationship between precision-recall and ROC curves. Proceedings of the 23rd International Conference on Machine LearningNew YorkAssociation for Computing Machinery (2006) 233ā240.
16
de BatlleJ.BenĆtezI. D.MoncusĆ-MoixA.AndroutsosO.Angles BarbastroR.AntoniniA.et al. (2023). GATEKEEPERās strategy for the multinational large-scale piloting of an eHealth platform: tutorial on how to identify relevant settings and use cases. J. Med. Internet Res.25:e42187. doi: 10.2196/42187,
17
de JongM. R.Van der ElstM.HartholtK. A. (2013). Drug-related falls in older patients: implicated drugs, consequences, and possible prevention strategies. Ther. Adv. Drug Saf.4, 147ā154. doi: 10.1177/2042098613486829,
18
de VriesM.SeppalaL. J.DaamsJ. G.van de GlindE. M. M.MasudT.van der VeldeN. (2018). Fall-risk-increasing drugs: a systematic review and meta-analysis: I. Cardiovascular drugs. J. Am. Med. Dir. Assoc.19, 371.e1ā371.e9. doi: 10.1016/j.jamda.2017.12.013,
19
DeandreaS.BraviF.TuratiF.LucenteforteE.La VecchiaC.NegriE. (2013). Risk factors for falls in older people in nursing homes and hospitals. A systematic review and meta-analysis. Arch. Gerontol. Geriatr.56, 407ā415. doi: 10.1016/j.archger.2012.12.006,
20
DeschepperM.EecklooK.VogelaersD.WaegemanW. (2019). A hospital wide predictive model for unplanned readmission using hierarchical ICD data. Comput. Methods Prog. Biomed.173, 177ā183. doi: 10.1016/j.cmpb.2019.02.007,
21
DomĆnguez-AlmendrosS.BenĆtez-ParejoN.Gonzalez-RamirezA. R. (2011). Logistic regression models. Allergol. Immunopathol.39, 295ā305. doi: 10.1016/j.aller.2011.05.002,
22
DormoshN.SchutM. C.HeymansM. W.van der VeldeN.Abu-HannaA. (2022). Development and internal validation of a risk prediction model for falls among older people using primary care electronic health records. J. Gerontol. A Biol. Sci. Med. Sci.77, 1438ā1445. doi: 10.1093/gerona/glab311,
23
DormoshN.van de LooB.HeymansM. W.SchutM. C.MedlockS.van SchoorN. M.et al. (2024). A systematic review of fall prediction models for community-dwelling older adults: comparison between models based on research cohorts and models based on routinely collected data. Age Ageing53:afae131. doi: 10.1093/ageing/afae131,
24
DuongD.VogelL. (2023). Overworked health workers are āpast the point of exhaustionā. CMAJ195, E309āE310. doi: 10.1503/cmaj.1096042,
25
FastbomJ.SchmidtI.Indikatorer fƶr god lƤkemedelsterapi hos Ƥldre. Sweden: National Board of Health and Welfare; (2010). (The Swedish National Bord for Health and Welfare 2010).
26
GarrisonL. P. J.NeumannP. J.EricksonP.GarrisonL. P.MarshallD.MullinsC. D. (2007). Using real-world data for coverage and payment decisions: the ISPOR real-world data task force report. Value Health10, 326ā335. doi: 10.1111/j.1524-4733.2007.00186.x,
27
GnjidicD.HilmerS. N.BlythF. M.NaganathanV.WaiteL.SeibelM. J.et al. (2012). Polypharmacy cutoff and outcomes: five or more medicines were used to identify community-dwelling older men at risk of different adverse outcomes. J. Clin. Epidemiol.65, 989ā995. doi: 10.1016/j.jclinepi.2012.02.018,
28
GonzĆ”lez-CastroA.BenĆtez-AndradesJ. A.GonzĆ”lez-GonzĆ”lezR.Prada-GarcĆaC.Leirós-RodrĆguezR. (2025). Predicting fall risk in older adults: a machine learning comparison of accelerometric and non-accelerometric factors. Digit. Health11:20552076251331752. doi: 10.1177/20552076251331752
29
GoodwinN.SmithJ.DaviesA.PerryC.RosenR.DixonA.et al. (2012). Integrated Care for Patients and Populations: Improving Outcomes by Working Together. London: The Kingās Fund.
30
HanlonJ. T.SchmaderK. E.KoronkowskiM. J.WeinbergerM.LandsmanP. B.SamsaG. P.et al. (1997). Adverse drug events in high risk older outpatients. J. Am. Geriatr. Soc.45, 945ā948. doi: 10.1111/j.1532-5415.1997.tb02964.x,
31
HoffmanG. J.HaJ.AlexanderN. B.LangaK. M.TinettiM.MinL. C. (2018). Underreporting of fall injuries of older adults: implications for wellness visit fall risk screening. J. Am. Geriatr. Soc.66, 1195ā1200. doi: 10.1111/jgs.15360,
32
HuangY.TalwarA.ChatterjeeS.AparasuR. R. (2021). Application of machine learning in predicting hospital readmissions: a scoping review of the literature. BMC Med. Res. Methodol.21:96. doi: 10.1186/s12874-021-01284-z,
33
INFAC (2019). Medicamentos relacionados con caĆdas. Vitoria-Gasteiz: Departamento de Salud.
34
JobM.DottorA.VicecontiA.TestaM. (2020). Ecological gait as a fall indicator in older adults: a systematic review. Gerontologist60, e395āe412. doi: 10.1093/geront/gnz113,
35
Konara MudiyanselageS. P.YaoC. T.MaithreepalaS. D.LeeB. O. (2025). Emerging digital technologies used for fall detection in older adults in aged care: a scoping review. J. Am. Med. Dir. Assoc.26:105330. doi: 10.1016/j.jamda.2024.105330,
36
KuspinarA.HirdesJ. P.BergK.McArthurC.MorrisJ. N. (2019). Development and validation of an algorithm to assess risk of first-time falling among home care clients. BMC Geriatr.19:264. doi: 10.1186/s12877-019-1300-2,
37
LabelEncoder. scikit-learn. Available online at: https://scikit-learn/stable/modules/generated/sklearn.preprocessing.LabelEncoder.html (Accessed August 25, 2026).
38
LageI.BragaF.AlmendraM.MenesesF.TeixeiraL.AraĆŗjoO. (2023). Older people living alone: a predictive model of fall risk. Int. J. Environ. Res. Public Health20:13. doi: 10.3390/ijerph20136284,
39
LaValleS. M.BranickyM. S.LindemannS. R. (2004). On the relationship between classical grid search and probabilistic roadmaps. Int. J. Robot. Res.23, 673ā692. doi: 10.1177/0278364904045481
40
LĆ”zaro del NogalM. (2009). CaĆdas en el anciano. Med. Clin.133, 147ā153. doi: 10.1016/j.medcli.2008.12.029,
41
LeeT. C.ShahN. U.HaackA.BaxterS. L. (2020). Clinical implementation of predictive models embedded within electronic health record systems: a systematic review. Informatics (MDPI)7:25. doi: 10.3390/informatics7030025,
42
LoY.LynchS.UrbanowiczR.OlsonR.RitterA.WhitehouseC.et al. (2019). Using machine learning on home health care assessments to predict fall risk. Stud. Health Technol. Inform.264, 684ā688. doi: 10.3233/SHTI190310,
43
LockhartT. E.SoangraR.YoonH.WuT.FramesC. W.WeaverR.et al. (2021). Prediction of fall risk among community-dwelling older adults using a wearable system. Sci. Rep.11:20976. doi: 10.1038/s41598-021-00458-5,
44
MarierA.OlshoL. E. W.RhodesW.SpectorW. D. (2016). Improving prediction of fall risk among nursing home residents using electronic medical records. J. Am. Med. Inform. Assoc.23, 276ā282. doi: 10.1093/jamia/ocv061
45
MartĆnezJ.MayordomoM.DiagoA.FernĆ”ndezG.GagoI.GarrastatxuM.et alProtocolo para la prevención de caĆdas. Vitoria-Gasteiz: Departamento de Salud; (2017). Report No.: SS-350-2017.
46
McDermottM. B.ZhangH.HansenL. H.AngelottiG.GallifantJ.A closer look at AUROC and AUPRC under class imbalance. Proceedings of the 38th International Conference on Neural Information Processing SystemsRed HookCurran Associates Inc. (2025).
47
MeloF. (2013). āArea under the ROC curve,ā in Encyclopedia of Systems Biology, eds. DubitzkyW.WolkenhauerO.ChoK. H.YokotaH. (New York: Springer), 38ā39.
48
MenvielleL.Audrain-PonteviaA.MenvielleW. (2017). The Digitization of Healthcare: New Challenges and Opportunities. London: Palgrave Macmillan, 496.
49
MilletA.MadridA.Alonso-WeberJ.RodrĆguez-MaƱasJ.PĆ©rez-RodrĆ”-GuezR. (2023). Machine learning techniques applied to the development of a fall risk index for older adults. IEEE Access11, 84795ā84809. doi: 10.1109/access.2023.3299489
50
MilosV.Bondessonà .MagnussonM.JakobssonU.WesterlundT.MidlövP. (2014). Fall risk-increasing drugs and falls: a cross-sectional study among elderly patients in primary care. BMC Geriatr.14:40. doi: 10.1186/1471-2318-14-40,
51
Ministerio de Sanidad, Servicios Sociales e Igualdad (2014). Documento de consenso sobre prevención de fragilidad y caĆdas en la persona mayor. Estrategia de promoción de la salud y prevención en el SNS. Madrid: Ministerio de Sanidad, Servicios Sociales e Igualdad.
52
MinMaxScaler. scikit-learn. Available online at: https://scikit-learn/stable/modules/generated/sklearn.preprocessing.MinMaxScaler.html (Accessed August 25, 2026).
53
MisraP.YadavA. S. (2020). Improving the classification accuracy using recursive feature elimination with cross-validation. Int. J. Emerg. Technol.11, 659ā665.
54
MuƱoz EsquivelK.GillespieJ.KellyD.CondellJ.DaviesR.McHughC.et al. (2023). Factors influencing continued wearable device use in older adult populations: quantitative study. JMIR Aging6:e36807. doi: 10.2196/36807,
55
National Institute of Health and Care Excellence (2013). Falls in Older People: Assessing Risk and Prevention. London: National Institute of Health and Care Excellence (NICE).
56
NoffsingerR.ChinS. (2000). Improving the delivery of care and reducing healthcare costs with the digitization of information. J. Healthc. Inf. Manag.14, 23ā30.
57
NohB.YoumC.GohE.LeeM.ParkH.JeonH.et al. (2021). XGBoost based machine learning approach to predict the risk of fall in older adults using gait outcomes. Sci. Rep.11:12183. doi: 10.1038/s41598-021-91797-w,
58
Núñez MontenegroA. J.Montiel LuqueA.MartĆn AuriolesE.Torres VerdĆŗB.Lara MorenoC.GonzĆ”lez CorreaJ. A. (2014). Adherencia al tratamiento en pacientes polimedicados mayores de 65 aƱos con prescripción por principio activo. Aten. Primaria46, 238ā245. doi: 10.1016/j.aprim.2013.10.003,
59
OECD (2009). OECD Factbook 2009: Economic, Environmental and Social Statistics. Paris: OECD Publishing.
60
OECD (2011). Health at a Glance 2011: OECD Indicators. Paris: OECD Publishing.
61
Orueta MendiaJ. F.GarcĆa-ĆlvarezA.Alonso-MorĆ”nE.NuƱo-SolinisR. (2014). Desarrollo de un modelo de predicción de riesgo de hospitalizaciones no programadas en el PaĆs Vasco. Rev. Esp. Salud Publica88, 251ā260. doi: 10.4321/S1135-57272014000200007,
62
PĆ©rezM.CastilloR.RodrĆguezJ.MartosE.MoralesA. (2002). Adecuación del tratamiento farmacológico en población anciana polimedicada. Med. Fam.3, 23ā28.
63
PluijmS. M. F.SmitJ. H.TrompE. A. M.StelV. S.DeegD. J. H.BouterL. M.et al. (2006). A risk profile for identifying community-dwelling elderly with a high risk of recurrent falling: results of a 3-year prospective study. Osteoporos. Int.17, 417ā425. doi: 10.1007/s00198-005-0002-0,
64
RajagopalanR.LitvanI.JungT. P. (2017). Fall prediction and prevention systems: recent trends, challenges, and future research directions. Sensors17:2509. doi: 10.3390/s17112509,
65
RamadanO. M. E.AlruwailiM. M.AlruwailiA. N.ElsehrawyM. G.AlanaziS. (2024). Facilitators and barriers to AI adoption in nursing practice: a qualitative study of registered nursesā perspectives. BMC Nurs.23:891. doi: 10.1186/s12912-024-02571-y,
66
RFECV. scikit-learn. Available online at: https://scikit-learn/stable/modules/generated/sklearn.feature_selection.RFECV.html (Accessed August 25, 2026).
67
RothA. E., (ed). (1988). The Shapley Value: Essays in Honor of Lloyd S. Shapley. Cambridge: Cambridge University Press.
68
RubensteinL. Z. (2006). Falls in older people: epidemiology, risk factors and strategies for prevention. Age Ageing35, ii37āii41. doi: 10.1093/ageing/afl084,
69
SchefferA. C.SchuurmansM. J.van DijkN.van der HooftT.de RooijS. E. (2008). Fear of falling: measurement strategy, prevalence, risk factors and consequences among older persons. Age Ageing37, 19ā24. doi: 10.1093/ageing/afm169,
70
SchoeneD.HellerC.AungY. N.SieberC. C.KemmlerW.FreibergerE. (2019). A systematic review on the influence of fear of falling on quality of life in older people: is there a role for falls?Clin. Interv. Aging14, 701ā719. doi: 10.2147/CIA.S197857,
71
Scikit-Learn. scikit-learn. Available online at: https://scikit-learn.org/stable/api/sklearn.metrics.html#module-sklearn.metrics (Accessed August 25, 2026).
72
SczukaK. S.SchneiderM.SchellenbachM.KerseN.BeckerC.KlenkJ. (2023). Evaluating the effect of activity and environment on fall risk in a paradigm-depending laboratory setting: protocol for an experimental pilot study. JMIR Res. Protoc.12:e46930. doi: 10.2196/46930,
73
SeamanK.LudlowK.WabeN.DoddsL.SietteJ.NguyenA.et al. (2022). The use of predictive fall models for older adults receiving aged care, using routinely collected electronic health record data: a systematic review. BMC Geriatr.22:210. doi: 10.1186/s12877-022-02901-2,
74
Shwartz-ZivR.ArmonA. (2022). Tabular data: deep learning is not all you need. Inf. Fusion81, 84ā90. doi: 10.1016/j.inffus.2021.11.011
75
SitompulL. R.NababanA. A.ManihurukM. L.PonsenW. A.SupriyandiS. (2025). Comparison of XGBoost, random forest, and logistic regression algorithms in stroke disease classification. SinkrOn9, 957ā968. doi: 10.33395/sinkron.v9i2.14794
76
StalenhoefP. A.DiederiksJ. P. M.KnottnerusJ. A.KesterA. D. M.CrebolderH. F. J. M. (2002). A risk model for the prediction of recurrent falls in community-dwelling elderly: a prospective cohort study. J. Clin. Epidemiol.55, 1088ā1094. doi: 10.1016/S0895-4356(02)00502-4,
77
SteinmanM. A.HanlonJ. T. (2010). Managing medications in clinically complex elders: āthereās got to be a happy mediumā. JAMA304, 1592ā1601. doi: 10.1001/jama.2010.1482,
78
van den GoorberghR.van SmedenM.TimmermanD.Van CalsterB. (2022). The harm of class imbalance corrections for risk prediction models: illustration and simulation using logistic regression. J. Am. Med. Inform. Assoc.29, 1525ā1534. doi: 10.1093/jamia/ocac093,
79
VestJ. R.KashB. A. (2016). Differing strategies to meet information-sharing needs: publicly supported community health information exchanges versus health systemsā enterprise health information exchanges. Milbank Q.94, 77ā108. doi: 10.1111/1468-0009.12180,
80
VillafainaA.GavilĆ”nE. (2011). Pacientes polimedicados frĆ”giles, un reto para el sistema sanitario. Inf. Ter. Sist. Nac. Salud.35, 114ā123.
Summary
Keywords
artificial intelligence, chronic diseases, falls, first-fall prediction, machine learning, polypharmacy, predictive modelling
Citation
LarraƱaga I, Rujas M, Alayo I, Erreguerena I, Capo M, Lopez-Perez L, Mediavilla MM, MartĆnez MD, ChuliĆ”-Gil H, Georga EI, Haleem MS, Pecchia L, Fotiadis DI, Fico G and Fullaondo A (2026) Predicting first falls among older adults with chronic conditions and polypharmacy using routinely collected health records. Front. Artif. Intell. 9:1852439. doi: 10.3389/frai.2026.1852439
Received
10 April 2026
Revised
04 August 2026
Accepted
13 August 2026
Published
04 September 2026
Volume
9 - 2026
Edited by
Jinghua Wang, Tianjin Neurological Institute, China
Reviewed by
Jake Luo, University of WisconsināMilwaukee, United States
Raquel Leirós-RodrĆguez, University of León, Spain
Updates
Copyright
Ā© 2026 LarraƱaga, Rujas, Alayo, Erreguerena, Capo, Lopez-Perez, Mediavilla, MartĆnez, ChuliĆ”-Gil, Georga, Haleem, Pecchia, Fotiadis, Fico, Fullaondo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Igor LarraƱaga, igor.larranagauribeetxebarria@bio-sistemak.eus
ā These authors have contributed equally to this work
ā” These authors share senior authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.