Uji hipotesis: apakah nilai b berbeda?
Lakukan uji hipotesis seperti yang digambarkan pada latihan sebelumnya. Variabel mags_pre dan mags_post sudah dimuat ke dalam namespace Anda, begitu juga mt = 3.
Latihan ini adalah bagian dari kursus
Studi Kasus dalam Pemikiran Statistik
Petunjuk latihan
- Ambil irisan magnitudo gempa sebelum 2010 yang memiliki magnitudo di atas (atau sama dengan) ambang kelengkapan dan timpa
mags_predengan hasilnya. Lakukan hal yang sama untukmags_post. - Hitung selisih mean magnitudo yang teramati, dengan mengurangkan magnitudo gempa pra-2010 dari gempa pasca-2010.
- Hasilkan 10.000 replikasi permutasi menggunakan
dcst.draw_perm_reps(). Gunakandcst.diff_of_meanssebagai argumen untukfunc. - Hitung dan cetak p-value dengan menafsirkan "setidaknya sama ekstremnya" sebagai statistik uji lebih kecil daripada yang teramati.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
# Only magnitudes above completeness threshold
mags_pre = ____[____ >= ____]
mags_post = ____[____ >= ____]
# Observed difference in mean magnitudes: diff_obs
diff_obs = ____ - ____
# Generate permutation replicates: perm_reps
perm_reps = ____
# Compute and print p-value
p_val = ____(____ < ____) / 10000
print('p =', p_val)