Regularized linear discriminant analysis (LDA) was used to classify Day 7 embryo stage (gamma, delta hyperparameters were estimated to minimize the 5-fold cross-validation loss). All classes were assumed to have the same pooled variance-covariance matrix, and LDA analysis was restricted to the 54 genes with greatest gene expression (to ensure that the pooled covariance matrix of the training data set was positive definite).