Article Text

Original research
Food insecurity status is of added value in explaining poor health: a cross-sectional study among parents living in disadvantaged neighbourhoods in the Netherlands
  1. Laura A van der Velde1,
  2. Ewout W Steyerberg2,
  3. Mattijs E Numans1,
  4. Jessica C Kiefte-de Jong1
  1. 1Department of Public Health and Primary Care/ LUMC-Campus The Hague, Leiden University Medical Centre, Leiden, The Netherlands
  2. 2Department of Biomedical Data Sciences, Leiden University Medical Centre, Leiden, The Netherlands
  1. Correspondence to Professor Jessica C Kiefte-de Jong; J.C.Kiefte{at}lumc.nl

Abstract

Objectives The aim of this study was to examine the added value of food insecurity in explaining poor physical and mental health beyond other socioeconomic risk factors.

Design, setting, participants and outcome measures Data for this cross-sectional study were collected using questionnaires with validated measures for food insecurity status and health status, including 199 adult participants with at least 1 child living at home, living in or near disadvantaged neighbourhoods in The Hague, the Netherlands. To assess the added value of food insecurity, optimism-corrected goodness-of-fit statistics of multivariate regression models with and without food insecurity status as a covariate were compared.

Results In the multivariable models explaining poor physical health (Physical Component Summary: PCS) and mental health (Mental Component Summary: MCS), from all included socioeconomic risk factors, food insecurity score was the most important covariate. Including food insecurity score in those models led to an improvement of explained variance from 6.3% to 9.2% for PCS, and from 5.8% to 11.0% for MCS, and a slightly lower root mean square error. Further analyses showed that including food insecurity score improved the discriminative ability between those individuals most at risk of poor health, reflected by an improvement in C-statistic from 0.64 (95% CI 0.59 to 0.71) to 0.69 (95% CI 0.62 to 0.73) for PCS and from 0.65 (95% CI 0.55 to 0.68) to 0.70 (95% CI 0.61 to 0.73) for MCS. Further, explained variance in these models improved with approximately one-half for PCS and doubled for MCS.

Conclusions From these results it follows that food insecurity score is of added value in explaining poor physical and mental health beyond traditionally used socioeconomic risk factors (ie, age, educational level, income, living situation, employment status and migration background) in disadvantaged communities. Therefore, routine food insecurity screening may be important for effective risk stratification to identify populations at increased risk of poor health and provide targeted interventions.

  • mental health
  • epidemiology
  • nutrition & dietetics
  • preventive medicine
  • social medicine

Data availability statement

Data are available on reasonable request.

http://creativecommons.org/licenses/by-nc/4.0/

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

Request Permissions

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.

Strengths and limitations of this study

  • Socioeconomic risk factors such as educational level, household income level, living situation, employment status and migration background are associated with poor health, but the ability to explain poor health with these traditional socioeconomic risk factors is limited.

  • Our study is among the first to investigate the value of assessing food insecurity and adding this to traditional social determinants of health for explaining poor physical and mental health.

  • Food insecurity is a relatively understudied area in the Netherlands, and the presented results can stimulate larger scale, routine screening for food insecurity in the Netherlands.

  • Our study population mainly included women living in a disadvantaged urban setting, and therefore, the results may not be generalisable to the general Dutch population.

  • Our study is strengthened by the use of validated measures of our main outcome and covariate and by accounting for statistical optimism in our multivariate models, however, future studies are warranted to externally validate our results to verify your findings, also in other populations and settings.

Introduction

It has been extensively shown that individuals of lower socioeconomic position (SEP) groups generally have poorer health outcomes.1 Therefore, improving health in these groups and being able to identify those that are most at risk of poor health has great potential for improving population health. An emerging concept in aiming to improve population health is population health management, which strives to simultaneously improve population health, improve experienced quality of care (by both the patient and healthcare provider), and reduce healthcare costs (referred to as the Quadruple Aim).2 A crucial element of effective population health management is risk stratification: identification of populations that are most at risk. In risk stratification, several biomedical and social characteristics of individuals can be combined to establish a risk profile towards poor health outcomes or healthcare utilisation. This can be used to proactively identify populations at increased risk of poor health and target prevention (or care) resources specifically to these populations in order to improve successfulness and cost-effectiveness of interventions.3 Predictive modelling is a method that can be used to identify populations at increased risk of poor health and can therefore be used for risk stratification.3

