Оцінки інтервалів між землетрусами для Паркіфілда

У цій вправі ви спершу обчислите найкращі оцінки параметрів для експоненційної та нормальної моделей інтервалів між землетрусами. Потім побудуєте теоретичні CDF для відповідних моделей разом з формальною ECDF фактичних інтервалів між землетрусами в Паркіфілді.

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

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

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

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

Обчисліть середній інтервал між землетрусами та збережіть його як mean_time_gap. Проміжки часу між сильними землетрусами, у роках, збережено в time_gap.
Обчисліть стандартне відхилення інтервалів між землетрусами та збережіть його як std_time_gap.
Використайте np.random.exponential() щоб згенерувати 10 000 вибірок з експоненційного розподілу з відповідним середнім. Збережіть їх у змінній time_gap_exp.
Використайте np.random.normal() щоб згенерувати 10 000 вибірок із нормального розподілу з відповідними середнім та стандартним відхиленням. Збережіть їх у змінній time_gap_norm.
Побудуйте теоретичні CDF по одному рядку кожну, використовуючи підхід *dcst.ecdf(), який було представлено раніше в цьому розділі.
Побудуйте ECDF, використовуючи ключові аргументи formal=True, min_x=-10 і max_x=50.

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

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

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

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