Test de permutare pe datele broaștelor
Forța medie de lovire a Broaștei A a fost de 0,71 Newtoni (N), iar cea a Broaștei B de 0,42 N, rezultând o diferență de 0,29 N. Este posibil ca broaștele să lovească cu aceeași forță și ca această diferență observată să fie întâmplătoare. Vei calcula probabilitatea de a obține o diferență de cel puțin 0,29 N între forțele medii de lovire, în ipoteza că distribuțiile forțelor de lovire ale celor două broaște sunt identice. Vom folosi un test de permutare cu statistica de test bazată pe diferența dintre medii pentru a testa această ipoteză.
Pentru comoditate, datele sunt stocate în tablourile force_a și force_b.
Acest exercițiu face parte din cursul
Gândire statistică în Python (Partea 2)
Instrucțiuni pentru exercițiu
- Definește o funcție cu signatura
diff_of_means(data_1, data_2)care returnează diferența dintre medii pentru două seturi de date: media luidata_1minus media luidata_2. - Folosește această funcție pentru a calcula diferența empirică dintre medii observată în cazul broaștelor.
- Generează 10.000 de replici de permutare ale diferenței dintre medii.
- Calculează p-valoarea.
- Afișează p-valoarea.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
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)