prob.py

import numpy as np 
import scipy.stats as stats
import collections
import matplotlib.pyplot as plt


x = [1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 4, 4, 4, 4, 5, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 9, 9]

c = collections.Counter(x)
count_sum = sum(c.values())

for k,v in c.items():
	print("The frequency of number " + str(k) + " is " + str(float(v) / count_sum))

plt.boxplot(x)
plt.savefig("boxplot.png")

plt.hist(x, histtype='bar')
plt.savefig("histogram.png")

plt.figure()
data = x
graph = stats.probplot(x, dist="norm", plot=plt)
plt.savefig("qqplot.png")