Uji hipotesis bootstrap satu-sampel
Seekor katak juvenil lain diteliti, Katak C, dan Anda ingin melihat apakah Katak B dan Katak C memiliki gaya tumbukan yang serupa. Sayangnya, Anda tidak memiliki data gaya tumbukan Katak C, tetapi Anda mengetahui rataan gaya tumbukannya sebesar 0,55 N. Karena Anda tidak memiliki data asli, Anda tidak dapat melakukan uji permutasi, dan Anda tidak dapat menilai hipotesis bahwa gaya dari Katak B dan Katak C berasal dari sebaran yang sama. Oleh karena itu, Anda akan menguji hipotesis lain yang kurang ketat: Rataan gaya pukulan Katak B sama dengan rataan gaya pukulan Katak C.
Untuk menyiapkan uji hipotesis bootstrap, Anda akan menggunakan rataan sebagai statistik uji. Ingat, tujuan Anda adalah menghitung probabilitas memperoleh rataan gaya tumbukan yang kurang dari atau sama dengan yang teramati untuk Katak B jika hipotesis bahwa rataan sebenarnya dari gaya tumbukan Katak B sama dengan rataan Katak C adalah benar. Pertama, Anda mentranslasi seluruh data Katak B sehingga rataan menjadi 0,55 N. Ini dilakukan dengan menambahkan rataan gaya Katak C dan mengurangkan rataan gaya Katak B dari setiap pengukuran Katak B. Cara ini mempertahankan properti lain dari sebaran Katak B, seperti varians, agar tetap tidak berubah.
Latihan ini adalah bagian dari kursus
Pemikiran Statistik dengan Python (Bagian 2)
Petunjuk latihan
- Translasi gaya tumbukan Katak B sehingga rataan (mean)‑nya menjadi 0,55 N.
- Gunakan fungsi
draw_bs_reps()Anda untuk mengambil 10.000 replikasi bootstrap dari rataan gaya yang sudah ditranslasi. - Hitung p-value dengan mencari fraksi replikasi bootstrap yang lebih kecil daripada rataan gaya tumbukan teramati milik Katak B. Perhatikan bahwa variabel yang digunakan di sini adalah
force_b. - Cetak p-value Anda.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Make an array of translated impact forces: translated_force_b
translated_force_b = ____
# Take bootstrap replicates of Frog B's translated impact forces: bs_replicates
bs_replicates = draw_bs_reps(____, ____, 10000)
# Compute fraction of replicates that are less than the observed Frog B force: p
p = np.sum(____ <= np.mean(____)) / 10000
# Print the p-value
print('p = ', ____)