Aan de slagBegin gratis

Variatie in steekproefstatistieken

Als we één steekproef met size=1000 maken door zoveel punten uit een populatie te trekken, en vervolgens een steekproefstatistiek berekenen, zoals het gemiddelde, dan krijgen we één waarde die de steekproef samenvat.

Als je dat steekproefproces num_samples=100 keer herhaalt, krijg je 100 steekproeven. Als je voor elke verschillende steekproef de steekproefstatistiek, zoals het gemiddelde, berekent, levert dat een verdeling van waarden van het gemiddelde op. Het doel is dan om het gemiddelde van de gemiddelden en de standaardafwijking van de gemiddelden te berekenen.

Hier gebruik je de vooraf ingeladen population, num_samples en num_pts, en merk op dat de arrays means en deviations zijn geïnitialiseerd op nul om je containers te geven voor gebruik in de for-lus.

Deze oefening maakt deel uit van de cursus

Introductie tot lineaire modellering in Python

Bekijk cursus

Oefeninstructies

  • Genereer voor elk van num_samples=100 een steekproef en bereken en sla de steekproefstatistieken op.
  • Maak per iteratie een sample door met np.random.choice() 1000 willekeurige punten uit de populatie te trekken.
  • Bereken en sla per iteratie de methoden sample.mean() en sample.std() op om het gemiddelde en de standaardafwijking van de steekproef te bepalen.
  • Bereken voor de array means en de array deviations zowel het gemiddelde als de standaardafwijking, en print de resultaten.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

# Initialize two arrays of zeros to be used as containers
means = np.zeros(num_samples)
stdevs = np.zeros(num_samples)

# For each iteration, compute and store the sample mean and sample stdev
for ns in range(num_samples):
    sample = np.____.choice(population, num_pts)
    means[ns] = sample.____()
    stdevs[ns] = sample.____()

# Compute and print the mean() and std() for the sample statistic distributions
print("Means:  center={:>6.2f}, spread={:>6.2f}".format(means.mean(), means.std()))
print("Stdevs: center={:>6.2f}, spread={:>6.2f}".format(stdevs.____(), stdevs.____()))
Code bewerken en uitvoeren