Üstel dağılım için K-S testi
Parkfield dizisindeki depremler arası sürelerin Üstel dağılıma uyduğu sıfır hipotezini test et. Yani, depremler son depremin ne zaman olduğuna dair hiçbir hafıza olmadan rastgele gerçekleşir. Not: Bu hesaplama hesaplama açısından yoğundur (108'den fazla rastgele sayı üreteceksin), bu yüzden yaklaşık 10 saniye sürecektir.
Bu egzersiz
İstatistiksel Düşünmede Vaka Çalışmaları
kursunun bir parçasıdırEgzersiz talimatları
np.random.exponential()kullanarak Üstel dağılımdan 10.000 tekrar çiz. Depremler arasındaki ortalama zaman aralığı, önceki bir egzersizde hesapladığın vemean_time_gapolarak saklanan değerdir. Sonucux_fiçinde sakla.- Bu örnekleri (
x_f), gerçek zaman aralıkları (time_gap) ile birlikte kullanarakdcst.ks_stat()ile Kolmogorov–Smirnov istatistiğini hesapla. - Son egzersizde yazdığın ve artık kullanışlı bir şekilde
dcst.draw_ks_reps()olarak saklanan fonksiyonu kullanarak Üstel dağılımdan 10.000 K-S tekrarı çiz. Hedef Üstel dağılımdan örneklemek içinsize=10000anahtar sözcük argümanını kullan. Tekrarlarırepsolarak sakla. - p-değerini hesapla ve yazdır. Unutma: Bu durumda "en az bunun kadar uç" ifadesi, sıfır hipotezi altında elde edilen test istatistiğinin gözlenen değerden büyük veya ona eşit olması olarak tanımlanır.
Uygulamalı interaktif egzersiz
Bu örnek kodu tamamlayarak bu egzersizi bitirin.
# Draw target distribution: x_f
x_f = ____
# Compute K-S stat: d
d = ____
# Draw K-S replicates: reps
reps = ____(len(____), ____,
args=(mean_time_gap,), size=____, n_reps=____)
# Compute and print p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)