Initially descriptive analysis to assess the balance between two groups will be conducted. The primary analysis will employ multilevel multivariable lineal regression to investigate differences in the PHQ-9 scores between groups at 3 and 6 months after randomization, adjusting for baseline outcome variable scores. Sensitivity analysis making different assumptions will be conducted to investigate the potential effects of missing data. Similar analysis will be done for the secondary outcomes measures.