تحليل البيانات / المتقدم: قرارات مبررة وقابلة للتكرار / الدرس 15
إعادة المعاينة والاختبارات وتعدد المقارنات
استخدم إعادة المعاينة لشرح عدم اليقين تحت افتراضات محددة، وافصل الاختيار الاستكشافي من الدليل التأكيدي.
في هذه الصفحة
أهداف الدرس
- تمييز سحب bootstrap من تبديل تسميات المجموعات.
- تفسير قيم p ومعدلات الخطأ دون ادعاء احتمال حقيقة الفرض.
- مراعاة عائلة مقارنات محددة ومعلنة.
لا تنشئ إعادة المعاينة معلومات مصدر جديدة
تعيد المعاينة استخدام البيانات المرصودة بإجراء محدد لتقريب توزيع مرجعي أو عده بدقة. تسحب bootstrap المشاهدات مع الإرجاع لدراسة تغير المقدر، بينما تعيد طرق التبديل توزيع التسميات تحت حجة قابلية تبادل أو تخصيص عشوائي. يختلف السؤال والافتراض بينهما. تقلل زيادة التكرارات ضجيج المحاكاة، ولا تصلح انحياز اختيار العينة الأصلية.
قد يكسر سحب الصفوف المستقلة في البيانات العنقودية أو الزمنية البنية التي تحدد عدم اليقين. ربما يلزم سحب عناقيد أو كتل بتصميم يناسب الهدف. وقد يكون توزيع bootstrap للعينة الضئيلة متقطعا جدا وغير مفيد لذيول لم ترصد. افصل سهولة الحساب من صحة الاستدلال.
عد bootstrap صغيرة بدقة
يعد المثال الأول كل سحب مرتب من عنصرين مع الإرجاع للقيمتين1 و3. توجد أربع عينات متساوية الوزن تحت هذا التوزيع التجريبي الصغير، ومتوسطاتها1 و2 و2 و3، ومتوسطها العام2. هذا عرض حسابي واضح، لا إثبات أن فترة bootstrap من مشاهدتين تملك تغطية موثوقة للمجتمع.
يسجل العمل الواقعي وحدة السحب والإحصائية وطريقة الفترة والبذرة وعدد السحوبات ونسخ الأدوات. لا تختر صيغة فترة لأنها سهلة فقط. افحص هل يدعم الاستقلال والحجم وسلوك المقدر تلك الطريقة، وقارنها بحساسية مناسبة للمجال.
from itertools import product
from statistics import mean
values = [1, 3]
means = [mean(sample) for sample in product(values, repeat=2)]
print("means", ",".join(map(str, means)))
print(f"resamples={len(means)} average={mean(means):.1f}")
means 1,2,2,3
resamples=4 average=2.0
تنتج السحوبات الأربع متوسطات1 و2 و2 و3 دون عشوائية مونت كارلو.
تحديد توزيع العشوائية تحت الفرض الصفري
يفحص المثال الثاني الطرق الست لتخصيص نتيجتين من أربع نتائج ثابتة للمجموعةA. يعد فروق المتوسطات المطلقة التي لا تقل عن الفرق المشاهد، وهو اثنان. تحقق طريقتان الشرط، فتعطيان نسبة تبديل دقيقة ثنائية الطرف مقدارها الثلث. يحتاج ذلك مرجعا صفريا مبررا متساوي التخصيص، لا مجرد إمكان خلط مجموعتين رصديتين.
تصف p التطرف تحت الإجراء الصفري المحدد، لا احتمال صحة الفرض ولا احتمال أن النتيجة خطأ ولا حجم الأثر العملي. حدد الإحصائية والطرف، وأدخل التعادلات وفق القاعدة المعلنة. تحتاج تقديرات مونت كارلو أيضا إلى سياسة واضحة لعدد المحاكاة المحدود.
from itertools import combinations
from statistics import mean
values = [1, 2, 3, 4]
observed = 2.0
differences = []
for indices in combinations(range(4), 2):
chosen = set(indices)
first = [value for index, value in enumerate(values) if index in chosen]
second = [value for index, value in enumerate(values) if index not in chosen]
differences.append(abs(mean(second) - mean(first)))
extreme = sum(value >= observed for value in differences)
print(f"assignments={len(differences)} extreme={extreme} p={extreme / len(differences):.4f}")
assignments=6 extreme=2 p=0.3333
لطريقتين من طرق التخصيص الست المتساوية الوزن فرق مطلق لا يقل عن اثنين.
تغير المحاولات الكثيرة فرص الخطأ
يصنع اختبار مخرجات وفئات وفترات وتحويلات كثيرة عائلة فرص لنتيجة تبدو مهمة بالمصادفة. لا يحفظ انتقاء أصغر p بعد الاستكشاف تفسير الخطأ نفسه لاختبار واحد معلن مسبقا. حدد العائلة والغرض قبل اختيار التصحيح؛ فضبط الخطأ العائلي وضبط نسبة الاكتشافات الكاذبة يستهدفان كميتين مختلفتين.
تقارن Bonferroni كل قيمة من m قيم p صالحة بـalpha/m، وتوفر حدا عائليا محافظا دون اشتراط الاستقلال. لا تصلح p غير صالحة أو انتقاء مخرجات مخفيا أو عائلة سيئة التعريف. كذلك يحتاج الرصد المتتابع خطة؛ فتكرار عتبة عينة ثابتة حتى تعبر يغير سلوكها.
ربط الاستدلال بالأثر والقرار
اعرض تقديرات الأثر وعدم اليقين وتصميم العينة وتفاصيل الاختبار وكل مقارنة مناسبة، بما فيها غير المثيرة. قد تعكس p الكبيرة معلومات ضعيفة لا تكافؤا، وقد تصاحب الصغيرة أثرا تافها. افصل الاستنتاج التأكيدي من فرضية نشأت لدراسة مستقبلية.
تشمل قابلية التكرار خوارزمية إعادة المعاينة وسياق السلسلة العشوائية، لا البذرة وحدها. استخدم العَد الحتمي للأمثلة الصغيرة ومولدات محلية للمحاكاة. افحص الاستقرار العددي وخطأ المحاكاة، لكن لا تجعل كثافة الحساب تخفي الافتراض الذي يجعل التوزيع المرجعي مناسبا.
الأخطاء الشائعة وإصلاحها
تعد زيادة سحوبات bootstrap إصلاحا لعينة منحازة.
- لماذا يحدث؟
- اختلطت دقة المحاكاة بتمثيل المصدر.
- كيف تصلحه؟
- عالج المعاينة والاختيار وافتراضاتهما مستقلا.
- كيف تتجنب تكراره؟
- اعرض وحدة إعادة السحب وقيود المصدر.
تسمى p احتمال عدم وجود أثر.
- لماذا يحدث؟
- عُكس تفسير احتمال مشروط بتوزيع صفري.
- كيف تصلحه؟
- اذكر الإجراء الصفري واعرض عدم يقين الأثر منفصلا.
- كيف تتجنب تكراره؟
- راجع نص التفسير بعناية مراجعة الشيفرة.
تعرض أفضل نتيجة لفئة فرعية فقط.
- لماذا يحدث؟
- يخفي الاختيار الاستكشافي عائلة المقارنات.
- كيف تصلحه؟
- اعرض البحث واستخدم خطة تعدد أو تحقق مبررة.
- كيف تتجنب تكراره؟
- أعلن المخرجات الأساسية مسبقا حيث يمكن.
اختبر فهمك
1. تشكل أربعة اختبارات صالحة عائلة معلنة عند alpha0.05. ما عتبة Bonferroni وما الذي لا تضمنه؟
اعرض الحل المشروح
تساوي العتبة0.05/4، أي0.0125. مع صلاحية قيم p الفردية تحد الطريقة احتمال رفض خاطئ واحد على الأقل في العائلة بـalpha المختارة. لا تثبت كبر الأثر ولا تصلح الإرباك ولا تبرر اختبارات مخفية إضافية. اذكر العائلة وتقديرات الأثر مع القرارات المعدلة.
2. لماذا لا تستخدم bootstrap لصفوف مستقلة تلقائيا لمشاهدات يومية قوية الارتباط المتسلسل؟
اعرض الحل المشروح
يدمر سحب الأيام المنفردة نمط الاعتماد وقد يجعل تغير المقدر أصغر ظاهريا أو خاطئا بطريقة أخرى. اختر إجراء زمنيا أو كتليا مبررا وتصميما مناسبا للكتل، أو نموذجا آخر يطابق العملية. لا تفعل زيادة سحب الوحدة الخاطئة إلا تقريب المرجع الخاطئ بدقة أكبر.
أهم ما تعلمته
- تحتاج bootstrap والتبديل حجتي صلاحية مختلفتين.
- ليست p احتمالا أن الفرض صحيح.
- يجب أن يظهر التعدد والاختيار وحجم الأثر في الاستدلال المعروض.