Many factors have been identified as risk factors in the association between lower SEP and poor health.4–8 Even though numerous studies have examined these associations with poor health, the ability to explain or predict poor health with traditional risk factors and social determinants of health (such as employment status, educational level and income9) often proves to be limited. Therefore, we hypothesise that less traditional social determinants of health, such as food insecurity, might be worthwhile to include in models aiming to explain poor health as a proxy to better identify risk groups and to be used for improving integration of social needs–informed care into medical care.10 11

Food insecurity can be defined as an insufficient physical and economic access to adequate food that meets dietary needs and food preferences.12 Food insecurity is a public health concern facing low-income, middle-income and high-income regions, including Europe: a large global study found a food insecurity prevalence of 25% across 39 European countries.13 Food insecurity can be considered as an adverse health outcome in itself, but also a determinant of poor health,11 14 and food insecurity is associated with increased healthcare utilisation and costs, even when socioeconomic factors are taken into account.15 To date, few studies have focused on food insecurity prevalence in the Netherlands. These studies indicate a food insecurity prevalence of approximately 25% among people living in an urban disadvantaged setting, and 70% among foodbank recipients.16 17 Also in the Netherlands, living on a low income is associated with poorer health. However, living on a low income is not one-on-one related to experiencing food insecurity, as the latter reflects not only a scarcity of financial means to acquire adequate food, but among others also induces psychosocial stress.14

Therefore, we hypothesise that it is worthwhile to include food insecurity for better explaining health outcomes in addition to traditional social determinants such as income, to better identify people most at risk of poor health. In the current study, we aim to explore the value of assessing food insecurity and adding this to traditional social determinants of health for better explaining poor physical and mental health.

Methods

Study design and population

Data for this cross-sectional study were collected between April 2017 and June 2018. This study was conducted among families living in highly urbanised disadvantaged neighbourhoods in the Dutch city The Hague. Participants were actively recruited at various public places, such as community centres, in four preselected disadvantaged neighbourhoods, based on criteria already in use by the Dutch Government to identify disadvantaged neighbourhoods.18 Participants were eligible for the study if they were living in or near one of the selected disadvantaged neighbourhoods; were aged ≥18 years; and had at least one child aged <18 years living at home. Only one parent per household could participate. A total of 199 participants were included in the current study.

Patient and public involvement

Participants were not involved in the design, or conduct, or reporting, or dissemination plans of our research.

Data collection

Data collection was done using paper-based or online questionnaires, available in the Dutch, English and Turkish language. Most participants completed the questionnaire and informed consent form at the site of recruitment immediately after being invited to the study. Participants were offered help completing the questionnaire if they had difficulty reading or writing. If participants provided contact information, they were contacted by phone or email to complement missing data from their questionnaire if applicable.

Primary outcome assessment: general health status

The primary outcome of our models is general health status, assessed using the 12-Item Short Form Health Survey (SF-12).19 The SF-12 consists of two summary scores: the Physical Component Summary (PCS) score and the Mental Component Summary (MCS) score. The SF-12 is a widely used, reliable and validated instrument with a relative validity ranging from 0.63 to 0.93 for the 12-item PCS, and 0.60 to 1.07 for the 12-item MCS compared with the best SF-36 scale in an adult population.19 The SF-12 assesses self-rated general health and therefore reflects the subjective perception of how physically (PCS) and mentally (MCS) healthy a person feels. In our analyses, we used the two continuous summary scores of general health status: the PCS and MCS. PCS and MCS scores were created according to the SF-12 scoring guide by Ware et al.20 The PCS and MCS scores range from 0 to 100, and these scores were reversed so that higher scores represent poorer health. The PCS and MCS are scored using norm-based methods. In both summary scores all SF-12 items are included, but different weights are assigned to each SF-12 item for the PCS and MCS score calculations. These item weights are chosen so that both scores have a mean of 50 and a SD of 10 in the general US population, as described in the SF-12 scoring guide by Ware et al.20 An advantage of using this norm-based scoring is that it enables comparison of our results and to interpret them in relation to scores in the general US population and across other studies using the same scoring weights.20 For instance, scores above 50 indicate a better health than the general US population and scores below 50 indicate a poorer health than the general US population.

