Pars-bootstrap för läs- och skrivkunnighet/fertilitet
Använd funktionen du precis skrev och utför pars-bootstrap för att rita ett histogram som beskriver skattningen av lutningen från data om analfabetism/fertilitet. Rapportera även 95-procentigt konfidensintervall för lutningen. Data finns tillgänglig i NumPy-arrayerna illiteracy och fertility.
Som påminnelse: draw_bs_pairs_linreg() har signaturen draw_bs_pairs_linreg(x, y, size=1) och returnerar två värden: bs_slope_reps och bs_intercept_reps.
Den här övningen är en del av kursen
Statistiskt tänkande i Python (del 2)
Övningsinstruktioner
- Använd din funktion
draw_bs_pairs_linreg()för att ta1000bootstrap-replikat av lutningen och interceptet. Data för x-axeln ärilliteracyoch för y-axelnfertility. - Beräkna och skriv ut 95-procentigt bootstrap-konfidensintervall för lutningen.
- Rita och visa ett histogram över replikaten av lutningen. Se till att märka ut dina axlar. Detta har redan gjorts åt dig, så klicka på Skicka in svar för att se ditt histogram!
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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()