In a secondary set of sensitivity analyses, we examined whether our pattern of results was robust to the effects of divorce and separation, college attendance, gene-by-covariate and covariate-by-environment interactions, and family ascertainment status. We observed the same pattern of results after removing individuals who reported divorce or separation from the sample and when controlling for college attendance. Likewise, we observed the same pattern of results after controlling for all gene-by-covariate and covariate-by-environment interactions. Lastly, we observed the same pattern of results when controlling for family ascertainment status. None of our effects significantly changed with the inclusion of these covariates or interaction terms, which guided our decision to report the results with the fewest parameters in this paper (full results available upon request from the first author).
