Skip to main content

Non-response bias in the analysis of the association between mental health and the urban environment: a cross-sectional study in Brussels, Belgium

Abstract

Background

This paper aims at analysing the impact of partial non-response in the association between urban environment and mental health in Brussels. The potential threats of the partial non-response are biases in survey estimates and statistics. The effect of non-response on statistical associations is often overlooked and evidence in the research literature is lacking.

Methods

Data from the Belgian Health Interview Survey 2008 and 2013 were used. The association between non-response and potential determinants was explored through logistic regressions.

Results

Participants with low income, low educational levels, lower or higher age or in households with children were less likely to respond. When adjusting for socio-economic variables, non-response was higher in areas which are less vegetated, more polluted or more urbanised. Because the determinants of non-response and depressive disorders were similar, it is reasonable to assume that there will be more people with mental health problems among the non-respondents. And because more non-responses were found in low vegetation areas, the protective association between green spaces and mental health may be underestimated.

Conclusion

Our capacity to measure the association between the urban environment and health is affected by non-response in surveys. The non-random spatial and socio-economic distribution of this bias affects the research findings.

Peer Review reports

Background

The issue of non-response and incomplete data, that often confronts population studies such as health surveys, is a well-known challenge in epidemiology [1, 2] that has already been extensively studied [3,4,5,6,7,8,9]. Non-response can occur at different stages of the survey and varies from the refusal to participate (initial non-response) to a non-response to one of the questions (partial non-response). While the determinants of initial non-response are difficult to study (via auxiliary data as in Ekholm et al., 2010 or via non-responder questionnaire as in Jaehn et al., 2020), partial non-response has the advantage that it can be studied on the basis of participants' partial responses [3, 4, 9].

Personal characteristics, such as economic status, are important in the variation of non-response [5,6,7,8]. In the case of the Belgian Health Interview Survey (BHIS), determinants of partial non-response have already been evaluated by Berete et al. [3]: non-response is more frequent among youngsters, non-Belgians, low educational levels, lower income, residents of Brussels and Wallonia, and people with poor perceived health. Partial non-response is also strongly associated with the interviewer: it explains almost half (46.6%) of the variability in partial non-response. Other studies [4, 9] on partial non-response have additionally shown that non-response is greater for men and for unskilled workers. The potential threats of this partial non-response are biases in survey estimates and statistics. The effect of non-response on statistical associations is often overlooked and evidence in the research literature is lacking [2].

Understanding the association between the urban environment and mental health is a challenging exercise. Cities are complex systems, with urban mental health outcomes affected by many interactions [10]. This topic has been attracting a lot of interest and the literature has grown substantially over the past few years [11,12,13,14,15,16,17,18]. A study conducted in Brussels concluded that exposure to traffic-related air pollution (black carbon, NO2, PM10) was positively associated with higher odds of depressive disorders but no association between mental health and green surroundings, noise, or building morphology could be demonstrated [19]. Concerning vegetation, these results differ from many studies that describe how higher levels of exposure to greenness can reduce levels of depressive symptoms [20]. This raises the question of why there is no association in the case of Brussels.

Personal characteristics, such as economic status, are also important in mental health while the urban environment appears to be a less important determinant [21]. Moreover, in Brussels, as in other cities, socio-spatial inequalities complicate matters: environmental factors are strongly correlated with socio-economic ones [22]. It is therefore difficult to disentangle the variance due to socio-economic status from the variance due to the urban environment.

If personal characteristics play an important role in explaining non-response and mental health status, we can question the magnitude of bias generated by non-response in the analysis of the association between urban environment and mental health. This paper therefore aims at analysing, characterising and discussing the impacts of non-responses in this association with the example of Brussels.

Methods

Study area

Brussels is the capital of Belgium, centrally located in the country. The city extends into three administrative regions: the Brussels-Capital Region (BCR), Flanders and Wallonia [23].We focus here on the urban area encompassed within the administrative border of the Brussels-Capital Region (BCR) because the Region was more intensively sampled compared to the rest of Belgium and because of data consistency issues for the urban environment indicators. The BCR is divided into 19 municipalities. It covers 161.38 km2 and has 1.2 million inhabitants (source: Statistics Belgium – 01/01/2018), which means an average density of 7438 inhabitants per square kilometre.

