Aan de slagGa gratis aan de slag

Permutatiereplicaten genereren

Zoals in de video besproken, is een permutatiereplicaat één enkele waarde van een statistiek die is berekend uit een permutatiesteekproef. Net zoals de functie draw_bs_reps() die je in hoofdstuk 2 schreef handig is om bootstrapreplicaten te genereren, is het nuttig om een vergelijkbare functie, draw_perm_reps(), te hebben om permutatiereplicaten te genereren. In deze oefening ga je deze handige functie schrijven.

De aanroep van de functie is draw_perm_reps(data_1, data_2, func, size=1). Belangrijk: func moet een functie zijn die twee arrays als argumenten accepteert. In de meeste gevallen is func een functie die je zelf schrijft.

Deze oefening maakt deel uit van de cursus

Statistical Thinking in Python (deel 2)

Cursus bekijken

Oefeninstructies

  • Definieer een functie met deze signatuur: draw_perm_reps(data_1, data_2, func, size=1).
    • Initialiseer een array om de permutatiereplicaten in op te slaan met np.empty().
    • Schrijf een for-lus om:
      • Een permutatiesteekproef te berekenen met je functie permutation_sample()
      • De steekproeven door te geven aan func() om het replicaat te berekenen en het resultaat op te slaan in je array met replicaten.
    • Retourneer de array met replicaten.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

def draw_perm_reps(data_1, data_2, func, size=1):
    """Generate multiple permutation replicates."""

    # Initialize array of replicates: perm_replicates
    perm_replicates = ____

    for i in ____:
        # Generate permutation sample
        perm_sample_1, perm_sample_2 = ____

        # Compute the test statistic
        perm_replicates[i] = ____

    return perm_replicates
Code bewerken en uitvoeren