تقدم الفصل 40%

فهرس الفصل

أدوات تفاعلية

دليل تحديد النماذج

AR(p): قطع PACF عند p
MA(q): قطع ACF عند q
ARMA: تناقص تدريجي

الفصل الثاني: منهجية بوكس جينكينز

منهجية بوكس جينكينز هي إطار منهجي لبناء نماذج السلاسل الزمنية ARIMA. تتكون هذه المنهجية من أربع مراحل متكررة تهدف إلى تحديد وتقدير وتشخيص نماذج السلاسل الزمنية وإجراء التنبؤات.

1

التعرف

تحديد النموذج

2

التقدير

تقدير المعلمات

3

التشخيص

فحص النموذج

4

التنبؤ

إجراء التنبؤات

الرسم التخطيطي لمنهجية بوكس جينكينز

1

مرحلة التعرف (التحديد)

هدف المرحلة:

تحديد النموذج المناسب للبيانات من خلال فحص خصائص السلسلة الزمنية، والتأكد من استقراريتها، وتحديد رتب النموذج (p, d, q).

1

فحص الاستقرارية

  • فحص الاتجاه والموسمية
  • اختبار ديكي-فولر للاستقرارية
  • التحويلات اللازمة (تفاضل، تحويل لوغاريتمي)
2

تحليل ACF و PACF

  • رسم دالة الارتباط الذاتي (ACF)
  • رسم دالة الارتباط الذاتي الجزئي (PACF)
  • تحديد أنماط التناقص والقطع
3

تحديد النموذج

  • استخدام قواعد تحديد النموذج
  • اختبار عدة نماذج محتملة
  • التعامل مع الموسمية (P, D, Q)

أداة تفاعلية: مرحلة التعرف

توليد سلسلة زمنية
100 300 1000
اختبار الاستقرارية
اختبار ديكي-فولر
إحصائية الاختبار: -
القيمة الحرجة (5%): -2.86
النتيجة: لم يتم الاختبار
تحليل السلسلة
السلسلة الزمنية
دالة ACF
دالة PACF
توصية النموذج

توليد سلسلة أولاً لرؤية التوصية...

دليل تحديد نموذج ARIMA من ACF و PACF

النموذج سلوك ACF سلوك PACF رتب النموذج
AR(p) تناقص تدريجي (أسي/جيبي) قطع بعد التأخر p p من PACF، q=0
MA(q) قطع بعد التأخر q تناقص تدريجي q من ACF، p=0
ARMA(p,q) تناقص تدريجي بعد التأخر q تناقص تدريجي بعد التأخر p p و q غير واضحين، تجريب
ARIMA(p,d,q) بطيء التناقص (عدم استقرار) بطيء التناقص تفاضل d مرات أولاً
2

مرحلة التقدير

هدف المرحلة:

تقدير معاملات النموذج المحدد باستخدام الطرق الإحصائية المناسبة مثل المربعات الصغرى أو طريقة الامكان الأعظم.

طرق التقدير

1. المربعات الصغرى (LS)

تُستخدم لنماذج AR. تعمل على تقليل مجموع مربعات البواقي: \[ \min \sum_{t=1}^n \varepsilon_t^2 \]

2. الامكان الأعظم (MLE)

الطريقة المفضلة لـ ARMA. تعظم دالة الامكان: \[ L(\theta) = \prod_{t=1}^n f(\varepsilon_t|\theta) \]

3. لحظات العينة

تُستخدم لـ MA. تحل معادلات بناءً على لحظات العينة (الارتباطات الذاتية).

معايير الجودة

دقة التقدير
  • • أخطاء معيارية صغيرة
  • • إحصائيات t كبيرة (قيم p صغيرة)
  • • فترات ثقة ضيقة
معايير المعلومات
  • • AIC: \( \text{AIC} = -2\log(L) + 2k \)
  • • BIC: \( \text{BIC} = -2\log(L) + k\log(n) \)
  • • HQIC: \( \text{HQIC} = -2\log(L) + 2k\log(\log(n)) \)

أداة تفاعلية: مقارنة طرق التقدير

تقدير معاملات النموذج

3

مرحلة التشخيص

هدف المرحلة:

فحص جودة النموذج المقدر والتحقق من كفاءته من خلال تحليل البواقي واختبار فروض النموذج.

1. تحليل البواقي
  • فحص استقلالية البواقي (اختبار Ljung-Box)
  • فحص الثبات (عدم وجود ارتباط ذاتي)
  • فحص التوزيع الطبيعي (اختبار Jarque-Bera)
2. اختبارات شكلية
  • اختبار عدم ثبات التباين (ARCH)
  • اختبار عدم الخطية
  • اختبار استقرار المعلمات
3. معايير كفاءة النموذج
  • مقارنة مع نماذج بديلة
  • تحليل القوة التنبؤية
  • اختبارات خارج العينة

