Usual physical activity and subsequent hospital usage over 20 years in a general population: the EPIC-Norfolk cohort

Background While physical activity interventions have been reported to reduce hospital stays, it is not clear if, in the general population, usual physical activity patterns may be associated with subsequent hospital use independently of other lifestyle factors. Objective We examined the relationship between reported usual physical activity and subsequent admissions to hospital and time spent in hospital for 11,228 men and 13,786 women aged 40–79 years in the general population. Methods Participants from a British prospective population-based cohort study were followed for 20 years (1999–2019) using record linkage to document hospital usage. Total physical activity was estimated by combining workplace and leisure time activity reported in a baseline lifestyle questionnaire and repeated in a subset at a second time point approximately 12 years later. Results Compared to those reporting no physical activity, participants who were the most active had a lower likelihood of spending more than 20 days in hospital odds ratio (OR) 0.88 (95% confidence interval (CI) 0.81–0.96) over the next 20 years after multivariable adjustment for age, sex, smoking status, education, social class and body mass index. Participants reporting any activity had a mean of 0.42 fewer hospital days per year between 1999 and 2009 compared to inactive participants, an estimated potential saving to the National Health Service (NHS) of £247 per person per year, or approximately 7% of UK health expenditure. Participants who remained physically active or became active 12 years later had lower risk of subsequent hospital usage than those who remained inactive or became inactive, p-trend < 0.001. Conclusion Usual physical activity in this middle-aged and older population predicts lower future hospitalisations - time spent in hospital and number of admissions independently of behavioural and sociodemographic factors. Small feasible differences in usual physical activity in the general population may potentially have a substantial impact on hospital usage and costs.


Introduction
Historically UK government spending on health has risen on average by 3.7% per year since 1948, outpacing economic growth over the period [1,2]. As a result, health expenditure as a proportion of UK Gross Domestic Product (GDP) has increased from 3.6 to 7.5% over the same period. Approximately a half of government health expenditure is used for hospitals [3]. There are many factors which may influence hospital usage, not all of which are related to ill health while increases in expenditure are only partly explained by demographic changes [4]. Changes in modifiable lifestyle factors have the potential to lower hospital length of stay. There is growing evidence of the effectiveness of preoperative exercise programmes and other pre-admission interventions in reducing hospital length of stay and readmission rates [5][6][7][8][9] but it is unclear whether in the general population, usual physical activity is related to hospital use. Long-term randomised controlled trials (RCTs) of physical activity interventions with health endpoints are not generally feasible, so evidence is largely based on observational studies.
Physical activity is associated with lower rates of mortality from all causes and cardiovascular disease [10][11][12]. It is also associated with a lower risk of many non-fatal diseases [13][14][15][16] but few studies have examined the relationship between usual physical activity in middle and later life and subsequent hospital usage the general population [17]. The measurement of usual physical activity is problematic. Objective measurements, such as accelerometry have only been developed relatively recently and hence studies based on large, free-living, community-based populations with long follow-up have used self-reported activity from questionnaires. Studies with longer follow-up time are less likely to be affected by reverse causality, which is a feature of studies with short duration of follow-up where individuals who report low physical activity at baseline are inactive by virtue of being affected by the outcome of interest. Selfreported physical activity is most often assessed by questions related to leisure-time activities [18,19]. Few studies capture both occupational and leisure-time activity.
Hospital usage can be measured by total admissions and length of stay over a fixed follow-up period. These non-disease specific outcome measures can be used to examine the overall level of health service usage [20]. Ageing populations put ever-increasing pressure on health care services and it is therefore important to establish if modest differences in modifiable lifestyle behaviours such as physical activity are related to hospitalisation [21][22][23][24].
This study examines the relationship between measures of physical activity using a validated physical activity scale, change in physical activity, and subsequent hospital usage, in older men and women living in the general community over a 10-year period, and a subsequent 10-year follow-up period, taking into account a range of demographic and lifestyle factors.

