Uji permutasi pada data katak
Gaya pukulan rata-rata Katak A adalah 0,71 Newton (N), dan Katak B adalah 0,42 N sehingga selisihnya 0,29 N. Ada kemungkinan kedua katak memukul dengan gaya yang sama dan perbedaan yang teramati ini terjadi karena kebetulan. Anda akan menghitung probabilitas memperoleh setidaknya selisih 0,29 N pada rata-rata gaya pukulan dengan hipotesis bahwa sebaran gaya pukulan kedua katak identik. Kita menggunakan uji permutasi dengan statistik uji berupa selisih rata-rata untuk menguji hipotesis ini.
Untuk kenyamanan Anda, data telah disimpan dalam larik force_a dan force_b.
Latihan ini adalah bagian dari kursus
Pemikiran Statistik dengan Python (Bagian 2)
Petunjuk latihan
- Definisikan fungsi dengan tanda panggil
diff_of_means(data_1, data_2)yang mengembalikan selisih rata-rata antara dua himpunan data, yaitu rata-ratadata_1dikurangi rata-ratadata_2. - Gunakan fungsi ini untuk menghitung selisih rata-rata empiris yang teramati pada data katak.
- Ambil 10.000 replikasi permutasi dari selisih rata-rata.
- Hitung p-value.
- Cetak p-value.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
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)