Previous literature clearly shows that poorer PCS and MCS scores are associated with higher healthcare costs.21 To enable evaluation of the discriminative performance of our models, we also dichotomised the PCS and MCS into scores below 50 and scores above 50, where scores above 50 reflect poorest physical and mental health and therefore highest expected healthcare use and costs.21 22

Food insecurity status assessment

Household food insecurity status was assessed using the 18-item United States Department of Agriculture Household Food Security Survey Module (USDA-HFSSM).23 The original USDA-HFSSM was translated from the English to the Dutch language based on the translation by Neter et al16 who applied the translation and back-translation technique.16 In the survey, conditions and behaviours that are characteristic for households having difficulty meeting basic food needs are addressed, with the past 12 months as reference period. Affirmative responses to these questions were summed, resulting in a continuum of food insecurity score ranging from 0 to 18, with higher scores reflecting a higher food insecurity. The food insecurity score was dichotomised into ‘food secure’ (FS: 0–2 affirmative responses), and ‘food insecure’ (FI: 3–18 affirmative responses), according to the USDA standards.23

Sociodemographic and lifestyle variables assessment

Sociodemographic and lifestyle information was collected, including age or date of birth, sex, height, weight, gross monthly household income, marital status, educational level, country of birth of the participant and their parents, employment status, smoking status and presence of common lifestyle-related diseases and medication use. Detailed information on how these data were used to calculate and categorise age, Body Mass Index (BMI), household income, educational level, employment status, living situation and migration background, is described elsewhere.17

Further, the presence of the following common health issues was assessed: high blood pressure, high cholesterol, surgery on the heart, heart attack, asthma, chronic obstructive pulmonary disease, diabetes mellitus (participants could additionally specify whether it was type 1 or 2) and anaemia (in the previous 12 months). Additionally, obesity status was included (ie, BMI >30 kg/m2). The total number of present health issues was calculated as a reflection of comorbid health issues.

Covariates explaining poor health

We selected age (in years, continuous), educational level (low/higher), household income level (below/above basic needs budget), living situation (partner/single), employment status (currently employed/not currently employed) and migration background (Western/non-Western) as covariates explaining poor health. These covariates were selected on the basis of variables routinely assessed in health monitors of the Netherlands.24 Food insecurity score and food insecurity status (FS/FI) were included as covariates to assess their added value in explaining poor health.

Statistical analysis

Power calculation

The current study describes secondary analyses of our study on food insecurity and obesity,17 for which a conservative power calculation was performed based on obesity prevalence. For the current study, we compared 150 FS to 49 FI participants. With an alpha of 0.05, the power was more than 90% to detect a difference in health outcomes of 5.8–7.6 points with SD of 8.3–11.3. For reliable explanatory and prediction modelling, we generally need at least 2 subjects per variable with a continuous outcome; with 199 participants, our number of subjects per variable was well over the minimum required number.25

Population description

Participant characteristics were described for the total population and separately for participants that reported their health being fair to poor and good to excellent. Continuous variables were reported as median and IQR. Categorical variables were described as frequencies and percentages.

Models explaining poor physical health (PCS) and mental health (MCS)

First, the crude associations between all separate covariates (age, educational level, household income level, living situation, employment status, migration background, food insecurity score and food insecurity status) and the individual outcome measures PCS and MCS were assessed using bivariate linear regression models. Second, two separate multinomial linear regression models were built with both PCS and MCS as individual outcome variables, including all selected covariates except food insecurity score. Third, the same methods as described above were repeated but now additionally including food insecurity score as a covariate.

For the multivariate models, besides the β-coefficients also the standardised β-coefficients were presented to enable a comparison of the relative importance of each covariate. The relative importance of the food insecurity score in explaining poor health would be reflected by a relatively high standardised β-Coefficient.

The potential added value of including food insecurity score in explaining poor health is reflected in an improvement in the goodness-of-fit statistics, namely R-squared (R2) and the root mean square error (RMSE). R2 presents the proportion of variance in the dependent variable that can be explained by the independent variables. R2 indicates the percentage of the total variation observed for PCS and MCS that can be explained by the model (a value of 0 indicates that the model explains none of the variation in PCS and MCS, while a value of 1 indicates that the model explains all of the variation). An increase in R2 and a decrease in RMSE after adding food insecurity score to the model, would imply that adding food insecurity score to the model improves its performance.