Materials and methods
The European Prospective Investigation into Cancer in Norfolk (EPIC-Norfolk) is a general population cohort study of men and women aged 40-79 years living in Norfolk recruited from general practices between 1993 and 1997. The response rate for recruitment was approximately 40%. The cohort has similar characteristics to national population surveys except for a lower prevalence of current smokers [25]. The study has ethics committee approval and all participants gave informed, signed consent for study participation including access to medical records. The cohort is flagged for mortality and hospital admissions from linkage to national databases held by NHS Digital and hence there is virtually no loss to follow-up.
At recruitment, participants completed a lifestyle questionnaire where they were asked about their occupational and leisure physical activity. Occupational activity was assessed using a four category question ("sedentary", "standing", "moderate physical work" and "heavy manual work") with examples such as office worker, shop assistant, plumber and construction worker respectively. Leisure activity in both summer and winter was assessed from the number of hours per week spent cycling, attending keep fit classes or aerobics and swimming or jogging. Estimated average hours of leisure activity was calculated as the mean of summer and winter activities and categorised using 0, (0,3.5], (3.5,7] and > 7. A combined score, divided into four ordered categories with individuals labelled as "inactive", "moderately inactive", "moderately active" and "active" was created combining leisure and occupational elements. Those who did not complete the activity question were placed in the inactive category. The score was validated against energy expenditure measured by free-living heart rate monitoring with individual calibration [26]. It has been reported to predict all-cause mortality and cardiovascular disease incidence [27]. Participants attending the baseline health examination had their height to the nearest 0.1 kg measured using a stadiometer (Chasemores, UK) and their weight to the nearest 100 g measured in light clothing without shoes (Salter, West Bromwich, UK). Body mass index (BMI) was calculated using measured weight in kilograms divided by the square of measured height in square metres. Two yes/no questions were used to derive smoking status: "Have you ever smoked as much as one cigarette a day for as long as a year?" and, where a positive response was given, "Do you smoke cigarettes now?" Participants also completed questions about their employment and that of their partner with details of both current and past employment recorded. Occupational social class was defined according to the Registrar General's classification [28,29]. A list of common UK qualifications was used to establish educational attainment and participants were asked to mark all relevant qualifications. These were then categorised using the highest qualification attained. Participants were asked at baseline "Has the doctor ever told you that you have any of the following?" followed by a list of common conditions including "Heart attack (myocardial infarction)", "Stroke" and "Cancer".
Surviving participants were invited to complete a lifestyle questionnaire and attend a health examination (second time-point, "TP2") between 2006 and 2011 [30] . Questions on physical activity and cigarette smoking, similar to those at baseline, were included in a postal questionnaire, completed by a subset of 9827 of the original cohort. Weight and height were measured on 8094 by clinic staff and body mass index calculated in the same way as at baseline described previously.

Ascertainment of hospital usage through record linkage
The National Health Service (NHS) in Britain treats residents without charge at the point of service so covers virtually all major health service usage. The EPIC-Norfolk cohort was regularly linked to hospital records from 1999 onwards as previously reported [20]. Briefly, NHS numbers were used to perform linkage to hospital databases between 1999 and 2019. Initially, up to 2009, linkage was made via the East Norfolk Primary Health Care Trust while later, national databases held by NHS Digital were used [31]. All hospital activity for EPIC-Norfolk participants was captured wherever they were treated in England and Wales. Hospital episode statistics (HES) records which included admission and discharge dates were used to calculate time in hospital and numbers of admissions. Contiguous admissions were merged and counted as a single admission.

