Because the mean score of depressive symptoms from the follow-up
survey displayed a positively skewed distribution, we constructed three
binary variables. The first binary variable indicated whether an adoles-
cent scored above the median of depressive symptoms or not