مخطط الارتباط

رسم بياني يوضح 100 رقم عشوائي مع دالة جيبية "مخفية" ، وارتباط ذاتي (مخطط الارتباط) للسلسلة في الأسفل.

في تحليل البيانات، يُعدّ مخطط الارتباط الذاتي رسمًا بيانيًا لإحصاءات الارتباط . على سبيل المثال، في تحليل السلاسل الزمنية ، يُستخدم رسم بياني للارتباطات الذاتية للعينة.رح{\displaystyle r_{h}}عكسح{\displaystyle h\,}(الفترات الزمنية) هو مخطط الارتباط الذاتي . إذا تم رسم الارتباط المتبادل ، فإن النتيجة تسمى مخطط الارتباط المتبادل .

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

بالإضافة إلى ذلك، تُستخدم مخططات الارتباط الذاتي في مرحلة تحديد النموذج لنماذج السلاسل الزمنية ذات المتوسط ​​المتحرك للانحدار الذاتي من نوع بوكس-جينكينز . يجب أن تكون معاملات الارتباط الذاتي قريبة من الصفر لضمان العشوائية؛ فإذا لم يتحقق المحلل من العشوائية، فإن صحة العديد من الاستنتاجات الإحصائية تصبح موضع شك. يُعد مخطط الارتباط الذاتي وسيلة ممتازة للتحقق من هذه العشوائية.

في التحليل متعدد المتغيرات ، قد تُسمى مصفوفات الارتباط المعروضة كصور ملونة "مخططات الارتباط" أو "مخططات الارتباط". [ 1 ] [ 2 ] [ 3 ]

التطبيقات

يمكن أن يساعد مخطط الارتباط في تقديم إجابات للأسئلة التالية: [ 4 ]

  • هل البيانات عشوائية؟
  • هل هناك علاقة بين إحدى الملاحظات وملاحظة مجاورة لها؟
  • هل ترتبط إحدى الملاحظات بملاحظة أخرى بعيدة عنها؟ (إلخ)
  • هل السلسلة الزمنية المرصودة عبارة عن ضوضاء بيضاء ؟
  • هل السلسلة الزمنية المرصودة جيبية؟
  • هل السلسلة الزمنية المرصودة ذاتية الانحدار؟
  • ما هو النموذج المناسب للسلسلة الزمنية المرصودة؟
  • هل النموذج
Y=ثابت+خطأ{\displaystyle Y={\text{constant}}+{\text{error}}}
صحيح وكافٍ؟
  • هل الصيغةsY¯=s/شمال{\displaystyle s_{\bar {Y}}=s/{\sqrt {N}}}صالح؟

أهمية

تُعدّ العشوائية (إلى جانب النموذج الثابت، والتباين الثابت، والتوزيع الثابت) أحد الافتراضات الأربعة التي تقوم عليها عادةً جميع عمليات القياس. ويُعتبر افتراض العشوائية بالغ الأهمية للأسباب الثلاثة التالية:

  • تعتمد معظم الاختبارات الإحصائية القياسية على العشوائية. وترتبط صحة نتائج الاختبار ارتباطًا مباشرًا بصحة فرضية العشوائية.
  • تعتمد العديد من الصيغ الإحصائية الشائعة الاستخدام على افتراض العشوائية، وأكثر الصيغ شيوعًا هي صيغة تحديد الخطأ المعياري لمتوسط ​​العينة:
sY¯=s/شمال{\displaystyle s_{\bar {Y}}=s/{\sqrt {N}}}

حيث يمثل s الانحراف المعياري للبيانات. على الرغم من شيوع استخدام هذه الصيغة، إلا أن نتائجها لا قيمة لها ما لم يتحقق افتراض العشوائية.

  • بالنسبة للبيانات أحادية المتغير، يكون النموذج الافتراضي هو
Y=ثابت+خطأ{\displaystyle Y={\text{constant}}+{\text{error}}}

إذا لم تكن البيانات عشوائية، فإن هذا النموذج غير صحيح وغير صالح، وتصبح تقديرات المعلمات (مثل الثابت) غير منطقية وغير صالحة.

تقدير الارتباطات الذاتية

معامل الارتباط الذاتي عند التأخر الزمني h يُعطى بالعلاقة التالية:

رح=جح/ج0{\displaystyle r_{h}=c_{h}/c_{0}\,}

