pituuksien_normaalijakauma.py

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

height = [186, 176, 158, 180, 186, 168, 168, 164, 178, 170, 189, 192, 172,
     187, 180, 186, 185, 168, 179, 178, 183, 179, 170, 175, 186, 159,
     161, 178, 175, 185, 175, 162, 173, 172, 177, 175, 172, 177, 180]
height.sort() # lajitellaan pistuusjärjestykseen
heightmean = np.mean(height) # pituuksien keskiarvo
heightmedian = np.median(height) # pituuksien mediaani
heightstd = np.std(height)   # pituuksien keskihajonta
pituuksia = str(round(len(height)))
keskipit = round(heightmean)
keskihajonta = round(heightstd)
mediaani = round(np.median(height))
pdf = stats.norm.pdf(height, heightmean, heightstd)
plt.title('Pituuksia '+pituuksia+ ', keskipituus '+str(keskipit)+ ' ,mediaani '+str(mediaani)+', keskihaj '+str(keskihajonta))
plt.xlabel('pituus')
plt.ylabel('todennäköisyys')
plt.plot(height, pdf)
plt.show()