To characterise the specificity of the allelic effects for our main findings, we examined the association evidence in three ways: we generated forest plots of the disorder beta coefficients with 95% CIs; we calculated a heterogeneity p value for the disorder-specific effects contributing to the overall statistics for meta-analytic association; and we undertook a multinomial logistic regression procedure with model selection19 for each main SNP for all five disorders to assess the pattern of phenotypic effects (appendix pp 8–11). To compare the fit of various models of genotype–phenotype associations, we applied established goodness-of-fit metrics (the Bayesian information criteria and the Akaike information criteria). We report the best-fitting model by Bayesian criteria and show results of both metrics for a range of models (appendix pp 38–45, 51–61).