अगला बड़ा Parkfield भूकंप कब आएगा?
Parkfield क्षेत्र में आखिरी बड़ा भूकंप 27 सितंबर, 2004 की शाम (लोकल टाइम) को आया था. आपका काम यह अनुमान लगाना है कि अगला Parkfield भूकंप कब आएगा, मानते हुए कि मॉडल Exponential है और साथ ही Gaussian भी. दोनों ही मामलों में, सबसे अच्छा अनुमान mean time gap से मिलता है, जिसे आपने पिछले अभ्यास में 24.62 वर्ष निकाला था, यानी अगला भूकंप 2029 में होगा. 95% confidence intervals निकालिए कि अगला भूकंप कब होगा, मानते हुए कि वितरण Exponential है और उसे आप पिछले अभ्यास में निकाले गए mean_time_gap से पैरामीटराइज़ कर रहे हैं. यही काम Normal वितरण के लिए भी कीजिए, जहाँ पैरामीटर mean_time_gap और std_time_gap हैं.
यह अभ्यास पाठ्यक्रम का हिस्सा है
Statistical Thinking के केस स्टडीज़
अभ्यास निर्देश
mean_time_gapद्वारा दिए गए mean के साथ Exponential वितरण से 100,000 सैंपल ड्रॉ कीजिए. परिणामexp_samplesमें स्टोर करें.mean_time_gapmean औरstd_time_gapstandard deviation के साथ Normal वितरण से 100,000 सैंपल ड्रॉ कीजिए. परिणामnorm_samplesमें स्टोर करें.- चूँकि आज तक Parkfield भूकंप नहीं आया है, वे सैंपल slice कर दीजिए जो
today - last_quakeसे बड़े हैं. यहाँ आज का दशमलव वर्षtodayमें स्टोर है, औरlast_quake = 2004.74(आखिरी Parkfield भूकंप का दशमलव वर्ष). इन sliced ऐरेज़ से संबंधितexp_samplesऔरnorm_samplesवैरिएबल्स को ओवरराइट कर दें. np.percentile()का उपयोग करके यह निकालिए कि अगला Parkfield भूकंप कब होगा, उसका 95% confidence interval क्या है. इसी फ़ंक्शन कॉल में आप 50th percentile शामिल करके median भी निकाल सकते हैं.
इंटरैक्टिव व्यावहारिक अभ्यास
इस अभ्यास को इस नमूना कोड को पूरा करके आज़माएँ।
# 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)