Study design

The study presented here consists of a secondary analysis of a population-based, cross-sectional study of the association between urban environment and mental health in Brussels, Belgium. Results from the original research could not demonstrate any association between urban environment (green surroundings, noise, building morphology) and mental health except for traffic-related air pollution (black carbon, NO2, PM10) exposure which was positively associated with higher odds of depressive disorders [19].

Study population and data

The Belgian Health Interview Surveys (BHIS) of 2008 and 2013 were used. We kept participants aged over 15 (minimum age to fill the self-administered questionnaire) and living at the place of residence for at least one year (n = 4,355).

In the BHIS, three types of non-response can be highlighted: (i) Initial non-response: refusal to participate in the survey [24], not addressed in this paper; (ii) Self-administered questionnaire (SAQ) non-response: for those over 15 years of age, the survey consists in two parts: a face-to-face (F2F) interview and a SAQ, the latter of which is not filled out; and (iii) Item non-response: at various points of the survey some questions were not answered by the participants. These last two types can be designated as partial non-response. Questions related to mental health are particularly sensitive. They are generally not asked in a F2F interview to avoid desirability bias and are included in a SAQ, as is the case for the BHIS.

Mental health and non-response definition

The response to six indicators related to mental health was extracted from the BHIS [25]: Energy level (SF-36 ‘vitality scale’) [based on 4 items in SAQ] [26], General Health Questionnaire (GHQ-12) [based on 13 items in SAQ] [27], Self-reported depression during the last year [based on 1 item in F2F interview], Depressive disorder [based on 13 items in SAQ], Anxiety disorder [based on 10 items in SAQ] and Sleeping disorder [based on 3 items in SAQ] (SCL-90-R subscales) [28].

Non-response to a mental health indicator means that there are not enough items answered to compute the indicator. The reason may be either that the SAQ was not completed or that some or all of the items related to an indicator were not answered. If at least 75% of the items related to an indicator have valid responses, the missing item(s) are substituted by the mean value of the valid responses. For example, if 1 to 3 items out of the 12 items of the GHQ scale have missing values, the GHQ is still computed; whereas if more than 3 items are missing, the indicator score is considered missing [29].

Urban environment

The urban environment was measured by eight indicators computed at the residence address of the BHIS participants. Urban greenness was assessed at three different geographical levels (residence, street, neighbourhood). At residence level, view of green was determined by the ratio of the total green area from Google Street View panorama (~ picture 360°) to the total area of the panorama. At street level, Linear tree density was defined as the ratio between the number of trees (point data from Urbis) [30] and the street length and Street visible vegetation coverage is the vegetation coverage (from Brussels Environment) [31] on the street, and 10 m on either side. At neighbourhood level, Vegetation coverage was defined as the ratio of vegetation coverage (from Brussels Environment) in a 1000 m circle around the respondent’s residence. The built-up environment was assessed by two indicators at the street level: the Street canyon effect which is the ratio between average building height and average open space width and the Street corridor effect which is the ratio between parallel facades length and street length. Noise was extracted from noise from multiple sources map (day – evening – night noise level, Lden) from Bruxelles Environnement for the years 2006 and 2011 at residence level. Black carbon exposition was defined as the annual average in the year of BHIS participation at residence level [32,33,34]. See Pelgrims et al. [19] for more detailed descriptions of the indicators. These indicators were missing when the health survey could not be coupled with the participant’s address (in Brussels). To address this, a ‘no data’ category was added to each indicator with missing value.

Socio-economic status

Five socio-economic indicators were used: Reported household income, Age, Gender, Family composition and Highest educational level in the household. These indicators were extracted from the BHIS. Non-responses were observed if the item has not been answered. To address this, a ‘no answer’ category was added to each indicator with non-responses.

Statistical analyses

The frequency of non-response was summarised for each mental health and socio-economic indicator. Descriptive statistics were calculated to account for the sampling strategy. ‘The weight for each sampled individual in the BHIS is the product of the reciprocal of the selection probability within a household and of a post stratification factor for each province according to age, gender, household size and quarter of the year in which the interview was done’ (Scientific Institute of Public Health, 2013).

