Hypothesis test: क्या भूकंप की आवृत्ति बदली?
स्पष्ट है कि जैसे ही वेस्टवॉटर इंजेक्शन शुरू हुआ, भूकंपों की आवृत्ति में भारी बढ़ोतरी हुई। फिर भी, अभ्यास के लिए आप एक परिकल्पना परीक्षण करेंगे। आप यह परिकल्पना नहीं जाँचेंगे कि 2010 से पहले और बाद के इंटर-भूकंप समयों का वितरण समान है, क्योंकि वेस्टवॉटर इंजेक्शन वितरण को प्रभावित कर सकता है। इसके बजाय, आप मानेंगे कि उनका औसत समान है। इसलिए, उस परिकल्पना से जुड़ी p-value निकालिए कि 2010 से पहले और बाद के इंटर-भूकंप समयों का औसत समान है, और टेस्ट स्टैटिस्टिक के रूप में 2010-पूर्व समयांतरों के औसत में से 2010-पश्चात समयांतरों के औसत को घटाइए।
यह अभ्यास पाठ्यक्रम का हिस्सा है
Statistical Thinking के केस स्टडीज़
अभ्यास निर्देश
- प्रेक्षित टेस्ट स्टैटिस्टिक निकालिए। पिछले अभ्यासों से
mean_dt_preऔरmean_dt_postवैरिएबल आपकी namespace में हैं. - 2010-पश्चात डेटा को 2010-पूर्व डेटा जितना ही औसत देने के लिए शिफ्ट करें। परिणाम
dt_post_shiftके रूप में सहेजें. dt_preऔरdt_post_shiftके माध्य के 10,000-10,000 बूटस्ट्रैप रिप्लिकेट्स ड्रॉ करें। क्रमशः परिणामbs_reps_preऔरbs_reps_postमें सहेजें.- माध्यों के अंतर के रिप्लिकेट्स निकालने के लिए
bs_reps_preमें सेbs_reps_postघटाएँ. - p-value निकालें और प्रिंट करें। "कम से कम इतना ही चरम" मानदंड के लिए मानें कि टेस्ट स्टैटिस्टिक प्रेक्षित मान के बराबर या उससे बड़ा है.
इंटरैक्टिव व्यावहारिक अभ्यास
इस अभ्यास को इस नमूना कोड को पूरा करके आज़माएँ।
# Compute the observed test statistic
mean_dt_diff = ____ - ____
# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____
# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____
# Get replicates of difference of means
bs_reps = ____ - ____
# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)