Pairs-bootstrap van analfabetisme/vruchtbaarheid-gegevens
Gebruik de functie die je zojuist hebt geschreven om een pairs-bootstrap uit te voeren en een histogram te plotten dat de schatting van de helling beschrijft op basis van de analfabetisme-/vruchtbaarheidsgegevens. Rapporteer ook het 95%-betrouwbaarheidsinterval van de helling. De gegevens staan voor je klaar in de NumPy-arrays illiteracy en fertility.
Als herinnering: draw_bs_pairs_linreg() heeft de functiehandtekening draw_bs_pairs_linreg(x, y, size=1) en geeft twee waarden terug: bs_slope_reps en bs_intercept_reps.
Deze oefening maakt deel uit van de cursus
Statistical Thinking in Python (deel 2)
Oefeninstructies
- Gebruik je
draw_bs_pairs_linreg()-functie om1000bootstrap-replicaties van de helling en de intercept te nemen. De x-asgegevens zijnilliteracyen de y-asgegevens zijnfertility. - Bereken en print het 95% bootstrap-betrouwbaarheidsinterval voor de helling.
- Plot en toon een histogram van de hellingreplicaties. Vergeet niet je assen te labelen. Dit is al voor je gedaan, dus klik op Verzenden om je histogram te zien!
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Generate replicates of slope and intercept using pairs bootstrap
bs_slope_reps, bs_intercept_reps = ____
# Compute and print 95% CI for slope
print(np.percentile(____, ____))
# Plot the histogram
_ = plt.hist(bs_slope_reps, bins=50, normed=True)
_ = plt.xlabel('slope')
_ = plt.ylabel('PDF')
plt.show()