In order to assess the association of socio-economic and urban environmental with non-response, logistic regressions with the non-response as dependent variable were computed: no missing data (0) vs. at least one non-response among six mental health indicators (1). Taking non-response to at least one mental health outcome as a proxy for non-response is a somewhat crude generalisation. However, it reflects a practice: if we want to analyse a complete case, i.e. take only those individuals who responded to all the questions we are interested in, we will then remove all individuals who did not respond to at least one question. Our purpose here was first to identify the characteristics of those excluded from the complete case by running models for socio-economic indicators. Secondly, we aimed at analysing whether these individuals were concentrated in specific environments by running single-exposure models for urban environment indicators (models A) and adjusted models (with gender, age, reported household income and year of the BHIS) for urban environment indicators (models B).

In order to assess the differences and similarities between the socio-economic characteristics of non-respondents and participants with depressive disorders, two groups of logistic regressions were performed for socio-economic variables: (i) with depressive disorders as dependent variable (models C) and (ii) with the non-response to depressive disorders related questions as dependent variable (models D).

Correct estimates and valid inferences of odds ratio (OR) were obtained by taking into account the survey weights, strata and clusters relative to the sample design. All analyses were performed using the statistical software R [35] using the survey package [36].

Results

Data description

A total of 4355 residents of the BCR were included in the study population. Table 1 describes all considered variables.

Table 1 Characteristics of the study population (weighted percentages and weighted mean)

Overview of non-responses among the study population is displayed in Table 2. The number of individuals refers to the number of respondents in the sample while the percentage is weighted to represent the study population. In our sample of 4,355 individuals, the SAQ was not available for 1,602 participants (35.48% of the population). For 18.97%, SAQ was required but not available and for 16.51%, SAQ was not required and not available (when the interview is done by a proxy such as a parent). Non-response to mental health items (excluding SAQ non-response) ranged from 1.03% (Self-reported depression) to 7.33% (Energy level). In the subset of participants for whom the SAQ was available and who answered all mental health and socio-economic items, data was available for 1,929 individuals (45.92% of the initial population).

Table 2 Non-response in the study population. No answer (NA) is due to an item non-response or to the self-administered questionnaire (SAQ) not filled out (rows in grey). Complete cases are individuals with all socio-economic (SE) and/or mental health (MH) indicators answered. Weighted percentages

Determinants of partial non-response

Figure 1 shows the association between socio-economic variables and non-responses through univariate models and fully adjusted model. These models are also displayed in Additional File 1. Considering that the coefficients between the univariate models and the fully adjusted model are very similar, only the coefficients of the univariate models will be discussed below. Participants in the lower socio-economic quartiles (Quartile 1 and Quartile 2) were less likely to respond than those in the higher quartile (Quartile 4) (p < 0.001, OR = 2.19, 95% CI = 1.66–2.90 and p < 0.001, OR = 2.15, 95% CI = 1.63–2.82 respectively). Compared to 25–44 yr, older (65 + yr) and younger (15–24 yr) participants were less likely to respond (p < 0.001, OR = 1.66, 95% CI = 1.36–2.02 and p < 0.001, OR = 1.9, 95% CI = 1.47–2.46 respectively). Compared to single participants, couples with child (ren) and one parent with child (ren) were less likely to respond (p < 0.001, OR = 1.61, 95% CI = 1.3–2 and p = 0.008, OR = 1.45, 95%CI = 1.1–1.91 respectively). Finally, participants from households with a lower educational level (higher secondary, lower secondary and no diploma or primary education) were less likely to respond than those with a higher education (after secondary) and the effect size increases as the level of education decrease (p = 0,004, OR = 1.36, 95% CI = 1.1–1.69 and p < 0.001, OR = 1.85, 95% CI = 1.43–2.39 and p < 0,001, OR = 3.6, 95% CI = 2.7–4.79 respectively).

Fig. 1
figure 1

Association between non-response and socio-economic indicators (univariate regression models and fully adjusted regression model). Levels of significance: ‘***’ = p < 0.0001, ‘**’ = p < 0.001, ‘*’ = p < 0.01, ‘·’ = p < 0.05