Discriminative performance

The power of the model to discriminate between those individuals most at risk of poor health and associated healthcare use and costs was evaluated by building additional models using logistic regression, including the same covariates as described above but with dichotomous outcome measures of PCS and MCS (ie, PCS and MCS scores below or above 50). The discriminative performance of the logistic regression models was presented by the C-statistic and Nagelkerke’s R2.26

The C-statistic is an indicator of how well the model can discriminate between the two groups and it ranges from 0.5 (no discrimination) to 1.0 (perfect discrimination). The C-statistic represents the area under the receiver operating characteristic curve. Herein, the sensitivity (percentage of persons that correctly is predicted to have poor health) is on the y-axis and one minus the specificity (percentage of persons that correctly is predicted not to have poor health) on the x-axis. Nagelkerke’s R2 is an adjusted version of the Cox and Snell R2 so that it ranges from 0 to 1. It can be interpreted similarly to the R2 as described above, that is, higher values indicate a larger proportion of variance in the dependent variable that can be explained by the independent variables. The added value of including food insecurity score to discriminate between those individuals most at risk of poor health is reflected by an improvement in the C-statistic and Nagelkerke’s R2.

Internal validation to estimate optimism-corrected model performance

We used the same dataset to fit the models and to assess the validity of the model, which can lead to optimistic estimates of the model performance (ie, statistical optimism).27 All performance measures (ie, R2, RMSE, the C-statistic and Nagelkerke’s R2) were, therefore, adjusted for statistical optimism by a bootstrap resampling and cross-validation procedure (n=1000). With this procedure, we estimate the loss in predictive accuracy of our model in a new sample and correct for this. Bootstrapping included resampling with replacement from the original sample.28 To correct for the statistical optimism, the performance measures of a model in a bootstrapped sample and the original sample was compared and the average difference between the performance measures of these samples was used as the optimism bias. This optimism was subtracted from the original performance measures to obtain the optimism-corrected performance measures.28 29

Multiple imputation

Multiple imputation was used to reduce potential bias associated with missing data in our study. Missing data were imputed and 10 independent datasets were created using fully conditional specification (Markov chain Monte Carlo method) with a maximum of 10 iterations. Predictive mean matching was used for non-normally distributed variables and logistic regression models for categorical variables. A more detailed description of the multiple imputation process including online supplemental material providing details of the multiple imputation process and participant characteristics in original and imputed data are provided elsewhere.17 Because results were similar in the imputed and unimputed data, pooled results after the multiple imputation were presented.

The bootstrap procedure to obtain optimism-corrected goodness-of-fit statistics was performed in one randomly selected imputed dataset using R-studio. All other statistical analyses were performed using SPSS V.25.0 (IBM). A two-sided p value of 0.05 was considered statistically significant.

Results

Population description

A total of 199 participants were included, of whom approximately one-quarter rated their health fair to poor (table 1). The median (IQR) PCS and MCS scores were 49.0 (45.2; 57.6) and 48.3 (42.1; 54.6), respectively, with higher scores indicating a poorer experienced health. Approximately one-quarter of the participants experienced food insecurity. Participants had a median (IQR) age of 38.0 (33.8–43.5) years. The majority of participants were women (84.9%), had an income below the basic needs budget (64.8%), had an upper secondary educational level or more (61.3%), were married or cohabiting (69.8%), and were currently unemployed (55.8%). Compared with participants who rated their health good to excellent, participants with fair to poor health more often experienced food insecurity (42.0% vs 18.8%), more often had an income below the basic needs budget (78.0% vs 60.4%), more often were lower educated (54.0% vs 32.9%), more often were single (50.0% vs 23.5%) and less often were currently employed (32.0% vs 48.3%). They further had a slightly higher BMI (table 1).

Table 1

General health status, food insecurity status and participant characteristics for the total population and split by general health status categories

Compared with FS participants, FI participants more often reported fair to poor health, and also had a higher median (IQR) PCS score (56.2 (46.4; 66.1) vs 47.4 (45.2; 54.8)) and MCS score (54.0 (46.3; 63.6) vs 46.3 (41.3; 52.9)), indicating poorer physical and mental health (online supplemental table 1).

