A method for biomarker discovery for substance use disorders wherein high dimensional data containing a plurality of variables based on a sample size and a number of variables exceeding that sample size are applied to an ensemble of statistical learning models whereby biomarkers of a substance use disorder are identified.