أداة تفاعلية: تشخيص النموذج

تفسير نتائج اختبارات التشخيص

الاختبار الفرضية الصفرية متى تقبل الفرضية الصفرية؟ ماذا يعني ذلك؟ ما العمل؟
Ljung-Box البواقي بيضاء (لا ارتباط ذاتي) قيمة p > 0.05 النموذج مناسب تابع للمرحلة التالية
Jarque-Bera البواقي تتبع التوزيع الطبيعي قيمة p > 0.05 التوزيع طبيعي (جيد لـ MLE) لا حاجة للتحويل
ARCH لا يوجد تأثير ARCH قيمة p > 0.05 تباين ثابت استمرار مع ARIMA
رامزي RESET النموذج صحيح الشكل قيمة p > 0.05 الشكل الخطي مناسب لا حاجة لنموذج غير خطي
4

مرحلة التنبؤ

هدف المرحلة:

استخدام النموذج المختار والمقدر لإجراء تنبؤات للمستقبل، مع حساب فترات الثقة للتنبؤات وتقييم دقة النموذج.

طرق التنبؤ

تنبؤ نقطي

تقدير القيمة المتوقعة المستقبلية: \[ \hat{X}_{t+h} = E[X_{t+h} | X_t, X_{t-1}, \dots] \]

فترات الثقة

حساب مجال الثقة للتنبؤ: \[ \hat{X}_{t+h} \pm z_{\alpha/2} \cdot \text{SE} \] حيث SE هو الخطأ المعياري للتنبؤ.

مقاييس الدقة
  • • MSE: متوسط مربعات الخطأ
  • • MAE: متوسط الخطأ المطلق
  • • MAPE: متوسط النسبة المئوية للخطأ
  • • RMSE: جذر متوسط مربعات الخطأ

خصائص تنبؤات ARIMA

الانحدار نحو المتوسط

تنبؤات نماذج ARIMA المستقرة تتجه نحو متوسط السلسلة مع زيادة أفق التنبؤ.

تزايد فترات الثقة

فترات الثقة تتسع مع زيادة أفق التنبؤ بسبب تراكم عدم اليقين.

تحديث التنبؤات

عندما تتوفر بيانات جديدة، يمكن تحديث التنبؤات بكفاءة باستخدام خوارزمية التحديث.

أداة تفاعلية: التنبؤ باستخدام ARIMA

إعداد النموذج
50% 80% 90%
إعدادات التنبؤ
5 20 50
توليد البيانات

مقارنة أداء النماذج في التنبؤ

5

مثال شامل: تطبيق منهجية بوكس جينكينز

المشكلة:

لدينا سلسلة زمنية لعدد الزوار الشهري لمتحف خلال 5 سنوات (60 مشاهدة). نريد بناء نموذج ARIMA للتنبؤ بعدد الزوار للشهور القادمة.

1

تحميل البيانات والفحص البصري

البيانات: عدد الزوار الشهري لمتحف (يناير 2018 - ديسمبر 2022)

2

مرحلة التعرف: فحص الاستقرارية والموسمية

السلسلة الأصلية
بعد الفروق (d=1)
اختبار ديكي-فولر

اضغط "تحليل الاستقرارية" لعرض النتائج

3

مرحلة التعرف: تحليل ACF و PACF

ACF للسلسلة المتفاضلة
PACF للسلسلة المتفاضلة
توصية النموذج بناءً على ACF و PACF:

سيظهر بعد تحليل ACF و PACF...

4

مرحلة التقدير: تقدير النموذج المختار

5

مرحلة التشخيص: فحص كفاءة النموذج

6

مرحلة التنبؤ: التنبؤ بعدد الزوار المستقبلي

7

الخلاصة والتوصيات

سيقوم بتشغيل جميع خطوات منهجية بوكس جينكينز تلقائياً

تمارين الفصل الثاني

التمرين 1: تحديد النموذج من ACF و PACF

بالنظر إلى الرسوم البيانية التالية لـ ACF و PACF، حدد النموذج المناسب:

ACF
PACF

التمرين 2: تفسير نتائج اختبار Ljung-Box

بالنظر إلى نتائج اختبار Ljung-Box التالية للبواقي، هل النموذج مناسب؟

نتائج اختبار Ljung-Box:
إحصائية الاختبار (Q): 15.32
درجات الحرية (df): 10
قيمة p: 0.12

التمرين 3: مقارنة نماذج ARIMA

بالنظر إلى مقاييس الأداء التالية، أي نموذج تختار ولماذا؟

النموذج AIC BIC RMSE MAPE (%)
ARIMA(1,0,0) 245.3 250.1 12.5 5.2
ARIMA(0,0,1) 248.7 253.5 13.1 5.5
ARIMA(1,0,1) 242.8 249.9 11.9 4.8
ARIMA(2,0,2) 241.5 251.0 11.7 4.6