Statistical analyses were performed in IBM SPSS Statistics 26 (SPSS, Chicago, IL, USA). The primary outcome measure was the dichotomous variable for having clinically significant depressive symptoms. Descriptive statistics were used to provide general characteristics of all participants and were used to look for differences in the experienced depressive symptoms between subgroups. After descriptive statistics were performed, missing data were handled using a multiple imputation technique in 10 different data sets.
52 This was done by using the Markov chain Monte Carlo method,
53 providing pooled estimates for all analyses. After assumptions were checked,
54 multivariable logistic regression models were used to investigate the associations of interest. The independent variable in the first analysis was season, which was divided into four categories (spring, summer, fall, and winter). Dummy variables were added to the model. The independent variable in the second analysis was the hours of sunlight, which was a continuous variable. The independent variable in the third analysis was sensitivity to bright light, which was dichotomous. All three analyses were performed unadjusted and adjusted for age, sex, level of education, having paid work, visual acuity, comorbidity, and household situation. These potential confounders were selected based on availability in the various data sets and the possibility that these variables could be unevenly distributed among the investigated independent variables and may possibly have an effect on depression in this population. All three analyses were expanded by exploring the role of the possible effect modifiers (i.e., severity of vision loss and eye condition). In all three analyses, interactions of people who were blind versus people who had low vision and people with different eye diseases were added to the adjusted models. Additionally, effect modification of sensitivity to bright light was explored in the first and second analyses (season and sunlight) by adding interactions to the adjusted models. During all analyses, a
P value of 0.05 (two-sided) was considered statistically significant.