Test d'ipotesi: i valori di b sono diversi?
Esegui il test d'ipotesi abbozzato nell'esercizio precedente. Le variabili mags_pre e mags_post sono già caricate nello spazio dei nomi, così come mt = 3.
Questo esercizio fa parte del corso
Casi di studio nel pensiero statistico
Istruzioni dell'esercizio
- Seleziona le magnitudini dei terremoti precedenti al 2010 che hanno una magnitudine superiore (o uguale) alla soglia di completezza e sovrascrivi
mags_precon il risultato. Fai lo stesso permags_post. - Calcola la differenza osservata tra le medie delle magnitudini, sottraendo le magnitudini dei terremoti pre-2010 da quelle dei terremoti post-2010.
- Genera 10.000 repliche per permutazione usando
dcst.draw_perm_reps(). Usadcst.diff_of_meanscome argomento perfunc. - Calcola e stampa il p-value intendendo per "almeno altrettanto estremo" che la statistica test sia più piccola di quella osservata.
esercizio interattivo pratico
Prova questo esercizio completando questo codice di esempio.
# 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)