Objective Identification of older patients at risk, among those accessing the emergency department (ED), may support clinical decision-making. To this purpose, we developed and validated the Dynamic Silver Code (DSC), a score based on real-time linkage of administrative data.
Design and setting The ‘Silver Code National Project (SCNP)’, a non-concurrent cohort study, was used for retrospective development and internal validation of the DSC. External validation was obtained in the ‘Anziani in DEA (AIDEA)’ concurrent cohort study, where the DSC was generated by the software routinely used in the ED.
Participants The SCNP contained 281 321 records of 180 079 residents aged 75+ years from Tuscany and Lazio, Italy, admitted via the ED to Internal Medicine or Geriatrics units. The AIDEA study enrolled 4425 subjects aged 75+ years (5217 records) accessing two EDs in the area of Florence, Italy.
Outcome measures Primary outcome: 1-year mortality. Secondary outcomes: 7 and 30-day mortality and 1-year recurrent ED visits.
Results Advancing age, male gender, previous hospital admission, discharge diagnosis, time from discharge and polypharmacy predicted 1-year mortality and contributed to the DSC in the development subsample of the SCNP cohort. Based on score quartiles, participants were classified into low, medium, high and very high-risk classes. In the SCNP validation sample, mortality increased progressively from 144 to 367 per 1000 person-years, across DSC classes, with HR (95% CI) of 1.92 (1.85 to 1.99), 2.71 (2.61 to 2.81) and 5.40 (5.21 to 5.59) in class II, III and IV, respectively versus class I (p<0.001). Findings were similar in AIDEA, where the DSC predicted also recurrent ED visits in 1 year. In both databases, the DSC predicted 7 and 30-day mortality.
Conclusions The DSC, based on administrative data available in real time, predicts prognosis of older patients and might improve their management in the ED.
- prognostic assessment
- emergency department
- dynamic silver code
- administrative data
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
Strengths and limitations of this study
The Dynamic Silver Code combines in real time administrative data available in all Italian regions, to automatically produce a score predicting 1-year mortality and other health outcomes in older persons accessing the emergency department (ED).
The tool was developed retrospectively and initially validated in a large, representative cohort of patients aged 75+ years; it was further validated in a new cohort of subjects prospectively recruited, where the tool was available in the software routinely used for clinical management of patients in the ED.
We did not include predictors such as data on outpatient services, functional status and census, nor outcomes such as functional impairment, cognitive decline and institutionalisation.
At the end of the disease era, when medical care is mostly directed towards management of chronic multimorbidity or non-disease-specific complaints in older persons,1 a gear shift in the scope of, and the approach to, prognostic assessment is necessary. Besides survival, other goals should be pursued as outcomes of disease prevention and treatment, such as relief of symptoms, maintenance of personal autonomy and preservation of quality of life, all of which may be targets of prognostic evaluation.2 Focusing on one single disease is of limited value, when other coexistent illnesses, age-related physiological changes and non-biological determinants of health may all affect prognosis. Patient’s preferences and future perspectives must be valued: many older patients at the end of their life receive futile therapies for minor conditions,3 whereas others with a reasonable disability-free life expectancy fail to receive appropriate treatments just because they are considered too old.4
Prognostic tools have the potential to help clinicians optimise the benefit/burden ratio of medical interventions2 and support their decision-making5; they may prove useful also in the health policy arena, to perform risk adjustment when comparing clinical interventions or healthcare models.2 5 6 Administrative data are increasingly used to this purpose, because of their variety, availability, low cost and accuracy.7–9 Such tools have been applied to screen inpatients,10 11 outpatients12 and free-living subjects5 12 13 for mortality, hospitalisation or disease-specific outcomes, aiming for personalised cures.14–16
We have previously described the Silver Code (SC), which combines administrative data into a score predicting 1-year survival in subjects aged 75+ years, admitted to the hospital via the emergency department (ED).17 18 The Dynamic SC (DSC) presented here is an evolution of the original score, developed from a large administrative data set and then applied prospectively in a new sample. Whereas the original SC was based on one single moment of observation, the DSC considers, for each individual, the dynamics of events occurring across time. The new tool has been implemented into the software routinely used in the ED of several hospitals in Tuscany, Italy, to provide automated, real-time risk stratification of older patients.
The DSC was obtained following the same general approach used for the SC, that is, combining into a score variables from healthcare administrative archives, to predict long-term survival of older persons admitted to the hospital. Compared with the SC, the cohort used to develop and validate the DSC was much larger and representative of the general older population. Moreover, important computational differences were introduced and external validation in a completely new cohort was obtained.
Study design and data source
Two different studies contributed to the present investigation. The first was the ‘Silver Code National Project (SCNP)’, which was sponsored by the Centre for Disease Control of the Italian Ministry of Health in 2008.19 The second was represented by the ‘Anziani in DEA (AIDEA)’ (standing for ‘Older Persons in the ED’) study, jointly sponsored by the Italian Ministry of Health and by the Tuscany Region in 2010 (RF-2010-2321801). The SCNP was a non-concurrent cohort study, whereas AIDEA followed a prospective cohort study design.
The SCNP database,19 representing the primary data source for the development and initial validation of the DSC, was obtained from the administrative archives of two regions in Italy, Tuscany and Lazio, which deliver healthcare services to a total population of more than 9.6 million persons. The archives included data on demographics, hospitalisations, drug prescriptions and deaths of beneficiaries aged 75+ years, who had been admitted via the ED to hospitals in the two regions between April 2004 and December 2009. Data were linked using a numeric unique identifier, which allowed records anonymisation before data processing to preserve beneficiaries’ confidentiality. Universal healthcare coverage in Italy allows completeness and comprehensiveness of the information contained in these databases, which have been used in previous epidemiological studies.17 20 The Italian Ministry of Health reported that Tuscany archives are 100% complete and 95% accurate.21 22
Further validation of the score developed from the SCNP was subsequently obtained in a different sample, assembled in the AIDEA study, which was conducted in the ED of two hospitals in Florence, Italy, the Azienda Ospedaliero-Universitaria Careggi (AOUC), an academic tertiary hospital, and the Ospedale S Maria Annunziata (OSMA), a large community hospital. Enrolment was consecutively conducted between June and August 2016 and again between February and March 2017 in the AOUC, and between August and September 2016 in the OSMA, for a total of 22 weeks. In AIDEA, an application generating the DSC was incorporated into the software routinely used in the ED of the hospitals of the Tuscany Region: as soon as an eligible patient was triaged, the repository of healthcare data of the Local Health Unit was queried to provide, thanks to on-demand linkage of the different archives involved, the information required to obtain the DSC. The score was then in real time calculated and shown, together with the corresponding risk class, onto the computer screen. The lag time between occurrence of events contributing to the DSC (hospitalisations and drug prescriptions) and their registration in the healthcare data repository was approximately 2 weeks.
The SCNP selected only records of subjects aged 75+ years, residing in the area where the study was conducted, who consecutively accessed the ED of the participating hospitals in the specified time windows and were eventually admitted to Internal Medicine or Geriatrics wards. Conversely, the AIDEA study enrolled all the residents aged 75+ years consecutively accessing the ED in the specified time periods, independent of the outcome (hospitalisation or discharge) of their access to the ED.
Patient and public involvement
Patients or the public were not involved in the design or conduct of our research, partially because of its retrospective nature. Patients’ associations will be involved in the upcoming, large-scale application of the DSC in hospitals of the Tuscany Region. Reporting was provided as requested by the funding institutions. Dissemination to the public was obtained through lay press.
Statistical analysis was performed with SPSS for Mac V.25 (IBM Corp., Armonk, NY, USA), STATA V.15.1 (StataCorp. 2017. Stata Statistical Software: Release 15. College Station, TX: StataCorp LLC) and R V.3.5.0 (R Core Team, 2018).
A total of 281 321 records were available in the SCNP database to create the DSC and test its validity. The χ2 was used to test differences in relative frequencies and to identify bivariate predictors of death, taking into account trends as appropriate.
The sample was randomly split into a development and a validation subsample. In the first subsample, a Cox proportional hazards model was fit to estimate the association, expressed as HR with 95% CI of demographics and other variables on 1-year all-cause mortality risk. Variables initially considered, by definition limited to those available prior to the index ED visit, included age, gender, number of drugs prescribed in the previous 3 months as resulting from pharmacy claims (categorised as 0–3, 4–5, 6–8 and 9+), days from previous hospital admission (no admission, 30–180 days, 0–30 days) and its associated main diagnostic group (cancer, respiratory disease, cardiovascular disease and other conditions), which were selected as the most frequently observed, in the study cohort, among those included in the International Classification of Diseases, Ninth Revision (ICD-9) classification. Compared with the original version of the SC, previous admission to a day hospital was not considered, because availability and utilisation of day hospital service are not consistent across regions in Italy. Marital status information was also omitted, because it was frequently missing in the discharge summary and in preliminary analyses it contributed marginally to the prediction of death. Different from the SC, repeated hospitalisations were taken into account to dynamically update patient’s information at each new hospital admission. To this purpose, data were arranged in order to have one observation per event or time interval, and the counting process approach proposed by Andersen and Gill23 was applied: this is a generalisation of Cox model, which assumes that the correlation between event times for a person can be explained by past events. Thus, we made the assumption that correlations among events for each individual were captured by appropriate time-dependent covariates.
We evaluated the adherence of predictors to the assumption of proportional hazards and tested for multicollinearity comparing different models in terms of Akaike information criteria.24 We then assigned each risk factor a score, calculated as the ratio, rounded to the nearest integer, between the regression b coefficient for that variable and the smallest significant b coefficient (b0) in the Cox model. We finally computed the DSC as a summary score, by adding the points individually assigned to each risk factor. Four prognostic classes were then created from DSC quartiles.
In the validation data set, Harrell’s C index of concordance was applied to evaluate model performance, as a measure of the predictive power of Cox regression model.25 Furthermore, the ability of DSC classes to predict 1-year mortality was analysed after adjusting for region of enrolment and ICD-9 coded principal diagnosis at discharge after the index hospitalisation. External validity of the DSC was performed in the AIDEA data set, using Cox regression models and Harrell’s C index of concordance to predict 1-year mortality and the risk of recurrent ED access after a first hospital admission.
Additional analyses were performed to verify the ability of the DSC to predict 7 and 30-day mortality, two outcomes that are more immediate and potentially of greater clinical interest in the perspective of application of the tool in the ED. Also these analyses were adjusted for main diagnostic group at discharge and region of enrolment.
A two-tailed p value less than 0.05 was considered statistically significant.
Participants in the SCNP
Out of the 281 321 records in the SCNP, 180 079 (64.01%) pertained to a single hospitalisation, the remaining to patients with multiple hospital admissions in the study period. The random split assigned 90 039 patients to the development and 90 040 to the validation subsample, with 140 716 and 140 605 records, respectively. As shown in table 1, baseline characteristics of participants on their first hospital admission were comparable between the two subsamples.
Predictors of 1-year mortality
From hospital admission through the following year, a total of 42 434 deaths were recorded; mortality was comparable in the development (21 250/90 039, 23.6%) and the validation (21 184/90 040, 23.5%) subsamples. In the development sample, bivariate predictors of death were older age (75–79 years: 16.4%, 80–84 years: 20.5%, 85+ years: 31.9%; p=0.001), gender (men: 24.8%, women: 22.7%; p<0.001), previous hospital admission with its corresponding main discharge diagnosis (no admission: 22.1%, cardiovascular disease: 23.5%, cancer: 49.1%, respiratory disease: 30.8%, other diagnoses: 30.2%; p<0.001), days from previous hospital admission (1–30 days: 34.5%, 31–180 days: 28.5%, p<0.001) or taking more drugs prior to the index ED admission (0–3: 22.6%, 4–5: 22.9%, 6–8: 23.3%, 9+: 26.2%; p<0.001).
Table 2 shows the results of Cox proportional hazards regression model in the SCNP development subsample. All the variables considered were significant predictors of 1-year mortality. Because the diagnostic groups of ‘cardiovascular disease’ and ‘others’ obtained similar HRs, they were combined into a single category, contrasted with ‘no previous hospital admission’, ‘cancer’ and ‘respiratory disease’. The predictors ‘days from the previous hospital admission’ and the associated ‘main diagnostic group’ showed some multicollinearity, but the models with only one of those variables performed worse in terms of Akaike information criteria: therefore, both variables were maintained in the final model. All predictors satisfied the assumption of proportional hazards.
The smallest significant b coefficient in the Cox model resulted to be associated with four to five drugs prescribed in the previous 3 months (b0=0.0364). Thus, each significant predictor was assigned a score by rounding the ratio between the corresponding b coefficient in the regression model and b0 up to the nearest integer (table 2). The DSC was finally calculated as summation of all scores.
Participants were classified into low, medium, high and very high-risk classes, based on DSC quartiles. One-year risk of death increased significantly from the first to the fourth DSC class, with almost identical rates and HRs in the two subsamples (table 3). In the validation subsample, the performance of the risk scoring system, assessed with Harrell’s C index, was 0.668 (95% CI 0.665 to 0.672). The ability of the DSC to predict 1-year mortality persisted controlling for main diagnostic group at discharge after the index hospitalisation and region of enrolment (Tuscany vs Lazio), as shown by HR (95% CI) of 1.93 (1.88 to 1.98), 2.71 (2.64 to 2.78) and 5.00 (4.88 to 5.13) for class II, III and IV, respectively, compared with class I (p<0.001) (figure 1).
External validity of the DSC was ascertained in the AIDEA sample, which included 5217 records for 4425 subjects, of whom 1292 died and 465 had a new ED admission within 1 year. Mortality increased progressively across DSC classes, with rates of 106, 197, 303 and 412 per 1000 person-years from class I through class IV, and HR (95% CI) of 2.06 (1.68 to 2.53), 3.32 (2.74 to 4.03) and 5.28 (4.37 to 6.39) in class II, III and IV, respectively versus class I (p<0.001). Harrell’s C index for mortality was 0.660 (95% CI 0.641 to 0.678). Also 1-year risk of recurrent ED access increased significantly across DSC classes, with rates of 69, 93, 89 and 139 per 1000 person-years from class I through class IV, and HR (95% CI) of 1.46 (1.11 to 1.93), 1.45 (1.10 to 1.92) and 2.60 (1.99 to 3.40) in class II, III and IV, respectively versus class I (p<0.001). For recurrent ED access, Harrell’s C index was 0.604 (95% CI 0.573 to 0.634).
Prediction of short-term mortality
The DSC predicted also 7 and 30-day mortality in both databases. In the SCNP database as a whole, 20 297 and 27 812 participants died within 7 and 30 days, respectively, out of the 180 079 enrolled. In AIDEA, the analysis on 7 and 30-day mortality included 4425 participants, of whom 102 and 345 died within 7 and 30 days, respectively. In both databases and for both follow-up times, short-term mortality increased significantly across DSC classes (p<0.001) (table 4). Harrell’s C index was 0.623 (95% CI 0.618 to 0.627) and 0.639 (95% CI 0.635 to 0.643) for 7 and 30-day mortality, respectively, in the SCNP and 0.690 (95% CI 0.619 to 0.761) and 0.683 (95% CI 0.645 to 0.721) in AIDEA.
Using the SCNP large, representative cohort of older persons accessing the ED and then admitted to hospitals in two Italian regions, we developed the DSC, which combines simple variables, extracted from administrative databases, into a score predicting 1-year mortality. The score was validated against mortality in a random subset of the same cohort, and then its external validity was proven in the newly enrolled sample of the AIDEA study. Remarkably, whereas in the SCNP cohort the DSC was obtained from offline processing of consolidated archives, in AIDEA it was generated in real time on arrival of each eligible patient, thanks to an application implemented in the software routinely used for patients’ clinical management in the EDs of Tuscany Region.
In older persons, because of coexisting chronic comorbidities and diminished life expectancy, fully informed clinical decision-making requires adequate knowledge of patient’s prognosis. In the presence of any given acute condition, comprehensive care planning should be carefully performed, by considering life expectancy, as well as premorbid functional and cognitive status.26 Nevertheless, accurate prognostic assessment of complex older patients is frequently overlooked in clinical practice, especially in the busy daily practice of an ED, where geriatric patients often arrive with non-specific complaints, such as worsening functional status, confusion, dizziness or fall. Their clinical presentation is frequently characterised not by a single, well-defined disease, but rather by an entangled combination of age-related changes, comorbidity, functional and cognitive impairment, polypharmacy and social problems. As a consequence, the risk of wrong triage, incorrect diagnosis and treatment, prolonged ED stay and inappropriate destination is substantial.27 Worsening disability, institutionalisation or death shortly after the ED access may be ultimate consequences.28
Instruments to screen older persons in the ED, based on comprehensive geriatric assessment, do exist,29 but are not used routinely in clinical practice, as they require some expertise and are time consuming. Thus, in spite of their inherent limitations, simple administrative data are increasingly explored as an attractive contribution to prognostic assessment, because they are accurate, objective and easily available at a low cost.7–9 Along this track, we have previously developed the SC, which was shown to predict long-term mortality17 18 and hospitalisations of older persons in the ED18 in two different cohorts, with a discrimination ability comparable to that of other tools requiring face-to-face interview.18 Moreover, the SC allowed risk stratification when comparing the effects of different therapeutic approaches.17 30 However, the algorithm to develop the SC relied only on the first ED admission, so that the score remained constant throughout the following observation time, independent of new hospitalisations or drug prescriptions. To overcome this limitation of the parent tool, the DSC was developed by taking into account each hospital admission as the unit of analysis. This would allow the score to reflect more closely the dynamics of risk status of older patients, changing after each hospitalisation or new drug prescription.
The variables contributing to the DSC can be all recognised as meaningful predictors of prognosis. In particular, the main driving forces in the DSC are advanced age and severe comorbidity, two well-known powerful prognostic factors. Yet, because short-term rehospitalisation usually indicates clinical instability, it is somewhat unexpected that, in the multivariable model reported in table 2, a previous hospital admission between 0 and 30 days had a b coefficient of zero, comparable to no admission. However, the Italian hospital discharge coding system does not allow distinguishing, among short-term rehospitalisations, those that might have been planned at the end of the index hospitalisation (eg, to complete diagnostics or treatments) from those occurring because of clinical instability.
The estimates we obtained are robust, because they were based on high-quality data from a large, randomly split population-based cohort, thus overcoming the problem of overestimation typical of small sample size data sets.31 The predictive ability of the tool persisted even after controlling for discharge diagnosis. As a further confirm, external validation was achieved in the new sample of the AIDEA study, whose participants might be hospitalised or not after the index ED access. It should be emphasised that the DSC predicted also short-term mortality and, at least in AIDEA, recurrent hospitalisations, outcomes for which it had not been created. As suggested by the lower short-term mortality rates reported in table 4, the AIDEA cohort was substantially healthier than the SCNP cohort: this is to be expected, as AIDEA enrolled also participants who were immediately discharged from the ED, without being hospitalised. Nevertheless, the incremental risk of death across DSC classes was comparable in the two studies. Moreover, short and long-term mortality rate figures in the AIDEA cohort compare well to those reported by Tanderup et al: in a Danish study on 3775 subjects aged 65+ years accessing the ED, 30-day and 1-year cumulative mortality increased progressively, from 2.2% to 10.6% and from 8% to 39%, respectively, with the number of geriatric conditions (from 0 to 4), identified on the basis of comprehensive geriatric assessment.27 The similarity between Tanderup’s and our estimates provides further, indirect support to the external validity of the DSC: the absolute mortality risk is indeed crucial, more than HR, to take appropriate decisions on treatment level.
In terms of Harrell’s C index, the discrimination of the DSC was only moderate, similar in the internal and external validation cohorts. Thus, although predictive at the group level, its ability to predict prognosis at an individual level is suboptimal. A possible partial explanation for this finding is that the DSC might misclassify subjects defined at low risk, if they reached the ED with acute life-threatening conditions. Also nursing home residents with severe disability or dementia might be erroneously considered as at low risk by the DSC, if they had no recent hospital admissions: other administrative archives collecting data on nursing home services, which are expected to become promptly available at least in Tuscany, should eliminate this source of misclassification. On the other hand, comparably suboptimal performances have been reported for other tools, proposed in the literature for prognostic assessment in the ED, such as the Identification of Seniors At Risk18 29 (C statistics of 0.65 for 6-month mortality) or the Hospital Frailty Risk Score32 (C statistics of 0.60 for 30-day mortality), or widely used to drive clinical decisions in acute coronary syndromes, such as the TIMI risk score (C statistics of 0.65 for a 14-day composite endpoint).33 This may reflect the fact that, as pointed out by other authors,32 individual outcomes are inherently unpredictable in acute settings, even when risk status is well characterised at a cohort level. This does not detract from the clinical impact of the DSC and similar scores, which help delivery targeted interventions to groups of patients with specific characteristics, whereas, at an individual level, should be used to support, not to substitute clinical judgement.
A 2012 systematic review on prognostic indices for older adults2 reported only one article—our original publication on the parent instrument17—that focused on ED triage of older persons using administrative data, but not in a real-time application. Although other more recent studies dealing with prognostic issues around older persons in an acute hospital setting can be found, none of them was developed for use in the ED nor proposed real-time linkage of administrative data. Thus, to the best of our knowledge, the DSC represents the first example in the literature of real-time utilisation of administrative data for prognostic purposes, at least in—but probably not limited to—older persons accessing the ED.
Several strengths of the DSC should be highlighted. First, it predicts short and long-term prognosis and risk of hospitalisation combining with a simple calculation a limited set of variables, obtained from administrative data: thus, it is objective and does not rely on patient’s collaboration, which is sometime difficult to obtain in the presence of cognitive decline or communication barriers. Second, as the AIDEA study shows, it can be calculated in real time, being therefore immediately available to the ED staff as soon as a target patient is triaged, even before any clinical assessment and with no need of dedicated human resources to gather information. To our knowledge, this represents a unique example of clinical utilisation of administrative data, which might be easily replicated in other hospitals, at least in Italy. Third, it is a non-disease-specific tool, an important characteristic allowing its universal application in older patients, in whom multimorbidity and atypical symptom presentation make clinical and prognostic assessment challenging. Fourth, similar to other prognostic tools based on administrative data, it can be applied also to obtain risk adjustment in the healthcare policy arena: to our knowledge, no previous prognostic tools have been developed to pursue both clinical and epidemiological purposes.
Study limitations are to be acknowledged. The original SCNP data set is relatively old: however, the validity of the score was confirmed in the definitively more recent AIDEA cohort. The score was developed based on data obtained only in subjects admitted to Internal Medicine or Geriatrics units: nonetheless, its validity was confirmed in a completely different series, which included patients visiting the ED for a variety of reasons and most of times discharged from the ED. The data so far available did not consent considering other predictors (eg, claims for specialised visits or other services, or census data), neither different outcomes besides mortality, such as functional or cognitive decline and institutionalisation, which are highly relevant in geriatric patients. However, we are confident that the increasing availability of other administrative databases, as suggested above in reference to nursing home services, will allow incorporating a wider set of predictors as well as of endpoints in the next future.
In conclusion, the DSC, a simple prognostic tool based on administrative data, available in real time in the software used in the ED of Italian hospitals, offers a valid prognostic assessment of older patients, virtually at no additional costs, once the system has been set. This might support clinical decision-making and improve the quality of the care provided. Future studies are needed to assess whether these expectations will be satisfied.
Contributors DB conceived the study, conducted data management and analysis, and contributed to drafting of the manuscript. GC conducted data analysis and contributed to drafting of the manuscript. FT contributed to obtaining the data and drafting of the manuscript. LDE and PM obtained the data and contributed to data management. AU and EB contributed to obtaining the data and provided institutional support. LG and GL provided institutional support. RB and NM provided institutional support and contributed to developing the research question. MDB conceived the study, wrote the protocol, supervised the analyses and drafted the manuscript. All authors read and approved the final manuscript. All the authors were independent of the funders, had access to all the data and take responsibility for the integrity of the data and the accuracy of the data analysis.
Funding The 'Silver Code National Project' was funded by the Centre for Disease Control of the Italian Ministry of Health. The 'Anziani in DEA' study was jointly funded by the Italian Ministry of Health and by the Tuscany Region (RF-2010-2321801).
Disclaimer The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
Competing interests None declared.
Patient consent for publication Not required.
Ethics approval Following the Italian legislation on observational studies using administrative databases in place at the time of data acquisition, ethics committee approval and participants’ subscription of informed consent was waived in the SCNP, because only anonymised administrative data were extracted. The protocol of the AIDEA study, which included also face-to-face interview (not reported in the present work), was approved by the Ethics Committee of the Azienda Ospedaliero-Universitaria Careggi (976/13_AOUC).
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available upon reasonable request.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.