Figure 2 shows the association between urban environment indicators and non-response through univariate models (Models A) and adjusted models (Models B) for gender, age, reported household income and year of the BHIS. These models are also displayed in Additional File 2. When adjusted for socio-economic variables, non-response was higher in areas with lower amounts of vegetation cover (vegetation at 1 km buffer, p = 0.001, OR = 1.44, 95% CI = 1.16–1.78 and street visible vegetation, p = 0.022, OR = 1.29, 95% CI = 1.04–1.59), lower in less polluted areas (black carbon, p = 0.003, OR = 0.71, 95% CI = 0.56–0.89) and in less urbanised areas (street canyon, p = 0.04, OR = 0.8, 95% CI = 0.65–0.99 and street corridor effect, p = 0.044, OR = 0.8, 95% CI = 0.65–0.99).

Fig. 2
figure 2

Association between non-response and urban environment indicators. Models A are univariate regression models. Models B are adjusted regression models for gender, age, reported household income and year of the BHIS. Levels of significance: ‘***’ = p < 0.0001, ‘**’ = p < 0.001, ‘*’ = p < 0.01, ‘·’ = p < 0.05

Comparison of partial non-response determinant and depressive disorders determinants

Figure 3 shows the associations between socio-economic variables and non-response to depressive disorders related questions (univariate regression models C and fully adjusted regression model E) and the associations between socio-economic variables and depressive disorders (univariate regression models D and fully adjusted regression model F). These models are also displayed in Additional File 3 and 4. It should be noted that the two populations studied are different: population of models D and F is limited to the respondents to depressive disorders related questions only. The comparison of the two types of models nevertheless allows us to have an overview of the similarities and differences even if it cannot be compared on a more formal basis. The determinants of non-response are quite similar to the determinants of depressive disorders. Odds ratios are similar for household income, age (except for 15–24 y) and education level, while there are substantial differences for family composition and gender.

Fig. 3
figure 3

Association between non-response to depressive disorders related questions and socio-economic indicators (univariate regression models C and fully adjusted regression model E) and between depressive disorders and socio-economic indicators (univariate regression models D and fully adjusted regression model F). Levels of significance: ‘***’ = p < 0.0001, ‘**’ = p < 0.001, ‘*’ = p < 0.01, ‘·’ = p < 0.05

Discussion

Our results confirm that personal characteristics play a major role in explaining both non-response and mental health status. The socio-economic variables associated with non-response found in this paper (Fig. 1) are similar to those found in the literature: non-response is more frequent among young [3, 6], low educated [3, 4, 7, 9] and low income [3] inhabitants. Family composition or type of household is not a common studied factor in the literature; Ekholm et al. [6] and Jaehn et al. [7] found that married people are more likely to respond to the survey than unmarried people but it is difficult to link this finding to our results because marital status has not been taken into account as such and being in a couple is mixed with having children. A higher non-response rate for men than for women is often found in the literature [4, 6, 9] but was not found in our study.

Based on the comparison between the personal characteristics associated with non-response and depressive disorders in the results (Fig. 3), it is reasonable to assume that there will be more people with mental health problems within the group of non-respondents compared to the group of respondents. This is coherent with the literature that examined poor subjective health [7, 9], chronic condition [3] or long-term health condition [4] as a determinant for non-response.

With these results, we can extrapolate the effect of non-response to the association between mental health and the urban environment. Considering the association between vegetation and mental health, several studies on the topic have shown that vegetation is a protective factor for mental health [37, 38]. However, in a previous study using the same data as this paper [19], no association was found. With our results, it can be assumed that there would be proportionally higher number of respondents with mental disorders and low vegetation cover (\({A}_{true}\)) than actually observed due to non-responses \(({A}_{observed})\) because both low vegetation cover and mental health problems are associated with a higher probability for non-responses. Indeed, results show that the personal characteristics associated with mental health have strong similarities with those of non-response and significantly more non-responses are found in places with lower vegetation cover.

Considering the equation of an odds ratio (\(OR\)) with \(A\) as the number of respondents with mental disorder (MD +) and exposed to low green space cover (LGS +), \(B\) as the number of respondents with no mental disorder (MD-) and exposed to low vegetation cover (LGS +), \(C\) as the number of respondents with mental disorder (MD +) and high vegetation cover (LGS-) and \(D\) as the number of respondents with no mental disorder (MD-) and high vegetation cover (LGS-) (Table 3):

