Vykreslení bootstrapových regresí
Skvělý způsob, jak vizualizovat variabilitu, kterou bychom mohli očekávat v lineární regresi, je vykreslit přímku odpovídající každému bootstrapovému replikátu sklonu a průsečíku. Proveď to pro prvních 100 bootstrapových replikátů sklonu a průsečíku (uložených jako bs_slope_reps a bs_intercept_reps).
Toto cvičení je součástí kurzu
Statistical Thinking in Python (Part 2)
Pokyny k cvičení
- Vytvoř pole \(x\)-hodnot obsahující
0a100pro vykreslení regresních přímek. Použij k tomu funkcinp.array(). - Napiš
forcyklus, ve kterém vykreslíš regresní přímku se sklonem a průsečíkem daným dvojicemi bootstrapových replikátů. Proveď to pro100přímek.- Při vykreslování regresních přímek v každé iteraci
forcyklu vycházej z rovnice regresey = a*x + b. Zde jearovnobs_slope_reps[i]abjebs_intercept_reps[i]. - Ve volání
plt.plot()zadej klíčové argumentylinewidth=0.5,alpha=0.2acolor='red'.
- Při vykreslování regresních přímek v každé iteraci
- Vytvoř bodový graf s
illiteracyna ose x afertilityna ose y. Nezapomeň zadat klíčové argumentymarker='.'alinestyle='none'. - Popisky os, 2% okraj a zobrazení grafu jsou již připraveny za tebe — stačí kliknout na Submit a bootstrapové regrese 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 array of x-values for bootstrap lines: x
x = ____
# Plot the bootstrap lines
for i in range(____):
_ = plt.plot(____,
____*x + ____,
____=0.5, ____=0.2, ____='red')
# Plot the data
_ = ____
# Label axes, set the margins, and show the plot
_ = plt.xlabel('illiteracy')
_ = plt.ylabel('fertility')
plt.margins(0.02)
plt.show()