Create pdstats.py