Table 3 Contingency table of disease vs. exposure to calculate the odds ratio
$$OR= \frac{A* D}{B* C}$$

When this OR is greater than 1, it indicates that the odds of exposure (to low green space) among case-patients (mental disorder, MD +) are greater than the odds of exposure among controls. The exposure (to low green space cover) is interpreted as a risk factor for the disease or condition (mental disorder).

If Atrue > Aobserved we may conclude that ORtrue > ORobserved. Therefore, the risk of having a mental disorder associated with exposure to low green space cover is underestimated (ORobserved < ORtrue), i.e. the odds ratio observed is smaller than the true odds ratio.

Despite the strength of the arguments presented and the valuable insights gained from this study, we have to acknowledge missing data cannot be recovered. Without complete data on all individuals, it is challenging to ascertain associations and their implications accurately. Future prospective studies can provide a more comprehensive understanding of the relationship between personal characteristics, non-response, mental health, and the urban environment.

Conclusions

In conclusion, we here show that the risk of low green space cover (or, conversely, the protective effect of high green space cover) is underestimated using the BHIS data. Indeed, the capacity to measure the association between urban environment and health is affected by non-response in surveys. The non-random spatial and socio-economic distribution of this bias affects the parameter estimates in the statistical models. In a sense, non-responses may be seen as a cause of a loss of statistical power, leading to a lower probability to statistically detect differences if such differences would exist. This may explain why when studying environment and human health, results are sometimes inconclusive or inconsistent with the literature. The environment and human health association is particularly challenging to study as we are dealing with multifactorial elements. More attention should be given to understand the biases that can be encountered due to non-responses and to the non-random character of this bias. A combination of methods or a mixed-method approach (quantitative and qualitative) can help to obtain a comprehensive understanding of these highly complex topics.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available due to privacy.

