ÎncepețiÎncepe gratuit

Când va avea loc următorul mare cutremur din Parkfield?

Ultimul cutremur major din regiunea Parkfield a avut loc în seara zilei de 27 septembrie 2004, ora locală. Sarcina ta este să estimezi când va fi următorul cutremur din Parkfield, presupunând atât modelul Exponențial, cât și modelul Gaussian. În ambele cazuri, cea mai bună estimare este dată de intervalul mediu de timp, pe care l-ai calculat în exercițiul anterior ca fiind 24,62 ani – ceea ce înseamnă că următorul cutremur ar urma să aibă loc în 2029. Calculează intervalele de încredere de 95% pentru momentul următorului cutremur, presupunând o distribuție Exponențială parametrizată cu mean_time_gap, calculat în exercițiul anterior. Fă același lucru presupunând o distribuție Normală parametrizată cu mean_time_gap și std_time_gap.

Acest exercițiu face parte din cursul

Studii de caz în gândirea statistică

Vezi cursul

Instrucțiuni pentru exercițiu

  • Generează 100.000 de eșantioane dintr-o distribuție Exponențială cu media dată de mean_time_gap. Stochează rezultatul în exp_samples.
  • Generează 100.000 de eșantioane dintr-o distribuție Normală cu media dată de mean_time_gap și deviația standard dată de std_time_gap. Stochează rezultatul în norm_samples.
  • Deoarece până în prezent nu a mai avut loc niciun cutremur în Parkfield, filtrează eșantioanele care sunt mai mari decât today - last_quake, unde am stocat anul zecimal al zilei de azi în today, iar last_quake = 2004.74 reprezintă anul zecimal al ultimului cutremur din Parkfield. Suprascrie variabilele exp_samples și norm_samples cu aceste arrays-uri filtrate.
  • Folosește np.percentile() pentru a calcula intervalul de încredere de 95% pentru momentul următorului cutremur din Parkfield. În același apel al funcției, poți calcula și mediana incluzând percentila 50.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# Draw samples from the Exponential distribution: exp_samples
exp_samples = ____

# Draw samples from the Normal distribution: norm_samples
norm_samples = ____

# No earthquake as of today, so only keep samples that are long enough
exp_samples = ____[____ > ____ - ____]
norm_samples = ____[____ > ____ - ____]

# Compute the confidence intervals with medians
conf_int_exp = ____(____, [____, ____, ____]) + last_quake
conf_int_norm = ____(____, [____, ____, ____]) + last_quake

# Print the results
print('Exponential:', conf_int_exp)
print('     Normal:', conf_int_norm)
Editează și rulează codul