Create sumstats.py

Some initial summary statistics