1. Učit se
    2. /
  2. Kurzy
    3. /
  3. Statistical Thinking in Python (Part 2)
Connected

cvičení

Bootstrapové replikáty jiných statistik

V předchozím cvičení jsme viděli, že průměr má normální rozdělení. To ale neplatí nutně pro jiné statistiky – žádný strach ale: jako hackeři si vždy poradíme s bootstrapem! V tomto cvičení vygeneruješ bootstrapové replikáty rozptylu ročního srážkového úhrnu na meteorologické stanici Sheffield a vykreslíš histogram těchto replikátů.

Využiješ přitom funkci draw_bs_reps(), kterou jsi definoval/a v jednom z předchozích cvičení. Pro připomenutí je uvedena níže:

def draw_bs_reps(data, func, size=1):
    """Draw bootstrap replicates."""
    # Initialize array of replicates
    bs_replicates = np.empty(size)
    # Generate replicates
    for i in range(size):
        bs_replicates[i] = bootstrap_replicate_1d(data, func)
    return bs_replicates

Pokyny

100 XP
  • Pomocí funkce draw_bs_reps() vygeneruj 10000 bootstrapových replikátů rozptylu ročního srážkového úhrnu z datasetu rainfall. Nápověda: Pro výpočet rozptylu předej np.var.
  • Vydělte replikáty rozptylu (bs_replicates) hodnotou 100, aby byl rozptyl vyjádřen v jednotkách čtverečních centimetrů.
  • Vykresli histogram hodnot bs_replicates s použitím argumentu normed=True a 50 košů.