References

  1. Johnson TP, Wislar JS. Response rates and nonresponse errors in surveys. JAMA. 2012;307(17):1805–6.

    Article  CAS  PubMed  Google Scholar 

  2. Peytchev A. Consequences of survey nonresponse. Ann Am Acad Pol Soc Sci. 2013;645(1):88–111.

    Article  Google Scholar 

  3. Berete F, der Heyden JV, Demarest S, Charafeddine R, Gisle L, Braekman E, et al. Determinants of unit nonresponse in multi-mode data collection: a multilevel analysis. PLoS ONE. 2019;14(4): e0215652.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Boshuizen HC, Viet AL, Picavet HSJ, Botterweck A, van Loon AJM. Non-response in a survey of cardiovascular risk factors in the Dutch population: Determinants and resulting biases. Public Health. 2006;120(4):297–308.

    Article  CAS  PubMed  Google Scholar 

  5. Burzykowski T, Molenberghs G, Tafforeau J, Van Oyen H, Demarest S, Bellamammer L. Missing data in the health interview survey 1997 in Belgium. Archives of Public Health. 1999;57(3):107–29.

    Google Scholar 

  6. Ekholm O, Gundgaard J, Rasmussen NKR, Hansen EH. The effect of health, socio-economic position, and mode of data collection on non-response in health interview surveys. Scand J Public Health. 2010;38(7):699–706.

    Article  PubMed  Google Scholar 

  7. Jaehn P, Mena E, Merz S, Hoffmann R, Gößwald A, Rommel A, et al. Non-response in a national health survey in Germany: An intersectionality-informed multilevel analysis of individual heterogeneity and discriminatory accuracy. PLoS ONE. 2020;15(8): e0237349.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Korkeila K, Suominen S, Ahvenainen J, Ojanlatva A, Rautava P, Helenius H, et al. Non-response and related factors in a nation-wide health survey. Eur J Epidemiol. 2001;17(11):991–9.

    Article  CAS  PubMed  Google Scholar 

  9. Volken T. Second-stage non-response in the Swiss health survey: determinants and bias in outcomes. BMC Public Health. 2013;13(1):167.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Evans GW. The built environment and mental health. J Urban Health. 2003;80(4):536–55.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Chi D, Aerts R, Van Nieuwenhuyse A, Bauwelinck M, Demoury C, Plusquin M, et al. Residential exposure to urban trees and medication sales for mood disorders and cardiovascular disease in Brussels, Belgium: an ecological study. Environ Health Perspect. 2022;130(5): 057003.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Gong Y, Palmer S, Gallacher J, Marsden T, Fone D. A systematic review of the relationship between objective measurements of the urban environment and psychological distress. Environ Int. 2016;1(96):48–57.

    Article  Google Scholar 

  13. Gruebner OA, Rapp M, Adli M, Kluge U, Galea S, Heinz A. Cities and mental health. Dtsch Arztebl Int. 2017;114(8):121–7.

    PubMed  PubMed Central  Google Scholar 

  14. Hankey S, Marshall JD. Urban form, air pollution, and health. Curr Envir Health Rpt. 2017;4(4):491–503.

    Article  CAS  Google Scholar 

  15. Okkels N, Kristiansen CB, Munk-Jørgensen P, Sartorius N. Urban mental health: challenges and perspectives. Curr Opin Psychiatry. 2018;31(3):258–64.

    Article  PubMed  Google Scholar 

  16. Rautio N, Filatova S, Lehtiniemi H, Miettunen J. Living environment and its relationship to depressive mood: A systematic review: International Journal of Social Psychiatry. 2017; Available from: https://journals.sagepub.com/doi/https://doi.org/10.1177/0020764017744582.Cited 2020 Jul 3.

  17. Ventriglio A, Torales J, Castaldelli-Maia JM, Berardis DD, Bhugra D. Urbanization and emerging mental health issues. CNS Spectr. 2021;26(1):43–50.

    Article  PubMed  Google Scholar 

  18. van der Wal JM, van Borkulo CD, Deserno MK, Breedvelt JJF, Lees M, Lokman JC, et al. Advancing urban mental health research: from complexity science to actionable targets for intervention. The Lancet Psychiatry. 2021;8(11):991–1000.

    Article  PubMed  Google Scholar 

  19. Pelgrims I, Devleesschauwer B, Guyot M, Keune H, Nawrot TS, Remmen R, et al. Association between urban environment and mental health in Brussels, Belgium. BMC Public Health. 2021;21(1):635.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Fong KC, Hart JE, James P. A review of epidemiologic studies on greenness and health: updated literature through 2017. Curr Environ Health Rep. 2018;5(1):77–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kabisch N. The Influence of Socio-economic and Socio-demographic Factors in the Association Between Urban Green Space and Health. In: Marselle MR, Stadler J, Korn H, Irvine KN, Bonn A, editors. Biodiversity and Health in the Face of Climate Change. Cham: Springer International Publishing; 2019 91–119. Available from: http://link.springer.com/https://doi.org/10.1007/978-3-030-02318-8_5.Cited 2019 Sep 17.

  22. Stessens P, Khan AZ, Huysmans M, Canters F. Analysing urban green space accessibility and quality: A GIS-based model as spatial decision support for urban ecosystem services in Brussels. Ecosystem Services. 2017; Available from: http://www.sciencedirect.com/science/article/pii/S2212041616303898.Cited 2017 Nov 27.

  23. Thomas I, Cotteels C, Jones J, Peeters D. Revisiting the extension of the Brussels urban agglomeration : new methods, new data … new results ? Belgeo Revue belge de géographie. 2012;(1–2). Available from: http://belgeo.revues.org/6074

  24. Van der Heyden J, Demarest S, Van Herck K, De Bacquer D, Tafforeau J, Van Oyen H. Association between variables used in the field substitution and post-stratification adjustment in the Belgian health interview survey and non-response. Int J Public Health. 2014;59(1):197–206.

    Article  PubMed  Google Scholar 

  25. Gisle L. Santé mentale. In: Van der Heyden J, Charafeddine R, editors. Enquête de santé 2013 Rapport 1 : Santé et Bien-être. Bruxelles,: WIV-ISP; 2014. p. 778–935. Available from: https://www.sciensano.be/fr/biblio/enquete-de-sante-2013-rapport-1-etat-de-sante-et-bien-etre

  26. Ware JE, Sherbourne CD. The MOS 36-item short-form health survey (SF-36). I. Conceptual framework and item selection. Med Care. 1992;30(6):473–83.

    Article  PubMed  Google Scholar 

  27. Goldberg DP, Williams P. A user’s guide to the General Health Questionnaire. Berkshire: NFER-NELSON; 1988. p. 129.

  28. Derogatis LR. SCL-90-R administration, scoring and procedures manual. Minneapolis: National Computer Systems; 1994. p. 123.

  29. Scientific Institute of Public Health. Manual for the users of the Health Interview Survey 2013 database. Brussels, Belgium: Scientific Institute of Public Health; 2013.

    Google Scholar 

  30. BRIC. UrbIS. Brussels Regional Informatics Centre (BRIC); 2017. Available from: https://datastore.brussels/web/urbis-download.Cited 2020 May 18.

  31. Van de Voorde T, Canters F, Chan JCW. Mapping update and analysis of the evolution of nonbuilt (green) spaces in the Brussels Capital Region. Bruxelles Environnement / Leefmilieu Brussel; 2010. Available from: http://document.environnement.brussels/opac_css/elecfile/Study_NonBuildSpaces_I_II_en.PDF.Cited 2019 Jul 19.

  32. Janssen S, Dumont G, Fierens F, Mensink C. Spatial interpolation of air pollution measurements using CORINE land cover data. Atmos Environ. 2008;42(20):4884–903.

    Article  CAS  Google Scholar 

  33. Lefebvre W, Vranckx S. Validation of the IFDM-model for use in urban applications. Study Accomplished in the Framework of the ATMOSYS-project, VITO report 2013/RMA/R/56. 2013. p. 208.

  34. Lefebvre W, Van Poppel M, Maiheu B, Janssen S, Dons E. Evaluation of the RIO-IFDM-street canyon model chain. Atmos Environ. 2013;1(77):325–37.

    Article  Google Scholar 

  35. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria; 2020. Available from: https://www.R-project.org/

  36. Lumley T. Analysis of complex survey samples. J Stat Softw. 2004;15(9):1–19.

    Google Scholar 

  37. Gascon M, Sánchez-Benavides G, Dadvand P, Martínez D, Gramunt N, Gotsens X, et al. Long-term exposure to residential green and blue spaces and anxiety and depression in adults: a cross-sectional study. Environ Res. 2018;162:231–9.

    Article  CAS  PubMed  Google Scholar 

  38. Lee ACK, Maheswaran R. The health benefits of urban green spaces: a review of the evidence. J Public Health (Oxf). 2011;33(2):212–22.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This project was partially funded by the ‘Nature Impact on Mental Health Distribution (NAMED)’ project (BR/175/A3/NAMED), a BRAIN-be (Belgian Research Action through Interdisciplinary Networks) project funded by BELSPO (Belgian Science Policy).

