The ordinal regression models predicting the outcomes of ClinDep, suicidality, and SelfRepDep were as follows: Logit(outcome)=β0+β1(sex)+β2(age)+β3(PLE)+β4(NLE)+β4(genotype)+β5(PLE×genotype)+β6(NLE×genotype) where the β weights were in log(odds) units and sex was coded males = 0, females = 1. Analyses were run for two genotype classifications, (1) 5HTTLPR and (2) 5HTTLPR + rs25531 (see above), where both additive (SL intermediate to SS and LL) and S-dominant (SS and SL combined) models were fitted. The interaction terms were centered for all analyses.