Variables explaining poor physical and mental health status

Crude associations with physical and mental health

The dichotomous food insecurity status was a strong individual covariate explaining both poorer PCS and MCS in the unadjusted models: FI participants had a 5.79 (95% CI 2.89 to 8.68) points higher PCS and a 7.61 (95% CI 4.67 to 10.54) points higher MCS compared with FS participants (table 2).

Table 2

Crude associations between selected covariates and the PCS and MCS

Multivariable models explaining poor physical and mental health

Adding the food insecurity score as a covariate to the model with PCS as outcome, this was the most important covariate (standardised β 0.21), followed by age (standardised β 0.16), household income (standardised β 0.14) and living situation (standardised β 0.13). With MCS as outcome, including food insecurity score as a covariate, again this was the most important covariate (standardised β 0.27), followed by employment status (standardised β 0.20) and age (standardised β 0.11) (table 3).

Table 3

Associations between selected covariates and the PCS and MCS, with and without including food insecurity status score as a covariate

The optimism-corrected R2 for the multivariable model with PCS as outcome improved from 6.3% to 9.2% when adding food insecurity score as a covariate, an improvement in explained variance of 2.9%. The optimism-corrected R2 for the multivariable model with MCS as outcome improved from 5.8% to 11.0% when food insecurity score was included as a covariate, an improvement in explained variance of 5.2%. The models including food insecurity score were a better fit compared with the models not including food insecurity score, as indicated by lower optimism-corrected RMSEs (table 3).

Discriminative performance

Including the food insecurity score as a covariate for the dichotomous PCS score improved the optimism-corrected C-statistic from 0.64 (95% CI 0.59 to 0.71) to 0.69 (95% CI 0.62 to 0.73) and Nagelkerke’s R2 from 9.6% to 14.0%, an improvement of 4.4%. Including the food insecurity score as a covariate for the dichotomous MCS score improved the C-statistic from 0.65 (95% CI 0.55 to 0.68) to 0.70 (95% CI 0.61 to 0.73) and Nagelkerke’s R2 from 5.4% to 11.0%, an improvement of 5.6% (table 4).

Table 4

Optimism-corrected C-statistic and Nagelkerke’s R2 for the multivariable models explaining dichotomous PCS and MCS scores, with and without including food insecurity status score as a covariate

Discussion

The results of our study indicate that food insecurity status was a strong covariate explaining both poorer physical and mental health in unadjusted models. In the multivariable models explaining PCS and MCS, from all included socioeconomic risk factors, the food insecurity score was the most important covariate. Including food insecurity score in those models led to an increase in explained variance of nearly one-half for PCS, an almost two-fold increase in explained variance for MCS, and a slightly better model fit. Further analyses showed that including food insecurity score improved the discriminative ability between those individuals most at risk of poor health (ie, the ability to distinguish between those having a score below 50 and those having a score above 50, which indicates poorest physical and mental health), reflected by an increased C-statistic and an improvement in explained variance for both PCS and MCS. From these results it follows that food insecurity status is of added value in explaining poor health, particularly mental health, beyond traditionally used socioeconomic and sociodemographic risk factors (ie, age, educational level, household income level, living situation, employment status and migration background). Therefore, including food insecurity status may be important for effective risk stratification to identify populations at increased risk of poor health.

In line with previous literature,11 14 our results show that experiencing food insecurity is associated with poorer physical and mental health. The differences between FS and FI participants in physical and mental health that were found in our study were well above the minimal ‘clinically important difference’ of 3–5 points proposed by Samsa et al.30 Food insecurity may be linked to poor health through multiple potential pathways such as shifting towards less expensive, lower-quality foods31 and elevated levels of depression and (chronic) stress.14 Also, impaired adherence to medical recommendations due to budgetary constraints may play a role, for example having to choose between food and medicine.32 Food insecurity is forecasted to increase due to the current COVID-19 pandemic, thereby further increasing the risk of poor health in the short term and long term through several pathways.33 For example, a recent study including over 2700 low-income Americans showed that food insecurity caused by the COVID-19 pandemic was highly associated with mental health issues.34

