PRS were generated using the PRS-cs software package70 consistent with our other (i.e., Yale-Penn 4, BioVU) PRS analyses described above. Associations between addiction-rf PRS and phenotypes were estimated using mixed-effects models in the lme475 package in R. PRS were scaled to unit variance. Family ID and site were included as random effects to account for non-independence of measurement associated with relatedness and scanner/site. We controlled for the first 10 ancestry principal components, age, sex, age by sex. We used a Bonferroni-corrected phenome-wide significance threshold of 0.05/1480= 3.38e-05; all results are presented in the Supplemental Table 21.