Statistical analysis
For the main analysis, 625 men and women who died before 1999 were excluded. Dichotomous variables were created for the socioeconomic status variables. Professional, managerial and technical and non-manual skilled occupations (codes I, II and IIIa respectively) were classed as non-manual while manual skilled, partly skilled and unskilled (codes IIIb, IV and V respectively) were classed as manual. Educational attainment was categorised into "Higher education level" (which includes those with qualifications at secondary level or above) and "Lower education level" (those with no qualifications). The numbers of individuals with missing values for covariables were: 53 BMI, 218 smoking status, 545 social class, 18 education level. Validation of the physical activity measures [26] suggested that participants with missing data be classified inactive.
Logistic regression was used to model hospitalisation outcomes on physical activity category, adjusting for covariables. Several dichotomous outcome categories were calculated based on total admissions and length of stay spanning two periods: 1999-2009 (10-year followup) and 1999-2019 (20-year follow-up). Total admissions from 10-year follow-up were used to define "any hospital admissions" and "7 or more admissions" while length of stay from 10-year follow-up was used to create "greater than 20 hospital days". These thresholds were chosen to represent those with higher levels of hospital usage and were consistent with previous work [20]. Dichotomous outcome categories based on 20-year follow-up and having approximately the same proportion of the population as their 10-year follow-up counterparts include "12 or more admissions" and "greater than 50 hospital days" while "7 or more admissions" and "greater than 20 hospital days" were also calculated for this period to serve as a comparison. Hospital days are defined as the sum of total bed days (overnight stays) and day-cases. Linear regression was used to calculate the absolute difference in adjusted mean bed days between inactive participants and participants reporting any activity.
To address change in physical activity, we also used physical activity measured at TP2 approximately 12 years later as a second baseline. We excluded 105 participants who died prior to 2009, leaving 9722. Multiple imputation was used to address missing values, in particular for body mass index at TP2 where data for 1733 were not available for participants who completed a TP2 questionnaire but did not attend a health examination. Predictive mean matching with 5 multiple imputations and 50 iterations was used with baseline variables BMI, occupational social class and education attainment and TP2 current smoking. Changed-activity categories use combinations of physical activity categories at the baseline and TP2. The category shown as "Inactive/Inactive" is the set of participants who reported being inactive at baseline and remained inactive when asked again at TP2. The group who initially reported any activity but became inactive later is shown as "Any-activity/Inactive" while the other two categories "Inactive/Any-activity" and "Any-activity/Any-activity" were similarly defined.
The cost to the NHS of one bed-day is £496, calculated using the Reference Costs for English Hospitals 2017/18 for elective (5.4 £bn) and non-elective (18 £bn) admissions [32] and the total available beds (approximately 129,200) [33]. The cost per hospital day (overnight stays and day-cases) is £587 when the cost of daycase activity is included (4.4 £bn per year). The reported OECD UK per capita expenditure on health in 2017, was £3375 (exchange rate at the time of writing) [34]. Perperson costs were calculated by multiplying the cost per hospital day and hospital days per person. Percentage of NHS per-capita health expenditure was calculated as the ratio of per-person cost and OECD UK per-capita expenditure.
Adjusted mean hospital days by physical activity category were determined first by calculating hospital days for each one year period restricted to participants surviving to the start of the given year. Linear regression of hospital days on physical activity adjusted for age, sex, occupational social class, educational attainment, current smoking and body mass index was then used. Adjusted means by category were obtained using estimated marginal means. The overall mean difference of days was calculated by taking the mean of the annual differences for each of two periods (1999-2009 and 2009-2019).
Sensitivity analyses were conducted in which the physical activity exposure was dichotomised into inactive and any-activity groups, using the outcome more than 20 hospital days over the period 1999-2019. Multivariable adjusted odds ratios were examined, stratified by sex, age < 65 and ≥ 65 years, manual and non-manual social class, lower (no qualifications) and higher level of education, former or never smoking and current smoking, BMI ≤30, > 30 kg/m 2 , chromic disease (heart attack, stroke or cancer) and no reported chronic disease, survival to the end of follow-up (March 2019) and died during follow-up period. A further multivariable model was performed using the narrower follow-up period of 2004-2019, a minimum of five years after participants reported their level of physical activity excluding participants who died prior to 2004.
All analyses were performed using the R statistical language (R Foundation for Statistical Computing, Vienna, Austria version 3.6.0 with packages ggeffects, knitr, Gmisc, tidyverse, intubate, mice).

Results
Characteristics of the study population according to the four categories of physical activity score are described in Table 1. Active participants tend to be younger, nonsmokers, without chronic disease and have higher educational attainment, however those with manual social class also tend to be more active.
Prevalent disease is self-reported heart attack, stroke or cancer at baseline. Higher education level represents those with qualifications to at least secondary level.
In Table 2 odds ratios are shown first age and sex adjusted and then additionally adjusted for social class, educational attainment, BMI and smoking status. For the 10-year follow-up period 1999-2009, outcomes of any hospital admission, 7 or more hospital admissions and more than 20 days stay in hospital are shown according to the baseline physical activity score. The multivariable-adjusted models indicate that participants with a physical activity score of at least moderately inactive had fewer hospital admissions and fewer days in hospital, than those who were inactive. The associations for inactive vs active were OR 0.73 (95% CI 0.65-0.82) p-trend < 0.001 across activity score for seven or more hospital admissions and OR 0.75 (95% CI 0.67-0.83) ptrend < 0.001 for more than 20 hospital days.
Physical activity category at TP2 baseline was determined in 9827 men and women. The associations for inactive vs active for 20 hospital days over the subsequent 10-year follow-up period (2009 to 2019) were stronger than those for the first 10-year follow-up period OR 0.60 (95% CI 0.50-0.72) p-trend < 0.001 and similar for 7 or more admissions OR 0.77 (95% CI 0.64-0.91) p-trend 0.001. Table 3 shows multivariable-adjusted odds ratios for outcome more than 20 hospital days during the 1999-2019 follow-up in participants who were inactive compared to those reporting any activity at baseline, stratified by key variables in subgroups. The directions of the associations did not differ by subgroup. Higher inverse associations were seen in women, in the under 65 s, in those with no chronic disease at baseline and those surviving to the end of follow-up although confidence intervals overlapped in each case. Table 3 also shows that the association for the period 2004-2019, excluding the first 5 years of the outcome period was OR 0.93 (95% CI 0.87-1.00). Table 4 shows odds ratios by all combinations of change in physical activity category between baseline and TP2 were determined using the TP2 baseline and subsequent 10-year follow-up. The multivariableadjusted odds ratios comparing "Inactive/Inactive" (the reference) and "Any-activity/Any-activity" were OR 0.66 (95% CI 0.57-0.77) p-trend < 0.001 across changedactivity categories for more than 20 hospital days and OR 0.91 (95% CI 0.78-1.07) p-trend 0.026 for seven or more hospital admissions. Participants who remained physically active or became active had lower risk of subsequent hospital usage than those who remained inactive or became inactive.
Supplementary Table S1 shows all terms in a series of multivariable logistic regression models for inactive physical activity (vs any-activity) and various dichotomous outcomes over the period 1999-2019 for all, men and women. Covariables age per 10 years, manual social class, lower education level, current smoking and BMI >         Table S2 shows the adjusted mean hospital days for inactive and any-activity participants by year, and the absolute difference in days between the categories. The mean of the differences was calculated for 1999-2009 using baseline physical activity and 2009-2019 using physical activity at TP2 and cumulative costs were determined assuming £587 per hospital day.
The difference in multivariable adjusted mean hospital days between inactive participants and participants reporting any activity was 0.42 days per year over the first 10 years of follow-up, an estimated potential saving to the NHS of £247 per person per year or approximately 7% of health expenditure. The difference in hospital days over the subsequent 10 years (2009-2019) was slightly higher, with any-activity participants having 0.46 fewer hospital days, an estimated potential saving of £268 or approximately 8% of health expenditure.

