Začněte nyníZačněte zdarma

Odhady doby mezi zemětřeseními v Parkfieldu

V tomto cvičení nejprve vypočítáš nejlepší odhady parametrů pro exponenciální a gaussovský model doby mezi zemětřeseními. Poté vykreslíš teoretické CDF pro jednotlivé modely spolu s empirickou ECDF skutečných časových intervalů mezi zemětřeseními v Parkfieldu.

Toto cvičení je součástí kurzu

Případové studie ve statistickém myšlení

Zobrazit kurz

Pokyny k cvičení

  • Vypočítej průměrný časový interval mezi zemětřeseními a ulož ho jako mean_time_gap. Časové intervaly mezi hlavními zemětřeseními v letech jsou uloženy v time_gap.
  • Vypočítej směrodatnou odchylku časových intervalů a ulož ji jako std_time_gap.
  • Pomocí np.random.exponential() vygeneruj 10 000 vzorků z exponenciálního rozdělení s odpovídajícím průměrem. Ulož je do proměnné time_gap_exp.
  • Pomocí np.random.normal() vygeneruj 10 000 vzorků z normálního rozdělení s odpovídajícím průměrem a směrodatnou odchylkou. Ulož je do proměnné time_gap_norm.
  • Vykresli teoretické CDF, každou zvlášť na jednom řádku, pomocí přístupu *dcst.ecdf() představeného dříve v této kapitole.
  • Vykresli ECDF s klíčovými argumenty formal=True, min_x=-10 a max_x=50.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

# Compute the mean time gap: mean_time_gap
mean_time_gap = ____

# Standard deviation of the time gap: std_time_gap
std_time_gap = ____

# Generate theoretical Exponential distribution of timings: time_gap_exp
time_gap_exp = ____

# Generate theoretical Normal distribution of timings: time_gap_norm
time_gap_norm = ____

# Plot theoretical CDFs
_ = plt.plot(*____)
_ = plt.plot(*____)

# Plot Parkfield ECDF
_ = plt.plot(*____(____, ____=____, ____=____, ____=____))

# Add legend
_ = plt.legend(('Exp.', 'Norm.'), loc='upper left')

# Label axes, set limits and show plot
_ = plt.xlabel('time gap (years)')
_ = plt.ylabel('ECDF')
_ = plt.xlim(-10, 50)
plt.show()
Upravit a spustit kód