สุ่มตัวอย่าง 100,000 ค่าจากการแจกแจงแบบ Exponential โดยมีค่าเฉลี่ยเท่ากับ mean_time_gap แล้วเก็บผลลัพธ์ไว้ในตัวแปร exp_samples
สุ่มตัวอย่าง 100,000 ค่าจากการแจกแจงแบบ Normal โดยมีค่าเฉลี่ยเท่ากับ mean_time_gap และส่วนเบี่ยงเบนมาตรฐานเท่ากับ std_time_gap แล้วเก็บผลลัพธ์ไว้ในตัวแปร norm_samples
เนื่องจากยังไม่มีแผ่นดินไหวที่ Parkfield จนถึงปัจจุบัน ให้ตัดเฉพาะตัวอย่างที่มีค่ามากกว่า today - last_quake โดยที่ today คือปีปัจจุบันในรูปแบบทศนิยม และ last_quake = 2004.74 คือปีของแผ่นดินไหวครั้งล่าสุดที่ Parkfield ในรูปแบบเดียวกัน จากนั้นนำอาร์เรย์ที่ตัดแล้วไปบันทึกทับตัวแปร exp_samples และ norm_samples ตามลำดับ
ใช้ np.percentile() คำนวณช่วงความเชื่อมั่น 95% สำหรับเวลาที่แผ่นดินไหวครั้งต่อไปที่ Parkfield จะเกิดขึ้น และในการเรียกฟังก์ชันเดียวกัน สามารถคำนวณค่ามัธยฐานได้พร้อมกันโดยรวม percentile ที่ 50 เข้าไปด้วย

จังหวะเวลาของแผ่นดินไหวขนาดใหญ่และลำดับเหตุการณ์ที่ Parkfield

Create Your Free Account