summary_statistics function

Extract summary statistics