حيث c h هي دالة التغاير الذاتي

جح=1شمالت=1شمال-ح(Yت-Y¯)(Yت+ح-Y¯){\displaystyle c_{h}={\frac {1}{N}}\sum _{t=1}^{Nh}\left(Y_{t}-{\bar {Y}}\right)\left(Y_{t+h}-{\bar {Y}}\right)}

و c 0 هي دالة التباين

ج0=1شمالت=1شمال(Yت-Y¯)2{\displaystyle c_{0}={\frac {1}{N}}\sum _{t=1}^{N}\left(Y_{t}-{\bar {Y}}\right)^{2}}

ستتراوح القيمة الناتجة لـ r h  بين -1 و +1.

تقدير بديل

قد تستخدم بعض المصادر الصيغة التالية لدالة التغاير الذاتي:

جح=1شمال-حت=1شمال-ح(Yت-Y¯)(Yت+ح-Y¯){\displaystyle c_{h}={\frac {1}{Nh}}\sum _{t=1}^{Nh}\left(Y_{t}-{\bar {Y}}\right)\left(Y_{t+h}-{\bar {Y}}\right)}

على الرغم من أن هذا التعريف أقل تحيزًا ، إلا أن صيغة (1/ N ) تتمتع ببعض الخصائص الإحصائية المرغوبة، وهي الصيغة الأكثر شيوعًا في أدبيات الإحصاء. انظر الصفحات 20 و49-50 في كتاب تشاتفيلد لمزيد من التفاصيل.

على عكس التعريف أعلاه، يسمح لنا هذا التعريف بحسابجح{\displaystyle c_{h}}بطريقة أكثر سهولة وبديهية. انظر إلى العينةY1،...،Yشمال{\displaystyle Y_{1},\dots ,Y_{N}}، أينYأناRن{\displaystyle Y_{i}\in \mathbb {R} ^{n}}لأنا=1،...،شمال{\displaystyle i=1,\dots ,N}ثم دع

X=[Y1-Y¯Yشمال-Y¯]Rن×شمال{\displaystyle X={\begin{bmatrix}Y_{1}-{\bar {Y}}&\cdots &Y_{N}-{\bar {Y}}\end{bmatrix}}\in \mathbb {R} ^{n\times N}}

ثم نقوم بحساب مصفوفة غرامسؤال=XX{\displaystyle Q=X^{\top }X}. أخيراً،جح{\displaystyle c_{h}}يتم حسابها كمتوسط ​​العينة لـح{\displaystyle h}القطر الرئيسي لـسؤال{\displaystyle Q}على سبيل المثال،0{\displaystyle 0}القطر الرئيسي (القطر الرئيسي)سؤال{\displaystyle Q}لديهشمال{\displaystyle N}العناصر، ومتوسط ​​العينة يتوافق معج0{\displaystyle c_{0}}. ال1{\displaystyle 1}القطر الأول (إلى يمين القطر الرئيسي) منسؤال{\displaystyle Q}لديهشمال-1{\displaystyle N-1}العناصر، ومتوسط ​​العينة يتوافق معج1{\displaystyle c_{1}}وهكذا دواليك.

الاستدلال الإحصائي باستخدام مخططات الارتباط

مثال على مخطط الارتباط الذاتي من عينة مكونة من 400 نقطة لعملية انحدار ذاتي من الدرجة الأولى بمعامل ارتباط 0.75 بين النقاط المتجاورة، بالإضافة إلى فترات الثقة 95% (مُرسمة حول تقديرات الارتباط باللون الأسود وحول الصفر باللون الأحمر)، كما تم حسابها بواسطة المعادلات الواردة في هذا القسم. يُظهر الخط الأزرق المتقطع دالة الارتباط الذاتي الفعلية للعملية المأخوذة منها العينة.
20 مخطط ارتباط من عينات مكونة من 400 نقطة لنفس العملية العشوائية كما في الشكل السابق.

في الرسم البياني نفسه، يمكن رسم حدود عليا وسفلى للارتباط الذاتي مع مستوى الدلالة.α{\displaystyle \alpha \,}:

ب=±z1-α/2Sهـ(رح){\displaystyle B=\pm z_{1-\alpha /2}SE(r_{h})\,}معرح{\displaystyle r_{h}\,}باعتبارها الارتباط الذاتي المقدر عند التأخيرح{\displaystyle h\,}.

