Discriminant analysis code