Začněte nyníZačněte zdarma

Lineární regrese

Proveď lineární regresi pro data z roku 1975 i 2012. Poté odhadni parametry regrese pomocí párového bootstrapu. Uveď 95% intervaly spolehlivosti pro sklon a průsečík regresní přímky.

Použ funkci draw_bs_pairs_linreg(), kterou jsi napsal/a ve 2. kapitole.

Pro připomenutí: její volání vypadá takto: draw_bs_pairs_linreg(x, y, size=1) a vrací bs_slope_reps a bs_intercept_reps. Data délky zobáku jsou uložena jako bl_1975 a bl_2012, data hloubky zobáku pak jako bd_1975 a bd_2012.

Toto cvičení je součástí kurzu

Statistical Thinking in Python (Part 2)

Zobrazit kurz

Pokyny k cvičení

  • Vypočítej sklon a průsečík pro datové sady z roku 1975 i 2012.
  • Pomocí funkce draw_bs_pairs_linreg() získej 1000 párových bootstrap vzorků pro lineární regrese.
  • Vypočítej 95% intervaly spolehlivosti pro sklony a průsečíky.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# Compute the linear regressions
slope_1975, intercept_1975 = ____
slope_2012, intercept_2012 = ____

# Perform pairs bootstrap for the linear regressions
bs_slope_reps_1975, bs_intercept_reps_1975 = \
        ____
bs_slope_reps_2012, bs_intercept_reps_2012 = \
        ____

# Compute confidence intervals of slopes
slope_conf_int_1975 = ____
slope_conf_int_2012 = ____
intercept_conf_int_1975 = ____

intercept_conf_int_2012 = ____


# Print the results
print('1975: slope =', slope_1975,
      'conf int =', slope_conf_int_1975)
print('1975: intercept =', intercept_1975,
      'conf int =', intercept_conf_int_1975)
print('2012: slope =', slope_2012,
      'conf int =', slope_conf_int_2012)
print('2012: intercept =', intercept_2012,
      'conf int =', intercept_conf_int_2012)
Upravit a spustit kód