Perform linear and quadratic discriminant function analysis with MASS package. Use lda() for linear discriminant analysis and assess prediction accuracy.