Discussion
Usual physical activity assessed at baseline survey in 1993-1997 was inversely associated with future hospital usage independently of sociodemographic and lifestyle factors in this middle-aged and older cohort of men and women over a 20-year follow-up period. Compared to study participants who were inactive, active participants had a lower likelihood of having more than 20 hospital days or more than 12 admissions. Stronger associations were seen over a 10-year follow-up period with moderate inactivity or greater being associated with lower risk of seven or more hospital admissions or more than 20 hospital days. There was a dose response over physical activity categories over both the 10-year and 20-year follow-up periods for both hospital duration and number of admissions. There are a number of possible explanations for these findings.

Strengths and limitations of study
Reverse causality may partly explain the associations we observed. Participants may be physically inactive occupationally or less able to take part in leisure time activity because of known or preclinical illness which may also predispose to increased later hospitalisation [35]. However, sensitivity analyses excluding those with a selfreported chronic disease at baseline (heart attack, stroke or cancer), who might have lower physical activity, did not differ materially from the main findings. Also, a sensitivity analysis excluding hospital admissions occurring in the first 5 years of follow-up (the period 2004-2019), that is, those who were more likely to have preclinical illness and lowered physical activity, again did not show materially different associations.
Confounding is a major issue in examining the relationship between lifestyle factors and health outcomes. Individuals who are more physically active are likely to differ from those who are less active with respect to other factors relating to the likelihood of future hospitalisation including age, sex, smoking, body mass index, social class and education. However, the associations were consistent after multivariable-adjustment for these factors and after stratification by these potential confounding variables. As we examined total hospital usage over long time periods, individuals who died during the follow-up period did not use hospital services for the full period. This may have affected the results if there was differential mortality by physical activity, whereby study participants who were inactive were more likely to have died earlier than the more active participants and hence less likely to use hospital services for the full follow-up period. Sensitivity analysis models restricted to those surviving to the end of 20-year follow-up showed stronger associations of physical activity with lower hospital use than models using the whole population, including those who died during the follow-up period, suggesting there was some attenuation due to selective follow-up.
This study has several strengths. Few studies have examined the physical activity of middle-aged and older men and women and their subsequent healthcare utilisation. The literature falls into two groups, studies based on exercise interventions and observational studies. While most intervention studies provide some evidence that a physically active lifestyle improves health, intervention protocols vary and differences in dropout rates between groups in RCTs limit generalisability [36]. Intervention studies may also typically have smaller study size and shorter follow-up time and while observational studies are generally larger, there are few studies comparable in size to the present study. Our study, being well characterised, allowed adjustment for a broad range of relevant factors. We also used linked hospital data and did not depend on self-reported outcome data. Many studies are based on particular population groups or particular disease outcomes and some rely on selfselection to exercise programs. Few studies examine free living community-based populations [37,38], however we used a prospective cohort design and were able to examine hospital usage over a long follow-up period with a reliable population-based denominator.
Our study was based on a free-living population of older men and women living in the general community in the United Kingdom where the NHS provides health care free at the point of delivery. Potential major confounders such as income, and ability to pay that might therefore affect and limit access and use of health services, are less likely to apply in this study. The NHS also enables record linkage for virtually complete follow-up of the population. Though admissions to private hospitals in Norfolk were not included in our data which only counts NHS hospitals, the use of private hospitals in Norfolk was minimal in comparison with the use of NHS facilities.
Measurement of usual occupational and leisure time physical activity was assessed using a self-reported questionnaire. Objective measures such as accelerometry and similar techniques were not available when the EPIC-Norfolk cohort was recruited. However, the physical activity score used was previously validated using heart rate monitoring with individual calibration and based on both occupational and leisure-based components of physical activity.
It is also clear that a single measurement of physical activity is insufficient to determine accurately usual levels of activity over the life course. Events such as retirement or illness or progressive ageing related conditions such as frailty may result in a change to the amount of physical activity undertaken [39]. While we are unable to establish the length of time over which consistent physical activity was maintained, we were able to examine longitudinal measurements of physical activity at two time-points in a subset of participants. The associations observed at the later time-point were comparable with (in fact stronger than) those observed at the first time-point, despite the cohort mean age being approximately 10 years older and having a much higher proportion of retirees. Change in behaviour over the 20-year follow-up period is a more likely explanation for the attenuated associations observed, rather than age or employment status. Participants who remained inactive or became inactive had the highest risk of subsequent hospitalisation. Additionally, random measurement error is likely only to attenuate associations, and therefore unlikely to explain any of the associations observed between physical activity and hospitalisation.

