Коли станеться наступний великий землетрус у Парктфілді?

Останній великий землетрус у регіоні Парктфілд стався ввечері 27 вересня 2004 року за місцевим часом. Ваше завдання — оцінити, коли відбудеться наступний землетрус у Парктфілді, припускаючи експоненційну модель, а також гауссову модель. В обох випадках найкращу оцінку дає середній інтервал між подіями, який ви обчислили в попередній вправі як 24,62 року. Це означає, що наступний землетрус очікується у 2029 році. Обчисліть 95% довірчі інтервали для часу наступного землетрусу, припускаючи експоненційний розподіл з параметром mean_time_gap, який ви знайшли у попередній вправі. Зробіть те саме для нормального розподілу з параметрами mean_time_gap і std_time_gap.

Ця вправа є частиною курсу

Кейс-стаді зі статистичного мислення

Переглянути курс

Інструкції до вправи

Згенеруйте 100 000 вибірок з експоненційного розподілу із середнім, заданим mean_time_gap. Збережіть результат у exp_samples.
Згенеруйте 100 000 вибірок з нормального розподілу із середнім, заданим mean_time_gap, і стандартним відхиленням, заданим std_time_gap. Збережіть результат у norm_samples.
Оскільки станом на сьогодні землетрусу в Парктфілді не було, відфільтруйте значення вибірок, що більші за today - last_quake, де десяткове значення поточного року збережено в today, а last_quake = 2004.74 — десятковий рік останнього землетрусу в Парктфілді. Перезапишіть відповідні змінні exp_samples і norm_samples цими зрізаними масивами.
Використайте np.percentile(), щоб обчислити 95% довірчий інтервал для часу наступного землетрусу в Парктфілді. У тому самому виклику функції можна також обчислити медіану, включивши 50-й перцентиль.

Інтерактивна практична вправа

Спробуйте виконати цю вправу, доповнивши цей зразок коду.

# 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)

Редагувати та запускати код