Add bonus_sample and stats_sample
Add multinomial.py
Add pmfs for stats and bonus
Rename x_sample to sample_x