Párovací bootstrap pro data o gramotnosti a porodnosti
Pomocí funkce, kterou jsi právě napsal/a, proveď párovací bootstrap a vykresli histogram zachycující odhad směrnice z dat o negramotnosti a porodnosti. Také uveď 95% interval spolehlivosti pro směrnici. Data máš k dispozici v NumPy polích illiteracy a fertility.
Pro připomenutí: draw_bs_pairs_linreg() má signaturu draw_bs_pairs_linreg(x, y, size=1) a vrací dvě hodnoty: bs_slope_reps a bs_intercept_reps.
Toto cvičení je součástí kurzu
Statistical Thinking in Python (Part 2)
Pokyny k cvičení
- Pomocí funkce
draw_bs_pairs_linreg()vygeneruj1000bootstrapových replikátů směrnice a interceptu. Data na ose x jsouilliteracy, data na ose y jsoufertility. - Vypočítej a vypiš 95% bootstrapový interval spolehlivosti pro směrnici.
- Vykresli a zobraz histogram replikátů směrnice. Nezapomeň popsat osy. Tato část je už hotová, takže stačí kliknout na Submit a histogram se zobrazí!
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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()