Histograms   [Back] When we have data we often has a given pattern, such as fitting a normal distribution - such as with IQs and height of a person - or could be random. One of the best ways to look at these is a probability distribution with a histogram. The following shows a normal distribution and a random distribution. For a uniform random generator, we should get a flat trend across all the values, whereas the normal distribution depends on the μ and σ values. The more samples we have, the nearer it will get to the expected distribution: Options Normal- μ: Normal- σ: Samples: Bins: Determine Try an example μ=5, σ=3, Samples=100 Calc μ=5, σ=3, Samples=1,000 Calc μ=5, σ=3, Samples=10,000 Calc μ=5, σ=3, Samples=100,000 Calc μ=5, σ=3, Samples=1000, Bins=10 Calc μ=5, σ=3, Samples=1000, Bins=50 Calc μ=5, σ=3, Samples=1000, Bins=100 Calc μ=5, σ=3, Samples=1000, Bins=200 Calc μ=5, σ=10 Calc μ=5, σ=1 Calc μ=5, σ=0.5 Calc μ=5, σ=5 Calc μ=55, σ=10 Calc. This is the mark distribution we often aim for in academia. Normal distribution: \( f(x)=\frac{1}{\sqrt{2\pi}\sigma}e^{-(x-\mu)^2/(2\sigma^2)}\) eg: μ=9, σ=2 \( f(9)=\frac{1}{\sqrt{2\pi}2}e^{-(9-9)^2/(2(2)^2)}= 0.2\). This should be the value for x=9: Calc Uni: [0.5558678447286911, 0.18052094994382673, 0.037407728030086185, 0.30710321162326415, 0.4982919013247802, 0.3859152551676841, 0.7086932744952285, 0.6804377718369481, 0.7951386473836901, 0.8491753366839513] Norm: [8.10827091295326, 10.08584388901063, 10.46418117415455, 9.237056449172107, 11.7364685963184, 10.298896600242724, 9.12908982779769, 8.14759571578831, 8.129015612106528, 9.46088951634887] Tri: [0.8881170438588167, 0.5322907710085042, 0.3043741298648745, 0.3970473104682769, 0.2590049079828238, 0.6828123933368226, 0.8027360305760114, 0.3383821544711181, 0.3037509783308496, 0.6651992963610167] Log: [0.5638071787664829, 11.650262011213039, 9.746257696215823, 6.225871489368107, 9.12331009334957, 6.811514086317282, 4.883323184511458, 9.480857710263919, 4.930737483180143, 9.597009874587535] Exp: [25.64228904795903, 6.887507638805563, 6.6666652851075945, 0.3284127976291044, 12.79360959188501, 8.36030812872179, 16.203400965331703, 7.259260747954427, 0.8139905375051338, 6.6054884234201285] [View]

Source code

The following outlines the Python code used:

import matplotlib.pyplot as plt
import numpy as np
import sys
import random

file ='1111'
mu=9.0
sig=2.0
samples=10000



if (len(sys.argv)>1):
        file=str(sys.argv[1])

if (len(sys.argv)>2):
        mu=float(sys.argv[2])

if (len(sys.argv)>3):
        sig=float(sys.argv[3])

if (len(sys.argv)>4):
        samples=int(sys.argv[4])



fig,myplot = plt.subplots(4, 1,figsize=(8,8))
plt.tight_layout(w_pad=1.5, h_pad=2.0)


uniSamples = [random.random() for i in xrange(samples)]


myplot[0].hist(uniSamples, bins=100, normed=True)
myplot[0].set_title("Uniform random number generator histogram")
myplot[0].set_xlabel("x")
myplot[0].set_ylabel("Frequency of occurrence")
print "Uni: ",uniSamples[0:10]  #Take a look at the first 10


normSamples = [random.normalvariate(mu,sig) for i in xrange(samples)]

myplot[1].hist(normSamples, bins=100, normed=True)
myplot[1].set_title(r"Normal Histogram RNG $\mu = "+str(mu)+"$ and $\sigma = "+str(sig)+"$")
myplot[1].set_xlabel("x")
myplot[1].set_ylabel("Frequency of occurrence")

print "Norm: ",normSamples[0:10]  #Take a look at the first 10


triSamples = [random.triangular(0,1,0.5) for i in xrange(samples)]

myplot[2].hist(triSamples, bins=100, normed=True)
myplot[2].set_title(r"Triangular Histogram RNG")
myplot[2].set_xlabel("x")
myplot[2].set_ylabel("Frequency of occurrence")

print "Tri: ",triSamples[0:10]  #Take a look at the first 10
    

logSamples = [random.weibullvariate(mu,sig) for i in xrange(samples)]

myplot[3].hist(logSamples, bins=100, normed=True)
myplot[3].set_title(r"Weibull Histogram RNG")
myplot[3].set_xlabel("x")
myplot[3].set_ylabel("Frequency of occurrence")
   
print "Log: ",logSamples[0:10]  #Take a look at the first 10

plt.show()