Comparison with other studies
Physical activity has been associated with many health benefits including protection against cardiovascular [10,12,40] and many other chronic diseases [41,42] so there are many plausible reasons why it might also be associated with lower hospitalisations from individual conditions. Chronic conditions such as cardiovascular disease remain leading causes of hospitalisation. However, in this study, we were able to examine total hospital usage in a general population irrespective of cause of admission.
Small increases in physical activity have been reported to obtain cost savings for health services by reducing hospital admissions [17,22,[43][44][45] with many studies reporting reductions of length of stay after preoperative physical activity interventions. Our study has observed a 12-13% lower risk of long stay and high numbers of admissions by physical activity category. The mean difference in bed days between inactive and any-activity participants in our study was 0.42 days per year over the first 10 years of follow-up. Assuming a cost of £587 per hospital day (inpatient bed-days and day-cases), the potential saving to the NHS is approximately £247 per person per year for every inactive person who starts to undertake at least some exercise, or about 7% of UK per capita health expenditure. Similar results were observed 10 years later when participants were aged 50-90 years. Calculations such as these are unavoidably crude but serve to illustrate the significant financial contribution, when scaled nationally, that modest changes in lifestyle can achieve quite apart from the obvious personal gain from the reduction in risk of being hospitalised.
While there is evidence suggesting that pre-admission physical activity programmes may lower duration of hospital stay [5][6][7][8]46], these are short term, requiring resources and targeted at only a limited number of individuals. Our data indicate that usual physical activity patterns in the general population predict hospital usage over the subsequent 2 decades.

Conclusions and policy implications
Usual physical activity in this middle-aged and older population predicts lower future hospitalisations -time spent in hospital and number of admissions independently of behavioural and sociodemographic factors. Small, feasible differences in usual physical activity in the general population may potentially have a substantial impact on hospital usage and costs.
Additional file 1: Table S1. Multivariable logistic regression of risk factors for any hospital admissions, ≥7 hospital admissions and > 20 days of hospital stay from 1999 to 2019 in 25,014 men and women. Table S2. Adjusted mean hospital days by physical activity category for two periods, mean difference in days and cumulative cost, 1999-2009 using baseline physical activity and 2009-2019 using physical activity at TP2. Figure S1. Flow diagram showing numbers of participants at various stages in the EPIC-Norfolk study including invitations, consents, attendance at health examinations and questionnaire completion.
Abbreviations BMI: Body mass index; CI: Confidence interval; EPIC-Norfolk: European Prospective Investigation into Cancer in Norfolk; GDP: Gross Domestic Product; HES: Hospital episode statistics; OECD: Organisation for Economic Co-operation and Development; OR: Odds ratio; TP2: Time-point two access to all the data. No conflicts of interest were declared by any author (apart from the two grants) and they have no other relationships or activities that could appear to have influenced the submitted work.