Add a stdev- calculation and output to the mean function