1. Nauka
    2. /
  2. Kursy
    3. /
  3. Myślenie statystyczne w Pythonie (część 2)
Connected

ćwiczenie

Repliki bootstrap innych statystyk

W poprzednim ćwiczeniu przekonałeś się, że średnia ma rozkład normalny. Nie musi to jednak dotyczyć innych statystyk – nie ma jednak powodu do obaw: jako haker statystyczny zawsze możesz sięgnąć po repliki bootstrap! W tym ćwiczeniu wygenerujesz repliki bootstrap wariancji rocznych opadów na Stacji Meteorologicznej w Sheffield i wykreślisz histogram tych replik.

Skorzystasz tutaj z funkcji draw_bs_reps(), którą zdefiniowano kilka ćwiczeń wcześniej. Poniżej znajdziesz jej kod do wglądu:

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

Instrukcje

100 XP
  • Wylosuj 10000 replik bootstrap wariancji rocznych opadów ze zbioru danych rainfall, korzystając z funkcji draw_bs_reps(). Podpowiedź: Przekaż np.var jako argument do obliczania wariancji.
  • Podziel repliki wariancji (bs_replicates) przez 100, aby wyrazić wariancję w wygodnych jednostkach – centymetrach kwadratowych.
  • Stwórz histogram bs_replicates, używając argumentu kluczowego normed=True oraz 50 przedziałów.