In this study, we differentiated homogeneous, highly heritable clinical subtypes of CD to identify subtype-specific genetic variants. To accomplish this, we designed a multivariate clustering method using clinical features and symptoms of CD to assign cocaine-using subjects to subgroups. Using traits derived from the CD subgroups, we examined a panel of genetic variants for association, including both main effects and pairwise epistatic interactions.