As described by Predmore et al,35 addressing social determinants of health within healthcare organisations contributes to achieving the Triple Aim.35 With regard to predictive risk modelling, one of their proposed applications is ‘social predictive modelling and case finding’ by incorporating social risk factors,35 as was done in our study. However, despite the large body of literature showing that incorporating social determinants of health improves the ability to identify people at risk for poor health,11 35 food insecurity status is barely used for the identification of populations at increased risk of poor health.

Elaborating on this knowledge, our results underline the importance of using food insecurity status data to identify populations at increased risk of poor health in a Dutch urban setting. Implementing this requires availability of data on food insecurity status, emphasising the urge to start routinely collecting data on food insecurity status in the Netherlands. Screening for food insecurity status has value beyond better identification of people at risk of poor health, because it also helps making healthcare providers aware of the existence of social risk factors such as food insecurity. Only when they are aware of these issues among their patients, they can address them and improve access to resources, if available.36 Multiple tools are currently available for screening for food insecurity, ranging from very short, one-item screening tools to more elaborate surveys.36 For example, short, validated screening tools are available that allow minimal additional time and costs associated with the screening, which helps to maintain acceptability for both the person being screened and the person performing the screening.37 In the Netherlands, screening among high-risk groups could be done in clinical settings such as the general practice (as most Dutch citizens regularly visit their primary care physician) and/or nonclinical settings such as community centres (as these centres are generally visited by disadvantaged people).35 Importantly, the identification of people at risk of food insecurity should ideally be followed by referral to effective interventions or resources, and options to integrate these into routine care in the Dutch context should be further explored. This may also call for referral to resources across domains, such as the social domain (ie, social prescribing), which is challenging in the current Dutch context due to different funding streams.

Our results suggest the need for screening high-risk groups for food insecurity and the development and implementation of interventions addressing food insecurity and its consequences (while incorporating the needs and preferences of this population and the healthcare provider that performs the screening). Together, these actions are expected to contribute to the Quadruple Aim by improving experienced quality of care (as underlying needs associated with food insecurity and its consequences can be addressed), reducing healthcare costs (which will follow from reduced food insecurity prevalence), improved provider experience (as also their needs and preferences are considered and they can offer better help to their patients in need), and ultimately improved population health.2 38

Our study is among the first to investigate the added value of food insecurity status in explaining poor health. Our study is strengthened by the use of validated measures of our main outcome and covariate. As a measure of poor health, we used the SF-12 which is a widely used, reliable and well-validated measure of general health,19 and strongly associated with both short-term and long-term mortality risk39 and higher healthcare use and costs.21 Previous research has indicated that the SF-12 is a suitable alternative for the more elaborate SF-36, also in the Dutch population.40

We assessed food insecurity status using the widely applied 18-item USDA-HFSSM, which is regarded as the golden standard for Western countries.41 Because being poor is not one-to-one related to experiencing food insecurity, it is important not to use indirect indicators such as income as a proxy for food insecurity status,42 as was done in the current study. Food insecurity is a complex phenomenon that encompasses many dimensions, reflecting a condition where there is unreliable (physical or economic) access to sufficient food. Food insecurity may for example include (anxiety and worries about) not having enough (healthy) foods, the inability to acquire food in socially acceptable ways, or (perceived) social exclusion because of the inability to participate in the social and cultural norms. One could argue that food insecurity interacts with adverse health outcomes, and therefore reflects a potential syndemic (ie, two or more mutually enhancing health conditions that cluster within a specific population, in light of socioecological inequality and inequity that enhances this adverse interaction43). Himmelgreen et al44 clearly describe this in their proposed dynamic model of the food insecurity-diet-related chronic diseases syndemic.44 In short, this model shows how socioecological inequality and inequity induce food insecurity and associated stress, which has an amplifying adverse effect on nutrition and health status (also depending on the life course stage), which can ultimately result in diet-related chronic disease(s). These diseases create a feedback loop that can create a vicious cycle, thereby amplifying adverse health outcomes.44 This theory helps explain the added value of food insecurity beyond traditional social determinants of health in explaining poor health, as food insecurity may also comprise this syndemic effect. It should be noted that our measure of food insecurity, based on the USDA-HFSSM, mostly focuses on economic access to food, and may still not fully capture other dimensions of food insecurity that are also important for explaining poor health. However, we found a strong association between the food insecurity status as assessed using the USDA-HFSSM and poor physical and mental health, indicating that this measure adequately captured the food insecurity dimensions important for health.

