Bootstrap güven aralıkları
İstatistiksel Düşünmede Vaka Çalışmaları
Justin Bois
Lecturer, Caltech
EDA ilk adımdır
"Keşifsel veri analizi tek başına tüm hikâyeyi anlatamaz, ama başka hiçbir şey temel taş, ilk adım olamaz."
--John Tukey
Aktif atak uzunluğu ECDF'leri
1 Veri: Avni Gandhi, Grigorios Oikonomou ve David Prober, Caltech
En iyi parametre değeri
- En iyi (optimal) parametre değeri: Veriyi en iyi tanımlayan olasılık dağılımı parametresi
- Exponential dağılım için en iyi parametre: Verinin ortalamasından hesaplanır
np.mean(nuclear_incident_times)
87.140350877192986
1 Veri kaynağı: Wheatley, Sovacool ve Sornette, Nuclear Events Database
Bootstrap örneği
Verinin yeniden örneklenmiş bir dizisi
# nuclear_incident_times dizisinden yerine koyarak yeniden örnekle
bs_sample = np.random.choice(
nuclear_incident_times,
replace=True,
size=len(inter_times)
)
Bootstrap tekrarları
1 Veri kaynağı: Wheatley, Sovacool ve Sornette, Nuclear Events Database
Bootstrap tekrarları
1 Veri kaynağı: Wheatley, Sovacool ve Sornette, Nuclear Events Database
Bootstrap tekrarları
1 Veri kaynağı: Wheatley, Sovacool ve Sornette, Nuclear Events Database
Bootstrap tekrarları
1 Veri kaynağı: Wheatley, Sovacool ve Sornette, Nuclear Events Database
Bootstrap tekrarları
Bootstrap tekrarı: Bir bootstrap örneğinden hesaplanan istatistik
dcst.draw_bs_reps()
Bir veri kümesinden bootstrap tekrarları çeken fonksiyon
# nuclear_incident_times için ortalamanın
# 10000 tekrarını çek
bs_reps = dcst.draw_bs_reps(
nuclear_incident_times, np.mean, size=10000
)
Bootstrap güven aralığı
1 Veri kaynağı: Wheatley, Sovacool ve Sornette, Nuclear Events Database
Bootstrap güven aralığı
Ölçümleri defalarca yineleseydik, gözlenen değerlerin %p'si %p güven aralığı içinde olurdu
Bootstrap güven aralığı
np.percentile(bs_reps, [2.5, 97.5])
array([ 73.31505848, 102.39181287])
