Considerable debate exists on whether exposure to vegetation cover is associated with better mental health outcomes. Past studies could not accurately capture people's exposure to surrounding vegetation and heavily relied on non-spatial models, where the spatial autocorrelation and latent covariates could not be adjusted. Therefore, a suite of five different vegetation measures was used to separately analyze the association between vegetation cover and the number of psychotic and non-psychotic disorder cases in the neighborhoods of Toronto, Canada. Three satellite-based and two area-based vegetation measures were used to analyze these associations using Poisson lognormal models under a Bayesian framework. Healthy vegetation cover was found to be negatively associated with both psychotic and non-psychotic disorders. Results suggest that the satellite-based indices, which can measure both the density and health of vegetation cover and are also adjusted for urban and environmental perturbations, could be better alternatives to simple ratio- and area-based measures for understanding the effect of vegetation on mental health. A strong dominance of spatially structured latent covariates was found in the models, highlighting the importance of adopting a spatial approach. This study can provide critical guidelines for selecting appropriate vegetation measures and developing spatial models for future population-based epidemiological research.
Keywords: Bayesian; enhanced vegetation index; mental health; non-psychotic; normalized difference vegetation index; psychotic; random forest; soil-adjusted vegetation index; spatial modeling; vegetation.