Another important consideration is that we treated food insecurity as a covariate explaining poor health and aiding risk-stratification, not as a health outcome on itself. Conceptualising health from a broader, multidimensional and positive perspective (eg, ‘positive health’), health can be seen as more than the mere absence of disease, as it also includes functioning/resilience, resources/supports and quality of life.45 From this perspective, one could argue that food insecurity is a health outcome on itself rather than a covariate explaining poor health. For treating food insecurity as an outcome, different analyses and models than the ones used in the current study would have been more appropriate. However, our approach using a social determinant such as food insecurity as a covariate for better identification of high-risk populations is better aligned with how the current Dutch healthcare system operates.

It should further be noted that, although including food insecurity in the models improved the explained variance in poor health, these models still explained only about 10% of health differences. As health is a multidimensional concept that is influenced by many factors, it is not uncommon to find a relatively low explained variance.46 This suggests that besides food insecurity, other factors such as lifestyle behaviours or chronic stress, or social factors such as social networks, are important for explaining poor health. For example, a large study among middle-aged and older adults in Norway showed that the association between SEP and health was mediated by loneliness, suggesting that this is an important factor contributing to poor health.46

Our study is strengthened by accounting for statistical optimism in our multivariate models explaining poor health. We used the same dataset to fit the models and to assess the validity of our model, whereas ideally we would have externally validated our results using a test dataset from the same population to verify your results, which was not possible in our study.27 This can lead to optimistic estimates of model performance (ie, the models built using the same dataset as the one that was used to fit the models performs better in explaining poor health than it would have if a different dataset was used). One solution to assess the model performance without having a test set is by using bootstrapping, as was done in our study.

An important methodological consideration is the use of cross-sectional data for our analyses, which is not suitable for a traditional clinical prediction models wherein a future outcome is predicted and temporality can be ensured. In addition, we assume that experiencing food insecurity precedes poor health which is plausible considering previous research, however, it is also possible that poor health leads to food insecurity (eg, through increased stress, or medical costs or job loss leading to reduced budgets for food). The issue of reverse causality cannot be ruled out using cross-sectional data. Our approach was, however, suitable for our main aim as it enabled us to show that including information on food insecurity and adding this to traditional social determinants of health seems to have value for better explaining poor health.

Further, our sample mainly included women living in a disadvantaged urban setting, and therefore the results may not be generalisable to the general Dutch population. Previous studies indicate that women are more at risk of food insecurity and its accompanying health consequences,47 but due to the small number of men in our study sample, we were unable to explore these gender differences further in the current study. Also, the sample size was relatively small, especially when compared with large-scale food insecurity screening surveys such as those annually conducted by the USDA. However, it should be noted that food insecurity is a relatively understudied area in the Netherlands, and the presented results can stimulate larger-scale, routine screening for food insecurity in the Netherlands as well. Future studies should validate our results in other populations and settings, ideally using longitudinal data to confirm the temporality assumption.

Conclusions

Food insecurity status is important for explaining poor health, particularly mental health, beyond other socioeconomic risk factors in disadvantaged communities. Our results need confirmation in other populations and settings. Food insecurity status hereto needs to be assessed in routine data collections. These data can be used to better identify people with increased risk of poor health and optimise the allocation of available resources to the people most in need.

Data availability statement

Data are available on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

The study was reviewed by the Medical Ethics Committee of Leiden University Medical Centre and confirmed not to be subject to the Medical Research Involving Human Subjects Act (WMO)(P17.164). Participants gave informed consent to participate in the study before taking part.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors JK-dJ and LAvdV designed the research project. JK-dJ and MEN supervised the overall study. LAvdV was involved in the data collection. EWS provided consultation regarding data analysis. LAvdV conducted data analyses with advice from JK-dJ and EWS. JK-dJ, EWS and MEN provided consultation regarding the interpretation of the data. LAvdV drafted the manuscript in close collaboration with JK-dJ. JK-dJ is the guarantor. All authors read, edited and approved the final version of the manuscript.

  • Funding This work was supported by the Municipally of The Hague. The Municipally of The Hague was not involved in the design of the study, collection, analysis and interpretation of data, or in writing the manuscript.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.