Finding significant combinations of features in the presence of categorical covariates