Kom igångKom igång gratis

Hypotestest: förändrades jordbävningsfrekvensen?

Det är uppenbart att jordbävningsfrekvensen ökade kraftigt när injektionen av avloppsvatten påbörjades. Du kommer ändå att genomföra ett hypotestest för träningens skull. Du testar inte hypotesen att tiderna mellan jordbävningar har samma fördelning före och efter 2010, eftersom injektionen av avloppsvatten kan påverka fördelningen. Däremot antar du att medelvärdena är desamma. Beräkna alltså p-värdet för hypotesen att tiderna mellan jordbävningar före och efter 2010 har samma medelvärde, med medelvärdet för tidsgapen före 2010 minus medelvärdet för tidsgapen efter 2010 som teststatistika.

Den här övningen är en del av kursen

Fallstudier i statistiskt tänkande

Visa kurs

Övningsinstruktioner

  • Beräkna den observerade teststatistikan. Variablerna mean_dt_pre och mean_dt_post från tidigare övningar finns i ditt namnrymd.
  • Förskjut data efter 2010 så att de får samma medelvärde som data före 2010. Lagra resultatet som dt_post_shift.
  • Dra 10 000 bootstrap-replikat vardera av medelvärdet för dt_pre och dt_post_shift. Lagra resultaten i bs_reps_pre respektive bs_reps_post.
  • Beräkna replikat av medelvärdesskillnaden genom att subtrahera bs_reps_post från bs_reps_pre.
  • Beräkna och skriv ut p-värdet. Tolka "minst lika extremt som" som att teststatistikan är större än eller lika med det observerade värdet.

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# Compute the observed test statistic
mean_dt_diff = ____ - ____

# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____

# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____

# Get replicates of difference of means
bs_reps = ____ - ____

# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)
Redigera och kör kod