Use MASS library to output SE of D50