إذا كان الارتباط الذاتي أعلى (أدنى) من هذا الحد الأعلى (الأدنى)، يتم رفض الفرضية الصفرية القائلة بعدم وجود ارتباط ذاتي عند فترة تأخير معينة وما بعدها عند مستوى دلالة قدره α{\displaystyle \alpha \,}هذا الاختبار تقريبي ويفترض أن السلسلة الزمنية غاوسية .

في ما سبق، يمثل z 1− α /2 الكمية المئوية للتوزيع الطبيعي ؛ وSE هو الخطأ المعياري، والذي يمكن حسابه بواسطة صيغة بارتليت لعمليات MA( ):

Sهـ(ر1)=1شمال{\displaystyle SE(r_{1})={\frac {1}{\sqrt {N}}}}
Sهـ(رح)=1+2أنا=1ح-1رأنا2شمال{\displaystyle SE(r_{h})={\sqrt {\frac {1+2\sum _{i=1}^{h-1}r_{i}^{2}}{N}}}}لح>1.{\displaystyle h>1.\,}

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

لاحظ أن هناك صيغتين مختلفتين لتوليد نطاقات الثقة:

1. إذا تم استخدام مخطط الارتباط لاختبار العشوائية (أي، لا يوجد اعتماد زمني في البيانات)، يُوصى باستخدام الصيغة التالية:

±z1-α/2شمال{\displaystyle \pm {\frac {z_{1-\alpha /2}}{\sqrt {N}}}}

حيث N هو حجم العينة ، و z هي دالة الكمية للتوزيع الطبيعي القياسي ، و α هو مستوى الدلالة . في هذه الحالة، تكون نطاقات الثقة ذات عرض ثابت يعتمد على حجم العينة.

2. تُستخدم مخططات الارتباط الذاتي أيضًا في مرحلة تحديد النموذج لتركيب نماذج ARIMA . في هذه الحالة، يُفترض نموذج متوسط ​​متحرك للبيانات، ويجب إنشاء نطاقات الثقة التالية:

±z1-α/21شمال(1+2أنا=1كرأنا2){\displaystyle \pm z_{1-\alpha /2}{\sqrt {{\frac {1}{N}}\left(1+2\sum _{i=1}^{k}r_{i}^{2}\right)}}}

حيث k هو التأخير. في هذه الحالة، تزداد نطاقات الثقة مع ازدياد التأخير.

برمجة

تتوفر مخططات الارتباط الذاتي في معظم المكتبات الإحصائية العامة.

مخططات الارتباط:

كورغرامز:

مراجع

  1. فريندلي، مايكل (19 أغسطس 2002). "مخططات الارتباط: عروض استكشافية لمصفوفات الارتباط" (ملف PDF) . الإحصائي الأمريكي . 56 (4). تايلور وفرانسيس : 316-324 . doi : 10.1198/000313002533 . تاريخ الاسترجاع: 19 يناير 2014 .
  2. 1 2 "CRAN – Package corrgram" . cran.r-project.org . 29 أغسطس 2013. تم الاطلاع عليه في 19 يناير 2014 .
  3. 1 2 "Quick-R: Correlograms" . statmethods.net . تم الاطلاع عليه بتاريخ 19 يناير 2014 .
  4. "1.3.3.1. مخطط الارتباط الذاتي" . www.itl.nist.gov . تم الاطلاع عليه بتاريخ 20 أغسطس 2018 .
  5. "التصور § مخطط الارتباط الذاتي" .

للمزيد من القراءة

  • هانكي، جون إي.؛ ريتش، آرثر جي.؛ ويشرن، دين دبليو. التنبؤ بالأعمال (  الطبعة السابعة). أبر سادل ريفر، نيوجيرسي: برنتيس هول.
  • بوكس، جي إي بي؛ جينكينز، جي. (1976). تحليل السلاسل الزمنية: التنبؤ والتحكم . هولدن-داي.
  • تشاتفيلد، سي. (1989). تحليل السلاسل الزمنية: مقدمة (  الطبعة الرابعة). نيويورك، نيويورك: تشابمان وهول.

المجال العام تتضمن هذه المقالة موادًا متاحة للعموم من المعهد الوطني للمعايير والتكنولوجيا