Genomic structural equation modeling (GenomicSEM) was performed in R using the GenomicSEM package84. Multivariable linkage disequilibrium matrices were created using the 1000 Genomes Project Phase 3 European reference. Exploratory factor analysis (EFA) was used to estimate the most appropriate number of latent factors represented by the psychiatric phenotypes and psychopathologies tested assuming a maximum number of latent factors equal to Ntraits-1. Confirmatory factor analysis (CFA) was used to calculate factor loadings onto each latent factor(s). Standardized loading values are reported.
