IniziaInizia gratis

Test di permutazione sui dati delle rane

La forza media del colpo della Rana A era 0,71 Newton (N) e quella della Rana B era 0,42 N, per una differenza di 0,29 N. È possibile che le rane colpiscano con la stessa forza e che questa differenza osservata sia dovuta al caso. Calcolerai la probabilità di ottenere una differenza di almeno 0,29 N nella forza media del colpo sotto l'ipotesi che le distribuzioni delle forze di colpo delle due rane siano identiche. Usiamo un test di permutazione con come statistica test la differenza delle medie per verificare questa ipotesi.

Per tua comodità, i dati sono stati salvati negli array force_a e force_b.

Questo esercizio fa parte del corso

Pensiero statistico in Python (Parte 2)

Visualizza il corso

Istruzioni dell'esercizio

  • Definisci una funzione con firma diff_of_means(data_1, data_2) che restituisca la differenza tra le medie di due insiemi di dati, media di data_1 meno media di data_2.
  • Usa questa funzione per calcolare la differenza empirica delle medie osservata nelle rane.
  • Genera 10.000 repliche per permutazione della differenza delle medie.
  • Calcola il p-value.
  • Stampa il p-value.

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

def diff_of_means(data_1, data_2):
    """Difference in means of two arrays."""

    # The difference of means of data_1, data_2: diff
    diff = ____

    return diff

# Compute difference of mean impact force from experiment: empirical_diff_means
empirical_diff_means = ____

# Draw 10,000 permutation replicates: perm_replicates
perm_replicates = draw_perm_reps(____, ____,
                                 ____, size=10000)

# Compute p-value: p
p = np.sum(____ >= ____) / len(____)

# Print the result
print('p-value =', p)
Modifica ed esegui il codice