Centrální limitní věta

Úvod do statistiky v Pythonu

Maggie Matsui

Content Developer, DataCamp

Hod kostkou 5krát

die = pd.Series([1, 2, 3, 4, 5, 6])

# Roll 5 times
samp_5 = die.sample(5, replace=True)
print(samp_5)

array([3, 1, 4, 1, 1])

np.mean(samp_5)

2.0

šestistěnná kostka

# Roll 5 times and take mean
samp_5 = die.sample(5, replace=True)
np.mean(samp_5)

4.4

samp_5 = die.sample(5, replace=True)
np.mean(samp_5)

3.8

Opakovat 10krát:

sample_means = []

for i in range(10):

  samp_5 = die.sample(5, replace=True)
  sample_means.append(np.mean(samp_5))

print(sample_means)

[3.8, 4.0, 3.8, 3.6, 3.2, 4.8, 2.6,
3.0, 2.6, 2.0]

Výběrové rozdělení výběrového průměru

histogram 10 výběrových průměrů

sample_means = []
for i in range(100):
    sample_means.append(np.mean(die.sample(5, replace=True)))

histogram 100 výběrových průměrů

sample_means = []
for i in range(1000):
    sample_means.append(np.mean(die.sample(5, replace=True)))

histogram 1000 výběrových průměrů

S rostoucím počtem pokusů se výběrové rozdělení statistiky přibližuje normálnímu rozdělení.

histogramy 10, 100 a 1000 výběrových průměrů – vyšší počet průměrů má tvar zvonové křivky

* Výběry by měly být náhodné a nezávislé

sample_sds = []
for i in range(1000):
  sample_sds.append(np.std(die.sample(5, replace=True)))

Rozdělení 1000 výběrových směrodatných odchylek z 5 hodů kostkou

sales_team = pd.Series(["Amir", "Brian", "Claire", "Damian"])

sales_team.sample(10, replace=True)

array(['Claire', 'Damian', 'Brian', 'Damian', 'Damian', 'Amir', 'Amir', 'Amir', 
      'Amir', 'Damian'], dtype=object)

sales_team.sample(10, replace=True)

array(['Brian', 'Amir', 'Brian', 'Claire', 'Brian', 'Damian', 'Claire', 'Brian', 
      'Claire', 'Claire'], dtype=object)

Rozdělení výběrových proporcí také vypadá normálně

# Estimate expected value of die
np.mean(sample_means)

3.48

# Estimate proportion of "Claire"s
np.mean(sample_props)

0.26

Výběrové rozdělení výběrových průměrů s přerušovanou čarou uprostřed

Úvod do statistiky v Pythonu