Author information

Authors and Affiliations

Authors

Contributions

MG, IP, EDC & SV contributed to conception and design of the study. MG developed the statistical analysis. IP performed the statistical analysis. EDC, IT & SV supervised. RA helped interpret the findings. MG wrote the first draft of the manuscript. HK coordinated the NAMED project in which this paper emerged. MG, IP, RA, HK, RR, EDC & IT were part of the NAMED team. All authors contributed to manuscript revision, read, and approved the submitted version.

Corresponding author

Correspondence to Sophie O. Vanwambeke.

Ethics declarations

Ethics approval and consent to participate

The data linkage, relying on the secondary use of the Belgian BHIS data, administrative, GIS and monitoring data did not require any approval from an Ethics committee. An application to the Sector Committee Statistics has been submitted and approved (see Decision STAT n°02/2018 from the 19/01/2018).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Association between non-response and socio-economic indicators.

Additional file 2.

Association between non-response and urban environment indicators. Models A are univariate regression models. Models B are adjusted regression models for gender, age, reported household income and year of the BHIS.

Additional file 3.

Association between non-response to depressive disorders related questions and socio-economic indicators (univariate regression models C) and between depressive disorders and socio-economic indicators (univariate regression models D).

Additional file 4.

Association between non-response to depressive disorders related questions and socio-economic indicators (fully adjusted regression model E) and between depressive disorders and socio-economic indicators (fully adjusted regression model F).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guyot, M., Pelgrims, I., Aerts, R. et al. Non-response bias in the analysis of the association between mental health and the urban environment: a cross-sectional study in Brussels, Belgium. Arch Public Health 81, 129 (2023). https://doi.org/10.1186/s13690-023-01118-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13690-023-01118-y

Keywords