summaryStats function

Calculate summary statistics