تحليل المكونات الرئيسية

تحليل المكونات الرئيسية لتوزيع غاوسي متعدد المتغيرات متمركز عند النقطة (1،  3) بانحراف معياري قدره 3 تقريبًا في  الاتجاه (0.866، 0.5) و1 في الاتجاه العمودي. المتجهات الموضحة هي المتجهات الذاتية لمصفوفة التغاير مضروبة في الجذر التربيعي للقيمة الذاتية المقابلة، ومزاحة بحيث تكون ذيولها عند المتوسط.

تحليل المكونات الرئيسية ( PCA ) هو أسلوب خطي لتقليل الأبعاد مع تطبيقات في تحليل البيانات الاستكشافي والتصور ومعالجة البيانات المسبقة .

يتم تحويل البيانات خطيًا إلى نظام إحداثيات جديد بحيث يمكن تحديد الاتجاهات (المكونات الرئيسية) التي تلتقط أكبر تباين في البيانات بسهولة.

المكونات الرئيسية لمجموعة من النقاط في فضاء إحداثيات حقيقي هي سلسلة منص{\displaystyle p}متجهات الوحدة ، حيثأنا{\displaystyle i}المتجه رقم - هو اتجاه الخط الذي يناسب البيانات على أفضل وجه مع كونه متعامدًا مع المتجه الأولأنا-1{\displaystyle i-1}المتجهات. هنا، يُعرَّف الخط الأنسب بأنه الخط الذي يُقلِّل متوسط ​​مربع المسافة العمودية من النقاط إلى الخط . تُشكِّل هذه الاتجاهات (أي المكونات الرئيسية) أساسًا متعامدًا تكون فيه الأبعاد الفردية المختلفة للبيانات غير مترابطة خطيًا . تستخدم العديد من الدراسات أول مكونين رئيسيين لرسم البيانات في بُعدين ولتحديد مجموعات نقاط البيانات المترابطة بصريًا. [ 1 ]

يُستخدم تحليل المكونات الرئيسية في العديد من المجالات مثل علم الوراثة السكانية ، ودراسات الميكروبيوم ، وعلوم الغلاف الجوي . [ 2 ]

ملخص

عند إجراء تحليل المكونات الرئيسية (PCA)، فإن المكون الرئيسي الأول لمجموعة منص{\displaystyle p}المتغير المشتق هو المتغير الناتج عن توليفة خطية من المتغيرات الأصلية، وهو الذي يفسر أكبر قدر من التباين. يفسر المكون الرئيسي الثاني أكبر قدر من التباين المتبقي بعد إزالة تأثير المكون الأول، ويمكننا المتابعة من خلالص{\displaystyle p}تُجرى عمليات التكرار حتى يتم تفسير جميع التباينات. يُستخدم تحليل المكونات الرئيسية (PCA) بشكل شائع عندما تكون العديد من المتغيرات مترابطة ارتباطًا وثيقًا، ويُفضّل تقليل عددها إلى مجموعة مستقلة . يمكن تعريف المكون الرئيسي الأول، بشكل مكافئ، بأنه الاتجاه الذي يزيد من تباين البيانات المُسقطة.أنا{\displaystyle i}يمكن اعتبار المكون الرئيسي رقم -th اتجاهًا متعامدًا مع الأولأنا-1{\displaystyle i-1}المكونات الرئيسية التي تزيد من تباين البيانات المتوقعة.

في كلا الهدفين، يمكن إثبات أن المكونات الرئيسية هي متجهات ذاتية لمصفوفة التغاير للبيانات . لذا، غالبًا ما تُحسب المكونات الرئيسية عن طريق تحليل القيم الذاتية لمصفوفة التغاير أو تحليل القيم المفردة لمصفوفة البيانات. يُعد تحليل المكونات الرئيسية (PCA) أبسط أنواع التحليلات متعددة المتغيرات القائمة على المتجهات الذاتية الحقيقية، ويرتبط ارتباطًا وثيقًا بتحليل العوامل . يتضمن تحليل العوامل عادةً افتراضات أكثر تحديدًا للمجال حول البنية الأساسية، ويحل المتجهات الذاتية لمصفوفة مختلفة قليلًا. يرتبط تحليل المكونات الرئيسية أيضًا بتحليل الارتباط الكنسي (CCA) . يُعرّف تحليل الارتباط الكنسي أنظمة إحداثيات تصف على النحو الأمثل التغاير المشترك بين مجموعتي بيانات، بينما يُعرّف تحليل المكونات الرئيسية نظام إحداثيات متعامدًا جديدًا يصف على النحو الأمثل التباين في مجموعة بيانات واحدة. [ 3 ] [ 4 ] [ 5 ] [ 6 ] كما تم اقتراح متغيرات قوية ومتغيرات قائمة على معيار L1 لتحليل المكونات الرئيسية القياسي. [ 7 ] [ 8 ] [ 9 ] [ 6 ]

تاريخ

تم اختراع PCA في عام 1901 بواسطة كارل بيرسون ، [ 10 ] كنظير لنظرية المحور الرئيسي في الميكانيكا؛ تم تطويره لاحقًا بشكل مستقل وتسميته بواسطة هارولد هوتلينج في الثلاثينيات. [ 11 ] اعتمادًا على مجال التطبيق، يُطلق عليه أيضًا اسم تحويل كارونين-لوف المنفصل (KLT) في معالجة الإشارات ، وتحويل هوتلينغ في مراقبة الجودة متعددة المتغيرات، والتحليل المتعامد المناسب (POD) في الهندسة الميكانيكية، وتحليل القيم المفردة (SVD) لـ X (الذي تم اختراعه في الربع الأخير من القرن التاسع عشر [ 12 ]وتحليل القيم الذاتية (EVD) لـ X في الجبر الخطي، وتحليل العوامل (لمناقشة الاختلافات بين تحليل المكونات الرئيسية (PCA) وتحليل العوامل، انظر الفصل 7 من كتاب جولييف " تحليل المكونات الرئيسية ")، [ 13 ] ونظرية إيكارت-يونغ (هارمان، 1960)، أو الدوال المتعامدة التجريبية (EOF) في علم الأرصاد الجوية (لورنز، 1956)، وتحليل الدوال الذاتية التجريبية (سيروفيتش، 1987)، والأنماط شبه التوافقية. (بروكس وآخرون، 1988)، التحلل الطيفي في الضوضاء والاهتزاز، والتحليل النمطي التجريبي في ديناميكيات الهياكل. 

حدس

مخطط التباين الذي يُساعد في تفسير تحليل المكونات الرئيسية (PCA) وتحديد عدد المكونات التي يجب الاحتفاظ بها. تُشير بداية انحناء الخط (نقطة الانعطاف أو "الركبة") إلى عدد المكونات التي يجب الاحتفاظ بها، وبالتالي في هذا المثال، يجب الاحتفاظ بثلاثة عوامل.

يمكن اعتبار تحليل المكونات الرئيسية (PCA) بمثابة مطابقة قطع ناقص ذي أبعاد p مع البيانات، حيث يمثل كل محور من محاور القطع الناقص مكونًا رئيسيًا. إذا كان أحد محاور القطع الناقص صغيرًا، فإن التباين على طول ذلك المحور يكون صغيرًا أيضًا.

لإيجاد محاور القطع الناقص، يجب أولًا توسيط قيم كل متغير في مجموعة البيانات حول الصفر، وذلك بطرح متوسط ​​القيم المرصودة لكل متغير من تلك القيم. تُستخدم هذه القيم المُحوَّلة بدلًا من القيم المرصودة الأصلية لكل متغير. بعد ذلك، نحسب مصفوفة التغاير للبيانات، ثم نحسب القيم الذاتية والمتجهات الذاتية المقابلة لها. بعد ذلك، يجب تطبيع كل متجه ذاتي متعامد لتحويله إلى متجهات وحدة. بمجرد إتمام ذلك، يمكن تفسير كل متجه ذاتي متعامد على أنه محور للقطع الناقص المُطابق للبيانات. سيؤدي اختيار الأساس هذا إلى تحويل مصفوفة التغاير إلى شكل قطري، حيث تمثل العناصر القطرية تباين كل محور. يمكن حساب نسبة التباين التي يمثلها كل متجه ذاتي عن طريق قسمة القيمة الذاتية المقابلة لهذا المتجه الذاتي على مجموع جميع القيم الذاتية.

تُستخدم الرسوم البيانية الثنائية ورسوم التباين (درجة التباين المُفسَّر ) لتفسير نتائج تحليل المكونات الرئيسية.

تفاصيل

يُعرَّف تحليل المكونات الرئيسية (PCA) بأنه تحويل خطي متعامد على فضاء جداء داخلي حقيقي يحول البيانات إلى نظام إحداثيات جديد بحيث يقع أكبر تباين بواسطة إسقاط قياسي للبيانات على الإحداثي الأول (يسمى المكون الرئيسي الأول)، وثاني أكبر تباين على الإحداثي الثاني، وهكذا. [ 13 ]

فكر فين×ص{\displaystyle n\times p}مصفوفة البيانات ، X ، مع متوسط ​​تجريبي صفري عموديًا (تم تحويل متوسط ​​العينة لكل عمود إلى الصفر)، حيث يمثل كل صف من الصفوف n تكرارًا مختلفًا للتجربة، ويعطي كل عمود من الأعمدة p نوعًا معينًا من الميزات (على سبيل المثال، النتائج من مستشعر معين).

رياضياً، يتم تعريف التحويل بواسطة مجموعة من الحجمل{\displaystyle l}(أينل{\displaystyle l}يتم اختيارها عادة لتكون أقل منص{\displaystyle p}لتقليل الأبعاد) منص{\displaystyle p}متجهات ذات أبعاد للأوزان أو المعاملاتw(ك)=(w1،...،wص)(ك){\displaystyle \mathbf {w} _{(k)}=(w_{1},\dots ,w_{p})_{(k)}}التي تحدد كل متجه صفx(أنا)=(x1،...،xص)(أنا){\displaystyle \mathbf {x} _{(i)}=(x_{1},\dots ,x_{p})_{(i)}}تحويل X إلى متجه جديد من درجات المكونات الرئيسيةت(أنا)=(ت1،...،تل)(أنا){\displaystyle \mathbf {t} _{(i)}=(t_{1},\dots ,t_{l})_{(i)}}، مقدمة من

تك(أنا)=x(أنا)w(ك)وoرأنا=1،...،نك=1،...،ل{\displaystyle {t_{k}}_{(i)}=\mathbf {x} _{(i)}\cdot \mathbf {w} _{(k)}\qquad \mathrm {for} \qquad i=1,\dots ,n\qquad k=1,\dots ,l}

بطريقة تجعل المتغيرات الفرديةت1،...،تل{\displaystyle t_{1},\dots ,t_{l}}يتم اعتبار t على مجموعة البيانات بشكل متتابع ، حيث يرثون أقصى تباين ممكن من X ، مع تقييد كل متجه معامل w ليكون متجه وحدة .

يمكن كتابة ما سبق بشكل مكافئ في شكل مصفوفة كما يلي:

تي=Xدبليو{\displaystyle \mathbf {T} =\mathbf {X} \mathbf {W} }

أين تيأناك=تك(أنا){\displaystyle {\mathbf {T} }_{ik}={t_{k}}_{(i)}}، Xأناج=xج(أنا){\displaystyle {\mathbf {X} }_{ij}={x_{j}}_{(i)}}، و دبليوجك=wج(ك){\displaystyle {\mathbf {W} } _ {jk} = {w_ {j}} _ {(ك)}}.

المكون الأول

لتحقيق أقصى قدر من التباين، يجب أن يحقق متجه الوزن الأول w (1) الشرط التالي:

w(1)=argالأعلىw=1{أنا(ت1)(أنا)2}=argالأعلىw=1{أنا(x(أنا)w)2}// _ {i}\left(\mathbf {x} _{(i)}\cdot \mathbf {w} \right)^{2}\right\}}

وبصورة مكافئة، فإن كتابة هذا في شكل مصفوفة يعطي

w(1)=argالأعلىw=1{Xw2}=argالأعلىw=1{wتيXتيXw}{\displaystyle \mathbf {w} _{(1)}=\arg \max _{\left\|\mathbf {w} \right\|=1}\left\{\left\|\mathbf {Xw} \right\|^{2}\right\}=\arg \max _{\left\|\mathbf {w} \right\|=1}\left\{\mathbf {w} ^{\mathsf {T}}\mathbf {X} ^{\mathsf {T}}\mathbf {Xw} \right\}}

بما أن w (1) قد تم تعريفه على أنه متجه وحدة، فإنه يحقق أيضاً بشكل مكافئ ما يلي:

w(1)=argالأعلى{wتيXتيXwwتيw}{\displaystyle \mathbf {w} _{(1)}=\arg \max \left\{{\frac {\mathbf {w} ^{\mathsf {T}}\mathbf {X} ^{\mathsf {T}}\mathbf {Xw} }{\mathbf {w} ^{\mathsf {T}}\mathbf {w} }}\right\}}

يمكن التعرف على الكمية المراد تعظيمها على أنها حاصل قسمة رايلي . ومن النتائج الشائعة لمصفوفة شبه موجبة مثل X T X أن القيمة القصوى الممكنة لحاصل القسمة هي أكبر قيمة ذاتية للمصفوفة، والتي تتحقق عندما يكون w هو المتجه الذاتي المقابل .

بعد إيجاد w (1) ، يمكن بعد ذلك إعطاء المكون الرئيسي الأول لمتجه البيانات x ( i ) كدرجة t 1( i ) = x ( i )w (1) في الإحداثيات المحولة، أو كمتجه مقابل في المتغيرات الأصلية، { x ( i )w (1) } w (1) .

مكونات إضافية

يمكن إيجاد المكون رقم k عن طريق طرح المكونات الرئيسية k  1 الأولى من X :

X^ك=X-s=1ك-1Xw(s)w(s)تي{\displaystyle \mathbf {\hat {X}} _{k}=\mathbf {X} -\sum _{s=1}^{k-1}\mathbf {X} \mathbf {w} _{(s)}\mathbf {w} _{(s)}^{\mathsf {T}}}

ثم إيجاد متجه الوزن الذي يستخرج أكبر قدر من التباين من مصفوفة البيانات الجديدة هذه

w(ك)=أرزمأxw=1{X^كw2}=argالأعلى{wتيX^كتيX^كwwتيw}{\displaystyle \mathbf {w} _{(k)}=\mathop {\operatorname {arg\,max} } _{\left\|\mathbf {w} \right\|=1}\left\{\left\|\mathbf {\hat {X}} _{k}\mathbf {w} \right\|^{2}\right\}=\arg \max \left\{{\tfrac {\mathbf {w} ^{\mathsf {T}}\mathbf {\hat {X}} _{k}^{\mathsf {T}}\mathbf {\hat {X}} _{k}\mathbf {w} }{\mathbf {w} ^{T}\mathbf {w} }}\right\}}

يتضح أن هذا يُعطي المتجهات الذاتية المتبقية للمصفوفة X T X ، حيث تُعطى القيم القصوى للكمية بين قوسين بقيمها الذاتية المقابلة. وبالتالي، فإن متجهات الوزن هي متجهات ذاتية للمصفوفة X T X.

وبالتالي يمكن إعطاء المكون الرئيسي رقم k لمتجه البيانات x ( i ) كدرجة t k ( i ) = x ( i )w ( k ) في الإحداثيات المحولة، أو كمتجه مقابل في فضاء المتغيرات الأصلية، { x ( i )w ( k ) } w ( k ) ، حيث w ( k ) هو المتجه الذاتي رقم k لـ X T X .

وبالتالي، يمكن التعبير عن تحليل المكونات الرئيسية الكامل لـ X على النحو التالي:

تي=Xدبليو{\displaystyle \mathbf {T} =\mathbf {X} \mathbf {W} }

حيث W مصفوفة أوزان من الرتبة p × p ، أعمدتها هي المتجهات الذاتية للمصفوفة X. يُطلق على منقولة W أحيانًا اسم تحويل التبييض أو تحويل التكوير . تُسمى أعمدة W مضروبة في الجذر التربيعي للقيم الذاتية المناظرة، أي المتجهات الذاتية مضروبة في التباينات، معاملات التحميل في تحليل المكونات الرئيسية (PCA) أو في تحليل العوامل.

التغايرات

يمكن اعتبار X T X نفسها متناسبة مع مصفوفة التغاير التجريبية لمجموعة البيانات X. [ 13 ] : 30-31

يُعطى التباين المشترك Q بين اثنين من المكونات الرئيسية المختلفة على مجموعة البيانات بالصيغة التالية:

سؤال(Pج(ج)،Pج(ك))(Xw(ج))تي(Xw(ك))=w(ج)تيXتيXw(ك)=w(ج)تيλ(ك)w(ك)=λ(ك)w(ج)تيw(ك){\displaystyle {\begin{aligned}Q(\mathrm {PC} _{(j)},\mathrm {PC} _{(k)})&\propto (\mathbf {X} \mathbf {w} _{(j)})^{\mathsf {T}}(\mathbf {X} \mathbf {w} _{(k)})\\&=\mathbf {w} _{(j)}^{\mathsf {T}}\mathbf {X} ^{\mathsf {T}}\mathbf {X} \mathbf {w} _{(k)}\\&=\mathbf {w} _{(j)}^{\mathsf {T}}\lambda _{(k)}\mathbf {w} _{(k)}\\&=\lambda _{(k)}\mathbf {w} _{(j)}^{\mathsf {T}}\mathbf {w} _{(k)}\end{aligned}}}

حيث استُخدمت خاصية القيمة الذاتية للمتجه w ( k ) للانتقال من السطر 2 إلى السطر 3. مع ذلك، فإن المتجهين الذاتيين w ( j ) و w ( k ) المناظرين للقيم الذاتية لمصفوفة متناظرة يكونان متعامدين (إذا كانت القيم الذاتية مختلفة)، أو يمكن جعلهما متعامدين (إذا تصادف أن يشترك المتجهان في قيمة متكررة متساوية). وبالتالي، فإن حاصل الضرب في السطر الأخير يساوي صفرًا؛ إذ لا يوجد تباين مشترك بين المكونات الرئيسية المختلفة في مجموعة البيانات.

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

يمكن كتابة مصفوفة التغاير التجريبية للمتغيرات الأصلية في شكل مصفوفة

سؤالXتيX=دبليوΛدبليوتي{\displaystyle \mathbf {Q} \propto \mathbf {X} ^{\mathsf {T}}\mathbf {X} =\mathbf {W} \mathbf {\Lambda } \mathbf {W} ^{\mathsf {T}}}

تصبح مصفوفة التغاير التجريبية بين المكونات الرئيسية

دبليوتيسؤالدبليودبليوتيدبليوΛدبليوتيدبليو=Λ{\displaystyle \mathbf {W} ^{\mathsf {T}}\mathbf {Q} \mathbf {W} \propto \mathbf {W} ^{\mathsf {T}}\mathbf {W} \,\mathbf {\Lambda } \,\mathbf {W} ^{\mathsf {T}}\mathbf {W} =\mathbf {\Lambda } }

حيث Λ هي المصفوفة القطرية للقيم الذاتية λ ( k ) لـ X T X . λ ( k ) يساوي مجموع المربعات على مجموعة البيانات المرتبطة بكل مكون k ، أي أن λ ( k ) = Σ i t k 2 ( i ) = Σ i ( x ( i )w ( k ) ) 2 .

تقليل الأبعاد

التحويل P = X W يُحوّل متجه البيانات x ( i ) من فضاء أصلي مكون من x متغيرات إلى فضاء جديد مكون من p متغيرات غير مترابطة عبر مجموعة البيانات. ولجعل البيانات المركزية غير بُعدية، لنفترض أن X<sub> c</sub> يُمثل القيم المميزة لمتجهات البيانات X<sub> i</sub> ، كما يلي:

  • X{\displaystyle \|X\|_{\infty }}(المعيار الأقصى)،
  • 1نX1{\displaystyle {\frac {1}{n}}\|X\|_{1}}(المتوسط ​​المطلق)، أو
  • 1نX2{\displaystyle {\frac {1}{\sqrt {n}}}\|X\|_{2}}(المعيار الإقليدي المعياري)،

بالنسبة لمجموعة بيانات بحجم n ، تُستخدم هذه المعايير لتحويل الفضاء الأصلي للمتغيرات x و y إلى فضاء جديد من المتغيرات غير المرتبطة p و q (مع إعطاء Y c بنفس المعنى)، بحيثصأنا=XأناXج،qأنا=YأناYج{\displaystyle p_{i}={\frac {X_{i}}{X_{c}}},\quad q_{i}={\frac {Y_{i}}{Y_{c}}}}والمتغيرات الجديدة مرتبطة خطيًا كما يلي:q=αص{\displaystyle q=\alpha p}لإيجاد العلاقة الخطية المثلى، نقوم بتقليل إجمالي مربع خطأ إعادة البناء: هـ(α)=11-α2أنا=1ن(αصأنا-qأنا)2{\displaystyle E(\alpha )={\frac {1}{1-\alpha ^{2}}}\sum _{i=1}^{n}(\alpha p_{i}-q_{i})^{2}}بحيث يؤدي ضبط مشتقة دالة الخطأ إلى الصفر(هـ(α)=0){\displaystyle (E'(\alpha )=0)}ينتج عنه:α=12(-λ±λ2+4){\displaystyle \alpha ={\frac {1}{2}}\left(-\lambda \pm {\sqrt {\lambda ^{2}+4}}\right)}أينλ=صص-qqصq{\displaystyle \lambda ={\frac {p\cdot p-q\cdot q}{p\cdot q}}}[ 14 ]

مخطط التشتت لتحليل المكونات الرئيسية لأنماط هابلوتايب Y-STR المحسوبة من قيم عدد التكرارات لـ 37 علامة STR على الكروموسوم Y من 354 فردًا. وقد نجح تحليل المكونات الرئيسية في إيجاد تراكيب خطية للعلامات تفصل بين مجموعات مختلفة تتوافق مع سلالات مختلفة من النسب الجيني للكروموسوم Y للأفراد.

يُعدّ تقليل الأبعاد خطوةً بالغة الأهمية لتصوّر ومعالجة مجموعات البيانات عالية الأبعاد، مع الحفاظ على أكبر قدر ممكن من التباين فيها. فعلى سبيل المثال، يُحدّد اختيار L  =  2 والاحتفاظ بالمكوّنين الرئيسيين الأولين فقط المستوى ثنائي الأبعاد الذي يمرّ عبر مجموعة البيانات عالية الأبعاد، حيث تكون البيانات أكثر تباعدًا. لذا، إذا احتوت البيانات على تجمعات ، فقد تكون هذه التجمعات أيضًا أكثر تباعدًا، وبالتالي أكثر وضوحًا عند تمثيلها بيانيًا في مخطط ثنائي الأبعاد. أما إذا تم اختيار اتجاهين عشوائيًا عبر البيانات (أو متغيرين من المتغيرات الأصلية)، فقد تكون التجمعات أقل تباعدًا عن بعضها، بل وقد تتداخل بشكل كبير، ما يجعل تمييزها صعبًا.

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

قد يكون تقليل الأبعاد مناسبًا أيضًا عندما تكون المتغيرات في مجموعة البيانات مشوشة. إذا احتوى كل عمود من أعمدة مجموعة البيانات على ضوضاء غاوسية مستقلة ومتطابقة التوزيع، فإن أعمدة المصفوفة T ستحتوي أيضًا على ضوضاء غاوسية متطابقة التوزيع (هذا التوزيع ثابت تحت تأثير المصفوفة W ، والتي يمكن اعتبارها دورانًا عالي الأبعاد لمحاور الإحداثيات). مع ذلك، نظرًا لتركز جزء أكبر من التباين الكلي في المكونات الرئيسية القليلة الأولى مقارنةً بتباين الضوضاء نفسه، فإن التأثير النسبي للضوضاء يكون أقل - إذ تحقق المكونات القليلة الأولى نسبة إشارة إلى ضوضاء أعلى . وبالتالي، يمكن أن يؤدي تحليل المكونات الرئيسية (PCA) إلى تركيز جزء كبير من الإشارة في المكونات الرئيسية القليلة الأولى، والتي يمكن استخلاصها بسهولة عن طريق تقليل الأبعاد؛ بينما قد تهيمن الضوضاء على المكونات الرئيسية اللاحقة، وبالتالي يمكن التخلص منها دون خسارة كبيرة. إذا لم تكن مجموعة البيانات كبيرة جدًا، فيمكن اختبار أهمية المكونات الرئيسية باستخدام إعادة التوزيع المعياري (parametric bootstrap )، كوسيلة مساعدة في تحديد عدد المكونات الرئيسية التي يجب الاحتفاظ بها. [ 15 ]

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

يمكن أيضًا ربط تحويل المكونات الرئيسية بتحليل مصفوفة آخر، وهو تحليل القيم المفردة (SVD) للمصفوفة X.

X=يوΣدبليوتي{\displaystyle \mathbf {X} =\mathbf {U} \mathbf {\Sigma } \mathbf {W} ^{T}}

هنا Σ هي مصفوفة قطرية مستطيلة من الرتبة n × p من الأعداد الموجبة σ ( k ) ، تسمى القيم المفردة لـ X ؛ U هي مصفوفة من الرتبة n × n ، أعمدتها عبارة عن متجهات وحدة متعامدة بطول n تسمى المتجهات المفردة اليسرى لـ X ؛ و W هي مصفوفة من الرتبة p × p أعمدتها عبارة عن متجهات وحدة متعامدة بطول p وتسمى المتجهات المفردة اليمنى لـ X.

من حيث هذا التحليل، يمكن كتابة المصفوفة X T X على النحو التالي

XتيX=دبليوΣتييوتييوΣدبليوتي=دبليوΣتيΣدبليوتي=دبليوΣ^2دبليوتي{\displaystyle {\begin{aligned}\mathbf {X} ^{T}\mathbf {X} &=\mathbf {W} \mathbf {\Sigma } ^{\mathsf {T}}\mathbf {U} ^{\mathsf {T}}\mathbf {U} \mathbf {\Sigma } \mathbf {W} ^{\mathsf {T}}\\&=\mathbf {W} \mathbf {\Sigma } ^{\mathsf {T}}\mathbf {\Sigma } \mathbf {W} ^{\mathsf {T}}\\&=\mathbf {W} \mathbf {\hat {\Sigma }} ^{2}\mathbf {W} ^{\mathsf {T}}\end{aligned}}}

أينΣ^{\displaystyle \mathbf {\hat {\Sigma }} }هي المصفوفة القطرية المربعة التي تحتوي على القيم المفردة لـ X والأصفار الزائدة التي تم حذفها والتي تحققΣ^2=ΣتيΣ{\displaystyle \mathbf {{\hat {\Sigma }}^{2}} =\mathbf {\Sigma } ^{\mathsf {T}}\mathbf {\Sigma } }تُظهر المقارنة مع تحليل المتجهات الذاتية لـ X T X أن المتجهات المفردة اليمنى W لـ X تُكافئ المتجهات الذاتية لـ X T X ، بينما القيم المفردة σ ( k ) لـX{\displaystyle \mathbf {X} }تساوي الجذر التربيعي للقيم الذاتية λ ( k ) لـ X T X .

باستخدام تحليل القيم المفردة، يمكن كتابة مصفوفة النتائج T

تي=Xدبليو=يوΣدبليوتيدبليو=يوΣ{\displaystyle {\begin{aligned}\mathbf {T} &=\mathbf {X} \mathbf {W} \\&=\mathbf {U} \mathbf {\Sigma } \mathbf {W} ^{\mathsf {T}}\mathbf {W} \\&=\mathbf {U} \mathbf {\Sigma } \end{aligned}}}

إذن ، يُعطى كل عمود من T بأحد المتجهات المفردة اليسرى لـ X مضروبًا في القيمة المفردة المناظرة. هذا الشكل هو أيضًا التحليل القطبي لـ T.

توجد خوارزميات فعالة لحساب SVD لـ X دون الحاجة إلى تكوين المصفوفة X T X ، لذا فإن حساب SVD هو الآن الطريقة القياسية لحساب تحليل المكونات الرئيسية من مصفوفة البيانات ، [ 16 ] ما لم تكن هناك حاجة إلى عدد قليل من المكونات فقط.

كما هو الحال مع تحليل القيم الذاتية، يمكن الحصول على مصفوفة درجات مقطوعة بحجم n × T من خلال النظر فقط في أول L قيمة مفردة كبيرة ومتجهاتها المفردة:

تيل=يولΣل=Xدبليول{\displaystyle \mathbf {T} _{L}=\mathbf {U} _{L}\mathbf {\Sigma } _{L}=\mathbf {X} \mathbf {W} _{L}}

يؤدي اقتطاع المصفوفة M أو T باستخدام تحليل القيمة المفردة المقتطعة بهذه الطريقة إلى إنتاج مصفوفة مقتطعة هي أقرب مصفوفة ممكنة من الرتبة L إلى المصفوفة الأصلية، بمعنى أن الفرق بين الاثنين له أصغر معيار فروبينيوس ممكن ، وهي نتيجة تُعرف باسم نظرية إيكارت-يونغ [1936].

نظرية (الملاءمة المثلى ذات الأبعاد k). لتكن P مصفوفة بيانات n×m تم توسيط أعمدتها وتوحيد مقياسها، ولتكن P=يوΣVتي{\displaystyle P=U\,\Sigma \,V^{T}} ليكن تحليلها إلى قيم مفردة. عندئذٍ يكون أفضل تقريب من الرتبة k لـ P بمعنى المربعات الصغرى (معيار فروبينيوس) هو Pك=يوكΣكVكتي{\displaystyle P_{k}=U_{k}\,\Sigma _{k}\,V_{k}^{T}}حيث تتكون V k من الأعمدة k الأولى من V. علاوة على ذلك، فإن التباين المتبقي النسبي هو R(ك)=ج=ك+1مσج2ج=1مσج2{\displaystyle R(k)={\frac {\sum _{j=k+1}^{m}\sigma _{j}^{2}}{\sum _{j=1}^{m}\sigma _{j}^{2}}}}.

[ 14 ]

اعتبارات إضافية

القيم المفردة (في Σ ) هي الجذور التربيعية للقيم الذاتية للمصفوفة X T X. تتناسب كل قيمة ذاتية مع جزء من "التباين" (أو بالأحرى، مجموع مربعات المسافات بين النقاط ومتوسطها متعدد الأبعاد) المرتبط بكل متجه ذاتي. مجموع جميع القيم الذاتية يساوي مجموع مربعات المسافات بين النقاط ومتوسطها متعدد الأبعاد. تقوم تقنية تحليل المكونات الرئيسية (PCA) أساسًا بتدوير مجموعة النقاط حول متوسطها لتتوافق مع المكونات الرئيسية. ينقل هذا أكبر قدر ممكن من التباين (باستخدام تحويل متعامد) إلى الأبعاد القليلة الأولى. لذلك، تميل القيم في الأبعاد المتبقية إلى أن تكون صغيرة، ويمكن حذفها مع الحد الأدنى من فقدان المعلومات (انظر أدناه ). غالبًا ما تُستخدم تقنية PCA بهذه الطريقة لتقليل الأبعاد . تتميز تقنية PCA بكونها التحويل المتعامد الأمثل للاحتفاظ بالفضاء الفرعي الذي يحتوي على أكبر "تباين" (كما هو مُعرّف أعلاه). إلا أن هذه الميزة تأتي على حساب متطلبات حسابية أكبر مقارنةً، على سبيل المثال، وعند الاقتضاء، بتحويل جيب التمام المنفصل ، وخاصةً تحويل جيب التمام المنفصل من النوع الثاني (DCT-II) المعروف ببساطة باسم "DCT". وتميل تقنيات تقليل الأبعاد غير الخطية إلى أن تكون أكثر تطلبًا من الناحية الحسابية من تحليل المكونات الرئيسية (PCA).

تتأثر تقنية تحليل المكونات الرئيسية (PCA) بتغيير مقياس المتغيرات. رياضياً، تنشأ هذه الحساسية من الطريقة التي يغير بها تغيير المقياس مصفوفة التباين المشترك للعينة التي تقوم تقنية PCA بقطريتها. [ 14 ]

يتركXج{\displaystyle \mathbf {X} _{\text{c}}}لتكن مصفوفة البيانات *المتمركزة* ( عدد الصفوف n ، عدد الأعمدة p ) وحدد التباين المشترك Σ=1نXجتيXج.{\displaystyle \Sigma ={\frac {1}{n}}\,\mathbf {X} _{\text{c}}^{\mathsf {T}}\mathbf {X} _{\text{c}}.} إذاج{\displaystyle j}يتم ضرب المتغير رقم ‑ بعاملαج{\displaystyle \alpha _{j}}نحصل على Xج(α)=Xجد،د=التشخيص(α1،...،αص).{\displaystyle \mathbf {X} _{\text{c}}^{(\alpha )}=\mathbf {X} _{\text{c}}D,\qquad D=\operatorname {diag} (\alpha _{1},\ldots ,\alpha _{p}).} وبالتالي فإن التغاير الجديد هو Σ(α)=دتيΣد.{\displaystyle \Sigma ^{(\alpha )}=D^{\mathsf {T}}\,\Sigma \,D.}

لأن القيم الذاتية والمتجهات الذاتية لـΣ(α){\displaystyle \Sigma ^{(\alpha )}}هي تلك الخاصة بـΣ{\displaystyle \Sigma }تم قياسه بواسطةد{\displaystyle D}، تدور المحاور الرئيسية باتجاه أي عمود تم تضخيم تباينه، تمامًا كما يوضح المثال ثنائي الأبعاد أدناه.

إذا كان لدينا متغيران فقط، ولهما نفس تباين العينة ، ومرتبطان ارتباطًا تامًا، فإن تحليل المكونات الرئيسية (PCA) يستلزم تدويرًا بزاوية 45 درجة، وستكون "الأوزان" (وهي جيوب تمام زاوية التدوير) للمتغيرين بالنسبة للمكون الرئيسي متساوية. ولكن إذا ضربنا جميع قيم المتغير الأول في 100، فسيكون المكون الرئيسي الأول مطابقًا تقريبًا لهذا المتغير، مع مساهمة طفيفة من المتغير الآخر، بينما سيكون المكون الثاني متوافقًا تقريبًا مع المتغير الأصلي الثاني. هذا يعني أنه عندما تكون وحدات المتغيرات المختلفة مختلفة (مثل درجة الحرارة والكتلة)، فإن تحليل المكونات الرئيسية (PCA) يُعد طريقة تحليل اعتباطية إلى حد ما. (ستُحصل على نتائج مختلفة إذا استخدمنا فهرنهايت بدلًا من مئوية، على سبيل المثال). كان عنوان ورقة بيرسون الأصلية "حول الخطوط والمستويات الأقرب ملاءمة لأنظمة النقاط في الفضاء" - و"في الفضاء" يشير إلى الفضاء الإقليدي الفيزيائي حيث لا تظهر مثل هذه الاعتبارات. إحدى طرق جعل تحليل المكونات الرئيسية أقل تعسفًا هي استخدام متغيرات مُقاسة بحيث يكون تباينها واحدًا، وذلك بتوحيد البيانات، وبالتالي استخدام مصفوفة الارتباط الذاتي بدلًا من مصفوفة التغاير الذاتي كأساس لتحليل المكونات الرئيسية. مع ذلك، فإن هذا يضغط (أو يوسع) التقلبات في جميع أبعاد فضاء الإشارة إلى تباين واحد.

يفترض تحليل المكونات الرئيسية الكلاسيكي أن سحابة النقاط قد تم نقلها بالفعل بحيث يكون مركزها عند نقطة الأصل. [ 14 ]

اكتب كل ملاحظة على النحو التالي: qأنا=μ+zأنا،μ=1نأنا=1نqأنا.{\displaystyle \mathbf {q} _{i}={\boldsymbol {\mu }}+\mathbf {z} _{i},\qquad {\boldsymbol {\mu }}={\tfrac {1}{n}}\sum _{i=1}^{n}\mathbf {q} _{i}.}

بدون طرحμ{\displaystyle {\boldsymbol {\mu }}}نحن في الواقع نقوم بالقطر

Σغير معروف=نμμتي+1نZتيZ،{\displaystyle \Sigma _{\text{unc}}\;=\;n\,{\boldsymbol {\mu }}{\boldsymbol {\mu }}^{\mathsf {T}}\;+\;{\tfrac {1}{n}}\,\mathbf {Z} ^{\mathsf {T}}\mathbf {Z} ,}

أينZ{\displaystyle \mathbf {Z} }هي المصفوفة المركزية. الحد ذو الرتبة الأولىنμμتي{\displaystyle n\,{\boldsymbol {\mu }}{\boldsymbol {\mu }}^{\mathsf {T}}}غالباً ما يهيمن، مما يجبر المتجه الذاتي الرئيسي على الإشارة بدقة تقريباً نحو المتوسط، ويمحو أي بنية في الجزء المركزي.Z{\displaystyle \mathbf {Z} }بعد طرح المتوسط، يختفي هذا الحد وتتوافق المحاور الرئيسية مع الاتجاهات الحقيقية للتباين الأقصى.

لا حاجة لتوسيط البيانات عند إجراء تحليل المكونات الرئيسية على مصفوفة الارتباط، لأن البيانات تكون مُوَسَّطة بالفعل بعد حساب الارتباطات. تُستخلص الارتباطات من حاصل ضرب قيمتين معياريتين (الدرجات Z) أو عزمين إحصائيين (ومن هنا جاء اسم: معامل ارتباط بيرسون ). انظر أيضًا مقالة كرومري وفوستر-جونسون (1998) بعنوان "توسيط البيانات في الانحدار المُعدَّل: ضجة بلا طائل" . بما أن التغايرات هي ارتباطات بين متغيرات مُعَيَّرة ( الدرجات Z أو الدرجات المعيارية ) ، فإن تحليل المكونات الرئيسية المُستند إلى مصفوفة ارتباط X يُعادل تحليل المكونات الرئيسية المُستند إلى مصفوفة التغاير لـ Z ، وهي النسخة المعيارية من X.

يُعد تحليل المكونات الرئيسية (PCA) تقنية أساسية شائعة في التعرف على الأنماط . مع ذلك، فهو غير مُحسَّن لفصل الفئات. [ 17 ] ومع ذلك، فقد استُخدم لتحديد المسافة بين فئتين أو أكثر من خلال حساب مركز الكتلة لكل فئة في فضاء المكونات الرئيسية، ثم الإبلاغ عن المسافة الإقليدية بين مراكز كتل الفئتين أو أكثر. [ 18 ] يُعد تحليل التمييز الخطي بديلاً مُحسَّنًا لفصل الفئات.

جدول الرموز والاختصارات

رمزمعنىأبعادالمؤشرات
X=[Xأناج]{\displaystyle \mathbf {X} =[X_{ij}]}مصفوفة البيانات، التي تتكون من مجموعة جميع متجهات البيانات، متجه واحد لكل صفن×ص{\displaystyle n\times p}أنا=1...ن{\displaystyle i=1\ldots n}ج=1...ص{\displaystyle j=1\ldots p}
ن{\displaystyle n}عدد متجهات الصفوف في مجموعة البيانات1×1{\displaystyle 1\times 1}كمية قياسية
ص{\displaystyle p}عدد العناصر في كل متجه صف (بعد)1×1{\displaystyle 1\times 1}كمية قياسية
ل{\displaystyle L}عدد الأبعاد في الفضاء الفرعي ذي الأبعاد المخفضة،1لص{\displaystyle 1\leq L\leq p}1×1{\displaystyle 1\times 1}كمية قياسية
u=[uج]{\displaystyle \mathbf {u} =[u_{j}]}متجه المتوسطات التجريبية ، متوسط ​​واحد لكل عمود j من مصفوفة البياناتص×1{\displaystyle p\times 1}ج=1...ص{\displaystyle j=1\ldots p}
s=[sج]{\displaystyle \mathbf {s} =[s_{j}]}متجه الانحرافات المعيارية التجريبية ، انحراف معياري واحد لكل عمود j من مصفوفة البياناتص×1{\displaystyle p\times 1}ج=1...ص{\displaystyle j=1\ldots p}
ح=[حأنا]{\displaystyle \mathbf {h} =[h_{i}]}متجه مكون من جميع القيم 11×ن{\displaystyle 1\times n}أنا=1...ن{\displaystyle i=1\ldots n}
ب=[بأناج]{\displaystyle \mathbf {B} =[B_{ij}]}الانحرافات عن متوسط ​​كل عمود j من مصفوفة البياناتن×ص{\displaystyle n\times p}أنا=1...ن{\displaystyle i=1\ldots n}ج=1...ص{\displaystyle j=1\ldots p}
Z=[Zأناج]{\displaystyle \mathbf {Z} =[Z_{ij}]}يتم حساب قيم z باستخدام المتوسط ​​والانحراف المعياري لكل عمود j من مصفوفة البيانات.ن×ص{\displaystyle n\times p}أنا=1...ن{\displaystyle i=1\ldots n}ج=1...ص{\displaystyle j=1\ldots p}
ج=[ججج]{\displaystyle \mathbf {C} =[C_{jj'}]}مصفوفة التغايرص×ص{\displaystyle p\times p}ج=1...ص{\displaystyle j=1\ldots p}ج=1...ص{\displaystyle j'=1\ldots p}
R=[Rجج]{\displaystyle \mathbf {R} =[R_{jj'}]}مصفوفة الارتباطص×ص{\displaystyle p\times p}ج=1...ص{\displaystyle j=1\ldots p}ج=1...ص{\displaystyle j'=1\ldots p}
V=[Vجج]{\displaystyle \mathbf {V} =[V_{jj'}]}مصفوفة تتكون من مجموعة جميع المتجهات الذاتية لـ C ، متجه ذاتي واحد لكل عمودص×ص{\displaystyle p\times p}ج=1...ص{\displaystyle j=1\ldots p}ج=1...ص{\displaystyle j'=1\ldots p}
د=[دجج]{\displaystyle \mathbf {D} =[D_{jj'}]}مصفوفة قطرية تتكون من مجموعة جميع القيم الذاتية للمصفوفة C على طول قطرها الرئيسي ، و0 لجميع العناصر الأخرى (ملاحظة).Λ{\displaystyle \mathbf {\Lambda } }(مستخدم أعلاه)ص×ص{\displaystyle p\times p}ج=1...ص{\displaystyle j=1\ldots p}ج=1...ص{\displaystyle j'=1\ldots p}
دبليو=[دبليوجل]{\displaystyle \mathbf {W} =[W_{jl}]}مصفوفة من متجهات الأساس، متجه واحد لكل عمود، حيث يكون كل متجه أساس أحد المتجهات الذاتية لـ C ، وحيث تكون المتجهات في W مجموعة جزئية من تلك الموجودة في Vص×ل{\displaystyle p\times L}ج=1...ص{\displaystyle j=1\ldots p}ل=1...ل{\displaystyle l=1\ldots L}
تي=[تيأنال]{\displaystyle \mathbf {T} =[T_{il}]}مصفوفة تتكون من n متجهات صفية، حيث يمثل كل متجه إسقاط متجه البيانات المقابل من المصفوفة X على متجهات الأساس الموجودة في أعمدة المصفوفة W.ن×ل{\displaystyle n\times L}أنا=1...ن{\displaystyle i=1\ldots n}ل=1...ل{\displaystyle l=1\ldots L}

الخصائص والقيود

ملكيات

تتضمن بعض خصائص تحليل المكونات الرئيسية ما يلي: [ 13 ]

الخاصية 1 : لأي عدد صحيحq، 1 ≤qpالتحويل الخطيالمتعامد
y=بx{\displaystyle y=\mathbf {B'} x}
أينy{\displaystyle y}هو متجه مكون من q عنصر وب{\displaystyle \mathbf {B'} }هي مصفوفة ( q × p )، ولتكنΣy=بΣب{\displaystyle \mathbf {\Sigma } _{y}=\mathbf {B'} \mathbf {\Sigma } \mathbf {B} }لتكن مصفوفة التباين والتباين المشترك لـy{\displaystyle y}ثم أثرΣy{\displaystyle \mathbf {\Sigma } _{y}}، المشار إليهtr(Σy){\displaystyle \operatorname {tr} (\mathbf {\Sigma } _{y})}، يتم تحقيق أقصى قدر من خلال أخذب=أq{\displaystyle \mathbf {B} =\mathbf {A} _{q}}، أينأq{\displaystyle \mathbf {A} _{q}}يتكون من أول q عمود منأ{\displaystyle \mathbf {A} }(ب{\displaystyle (\mathbf {B'} }هو منقولب){\displaystyle \mathbf {B} )}. (أ{\displaystyle \mathbf {A} }(غير مُعرَّف هنا)
الخاصية 2 : لننظر مرة أخرى فيالتحويل المتعامد المعياري
y=بx{\displaystyle y=\mathbf {B'} x}
معx،ب،أ{\displaystyle x,\mathbf {B} ,\mathbf {A} }وΣy{\displaystyle \mathbf {\Sigma } _{y}}كما هو موضح سابقاً. ثمtr(Σy){\displaystyle \operatorname {tr} (\mathbf {\Sigma } _{y})}يتم تقليله عن طريق الأخذب=أq*،{\displaystyle \mathbf {B} =\mathbf {A} _{q}^{*},}أينأq*{\displaystyle \mathbf {A} _{q}^{*}}يتكون من آخر q عمود منأ{\displaystyle \mathbf {A} }.

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

الخاصية 3 : (التحليل الطيفي لـ Σ )
Σ=λ1α1α1++λصαصαص{\displaystyle \mathbf {\Sigma } =\lambda _{1}\alpha _{1}\alpha _{1}'+\cdots +\lambda _{p}\alpha _{p}\alpha _{p}'}

قبل أن نتناول استخدامها، سنلقي نظرة أولاً على العناصر القطرية ،

متغير(xج)=ك=1Pλكαكج2{\displaystyle \operatorname {Var} (x_{j})=\sum _{k=1}^{P}\lambda _{k}\alpha _{kj}^{2}}

لعلّ أهمّ دلالة إحصائية لهذه النتيجة هي أنه لا يمكننا فقط تحليل التباينات المُجمّعة لجميع عناصر x إلى مساهمات متناقصة ناتجة عن كل مُكوّن رئيسي، بل يمكننا أيضًا تحليل مصفوفة التغاير بأكملها إلى مساهمات.λكαكαك{\displaystyle \lambda _{k}\alpha _{k}\alpha _{k}'}من كل جهاز كمبيوتر. على الرغم من أنها ليست متناقصة تمامًا، إلا أن عناصرλكαكαك{\displaystyle \lambda _{k}\alpha _{k}\alpha _{k}'}سيميل إلى أن يصبح أصغر كلماك{\displaystyle k}يزداد، معλكαكαك{\displaystyle \lambda _{k}\alpha _{k}\alpha _{k}'}غير متزايدة بالنسبة للزيادةك{\displaystyle k}بينما عناصرαك{\displaystyle \alpha _{k}}تميل إلى البقاء بنفس الحجم تقريبًا بسبب قيود التطبيع:αكαك=1،ك=1،...،ص{\displaystyle \alpha _{k}'\alpha _{k}=1,k=1,\dots ,p}.

القيود

كما ذُكر أعلاه، تعتمد نتائج تحليل المكونات الرئيسية على مقياس المتغيرات. ويمكن معالجة ذلك عن طريق قياس كل خاصية بانحرافها المعياري، بحيث نحصل في النهاية على خصائص عديمة الأبعاد ذات تباين موحد. [ 19 ]

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

من القيود الأخرى عملية إزالة المتوسط ​​قبل بناء مصفوفة التغاير لتحليل المكونات الرئيسية (PCA). في مجالات مثل علم الفلك، تكون جميع الإشارات غير سالبة، وستؤدي عملية إزالة المتوسط ​​إلى جعل متوسط ​​بعض التعريضات الفلكية يساوي صفرًا، مما ينتج عنه تدفقات سالبة غير منطقية [ 21 ] ، ولذا يجب إجراء نمذجة أمامية لاستعادة القيمة الحقيقية للإشارات [ 22 ] . كبديل، يُعد تحليل المصفوفات غير السالبة ، الذي يركز فقط على العناصر غير السالبة في المصفوفات ، مناسبًا تمامًا للرصد الفلكي [23] [24] [25]. للمزيد ، انظر العلاقة بين تحليل المكونات الرئيسية وتحليل المصفوفات غير السالبة .

تُعاني تقنية تحليل المكونات الرئيسية (PCA) من قصورٍ إذا لم يتم توحيد البيانات قبل تطبيق الخوارزمية عليها. إذ تُحوّل PCA البيانات الأصلية إلى بياناتٍ ذات صلةٍ بالمكونات الرئيسية لتلك البيانات، مما يعني أن متغيرات البيانات الجديدة لا يُمكن تفسيرها بنفس طريقة تفسير المتغيرات الأصلية، بل هي تفسيراتٌ خطيةٌ لها. كما أن عدم تطبيق PCA بشكلٍ صحيحٍ يزيد من احتمالية فقدان المعلومات. [ 26 ]

تعتمد تقنية تحليل المكونات الرئيسية (PCA) على نموذج خطي. إذا احتوت مجموعة البيانات على نمط غير خطي كامن، فقد تُؤدي هذه التقنية إلى توجيه التحليل في الاتجاه المعاكس تمامًا. [ 27 ] اكتشف باحثون في جامعة ولاية كانساس أن خطأ المعاينة في تجاربهم يؤثر على تحيز نتائج تحليل المكونات الرئيسية. "إذا كان عدد العينات أو المجموعات أقل من 30، و/أو كان الباحث مهتمًا بالمكونات الرئيسية التي تتجاوز المكون الأول، فقد يكون من الأفضل تصحيح الارتباط التسلسلي قبل إجراء تحليل المكونات الرئيسية". [ 28 ] كما وجد الباحثون في جامعة ولاية كانساس أن تحليل المكونات الرئيسية قد يكون "متحيزًا بشكل كبير إذا لم يتم التعامل مع بنية الارتباط الذاتي للبيانات بشكل صحيح". [ 28 ]

تحليل المكونات الرئيسية ونظرية المعلومات

يؤدي تقليل الأبعاد إلى فقدان المعلومات بشكل عام. ويميل تقليل الأبعاد القائم على تحليل المكونات الرئيسية إلى تقليل فقدان المعلومات هذا، في ظل نماذج إشارة وضوضاء معينة.

بافتراض أن

x=s+ن،{\displaystyle \mathbf {x} =\mathbf {s} +\mathbf {n} ,}

أي أن متجه البياناتx{\displaystyle \mathbf {x} }هو مجموع الإشارة المطلوبة الحاملة للمعلوماتs{\displaystyle \mathbf {s} }وإشارة ضوضاءن{\displaystyle \mathbf {n} }يمكن للمرء أن يثبت أن تحليل المكونات الرئيسية (PCA) يمكن أن يكون الأمثل لتقليل الأبعاد، من وجهة نظر نظرية المعلومات.

وعلى وجه الخصوص، أظهر لينسكر أنه إذاs{\displaystyle \mathbf {s} }هو توزيع غاوسي ون{\displaystyle \mathbf {n} }إذا كان الضجيج غاوسيًا بمصفوفة تباين تتناسب مع مصفوفة الوحدة، فإن تحليل المكونات الرئيسية (PCA) يزيد من المعلومات المتبادلة.أنا(y؛s){\displaystyle I(\mathbf {y} ;\mathbf {s} )} بين المعلومات المطلوبةs{\displaystyle \mathbf {s} }والناتج ذو الأبعاد المخفضةy=دبليولتيx{\displaystyle \mathbf {y} =\mathbf {W} _{L}^{T}\mathbf {x} }[ 29 ]

إذا كانت الضوضاء لا تزال غاوسية ولها مصفوفة تباين تتناسب مع مصفوفة الوحدة (أي مكونات المتجه)ن{\displaystyle \mathbf {n} }هي مستقلة ومتطابقة التوزيع)، ولكن الإشارة الحاملة للمعلوماتs{\displaystyle \mathbf {s} }إذا كان غير غاوسي (وهو سيناريو شائع)، فإن PCA على الأقل يقلل من الحد الأعلى لفقدان المعلومات ، والذي يتم تعريفه على النحو التالي [ 30 ] [ 31 ]

أنا(x؛s)-أنا(y؛s).{\displaystyle I(\mathbf {x} ;\mathbf {s} )-I(\mathbf {y}  ;\mathbf {s} ).}

يتم الحفاظ على مثالية تحليل المكونات الرئيسية أيضًا إذا كان التشويشن{\displaystyle \mathbf {n} }هي مستقلة ومتطابقة التوزيع، وعلى الأقل أكثر غاوسية (من حيث تباعد كولباك-لايبير ) من الإشارة الحاملة للمعلوماتs{\displaystyle \mathbf {s} }[ 32 ] بشكل عام ، حتى لو كان نموذج الإشارة المذكور أعلاه صحيحًا، فإن تحليل المكونات الرئيسية يفقد مثاليته من الناحية النظرية للمعلومات بمجرد ظهور الضوضاءن{\displaystyle \mathbf {n} }يصبح معتمداً على غيره.

الحساب باستخدام طريقة التغاير

فيما يلي وصف تفصيلي لتحليل المكونات الرئيسية باستخدام طريقة التغاير [ 33 ] بدلاً من طريقة الارتباط [ 34 ] .

الهدف هو تحويل مجموعة بيانات معينة X ذات بُعد p إلى مجموعة بيانات بديلة Y ذات بُعد أصغر L. وبصورة مكافئة، نسعى لإيجاد المصفوفة Y ، حيث Y هي تحويل كارونين-لوف (KLT ) للمصفوفة X.

Y=كلتي{X}{\displaystyle \mathbf {Y} =\mathbb {KLT} \{\mathbf {X} \}}

  1. قم بتنظيم مجموعة البيانات

    لنفترض أن لديك بيانات تتألف من مجموعة من المشاهدات لـ p متغير، وتريد تقليل البيانات بحيث يمكن وصف كل مشاهدة باستخدام L متغير فقط، حيث L < p . ولنفترض كذلك أن البيانات مرتبة كمجموعة من n متجه بيانات.x1...xن{\displaystyle \mathbf {x} _{1}\ldots \mathbf {x} _{n}}مع كلxأنا{\displaystyle \mathbf {x} _{i}}يمثل هذا ملاحظة واحدة مجمعة لمتغيرات p .

    • يكتبx1...xن{\displaystyle \mathbf {x} _{1}\ldots \mathbf {x} _{n}}كمتجهات صفية، كل منها يحتوي على p عنصر.
    • ضع متجهات الصف في مصفوفة واحدة X ذات أبعاد n × p .
  2. احسب المتوسط ​​التجريبي
    • أوجد المتوسط ​​التجريبي على طول كل عمود j = 1،  ...، p . 
    • ضع القيم المتوسطة المحسوبة في متجه متوسط ​​تجريبي u ذي أبعاد p × 1.uج=1نأنا=1نXأناج{\displaystyle u_{j}={\frac {1}{n}}\sum _{i=1}^{n}X_{ij}}
  3. احسب الانحرافات عن المتوسط

    يُعدّ طرح المتوسط ​​جزءًا لا يتجزأ من حلّ إيجاد أساس للمكونات الرئيسية يُقلّل من متوسط ​​مربع الخطأ في تقريب البيانات. [ 35 ] لذا، نتابع بتوسيط البيانات على النحو التالي:

    • اطرح متجه المتوسط ​​التجريبيuتي{\displaystyle \mathbf {u} ^{T}}من كل صف من صفوف مصفوفة البيانات X.
    • قم بتخزين البيانات بعد طرح المتوسط ​​في المصفوفة B ذات الأبعاد n × p .ب=X-حuتي{\displaystyle \mathbf {B} =\mathbf {X} -\mathbf {h} \mathbf {u} ^{T}}حيث h عبارة عن متجه عمودي n × 1 جميع  قيمه 1:حأنا=1ل أنا=1،...،ن{\displaystyle h_{i}=1\,\qquad \qquad {\text{for }}i=1,\ldots ,n}

    في بعض التطبيقات، يمكن أيضًا تعديل مقياس كل متغير (عمود من B ) بحيث يكون تباينه مساويًا لـ 1 (انظر درجة Z ). [ 36 ] تؤثر هذه الخطوة على المكونات الرئيسية المحسوبة، لكنها تجعلها مستقلة عن الوحدات المستخدمة لقياس المتغيرات المختلفة.

  4. أوجد مصفوفة التغاير
    • أوجد مصفوفة التغاير التجريبية C من الرتبة p × p من المصفوفة B :ج=1ن-1ب*ب{\displaystyle \mathbf {C} ={1 \over {n-1}}\mathbf {B} ^{*}\mathbf {B} }أين*{\displaystyle *}هو عامل النقل المترافق . إذا كانت المجموعة B تتكون بالكامل من أعداد حقيقية، وهو الحال في العديد من التطبيقات، فإن "النقل المترافق" هو ​​نفسه النقل العادي .
    • إن المنطق وراء استخدام n − 1 بدلاً من n لحساب التغاير هو تصحيح بيسل .
  5. أوجد المتجهات الذاتية والقيم الذاتية لمصفوفة التغاير
    • احسب مصفوفة V للمتجهات الذاتية التي تُقطّر مصفوفة التغاير C :V-1جV=د{\displaystyle \mathbf {V} ^{-1}\mathbf {C} \mathbf {V} =\mathbf {D} }حيث D هي المصفوفة القطرية للقيم الذاتية للمصفوفة C. تتضمن هذه الخطوة عادةً استخدام خوارزمية حاسوبية لحساب المتجهات الذاتية والقيم الذاتية . تتوفر هذه الخوارزميات بسهولة كمكونات فرعية لمعظم أنظمة جبر المصفوفات ، مثل SAS [ 37 ] ، وR ، و MATLAB [ 38 ] [ 39 ] ، وMathematica [ 40 ] ، و SciPy ، و IDL ( لغة البيانات التفاعلية )، و GNU Octave ، بالإضافة إلى OpenCV .
    • ستأخذ المصفوفة D شكل مصفوفة قطرية من الرتبة p × p ، حيثدك=λكل ك={\displaystyle D_{k\ell }=\lambda _{k}\qquad {\text{for }}k=\ell }يمثل القيمة الذاتية رقم j لمصفوفة التغاير C ، ودك=0ل ك.{\displaystyle D_{k\ell }=0\qquad {\text{for }}k\neq \ell .}
    • تحتوي المصفوفة V ، ذات البعد p × p أيضًا ، على p متجهات عمودية، كل منها بطول p ، والتي تمثل المتجهات الذاتية p لمصفوفة التغاير C.
    • يتم ترتيب القيم الذاتية والمتجهات الذاتية وتكوين أزواج منها. القيمة الذاتية رقم j تتوافق مع المتجه الذاتي رقم j .
    • يرمز المصفوفة V إلى مصفوفة المتجهات الذاتية اليمنى (مقابل المتجهات الذاتية اليسرى ). وبشكل عام، لا يشترط أن تكون مصفوفة المتجهات الذاتية اليمنى هي منقولة (مرافقة) مصفوفة المتجهات الذاتية اليسرى.
  6. أعد ترتيب المتجهات الذاتية والقيم الذاتية
    • قم بترتيب أعمدة مصفوفة المتجهات الذاتية V ومصفوفة القيم الذاتية D بترتيب تنازلي للقيمة الذاتية.
    • تأكد من الحفاظ على الاقترانات الصحيحة بين الأعمدة في كل مصفوفة.
  7. احسب محتوى الطاقة التراكمي لكل متجه ذاتي
    • تمثل القيم الذاتية توزيع طاقة بيانات المصدر بين كل متجه ذاتي، حيث تشكل المتجهات الذاتية أساسًا للبيانات. محتوى الطاقة التراكمي g للمتجه الذاتي j هو مجموع محتوى الطاقة عبر جميع القيم الذاتية من 1 إلى j مقسومًا على مجموع محتوى الطاقة عبر جميع القيم الذاتية (كما هو موضح في الخطوة 8).زج=ك=1جدككل ج=1،...،ص{\displaystyle g_{j}=\sum _{k=1}^{j}D_{kk}\qquad {\text{for }}j=1,\dots ,p}
  8. اختر مجموعة فرعية من المتجهات الذاتية كمتجهات أساسية
    • احفظ أول L عمود من المصفوفة V كمصفوفة W من الرتبة p × L :دبليوكل=Vكل ك=1،...،ص=1،...،ل{\displaystyle W_{kl}=V_{k\ell }\qquad {\text{for }}k=1,\dots ,p\qquad \ell =1,\dots ,L}أين1لص.{\displaystyle 1\leq L\leq p.}
    • استخدم المتجه g كدليل لاختيار قيمة مناسبة لـ L. الهدف هو اختيار قيمة L أصغر ما يمكن مع تحقيق قيمة عالية نسبيًا لـ g كنسبة مئوية. على سبيل المثال، قد ترغب في اختيار L بحيث تكون الطاقة التراكمية g أعلى من عتبة معينة، مثل 90%. في هذه الحالة، اختر أصغر قيمة لـ L بحيثزلزص0.9{\displaystyle {\frac {g_{L}}{g_{p}}}\geq 0.9}
  9. قم بإسقاط البيانات على الأساس الجديد
    • تمثل نقاط البيانات المسقطة صفوف المصفوفةتي=بدبليو{\displaystyle \mathbf {T} =\mathbf {B} \cdot \mathbf {W} }

    أي العمود الأول منتي{\displaystyle \mathbf {T} }يمثل العمود الأول إسقاط نقاط البيانات على المكون الرئيسي الأول، ويمثل العمود الثاني الإسقاط على المكون الرئيسي الثاني، وهكذا.

الاشتقاق باستخدام طريقة التغاير

ليكن X متجهًا عشوائيًا ذو بُعد d معبرًا عنه كمتجه عمودي. وبدون فقدان للعمومية، نفترض أن متوسط ​​X يساوي صفرًا.

نريد أن نجد(*){\displaystyle (\ast )}مصفوفة تحويل متعامدة من الرتبة d × d P بحيث يكون لـ PX مصفوفة تباين قطرية (أي أن PX عبارة عن متجه عشوائي مع جميع مكوناته المتميزة غير المرتبطة بشكل زوجي).

حساب سريع بافتراضP{\displaystyle P}كانت العوائد أحادية:

كوف(PX)=هـ[PX (PX)*]=هـ[PX X*P*]=Pهـ[XX*]P*=Pكوف(X)P-1{\displaystyle {\begin{aligned}\operatorname {cov} (PX)&=\operatorname {E} [PX~(PX)^{*}]\\&=\operatorname {E} [PX~X^{*}P^{*}]\\&=P\operatorname {E} [XX^{*}]P^{*}\\&=P\operatorname {cov} (X)P^{-1}\\\end{aligned}}}

لذلك(*){\displaystyle (\ast )}يتحقق إذا وفقط إذاكوف(X){\displaystyle \operatorname {cov} (X)}كانت قابلة للقطر بواسطةP{\displaystyle P}.

هذا أمر بنّاء للغاية، حيث أن cov( X ) مضمون أن يكون مصفوفة غير سالبة محددة، وبالتالي مضمون أن يكون قابلاً للتقطير بواسطة مصفوفة وحدوية ما.

الحساب بدون تباين مشترك

في التطبيقات العملية، وخاصةً مع البيانات عالية الأبعاد ( p كبيرة )، نادرًا ما تُستخدم طريقة التغاير البسيطة لعدم كفاءتها نظرًا لتكاليفها الحسابية والذاكرية العالية لتحديد مصفوفة التغاير بشكل صريح. يتجنب النهج غير المعتمد على التغاير np² عملية حساب وتخزين مصفوفة التغاير XᵀX ، ويستخدم بدلاً من ذلك إحدى الطرق غير المعتمدة على المصفوفات ، على سبيل المثال، تلك القائمة على دالة تقييم حاصل ضرب Xᵀ ( Xᵀ ) بتكلفة 2np عملية .

الحساب التكراري

إحدى طرق حساب المكون الرئيسي الأول بكفاءة [ 41 ] موضحة في الكود الزائف التالي، لمصفوفة بيانات X ذات متوسط ​​صفري، دون حساب مصفوفة التباين الخاصة بها.

r = متجه عشوائي طوله p r = r / norm( r ) كرر العملية c مرة: s = 0 (متجه طوله p ) لكل صف x في X s = s + ( xr ) x λ = r T s // λ هي القيمة الذاتية الخطأ = |λ ⋅ rs | r = s / norm( s ) اخرج إذا كان الخطأ < التسامح أرجع λ، r

تقوم خوارزمية التكرار الأسي هذه ببساطة بحساب المتجه X T (X r) ، ثم تُعَيِّره، وتُعيد النتيجة إلى r . تُقَرِّب القيمة الذاتية بواسطة r T (X T X) r ، وهي حاصل قسمة رايلي على متجه الوحدة r لمصفوفة التغاير X T X. إذا كانت أكبر قيمة مفردة منفصلة جيدًا عن القيمة التي تليها في الحجم، فإن المتجه r يقترب من المكون الرئيسي الأول لـ X خلال عدد التكرارات c ، وهو عدد صغير نسبيًا مقارنةً بـ p ، بتكلفة إجمالية قدرها 2cnp . يمكن تسريع تقارب التكرار الأسي دون التضحية بشكل ملحوظ بالتكلفة المنخفضة لكل تكرار باستخدام طرق أكثر تطورًا لا تعتمد على المصفوفات ، مثل خوارزمية لانكزوس أو طريقة التدرج المترافق المُهيأ مسبقًا الأمثل محليًا ( LOBPCG ).

يمكن حساب المكونات الرئيسية اللاحقة بشكل فردي باستخدام تقنية الانكماش، أو بشكل متزامن كوحدة واحدة. في الطريقة الأولى، تؤثر عدم الدقة في المكونات الرئيسية التقريبية المحسوبة مسبقًا بشكل تراكمي على دقة المكونات الرئيسية المحسوبة لاحقًا، مما يزيد الخطأ مع كل عملية حساب جديدة. أما الطريقة الثانية، وهي طريقة القوة الكتلية، فتستبدل المتجهين المفردين r و s بمتجهات كتلية، وهي المصفوفات R و S. يُقارب كل عمود من R أحد المكونات الرئيسية الرائدة، بينما تُكرر جميع الأعمدة في آن واحد. وتتمثل العملية الحسابية الرئيسية في حساب حاصل ضرب X T (XR) . عند تطبيق هذه الطريقة، كما في LOBPCG ، يُزيل التجميع الفعال تراكم الأخطاء، ويسمح باستخدام دوال ضرب المصفوفات عالية المستوى من مكتبة BLAS ، ويؤدي عادةً إلى تقارب أسرع مقارنةً بتقنية المتجه المفرد.

طريقة NIPALS

تُعدّ طريقة المربعات الصغرى الجزئية التكرارية غير الخطية (NIPALS) صيغةً مُعدّلة من طريقة التكرار الأسي الكلاسيكية مع تقليص المصفوفة بالطرح، وتُستخدم لحساب المكونات القليلة الأولى في تحليل المكونات الرئيسية أو تحليل المربعات الصغرى الجزئية . بالنسبة لمجموعات البيانات عالية الأبعاد، مثل تلك المُولّدة في علوم *omics (على سبيل المثال، علم الجينوم ، علم الأيض )، عادةً ما يكون من الضروري حساب المكونات الرئيسية القليلة الأولى فقط. تُحدّث خوارزمية المربعات الصغرى الجزئية التكرارية غير الخطية (NIPALS) التقريبات التكرارية للدرجات والأوزان الرئيسية t₁ و r₁T من خلال تكرار أسي، حيث يتم الضرب في كل تكرار بالمتغير X من اليسار واليمين، أي يتم تجنب حساب مصفوفة التغاير، تمامًا كما هو الحال في تطبيق التكرار الأسي الخالي من المصفوفات للمتغير XᵀX ، استنادًا إلى الدالة التي تُقيّم حاصل الضرب Xᵀ (Xᵀr) = ( (Xᵀr) ᵀX ) .

يتم إجراء عملية تقليص المصفوفة بالطرح عن طريق طرح حاصل الضرب الخارجي، t 1 r 1 من مما ينتج عنه مصفوفة البواقي المُقلّصة التي تُستخدم لحساب المكونات الرئيسية الرئيسية اللاحقة. [ 42 ] بالنسبة لمصفوفات البيانات الكبيرة، أو المصفوفات التي تتميز بدرجة عالية من الارتباط الخطي بين الأعمدة، يعاني برنامج NIPALS من فقدان تعامد المكونات الرئيسية نتيجة لأخطاء التقريب الناتجة عن دقة الآلة والمتراكمة في كل تكرار، بالإضافة إلى تقليص المصفوفة بالطرح. [ 43 ] وللتخلص من هذا الفقدان للتعامد، يتم تطبيق خوارزمية غرام -شميدت لإعادة التعامد على كل من الدرجات والأوزان في كل خطوة تكرارية. [ 44 ] لا يمكن لـ NIPALS الاعتماد على عمليات الضرب أحادية المتجه للاستفادة من BLAS عالي المستوى ويؤدي إلى تقارب بطيء للقيم المفردة الرائدة المجمعة - يتم حل كلا هذين القصورين في حلول الكتل الأكثر تطورا الخالية من المصفوفات، مثل طريقة التدرج المترافق المشروط مسبقًا للكتل الأمثل محليًا ( LOBPCG ).

التقدير عبر الإنترنت/التسلسلي

في حالة "البث المباشر" أو "التدفق المستمر" حيث تصل البيانات على دفعات بدلاً من تخزينها في دفعة واحدة، من المفيد تقدير إسقاط تحليل المكونات الرئيسية (PCA) الذي يمكن تحديثه بشكل متسلسل. يمكن القيام بذلك بكفاءة، ولكنه يتطلب خوارزميات مختلفة. [ 45 ]

المتغيرات النوعية

في تحليل البيانات الاستكشافي، يمكن إسقاط المتغيرات النوعية (الفئوية) غير المستخدمة في بناء النموذج على محاور المكونات الرئيسية كعناصر تكميلية. [ 46 ] [ 47 ] [ 48 ] وكان استخدام المتغيرات النوعية التكميلية سمة مميزة للمدرسة الفرنسية لتحليل البيانات التي تطورت في الستينيات والسبعينيات. [ 49 ] واستنادًا إلى عمل لودوفيك ليبارت ، تتوفر هذه الطريقة في بيئة R من خلال حزم برمجية مثل FactoMineR. [ 50 ] [ 51 ] [ 52 ]

التطبيقات

ذكاء

كان أول تطبيق لتحليل العوامل في تحديد وقياس مكونات الذكاء البشري. ساد الاعتقاد بأن الذكاء يتألف من مكونات متعددة غير مترابطة، مثل الذكاء المكاني، والذكاء اللفظي، والاستقراء، والاستنتاج، وغيرها، وأن درجات هذه المكونات يمكن استخلاصها من نتائج اختبارات مختلفة باستخدام تحليل العوامل، للحصول على مؤشر واحد يُعرف باسم حاصل الذكاء (IQ). وقد طوّر عالم النفس الإحصائي الرائد، سبيرمان، تحليل العوامل في عام 1904 لنظريته الثنائية للذكاء، مضيفًا بذلك أسلوبًا رسميًا إلى علم القياس النفسي . وفي عام 1924، بحث ثورستون عن 56 عاملًا للذكاء، مطورًا مفهوم العمر العقلي. وتستند اختبارات حاصل الذكاء القياسية اليوم إلى هذا العمل المبكر. [ 53 ]

التمايز السكني

في عام ١٩٤٩، قدّم شيفكي وويليامز نظرية البيئة العاملية ، التي هيمنت على دراسات التمايز السكني من خمسينيات القرن العشرين إلى سبعينياته. [ ٥٤ ] كانت الأحياء في المدينة قابلة للتمييز عن بعضها البعض من خلال خصائص مختلفة يمكن اختزالها إلى ثلاثة باستخدام تحليل العوامل. عُرفت هذه الخصائص باسم "الرتبة الاجتماعية" (مؤشر على الوضع الوظيفي)، و"الترابط الأسري" أو حجم الأسرة، و"العرق". بعد ذلك، يمكن تطبيق تحليل التجميع لتقسيم المدينة إلى مجموعات أو مناطق وفقًا لقيم متغيرات العوامل الرئيسية الثلاثة. وقد تطورت أدبيات واسعة النطاق حول البيئة العاملية في الجغرافيا الحضرية، لكن هذا النهج تراجع بعد عام ١٩٨٠ لكونه بدائيًا من الناحية المنهجية ولا مكان له في النماذج الجغرافية ما بعد الحداثية.

لطالما تمثّلت إحدى مشكلات تحليل العوامل في إيجاد أسماء مقنعة للعوامل الاصطناعية المختلفة. في عام 2000، أعاد فلود إحياء منهج علم البيئة العاملي ليُبيّن أن تحليل المكونات الرئيسية يُقدّم إجابات ذات مغزى بشكل مباشر، دون اللجوء إلى تدوير العوامل. كانت المكونات الرئيسية في الواقع متغيرات ثنائية أو أسعارًا ضمنية لـ"قوى" تدفع الناس إلى التجمع أو التباعد في المدن. كان المكون الأول هو "إمكانية الوصول"، وهو المفاضلة الكلاسيكية بين الطلب على السفر والطلب على المساحة، والتي يقوم عليها الاقتصاد الحضري الكلاسيكي. أما المكونان التاليان فهما "الحرمان"، الذي يُبقي الأشخاص ذوي الوضع الاجتماعي المماثل في أحياء منفصلة (بوساطة التخطيط)، والعرق، حيث يسعى الأشخاص ذوو الخلفيات العرقية المتشابهة إلى التواجد في نفس المكان. [ 55 ]

في نفس الفترة تقريبًا، حدد مكتب الإحصاء الأسترالي مؤشرات متميزة للميزة والعيب، معتمدًا على المكون الرئيسي الأول لمجموعات من المتغيرات الأساسية التي اعتُبرت مهمة. تُنشر هذه المؤشرات (SEIFA) بانتظام لمختلف المناطق القضائية، وتُستخدم بكثرة في التحليل المكاني. [ 56 ]

مؤشرات التنمية

يمكن استخدام تحليل المكونات الرئيسية (PCA) كطريقة رسمية لتطوير المؤشرات. وكبديل لذلك، تم اقتراح تحليل المكونات التأكيدي لتطوير المؤشرات وتقييمها. [ 57 ]

طُوِّر مؤشر تنمية المدن من قِبَل شركة PCA استنادًا إلى حوالي 200 مؤشر لنتائج المدن، وذلك في مسح أُجري عام 1996 وشمل 254 مدينة عالمية. خضع المكون الرئيسي الأول لتحليل الانحدار التكراري، حيث أُضيفت المتغيرات الأصلية بشكل منفرد حتى تم تفسير حوالي 90% من تباينه. استخدم المؤشر في نهاية المطاف حوالي 15 مؤشرًا، ولكنه كان مؤشرًا جيدًا للعديد من المتغيرات الأخرى. توافقت قيمته المقارنة بشكل كبير مع التقييم الذاتي لحالة كل مدينة. كانت معاملات بنود البنية التحتية متناسبة تقريبًا مع متوسط ​​تكاليف توفير الخدمات الأساسية، مما يشير إلى أن المؤشر كان في الواقع مقياسًا للاستثمار المادي والاجتماعي الفعال في المدينة.

إن مؤشر التنمية البشرية على مستوى الدولة (HDI) من برنامج الأمم المتحدة الإنمائي ، والذي تم نشره منذ عام 1990 ويستخدم على نطاق واسع في دراسات التنمية، [ 58 ] له معاملات متشابهة للغاية على مؤشرات مماثلة، مما يشير بقوة إلى أنه تم إنشاؤه في الأصل باستخدام PCA.

علم الوراثة السكانية

في عام 1978، كان كافالي-سفورزا وآخرون من الرواد في استخدام تحليل المكونات الرئيسية (PCA) لتلخيص البيانات المتعلقة بتغير ترددات الجينات البشرية عبر المناطق. وقد أظهرت المكونات أنماطًا مميزة، بما في ذلك التدرجات والموجات الجيبية. وفسروا هذه الأنماط على أنها ناتجة عن أحداث هجرة قديمة محددة.

منذ ذلك الحين، أصبح تحليل المكونات الرئيسية (PCA) شائع الاستخدام في علم الوراثة السكانية، حيث استخدمته آلاف الأبحاث كآلية عرض. تختلف الوراثة بشكل كبير تبعًا للقرب الجغرافي، لذا فإن أول مكونين رئيسيين يُظهران التوزيع المكاني، ويمكن استخدامهما لرسم خريطة الموقع الجغرافي النسبي لمجموعات سكانية مختلفة، مما يُظهر الأفراد الذين هاجروا من مواقعهم الأصلية. [ 59 ]

لطالما أثارت تقنية تحليل المكونات الرئيسية (PCA) في علم الوراثة جدلاً واسعاً من الناحية التقنية، إذ طُبقت على متغيرات منفصلة غير طبيعية، وغالباً على علامات أليلية ثنائية. كما يُعدّ غياب أي مقاييس للخطأ المعياري في هذه التقنية عائقاً أمام استخدامها بشكل أكثر اتساقاً. في أغسطس 2022، نشر عالم الأحياء الجزيئية إران إلحايك ورقة بحثية نظرية في مجلة "ساينتيفيك ريبورتس" تُحلل 12 تطبيقاً لتقنية PCA. وخلص إلى سهولة التلاعب بهذه الطريقة، مما أدى، في رأيه، إلى نتائج "خاطئة ومتناقضة وغير منطقية". وعلى وجه التحديد، جادل بأن النتائج التي تم التوصل إليها في علم الوراثة السكانية اتسمت بانتقاء البيانات بشكل انتقائي والاستدلال الدائري . [ 60 ]

أبحاث السوق ومؤشرات المواقف

لطالما كان البحث التسويقي من أكثر مستخدمي تحليل المكونات الرئيسية (PCA). يُستخدم هذا التحليل لتطوير مؤشرات رضا العملاء أو ولائهم للمنتجات، وبالتزامن مع التجميع، يُستخدم لتحديد شرائح السوق التي يمكن استهدافها بحملات إعلانية، على غرار ما يفعله علم البيئة العاملي لتحديد المناطق الجغرافية ذات الخصائص المتشابهة. [ 61 ]

تُحوّل تقنية تحليل المكونات الرئيسية (PCA) كميات هائلة من البيانات بسرعة إلى متغيرات أصغر وأسهل فهمًا، ما يُتيح تحليلها بشكل أسرع وأسهل. في أي استبيان للمستهلكين، توجد سلسلة من الأسئلة المصممة لاستخلاص آراء المستهلكين، وتبحث المكونات الرئيسية عن المتغيرات الكامنة وراء هذه الآراء. على سبيل المثال، استطلعت دراسة أكسفورد للإنترنت عام 2013 آراء 2000 شخص حول آرائهم ومعتقداتهم، ومن خلال هذه الأسئلة، استخلص المحللون أربعة أبعاد للمكونات الرئيسية، حددوها على أنها "الهروب"، و"التواصل الاجتماعي"، و"الكفاءة"، و"خلق المشكلات". [ 62 ]

في مثال آخر، استخلص جو فلود في عام 2008 مؤشرًا للاتجاهات نحو السكن من 28 سؤالًا حول الاتجاهات في مسح وطني شمل 2697 أسرة في أستراليا. مثّل المكون الرئيسي الأول اتجاهًا عامًا نحو العقارات وامتلاك المنازل. ويمكن إدخال هذا المؤشر، أو أسئلة الاتجاهات التي يتضمنها، في نموذج خطي عام لاختيار نوع السكن. وكان المؤشر هو المحدد الأقوى للاستئجار الخاص، متجاوزًا الدخل أو الحالة الاجتماعية أو نوع الأسرة. [ 63 ]

التمويل الكمي

في التمويل الكمي ، يتم استخدام PCA [ 64 ] في إدارة المخاطر المالية ، وقد تم تطبيقه على مشاكل أخرى مثل تحسين المحفظة .

يُستخدم تحليل المكونات الرئيسية (PCA) بشكل شائع في المشكلات المتعلقة بالأوراق المالية ذات الدخل الثابت والمحافظ الاستثمارية ، ومشتقات أسعار الفائدة . تعتمد التقييمات هنا على منحنى العائد بأكمله ، والذي يتضمن العديد من الأدوات المترابطة ارتباطًا وثيقًا، ويُستخدم تحليل المكونات الرئيسية لتحديد مجموعة من المكونات أو العوامل التي تُفسر تحركات أسعار الفائدة، [ 65 ] مما يُسهل عملية النمذجة. أحد التطبيقات الشائعة لإدارة المخاطر هو حساب القيمة المعرضة للخطر (VaR)، وذلك بتطبيق تحليل المكونات الرئيسية على محاكاة مونت كارلو . [ 66 ] هنا، لكل عينة محاكاة، يتم اختبار المكونات، ثم تُعاد بناء أسعار الفائدة، وبالتالي قيم الخيارات ؛ مع حساب القيمة المعرضة للخطر في النهاية على مدار عملية المحاكاة بأكملها. يُستخدم تحليل المكونات الرئيسية أيضًا في التحوط من التعرض لمخاطر أسعار الفائدة ، مع الأخذ في الاعتبار المدد الجزئية والحساسيات الأخرى. [ 65 ] في كلتا الحالتين، تُعد المكونات الرئيسية الثلاثة الأولى للنظام، والتي تُمثل عادةً "التحول" و"الالتواء" و"التقوس"، ذات أهمية. يتم اشتقاق هذه المكونات الرئيسية من تحليل القيم الذاتية لمصفوفة التغاير للعائد عند آجال استحقاق محددة مسبقًا؛ [ 67 ] حيث يكون تباين كل مكون هو قيمته الذاتية (وبما أن المكونات متعامدة ، فلا حاجة لإدراج أي ارتباط في النمذجة اللاحقة).

في مجال الأسهم ، تُعتبر المحفظة المثلى تلك التي تُحقق أقصى عائد متوقع لمستوى معين من المخاطر، أو تلك التي تُقلل المخاطر إلى أدنى حد ممكن لعائد معين؛ انظر نموذج ماركويتز للمزيد من التفاصيل. وبالتالي، يتمثل أحد الأساليب في تقليل مخاطر المحفظة، حيث تُطبق استراتيجيات التخصيص على "المحافظ الرئيسية" بدلاً من الأسهم الأساسية . أما الأسلوب الثاني فيتمثل في تعزيز عائد المحفظة، باستخدام المكونات الرئيسية لاختيار أسهم الشركات ذات إمكانات النمو. [ 68 ] [ 69 ] كما استُخدم تحليل المكونات الرئيسية (PCA) لفهم العلاقات [ 64 ] بين أسواق الأسهم الدولية ، وداخل الأسواق بين مجموعات الشركات في الصناعات أو القطاعات .

يمكن تطبيق تحليل المكونات الرئيسية (PCA) أيضًا على اختبارات الضغط ، [ 70 ] وهي في جوهرها تحليل لقدرة البنك على تحمل سيناريو اقتصادي سلبي افتراضي . تكمن فائدتها في "تبسيط المعلومات الواردة في [عدة] متغيرات اقتصادية كلية إلى مجموعة بيانات أكثر قابلية للإدارة، والتي يمكن استخدامها بعد ذلك في التحليل". [ 70 ] هنا، يتم ربط العوامل الناتجة، على سبيل المثال، بأسعار الفائدة - بناءً على أكبر عناصر المتجه الذاتي للعامل - ثم تتم ملاحظة كيف تؤثر "صدمة" لكل عامل على الأصول الضمنية لكل بنك.

علم الأعصاب

يُستخدم أحد أشكال تحليل المكونات الرئيسية في علم الأعصاب لتحديد الخصائص المحددة للمنبه التي تزيد من احتمالية توليد جهد فعل في العصبون . [ 71 ] [ 72 ] تُعرف هذه التقنية بتحليل التغاير المُحفَّز بالنبضات . في تطبيق نموذجي، يُقدِّم الباحث عملية ضوضاء بيضاء كمنبه (عادةً إما كمدخل حسي لعينة الاختبار، أو كتيار يُحقن مباشرةً في العصبون) ويسجل سلسلة من جهود الفعل، أو النبضات، التي ينتجها العصبون نتيجةً لذلك. من المفترض أن بعض خصائص المنبه تجعل العصبون أكثر عرضةً لتوليد النبضات. لاستخراج هذه الخصائص، يحسب الباحث مصفوفة التغاير لمجموعة النبضات المُحفَّزة ، وهي مجموعة جميع المنبهات (المُعرَّفة والمُجزَّأة على مدى فترة زمنية محدودة، عادةً في حدود 100 مللي ثانية) التي سبقت النبضة مباشرةً. تشير المتجهات الذاتية للفرق بين مصفوفة التغاير المُحفَّزة بالنبضات ومصفوفة التغاير لمجموعة المحفزات السابقة (مجموعة جميع المحفزات، المُعرَّفة على نفس الفترة الزمنية) إلى الاتجاهات في فضاء المحفزات التي اختلف فيها تباين مجموعة المحفزات المُحفَّزة بالنبضات اختلافًا كبيرًا عن تباين مجموعة المحفزات السابقة. وبالتحديد، تتوافق المتجهات الذاتية ذات القيم الذاتية الموجبة الأكبر مع الاتجاهات التي أظهر فيها تباين مجموعة المحفزات المُحفَّزة بالنبضات أكبر تغير إيجابي مقارنةً بتباين المجموعة السابقة. وبما أن هذه هي الاتجاهات التي أدى فيها تغيير المحفز إلى حدوث نبضة، فإنها غالبًا ما تُمثل تقريبات جيدة لخصائص المحفزات المطلوبة.

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

تُعدّ تقنية تحليل المكونات الرئيسية (PCA) لتقليل الأبعاد مناسبة بشكل خاص للكشف عن الأنشطة المنسقة لمجموعات الخلايا العصبية الكبيرة. وقد استُخدمت في تحديد المتغيرات الجماعية، أي معلمات الترتيب ، أثناء التحولات الطورية في الدماغ. [ 73 ]

العلاقة مع الطرق الأخرى

تحليل المراسلات

طُوِّر تحليل التوافق (CA) على يد جان بول بنزكري [ 74 ] ، وهو مشابهٌ من حيث المفهوم لتحليل المكونات الرئيسية (PCA)، ولكنه يُعيد تحجيم البيانات (التي ينبغي أن تكون غير سالبة) بحيث تُعامل الصفوف والأعمدة معاملةً متساوية. ويُطبَّق تقليديًا على جداول التوافق . يُحلِّل تحليل التوافق إحصائية مربع كاي المرتبطة بهذا الجدول إلى عوامل متعامدة. [ 75 ] ولأن تحليل التوافق أسلوبٌ وصفي، يُمكن تطبيقه على الجداول التي تكون إحصائية مربع كاي مناسبةً لها أو غير مناسبة. تتوفر عدة أنواع من تحليل التوافق، بما في ذلك تحليل التوافق المُعدَّل وتحليل التوافق الكنسي . ومن الامتدادات الخاصة تحليل التوافق المتعدد ، الذي يُمكن اعتباره نظيرًا لتحليل المكونات الرئيسية للبيانات الفئوية. [ 76 ]

تحليل العوامل

الصورة أعلاه مثال على الفرق بين تحليل المكونات الرئيسية (PCA) وتحليل العوامل. في الرسم البياني العلوي، يمثل "العامل" (مثل المسار الوظيفي) المتغيرات الثلاثة الملاحظة (مثل طبيب، محامٍ، معلم)، بينما في الرسم البياني السفلي، يتم اختزال المتغيرات الملاحظة (مثل معلم رياض أطفال، معلم مدرسة متوسطة، معلم مدرسة ثانوية) إلى المكون محل الاهتمام (مثل معلم).

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

يُشبه تحليل العوامل تحليل المكونات الرئيسية في كونه يتضمن أيضًا تراكيب خطية للمتغيرات. وعلى عكس تحليل المكونات الرئيسية، يركز تحليل العوامل على الارتباطات، ساعيًا إلى إعادة إنتاج الارتباطات المتبادلة بين المتغيرات، حيث تمثل العوامل "التباين المشترك للمتغيرات، باستثناء التباين الفريد". [ 77 ] وفيما يتعلق بمصفوفة الارتباط، يتوافق هذا مع التركيز على تفسير العناصر غير القطرية (أي التباين المشترك)، بينما يركز تحليل المكونات الرئيسية على تفسير العناصر القطرية. ومع ذلك، وكنتيجة جانبية، عند محاولة إعادة إنتاج العناصر القطرية، يميل تحليل المكونات الرئيسية أيضًا إلى مطابقة الارتباطات غير القطرية بشكل جيد نسبيًا. [ 13 ] : 158 تتشابه نتائج تحليل المكونات الرئيسية وتحليل العوامل في معظم الحالات، ولكن هذا ليس هو الحال دائمًا، وهناك بعض المشكلات التي تختلف فيها النتائج اختلافًا كبيرًا. يُستخدم تحليل العوامل عمومًا عندما يكون الغرض من البحث هو الكشف عن بنية البيانات (أي البنى الكامنة أو العوامل) أو النمذجة السببية . إذا تمت صياغة نموذج العوامل بشكل غير صحيح أو لم يتم استيفاء الافتراضات، فإن تحليل العوامل سيعطي نتائج خاطئة. [ 78 ]

التجميع باستخدام خوارزمية K -means

لقد زُعم أن الحل المُبسّط لتجميع البيانات باستخدام خوارزمية k -means ، والمُحدد بمؤشرات التجميع، يُعطى بالمكونات الرئيسية، وأن فضاء تحليل المكونات الرئيسية (PCA) المُمتد بواسطة الاتجاهات الرئيسية يُطابق فضاء مركز التجميع. [ 79 ] [ 80 ] ومع ذلك، فإن كون تحليل المكونات الرئيسية (PCA) تبسيطًا مفيدًا لتجميع البيانات باستخدام خوارزمية k -means لم يكن نتيجة جديدة، [ 81 ] ومن السهل العثور على أمثلة مضادة لعبارة أن فضاء مركز التجميع مُمتد بواسطة الاتجاهات الرئيسية. [ 82 ]

تحليل المصفوفات غير السالبة

مخططات التباين الجزئي المتبقي (FRV) لتحليل المكونات الرئيسية (PCA) وتحليل المصفوفات غير السالبة (NMF)؛ [ 25 ] بالنسبة لتحليل المكونات الرئيسية، تمثل القيم النظرية مساهمة القيم الذاتية المتبقية. بالمقارنة، تصل منحنيات التباين الجزئي المتبقي لتحليل المكونات الرئيسية إلى مستوى ثابت حيث لا يتم التقاط أي إشارة بفعالية؛ بينما تنخفض منحنيات التباين الجزئي المتبقي لتحليل المصفوفات غير السالبة باستمرار، مما يشير إلى قدرة أفضل على التقاط الإشارة. كما تتقارب منحنيات التباين الجزئي المتبقي لتحليل المصفوفات غير السالبة إلى مستويات أعلى من تحليل المكونات الرئيسية، مما يدل على انخفاض خاصية التوفيق الزائد لتحليل المصفوفات غير السالبة.

يُعدّ تحليل المصفوفات غير السالبة (NMF) أسلوبًا لتقليل الأبعاد، حيث تُستخدم فيه العناصر غير السالبة فقط في المصفوفات، مما يجعله أسلوبًا واعدًا في علم الفلك، [ 23 ] [ 24 ] [ 25 ] نظرًا لأن الإشارات الفيزيائية الفلكية غير سالبة. تكون مكونات تحليل المكونات الرئيسية (PCA) متعامدة فيما بينها، بينما تكون جميع مكونات تحليل المصفوفات غير السالبة غير سالبة، وبالتالي تُنشئ أساسًا غير متعامد.

في تحليل المكونات الرئيسية (PCA)، يُرتب إسهام كل مكون بناءً على قيمة قيمته الذاتية المقابلة، والتي تُعادل التباين الجزئي المتبقي (FRV) في تحليل البيانات التجريبية. [ 21 ] أما في تحليل المصفوفات غير السالبة (NMF)، فتُرتب مكوناته بناءً على منحنيات التباين الجزئي المتبقي التجريبية فقط. [ 25 ] وتُظهر مخططات القيم الذاتية الجزئية المتبقية، أي1-أنا=1كλأنا/ج=1نλج{\displaystyle 1-\sum _{i=1}^{k}\lambda _{i}{\Big /}\sum _{j=1}^{n}\lambda _{j}}كدالة لعدد المكوناتك{\displaystyle k}بالنظر إلى إجمالين{\displaystyle n}تتميز مكونات تحليل المكونات الرئيسية (PCA) بثباتها عند مستوى ثابت، حيث لا يتم جمع أي بيانات لإزالة الضوضاء شبه الساكنة، ثم تنخفض المنحنيات بسرعة كمؤشر على فرط التخصيص (الضوضاء العشوائية). [ 21 ] أما منحنيات FRV لتحليل المصفوفات غير السالبة (NMF) فتنخفض باستمرار [ 25 ] عند إنشاء مكونات NMF بالتتابع ، [ 24 ] مما يشير إلى استمرار جمع الضوضاء شبه الساكنة؛ ثم تتقارب إلى مستويات أعلى من PCA، [ 25 ] مما يشير إلى انخفاض خاصية فرط التخصيص في NMF.

أيقونات العلاقات

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

أيقونات الارتباطات – الجيوكيمياء للهباء الجوي البحري

أما أيقونات الارتباطات ، على النقيض من ذلك، والتي لا تُمثل إسقاطًا على نظام محاور، فلا تعاني من هذه العيوب. ولذلك، يمكننا الاحتفاظ بجميع المتغيرات.

يتمثل مبدأ الرسم البياني في تسليط الضوء على الارتباطات "الرائعة" لمصفوفة الارتباط، بخط متصل (ارتباط إيجابي) أو خط منقط (ارتباط سلبي).

لا يُعدّ الارتباط القوي "ملفتًا للنظر" إذا لم يكن مباشرًا، بل ناتجًا عن تأثير متغير ثالث. وعلى العكس، قد تكون الارتباطات الضعيفة "ملفتة للنظر". فعلى سبيل المثال، إذا كان المتغير Y يعتمد على عدة متغيرات مستقلة، فإن ارتباطات Y بكل منها تكون ضعيفة، ومع ذلك تُعدّ "ملفتة للنظر".

التعميمات

تحليل المكونات الرئيسية المتفرقة

من أبرز عيوب تحليل المكونات الرئيسية (PCA) أن المكونات الرئيسية عادةً ما تكون تراكيب خطية لجميع متغيرات الإدخال. يتغلب تحليل المكونات الرئيسية المتفرق ( Sparse PCA ) على هذا العيب من خلال إيجاد تراكيب خطية تحتوي على عدد قليل فقط من متغيرات الإدخال. وهو يوسع نطاق طريقة تحليل المكونات الرئيسية الكلاسيكية (PCA) لتقليل أبعاد البيانات بإضافة قيد التفرق على متغيرات الإدخال. وقد تم اقتراح العديد من المناهج، بما في ذلك

  • إطار الانحدار، [ 83 ]
  • إطار عمل للبرمجة شبه المحددة/الاسترخاء المحدب، [ 84 ]
  • إطار عمل عام لطريقة القوة [ 85 ]
  • إطار عمل التحسين المتناوب [ 86 ]
  • البحث الجشع الأمامي والخلفي والأساليب الدقيقة باستخدام تقنيات التفرع والتقييد، [ 87 ]
  • إطار الصياغة البايزية. [ 88 ]

تم استعراض التطورات المنهجية والنظرية لتقنية تحليل المكونات الرئيسية المتفرقة (Sparse PCA) بالإضافة إلى تطبيقاتها في الدراسات العلمية مؤخرًا في ورقة بحثية. [ 89 ]

تحليل المكونات الرئيسية غير الخطي

مقارنة بين تحليل المكونات الرئيسية الخطي والمشعبات الرئيسية غير الخطية [ 90 ] لتصوير بيانات المصفوفات الدقيقة لسرطان الثدي : أ) تكوين العقد والسطح الرئيسي ثنائي الأبعاد في مشعب PCA الخطي ثلاثي الأبعاد. مجموعة البيانات منحنية ولا يمكن تمثيلها بشكل كافٍ على مستوى رئيسي ثنائي الأبعاد؛ ب) التوزيع في إحداثيات السطح الرئيسي غير الخطي ثنائي الأبعاد الداخلي (ELMap2D) مع تقدير لكثافة النقاط؛ ج) نفس ما ورد في (ب)، ولكن لمشعب PCA الخطي ثنائي الأبعاد (PCA2D). يتم تصوير النوع الفرعي "القاعدي" من سرطان الثدي بشكل أكثر ملاءمة باستخدام ELMap2D، وتصبح بعض ميزات التوزيع أكثر وضوحًا مقارنةً بـ PCA2D. يتم إنتاج المشعبات الرئيسية بواسطة خوارزمية الخرائط المرنة . البيانات متاحة للمسابقات العامة. [ 91 ] البرنامج متاح للاستخدام غير التجاري مجانًا. [ 92 ]

تستند معظم الطرق الحديثة لتقليل الأبعاد غير الخطية إلى جذورها النظرية والخوارزمية في تحليل المكونات الرئيسية (PCA) أو خوارزمية K-means. كانت فكرة بيرسون الأصلية هي اختيار خط مستقيم (أو مستوى) يمثل "أفضل تطابق" لمجموعة من نقاط البيانات. وقد طوّر تريفور هاستي هذا المفهوم باقتراح المنحنيات الرئيسية [ 93 ] كامتداد طبيعي للتفسير الهندسي لتحليل المكونات الرئيسية، والذي يقوم بإنشاء متعدد الشعب لتقريب البيانات ، ثم إسقاط النقاط عليه. انظر أيضًا خوارزمية الخريطة المرنة وتحليل الجيوديسية الرئيسية [ 94 ] . ومن التعميمات الشائعة الأخرى تحليل المكونات الرئيسية باستخدام النواة ، والذي يتوافق مع تحليل المكونات الرئيسية الذي يتم إجراؤه في فضاء هيلبرت ذي النواة المُستنسخة المرتبط بنواة موجبة التحديد.

في مجال تعلم الفضاءات الفرعية متعددة الخطوط ، [ 95 ] [ 96 ] [ 97 ] تم تعميم تحليل المكونات الرئيسية (PCA) إلى تحليل المكونات الرئيسية متعدد الخطوط (MPCA) الذي يستخلص الميزات مباشرةً من تمثيلات الموترات. يتم حل MPCA من خلال تطبيق تحليل المكونات الرئيسية على كل نمط من أنماط الموترات بشكل تكراري. وقد طُبقت MPCA في مجالات مثل التعرف على الوجوه، والتعرف على المشية، وغيرها. كما تم توسيع نطاق MPCA ليشمل تحليل المكونات الرئيسية غير المرتبط، وتحليل المكونات الرئيسية غير السالب، وتحليل المكونات الرئيسية القوي.

يمكن إجراء تحليل المكونات الرئيسية N -way باستخدام نماذج مثل تحليل Tucker ، و PARAFAC ، وتحليل العوامل المتعددة، وتحليل التماثل المشترك، و STATIS، و DISTATIS.

تحليل المكونات الرئيسية القوي

على الرغم من أن تحليل المكونات الرئيسية (PCA) يجد الطريقة المثلى رياضيًا (أي تقليل مربع الخطأ)، إلا أنه لا يزال حساسًا للقيم الشاذة في البيانات التي تُنتج أخطاءً كبيرة، وهو أمرٌ يسعى التحليل لتجنبه في المقام الأول. لذلك، من الممارسات الشائعة إزالة القيم الشاذة قبل حساب تحليل المكونات الرئيسية. مع ذلك، في بعض السياقات، قد يصعب تحديد القيم الشاذة. [ 98 ] على سبيل المثال، في خوارزميات استخراج البيانات مثل التجميع الارتباطي ، لا يكون تخصيص النقاط للمجموعات والقيم الشاذة معروفًا مسبقًا. وقد ساهم تعميمٌ مُقترحٌ حديثًا لتحليل المكونات الرئيسية [ 99 ] ، قائمٌ على تحليل المكونات الرئيسية الموزون، في زيادة المتانة من خلال تخصيص أوزان مختلفة لعناصر البيانات بناءً على مدى ملاءمتها المُقدَّرة.

كما تم اقتراح متغيرات مقاومة للقيم الشاذة من تحليل المكونات الرئيسية، استنادًا إلى صيغ معيار L1 ( L1-PCA ). [ 7 ] [ 5 ]

يُعد تحليل المكونات الرئيسية القوي (RPCA) عبر التفكيك في مصفوفات منخفضة الرتبة ومتفرقة تعديلًا لتحليل المكونات الرئيسية (PCA) يعمل بكفاءة مع البيانات المشوهة بشدة. [ 100 ] [ 101 ] [ 102 ]

تقنيات مماثلة

تحليل المكونات المستقلة

يتم توجيه تحليل المكونات المستقلة (ICA) إلى مشاكل مماثلة لتحليل المكونات الرئيسية، ولكنه يجد مكونات قابلة للفصل بشكل إضافي بدلاً من التقريبات المتتالية.

تحليل مكونات الشبكة

بالنظر إلى مصفوفةهـ{\displaystyle E}، يحاول تحليلها إلى مصفوفتين بحيثهـ=أP{\displaystyle E=AP}يتمثل أحد الاختلافات الرئيسية عن تقنيات مثل PCA وICA في أن بعض مدخلاتأ{\displaystyle A}يتم تقييدها لتكون صفرًا. هناP{\displaystyle P}يُطلق عليها اسم الطبقة التنظيمية. في حين أن هذا التفكيك يمكن أن يكون له حلول متعددة بشكل عام، إلا أنها تثبت أنه إذا تحققت الشروط التالية  :

  1. أ{\displaystyle A}يتمتع بترتيب كامل في العمود
  2. كل عمود منأ{\displaystyle A}يجب أن يكون لديه على الأقلل-1{\displaystyle L-1}أصفار حيثل{\displaystyle L}عدد أعمدةأ{\displaystyle A}(أو بدلاً من ذلك عدد صفوفP{\displaystyle P}). إن مبرر هذا المعيار هو أنه إذا تمت إزالة عقدة من الطبقة التنظيمية مع جميع عقد الإخراج المتصلة بها، فيجب أن تظل النتيجة مميزة بمصفوفة اتصال ذات رتبة عمودية كاملة.
  3. P{\displaystyle P}يجب أن يكون الصف كاملاً.

عندئذٍ يكون التفكيك فريدًا حتى الضرب في عدد قياسي. [ 103 ]

تحليل التمييز للمكونات الرئيسية

يُعدّ تحليل التمييز للمكونات الرئيسية (DAPC) أسلوبًا متعدد المتغيرات يُستخدم لتحديد ووصف مجموعات الأفراد ذوي الصلة الجينية. يُقسّم التباين الجيني إلى مُكوّنين: التباين بين المجموعات والتباين داخل المجموعات، ويُركّز هذا الأسلوب على الأول. تُمثّل المُميّزات الخطية توليفات خطية من الأليلات التي تُفرّق بين المجموعات على أفضل وجه. وبالتالي، فإنّ الأليلات التي تُساهم بشكلٍ أكبر في هذا التمييز هي تلك التي تختلف بشكلٍ ملحوظ بين المجموعات. يُمكن أن تُتيح مُساهمات الأليلات في التجميعات التي يُحدّدها تحليل التمييز للمكونات الرئيسية (DAPC) تحديد مناطق الجينوم التي تُؤدّي إلى التباين الجيني بين المجموعات [ 104 ] . في تحليل التمييز للمكونات الرئيسية (DAPC)، تُحوّل البيانات أولًا باستخدام تحليل المكونات الرئيسية (PCA)، ثم تُحدّد المجموعات باستخدام تحليل التمييز (DA).

يمكن تنفيذ تحليل DAPC في لغة R باستخدام حزمة Adegenet. (للمزيد من المعلومات: adegenet على الإنترنت )

تحليل المكونات الاتجاهية

يُعد تحليل المكونات الاتجاهية (DCA) أسلوبًا يُستخدم في علوم الغلاف الجوي لتحليل مجموعات البيانات متعددة المتغيرات. [ 105 ] وكما هو الحال في تحليل المكونات الرئيسية (PCA)، يُتيح هذا الأسلوب تقليل الأبعاد، وتحسين التصور، وتعزيز قابلية تفسير مجموعات البيانات الكبيرة. ومثل تحليل المكونات الرئيسية، يعتمد تحليل المكونات الاتجاهية على مصفوفة التغاير المُستمدة من مجموعة البيانات المُدخلة. ويكمن الفرق بين تحليل المكونات الرئيسية وتحليل المكونات الاتجاهية في أن الأخير يتطلب إدخال اتجاه متجه، يُشار إليه بالتأثير. فبينما يُعظّم تحليل المكونات الرئيسية التباين المُفسَّر، يُعظّم تحليل المكونات الاتجاهية كثافة الاحتمال بالنظر إلى التأثير. ويكمن الدافع وراء استخدام تحليل المكونات الاتجاهية في إيجاد مكونات مجموعة البيانات متعددة المتغيرات التي تكون مُحتملة (تُقاس باستخدام كثافة الاحتمال) وهامة (تُقاس باستخدام التأثير). وقد استُخدم تحليل المكونات الاتجاهية لإيجاد أنماط موجات الحر الأكثر احتمالًا والأكثر خطورة في مجموعات التنبؤات الجوية، [ 106 ] والتغيرات الأكثر احتمالًا والأكثر تأثيرًا في هطول الأمطار نتيجة لتغير المناخ. [ 107 ]

البرنامج/شفرة المصدر

  • ALGLIB – مكتبة C++ و C# تُنفذ تحليل المكونات الرئيسية (PCA) وتحليل المكونات الرئيسية المقتطع (truncated PCA).
  • Analytica – تقوم وظيفة EigenDecomp المدمجة بحساب المكونات الرئيسية.
  • ELKI – يتضمن PCA للإسقاط، بما في ذلك المتغيرات القوية لـ PCA، بالإضافة إلى خوارزميات التجميع القائمة على PCA .
  • يمكن إجراء تحليل المكونات الرئيسية باستخدام Gretlpca إما عن طريق الأمر أو عن طريق princomp()الوظيفة.
  • جوليا – تدعم تحليل المكونات الرئيسية (PCA) باستخدام pcaالدالة الموجودة في حزمة MultivariateStats.
  • KNIME – برنامج قائم على لغة جافا لترتيب العقد لأغراض التحليل، حيث تُسمى العقد فيه PCA و PCA compute و PCA Apply و PCA inverse، مما يسهل الأمر.
  • برنامج Maple – يتم استخدام أمر PCA لإجراء تحليل المكونات الرئيسية على مجموعة من البيانات.
  • Mathematica – ينفذ تحليل المكونات الرئيسية باستخدام الأمر PrincipalComponents باستخدام كل من طرق التغاير والارتباط.
  • MathPHP – مكتبة رياضية بلغة PHP تدعم تحليل المكونات الرئيسية (PCA).
  • MATLAB – تُعد دالة SVD جزءًا من النظام الأساسي. في حزمة أدوات الإحصاء، تُعطي الدالتان princomp( pcaR2012b) المكونات الرئيسية، بينما pcaresتُعطي الدالة البواقي والمصفوفة المُعاد بناؤها لتقريب PCA منخفض الرتبة.
  • Matplotlib – مكتبة بايثون تحتوي على حزمة PCA في وحدة .mlab.
  • mlpack – يوفر تطبيقًا لتحليل المكونات الرئيسية في لغة C++ .
  • mrmath – مكتبة رياضية عالية الأداء لـ Delphi و FreePascal يمكنها إجراء تحليل المكونات الرئيسية (PCA)؛ بما في ذلك المتغيرات القوية.
  • مكتبة NAG - يتم تنفيذ تحليل المكونات الرئيسية عبر g03aaالروتين (المتوفر في كلا إصداري Fortran من المكتبة).
  • NMath – مكتبة عددية خاصة تحتوي على PCA لإطار عمل .NET .
  • برنامج GNU Octave – بيئة حسابية برمجية مجانية متوافقة إلى حد كبير مع برنامج MATLAB، princompوتعطي الدالة المكون الرئيسي.
  • OpenCV
  • قاعدة بيانات أوراكل 12c - يتم تنفيذها DBMS_DATA_MINING.SVDS_SCORING_MODEعن طريق تحديد قيمة الإعداد SVDS_SCORING_PCA.
  • برنامج Orange – يدمج تحليل المكونات الرئيسية (PCA) في بيئة البرمجة المرئية الخاصة به. يعرض تحليل المكونات الرئيسية مخططًا بيانيًا (درجة التباين المُفسَّر) حيث يمكن للمستخدم تحديد عدد المكونات الرئيسية بشكل تفاعلي.
  • أوريجين – يحتوي على PCA في نسخته الاحترافية.
  • Qlucore – برنامج تجاري لتحليل البيانات متعددة المتغيرات باستجابة فورية باستخدام PCA.
  • R – حزمة إحصائية مجانيةprincomp ، يمكن استخدام وظائفها prcompلتحليل المكونات الرئيسية؛ prcompوتستخدم تحليل القيم المفردة الذي يوفر عمومًا دقة عددية أفضل. تتضمن بعض الحزم التي تُنفذ تحليل المكونات الرئيسية في R، على سبيل المثال لا الحصر: ade4، vegan، ExPosition، dimRed، و FactoMineR.
  • SAS – برنامج احتكاري؛ على سبيل المثال، انظر. [ 108 ]
  • scikit-learn – مكتبة بايثون للتعلم الآلي تحتوي على PCA و Probabilistic PCA و Kernel PCA و Sparse PCA وتقنيات أخرى في وحدة التفكيك.
  • Scilab – حزمة حسابية عددية مجانية ومفتوحة المصدر ومتعددة المنصات، تقوم الدالة princompبحساب تحليل المكونات الرئيسية، وتقوم الدالة pcaبحساب تحليل المكونات الرئيسية باستخدام متغيرات موحدة.
  • SPSS – برنامج احتكاري يستخدمه علماء الاجتماع بشكل شائع لتحليل المكونات الرئيسية، وتحليل العوامل، وتحليل التجميع المرتبط بها.
  • Weka – مكتبة جافا للتعلم الآلي تحتوي على وحدات لحساب المكونات الرئيسية.

انظر أيضاً

مراجع

  1. ^ جيورز، فيليبي إل. فيريرا، غوستافو ر.؛ أرودا، هنريك إف دي؛ سيلفا، فيليبي ن.؛ كومين، سيزار H.؛ أمانسيو، دييغو ر.؛ كوستا، لوتشيانو دا ف. (24 مايو 2021). "تحليل المكونات الرئيسية: نهج طبيعي لاستكشاف البيانات" . ايه سي ام للحوسبة. Surv . 54 (4): 70:1-70:34. أرخايف : 1804.02502 . دوى : 10.1145/3447755 .
  2. جوليف، إيان ت.؛ كاديما، خورخي (13 أبريل 2016). "تحليل المكونات الرئيسية: مراجعة وتطورات حديثة" . المعاملات الفلسفية للجمعية الملكية أ: العلوم الرياضية والفيزيائية والهندسية . 374 (2065) 20150202. رمز Bibcode : 2016RSPTA.37450202J . doi : 10.1098/rsta.2015.0202 . PMC 4792409. PMID 26953178 .  
  3. بارنيت، تي بي و ر. بريسندورفر. (1987). "أصول ومستويات مهارة التنبؤ الشهري والموسمي لدرجات حرارة الهواء السطحي في الولايات المتحدة، كما تم تحديدها بواسطة تحليل الارتباط الكنسي" . مجلة الطقس الشهرية . 115 (9): 1825. Bibcode : 1987MWRv..115.1825B . doi : 10.1175/1520-0493(1987)115 < 1825:oaloma > 2.0.co ; 2 .
  4. هسو، دانيال؛ كاكادي، شام م.؛ تشانغ، تونغ (2008). خوارزمية طيفية لتعلم نماذج ماركوف المخفية . arXiv : 0811.4413 . Bibcode : 2008arXiv0811.4413H .
  5. 1 2 ماركوبولوس، بانوس ب.؛ كوندو، سانديبان؛ تشاماديا، شوبهام؛ بادوس، ديميتريس أ. (15 أغسطس 2017). "تحليل المكونات الرئيسية الفعال باستخدام معيار L1 عبر قلب البتات". معاملات IEEE في معالجة الإشارات . 65 (16): 4252-4264 . arXiv : 1610.01959 . Bibcode : 2017ITSP...65.4252M . doi : 10.1109/TSP.2017.2708023 . S2CID 7931130 . 
  6. 1 2 تشاشلاكيس، ديميتريس ج.؛ براتر-بينيت، آشلي؛ ماركوبولوس، بانوس ب. (22 نوفمبر 2019). "تحليل موتر تاكر بمعيار L1" . IEEE Access . 7 : 178454–178465 . arXiv : 1904.06455 . Bibcode : 2019IEEEA...7q8454C . doi : 10.1109/ACCESS.2019.2955134 .
  7. 1 2 ماركوبولوس، بانوس ب.؛ كاريستينوس، جورج ن.؛ بادوس، ديميتريس أ. (أكتوبر 2014). "الخوارزميات المثلى لمعالجة الإشارات في فضاء L1 الفرعي". معاملات IEEE في معالجة الإشارات . 62 (19): 5046-5058 . arXiv : 1405.6785 . Bibcode : 2014ITSP...62.5046M . doi : 10.1109/TSP.2014.2338077 . S2CID 1494171 . 
  8. زان، ج.؛ فاسواني، ن. (2015). "تحليل المكونات الرئيسية القوي مع معرفة جزئية بالفضاء الفرعي". معاملات IEEE في معالجة الإشارات . 63 (13): 3332-3347 . arXiv : 1403.1591 . Bibcode : 2015ITSP...63.3332Z . doi : 10.1109/tsp.2015.2421485 . S2CID 1516440 . 
  9. كانادي، ت.؛ كي، قيفا (يونيو 2005). "تحليل معيار L₁ القوي في وجود القيم الشاذة والبيانات المفقودة باستخدام البرمجة المحدبة البديلة". مؤتمر جمعية مهندسي الكهرباء والإلكترونيات (IEEE) لعام 2005 حول رؤية الحاسوب والتعرف على الأنماط (CVPR'05) . المجلد 1. IEEE. الصفحات 739-746 . CiteSeerX 10.1.1.63.4605 . doi : 10.1109/CVPR.2005.309 . ISBN    978-0-7695-2372-9. S2CID 17144854 . 
  10. بيرسون، ك. (1901). "حول الخطوط والمستويات الأقرب ملاءمة لأنظمة النقاط في الفضاء" . المجلة الفلسفية . 2 (11): 559-572 . doi : 10.1080/14786440109462720 . S2CID 125037489 . 
  11. هوتلينج، هـ. (1933). تحليل مجموعة من المتغيرات الإحصائية إلى مكونات رئيسية. مجلة علم النفس التربوي ، 24 ، 417-441، و498-520. هوتلينج، هـ. (1936). "العلاقات بين مجموعتين من المتغيرات". بيومتريكا . 28 (3/4): 321-377 . doi : 10.2307/2333955 . JSTOR 2333955 . 
  12. ستيوارت، جي دبليو (1993). "حول التاريخ المبكر لتحليل القيم المفردة" . مجلة SIAM . 35 (4): 551-566 . Bibcode : 1993SIAMR..35..551S . doi : 10.1137/1035134 . hdl : 1903/566 .
  13. 1 2 3 4 5 جولييف، آي تي ​​(2002). تحليل المكونات الرئيسية . سلسلة سبرينغر في الإحصاء. نيويورك: سبرينغر-فيرلاغ. doi : 10.1007/b98835 . ISBN 978-0-387-95442-4.
  14. 1 2 3 4 هولمز، مارك هـ. (2023). مقدمة في الحوسبة العلمية وتحليل البيانات . نصوص في علوم وهندسة الحوسبة ( الطبعة الثانية). سبرينغر. ص 475-490 . ISBN   978-3-031-22429-4.
  15. فوركمان، ج.، جوس، ج.، بيفو، هـ.ب. (2019). "اختبارات الفرضيات لتحليل المكونات الرئيسية عند توحيد المتغيرات" . مجلة الإحصاءات الزراعية والبيولوجية والبيئية . 24 (2): 289-308 . Bibcode : 2019JABES..24..289F . doi : 10.1007/s13253-019-00355-5 .{{cite journal}}: صيانة CS1: أسماء متعددة: قائمة المؤلفين ( رابط )
  16. بويد، ستيفن؛ فاندنبيرغ، ليفين (8 مارس 2004). التحسين المحدب . مطبعة جامعة كامبريدج. doi : 10.1017/cbo9780511804441 . ISBN 978-0-521-83378-3.
  17. فوكوناغا، كينوسوكي (1990). مقدمة في التعرف الإحصائي على الأنماط . إلسيفير. ISBN 978-0-12-269851-4.
  18. علي زاده، إلاهه؛ ليونز، سامانثا م؛ كاسل، جوردان م؛ براساد، أشوك (2016). "قياس التغيرات المنهجية في شكل خلايا السرطان الغازية باستخدام عزم زيرنيك" . علم الأحياء التكاملي . 8 (11): 1183-1193 . doi : 10.1039/C6IB00100A . PMID 27735002 . 
  19. ليزنيك، م؛ توفاليس، ج. 2005 تقدير المكونات الرئيسية الثابتة باستخدام الانحدار القطري.
  20. جوناثان شلينز، درس تعليمي حول تحليل المكونات الرئيسية.
  21. 1 2 3 سومر، ريمي؛ بويو، لوران؛ لاركين، جيمس (2012). "الكشف عن الكواكب الخارجية والأقراص وتوصيفها باستخدام إسقاطات على صور كارونين-لوف الذاتية". رسائل المجلة الفيزيائية الفلكية . 755 (2): L28. arXiv : 1207.4197 . Bibcode : 2012ApJ...755L..28S . doi : 10.1088/2041-8205/755/2/L28 . S2CID 51088743 . 
  22. بويو، لوران (2016). "الكشف عن الكواكب الخارجية وتوصيفها باستخدام إسقاطات على صور كارونين لوف الذاتية: النمذجة الأمامية" . المجلة الفيزيائية الفلكية . 824 (2): 117. arXiv : 1604.06097 . Bibcode : 2016ApJ...824..117P . doi : 10.3847/0004-637X/824/2/117 . S2CID 118349503 . 
  23. 1 2 بلانتون، مايكل ر.؛ رويس، سام (2007). "تصحيحات K وتحويلات المرشحات في نطاق الأشعة فوق البنفسجية والضوئية والأشعة تحت الحمراء القريبة". المجلة الفلكية . 133 (2): 734-754 . arXiv : astro-ph/0606170 . Bibcode : 2007AJ....133..734B . doi : 10.1086/510127 . S2CID 18561804 . 
  24. 1 2 3 تشو، غوانغتون ب. (19-12-2016). "تحليل المصفوفات غير السالبة (NMF) مع عدم اليقين غير المتجانس والبيانات المفقودة". arXiv : 1612.06037 [ astro-ph.IM ].
  25. 1 2 3 4 5 6 رين، بين؛ بويو، لوران؛ تشو، غوانغتون ب.؛ دوشين، غاسبار (2018). "تحليل المصفوفات غير السالبة: استخلاص قوي للهياكل الممتدة" . المجلة الفيزيائية الفلكية . 852 (2): 104. arXiv : 1712.10317 . Bibcode : 2018ApJ...852..104R . doi : 10.3847/1538-4357/aaa1f2 . S2CID 3966513 . 
  26. "ما هي مزايا وعيوب تحليل المكونات الرئيسية؟" . i2tutorials . 1 سبتمبر 2019. تم الاطلاع عليه في 4 يونيو 2021 .
  27. أبوت، دين (مايو 2014). التحليلات التنبؤية التطبيقية . وايلي. ISBN 978-1-118-72796-6.
  28. 1 2 جيانغ، هونغ؛ إسكردج، كينت م. (2000). "التحيز في تحليل المكونات الرئيسية بسبب الملاحظات المترابطة" . مؤتمر الإحصاء التطبيقي في الزراعة . doi : 10.4148/2475-7772.1247 . ISSN 2475-7772 . 
  29. لينسكر، رالف (مارس 1988). "التنظيم الذاتي في شبكة إدراكية". مجلة IEEE Computer . 21 (3): 105-117 . Bibcode : 1988Compr..21c.105L . doi : 10.1109/2.36 . S2CID 1527671 . 
  30. ديكو وأوبرادوفيتش (1996). مدخل نظري للمعلومات في الحوسبة العصبية . نيويورك، نيويورك: سبرينغر. ISBN 978-1-4612-4016-7.
  31. بلومبلي، مارك (1991). نظرية المعلومات والشبكات العصبية غير الخاضعة للإشراف .ملاحظة فنية
  32. جيجر، برنارد؛ كوبين، جيرنوت (يناير 2013). "تحسين الإشارة كتقليل لفقدان المعلومات ذات الصلة". وقائع مؤتمر ITG حول الأنظمة والاتصالات والترميز . arXiv : 1205.6935 . Bibcode : 2012arXiv1205.6935G .
  33. انظر أيضًا إلى البرنامج التعليمي هنا
  34. "دليل الإحصاء الهندسي، القسم 6.5.5.2" . تم الاطلاع عليه بتاريخ 19 يناير 2015 .
  35. أ. أ. ميراندا، ي. أ. لو بورن، وج. بونتيمبي. مسارات جديدة من الحد الأدنى لخطأ التقريب إلى المكونات الرئيسية ، المجلد 27، العدد 3 / يونيو 2008، رسائل المعالجة العصبية، سبرينغر
  36. عبدي، ح. وويليامز، ل. ج. (2010). "تحليل المكونات الرئيسية". مراجعات وايلي متعددة التخصصات: الإحصاءات الحاسوبية . 2 (4): 433-459 . arXiv : 1108.4372 . doi : 10.1002/wics.101 . S2CID 122379222 . 
  37. "دليل مستخدم SAS/STAT(R) 9.3" .
  38. توثيق دالة eig في برنامج Matlab
  39. "نظام التعرف على الوجوه - قائم على تحليل المكونات الرئيسية" . www.mathworks.com . ١٩ يونيو ٢٠٢٣.
  40. توثيق دالة القيم الذاتية في برنامج Mathematica
  41. رويس، سام. "خوارزميات EM لتحليل المكونات الرئيسية وتحليل المكونات الرئيسية المتسلسل". التطورات في أنظمة معالجة المعلومات العصبية. تحرير مايكل آي. جوردان، ومايكل جيه. كيرنز، وسارة أ. سولا. مطبعة معهد ماساتشوستس للتكنولوجيا، 1998.
  42. جيلادي، بول؛ كوالسكي، بروس (1986). "انحدار المربعات الصغرى الجزئية: دليل تعليمي". مجلة التحليل الكيميائي . 185 : 1-17 . رمز Bibcode : 1986AcAC..185....1G . doi : 10.1016/0003-2670(86)80028-9 .
  43. كريمر، ر. (1998). تقنيات القياس الكيميائي للتحليل الكمي . نيويورك: مطبعة سي آر سي. رقم ISBN 978-0-203-90980-5.
  44. أندريكوت، م. (2009). "تنفيذ متوازي لخوارزميات تحليل المكونات الرئيسية التكرارية باستخدام وحدة معالجة الرسومات". مجلة علم الأحياء الحاسوبي . 16 (11): 1593-1599 . arXiv : 0811.1081 . doi : 10.1089/cmb.2008.0221 . PMID 19772385. S2CID 1362603 .  
  45. وارموث، إم كيه؛ كوزمين، دي. (2008). "خوارزميات تحليل المكونات الرئيسية العشوائية عبر الإنترنت مع حدود ندم لوغاريتمية في البعد" (ملف PDF) . مجلة أبحاث تعلم الآلة . 9 : 2287-2320 .
  46. ^ هوسون ، فرانسوا. لو سيباستيان. باجيس ، جيروم (2017/04/25). التحليل الاستكشافي متعدد المتغيرات بالمثال باستخدام R (0 ed.). تشابمان وهول/CRC. دوى : 10.1201/ب21874 . رقم ISBN  978-0-429-22543-7.
  47. عبدي، هيرفيه؛ ويليامز، لين جيه. (15 يوليو 2010). "تحليل المكونات الرئيسية" . إحصاءات WIRE الحاسوبية . 2 (4): 433-459 . doi : 10.1002/wics.101 . ISSN 1939-5108 . 
  48. أتكينسون، ويل (2023-04-28). "رسم خرائط الحقول والفضاءات كميًا: من تحليل التطابق المتعدد إلى تحليل المكونات الرئيسية الفئوية" . الجودة والكمية . 58 (1): 829-848 . doi : 10.1007/s11135-023-01669-w . hdl : 1983/7cbbf631-d50e-477e-8904-4e30f85b4bed . ISSN 0033-5177 . 
  49. ^ هوسون ، فرانسوا. جوسيه، جولي؛ سابورتا ، جيلبرت (2016/09/12). "جان دي ليو والمدرسة الفرنسية لتحليل البيانات" . مجلة البرامج الإحصائية . 73 : 1– 18. دوى : 10.18637/jss.v073.i06 . ردمك 1548-7660 . 
  50. ليبارت، لودوفيك؛ مورينو، آلان؛ وارويك، كينيث م. (1984). التحليل الإحصائي الوصفي متعدد المتغيرات: تحليل التوافق والتقنيات ذات الصلة للمصفوفات الكبيرة . وايلي. ISBN 978-0-471-86743-2.
  51. ^ ليبارت، لودوفيك. مورينو، آلان؛ بيرون، ماري (1995). إحصائية استكشافية متعددة الأبعاد . باريس: دونود. رقم ISBN 978-2-10-002886-3.
  52. لي، سيباستيان؛ جوسيه، جولي؛ هوسون ، فرانسوا (2008-03-18). "FactoMineR: حزمة R للتحليل متعدد المتغيرات" . مجلة البرامج الإحصائية . 25 : 1– 18. دوى : 10.18637/jss.v025.i01 . ردمك 1548-7660 . 
  53. كابلان، آر إم، وساكوزو، دي بي (2010). الاختبارات النفسية: المبادئ والتطبيقات والقضايا. (الطبعة الثامنة). بيلمونت، كاليفورنيا: وادزورث، سينجج ليرنينج.
  54. شيفكي، إشرف؛ ويليامز، مارلين (1949). المناطق الاجتماعية في لوس أنجلوس: التحليل والتصنيف . مطبعة جامعة كاليفورنيا.
  55. فلود، ج (2000). سيدني المنقسمة: إعادة النظر في علم البيئة العاملي. ورقة بحثية مقدمة إلى مؤتمر الجمعية الأمريكية لعلم النفس 2000، ملبورن، نوفمبر، وإلى المؤتمر الرابع والعشرين للجمعية الأسترالية والنيوزيلندية لعلم النفس التطبيقي، هوبارت، ديسمبر 2000.
  56. "المؤشرات الاجتماعية والاقتصادية للمناطق" . مكتب الإحصاءات الأسترالي . 2011. تم الاطلاع عليه بتاريخ 5 مايو 2022 .
  57. ^ شامبيرجر، تمارا. شوبيرث، فلوريان. هينسلر، يورغ (2023). “التحليل المركب التأكيدي في بحوث التنمية البشرية”. المجلة الدولية للتنمية السلوكية . 47 (1): 88-100 . دوى : 10.1177/01650254221117506 . اتش دي ال : 10362/143639 .
  58. تقارير التنمية البشرية. "مؤشر التنمية البشرية" . برنامج الأمم المتحدة الإنمائي . تم الاطلاع عليه بتاريخ 6 مايو 2022 .
  59. نوفمبر، جون؛ ستيفنز، ماثيو (2008). " تفسير تحليلات المكونات الرئيسية للتغير الوراثي المكاني للسكان" . نات جينيت . 40 (5): 646-49 . doi : 10.1038/ng.139 . PMC 3989108. PMID 18425127 .  
  60. إلهايك، إران (2022). "نتائج تحليل المكونات الرئيسية (PCA) في دراسات علم الوراثة السكانية متحيزة للغاية ويجب إعادة تقييمها" . التقارير العلمية . 12 (1) 14683. Bibcode : 2022NatSR..1214683E . doi : 10.1038/ s41598-022-14395-4 . PMC 9424212. PMID 36038559. S2CID 251932226 .   
  61. دي ساربو، واين؛ هاوسمان، روبرت؛ كوكيتز، جيفري (2007). " تحليل المكونات الرئيسية المقيدة لأبحاث التسويق" . مجلة التسويق في الإدارة . 2 : 305-328 عبر ResearchGate.
  62. داتون، ويليام هـ؛ بلانك، غرانت (2013). ثقافات الإنترنت: الإنترنت في بريطانيا (ملف PDF) . معهد أكسفورد للإنترنت. ص 6. 
  63. فلود، جو (2008). "التحليل متعدد الحدود لمسح الوظائف في مجال الإسكان" . ورقة بحثية مقدمة إلى مؤتمر الشبكة الأوروبية لأبحاث الإسكان، دبلن . تم الاطلاع عليها بتاريخ 6 مايو 2022 .
  64. ١ ٢ انظر الفصل ٩ في كتاب مايكل ب. ميلر (٢٠١٣). الرياضيات والإحصاء لإدارة المخاطر المالية ، الطبعة الثانية. وايلي ISBN 978-1-118-75029-2
  65. ١ ٢ §٩.٧ في جون هول (٢٠١٨). إدارة المخاطر والمؤسسات المالية، الطبعة الخامسة. وايلي. ISBN 1119448115
  66. §III.A.3.7.2 في كارول ألكسندر وإليزابيث شيدي، محرران (2004). دليل مديري المخاطر المحترفين . PRMIA . ISBN 978-0976609704
  67. مثال على التفكيك ، جون هول
  68. ليبين يانغ. تطبيق تحليل المكونات الرئيسية على إدارة محافظ الأسهم . قسم الاقتصاد والمالية، جامعة كانتربري ، يناير 2015.
  69. جورجيا باسيني (2017)؛ تحليل المكونات الرئيسية لإدارة محافظ الأسهم . المجلة الدولية للرياضيات البحتة والتطبيقية . المجلد 115، العدد 1، 2017، 153-167
  70. ١ ٢ انظر الفصل ٢٥ § "اختبار السيناريوهات باستخدام تحليل المكونات الرئيسية" في لي أونغ (٢٠١٤). "دليل أساليب ونماذج اختبار الإجهاد لصندوق النقد الدولي" ، صندوق النقد الدولي
  71. تشابين، جون؛ نيكوليس، ميغيل (1999). "تحليل المكونات الرئيسية لنشاط المجموعة العصبية يكشف عن تمثيلات حسية جسدية متعددة الأبعاد". مجلة أساليب علم الأعصاب . 94 (1): 121-140 . doi : 10.1016/S0165-0270(99)00130-2 . PMID 10638820. S2CID 17786731 .  
  72. ^ برينر، ن.، بياليك، دبليو، ودي رويتر فان ستيفنينك، ر.ر (2000).
  73. جيرسا، فيكتور؛ فريدريش، ر؛ هاكن، هيرمان؛ كيلسو، سكوت (1994). "نموذج نظري للتحولات الطورية في الدماغ البشري". علم التحكم الحيوي . 71 (1): 27-35 . doi : 10.1007/bf00198909 . PMID 8054384. S2CID 5155075 .  
  74. ^ بنزكري، ج.-ب. (1973). تحليل البيانات. المجلد الثاني. تحليل المراسلات . باريس، فرنسا: دونود.
  75. غرينيكر، مايكل (1983). نظرية وتطبيقات تحليل التوافق . لندن: أكاديميك برس. ISBN 978-0-12-299050-2.
  76. لو رو؛ بريجيت وهنري روانيه (2004). تحليل البيانات الهندسية، من تحليل التطابق إلى تحليل البيانات المهيكلة . دوردريخت: كلوير. ISBN 978-1-4020-2235-7.
  77. تيموثي أ. براون. التحليل العاملي التأكيدي لمنهجية البحث التطبيقي في العلوم الاجتماعية . مطبعة جيلفورد، 2006
  78. ميغلين، ر. ر. (1991). "فحص قواعد البيانات الكبيرة: منهج كيميائي إحصائي باستخدام تحليل المكونات الرئيسية". مجلة الكيمياء الإحصائية . 5 (3): 163-179 . doi : 10.1002/cem.1180050305 . S2CID 120886184 . 
  79. هـ. تشا؛ ج. دينغ؛ م. غو؛ إكس. هي؛ إتش دي سيمون (ديسمبر 2001). "الاسترخاء الطيفي لتجميع البيانات باستخدام خوارزمية K-means" (ملف PDF) . أنظمة معالجة المعلومات العصبية، المجلد 14 (NIPS 2001) : 1057-1064 .
  80. كريس دينغ؛ شياوفينغ هي (يوليو 2004). "تجميع البيانات باستخدام خوارزمية K-means عبر تحليل المكونات الرئيسية" (ملف PDF) . وقائع المؤتمر الدولي للتعلم الآلي (ICML 2004) : 225-232 .
  81. درينياس، ب.؛ أ. فريز؛ ر. كانان؛ س. فيمبالا؛ ف. فيناي (2004). "تجميع الرسوم البيانية الكبيرة باستخدام تحليل القيم المفردة" (ملف PDF) . تعلم الآلة . 56 ( 1-3 ): 9-33 . Bibcode : 2004MLear..56....9D . doi : 10.1023/b:mach.0000033113.59016.96 . S2CID 5892850. تاريخ الاسترجاع: 2012-08-02 . 
  82. كوهين، م.؛ إس. إلدر؛ سي. موسكو؛ سي. موسكو؛ إم. بيرسو (2014). تقليل الأبعاد لتجميع البيانات باستخدام خوارزمية k-means وتقريب الرتبة المنخفضة (الملحق ب) . arXiv : 1410.6801 . Bibcode : 2014arXiv1410.6801C .
  83. هوي زو؛ تريفور هاستي؛ روبرت تيبشيراني (2006). "تحليل المكونات الرئيسية المتفرقة" (ملف PDF) . مجلة الإحصاءات الحاسوبية والرسومية . 15 (2): 262-286 . CiteSeerX 10.1.1.62.580 . doi : 10.1198/106186006x113430 . S2CID 5730904 .  
  84. ألكسندر داسبريمون؛ لوران الغاوي؛ مايكل آي. جوردان؛ جيرت آر جي لانكريت (2007). "صياغة مباشرة لتحليل المكونات الرئيسية المتفرقة باستخدام البرمجة شبه المحددة" (ملف PDF) . مجلة SIAM Review . 49 (3): 434-448 . arXiv : cs/0406021 . Bibcode : 2007SIAMR..49..434D . doi : 10.1137/050645506 . S2CID 5490061 . 
  85. ميشيل جورني؛ يوري نيستيروف؛ بيتر ريشتاريك؛ رودولف سيبولكر (2010). "طريقة القوة المعممة لتحليل المكونات الرئيسية المتفرقة" (ملف PDF) . مجلة أبحاث تعلم الآلة . 11 : 517-553 . arXiv : 0811.4724 . Bibcode : 2008arXiv0811.4724J . ورقة نقاش CORE 2008/70.
  86. بيتر ريشتاريك؛ مارتن تاكاك؛ إس. داملا أهيباساوغلو (2012). "التعظيم المتناوب: إطار موحد لثمانية صيغ PCA متفرقة ورموز متوازية فعالة". arXiv : 1212.4137 [ stat.ML ].
  87. باباك مقدم؛ يائير فايس؛ شاي أفيدان (2005). "الحدود الطيفية لتحليل المكونات الرئيسية المتفرقة: خوارزميات دقيقة وجشعة" (ملف PDF) . التطورات في أنظمة معالجة المعلومات العصبية . المجلد 18. مطبعة معهد ماساتشوستس للتكنولوجيا. 
  88. يو غوان؛ جينيفر داي (2009). "تحليل المكونات الرئيسية الاحتمالية المتفرقة" (ملف PDF) . وقائع ورشة عمل ومؤتمر مجلة أبحاث تعلم الآلة . 5 : 185.
  89. هوي زو؛ لينغتشو شيو (2018). "نظرة عامة انتقائية على تحليل المكونات الرئيسية المتفرقة" . وقائع معهد مهندسي الكهرباء والإلكترونيات . 106 (8): 1311-1320 . doi : 10.1109/JPROC.2018.2846588 .
  90. AN Gorban , AY Zinovyev, "الرسوم البيانية الرئيسية والمتشعبات" ، في: دليل البحث في تطبيقات واتجاهات التعلم الآلي: الخوارزميات والأساليب والتقنيات ، Olivas ES et al (محررون). مرجع علوم المعلومات، IGI Global: Hershey، PA، الولايات المتحدة الأمريكية، 2009. 28–59.
  91. ^ وانغ، واي. كلين، ج.ج. تشانغ، Y.؛ سيويرتس، صباحا؛ انظر أيها النائب؛ يانغ، ف؛ تالانتوف، د.؛ تيمرمانز، م. ماير فان جيلدر، أنا؛ يو، J.؛ وآخرون . (2005). “ملامح التعبير الجيني للتنبؤ بالورم الخبيث البعيد لسرطان الثدي الأولي السلبي للعقدة الليمفاوية”. لانسيت . 365 (9460): 671–679 . دوى : 10.1016 / S0140-6736(05)17947-1 . بميد 15721472 . S2CID 16358549 .   البيانات متاحة عبر الإنترنت
  92. زينوفييف، أ. "ViDaExpert - أداة تصور البيانات متعددة الأبعاد" . معهد كوري . باريس.(مجاني للاستخدام غير التجاري)
  93. هاستي، ت .؛ ستوتزل، و. (يونيو 1989). "المنحنيات الرئيسية" (ملف PDF) . مجلة الجمعية الإحصائية الأمريكية . 84 (406): 502-506 . Bibcode : 1989JASA...84..502H . doi : 10.1080/01621459.1989.10478797 .
  94. أ. ن. غوربان، ب. كيغل، د. س. وونش، أ. زينوفييف (محررون)، المتشعبات الرئيسية لتصور البيانات وتقليل الأبعاد ، سلسلة محاضرات في علوم وهندسة الحاسوب 58، سبرينغر، برلين - هايدلبرغ - نيويورك، 2007. رقم ISBN 978-3-540-73749-0
  95. فاسيلسكو، ماو؛ تيرزوبولوس، د. (2003). تحليل الفضاء الفرعي متعدد الخطوط لمجموعات الصور (ملف PDF) . وقائع مؤتمر IEEE حول رؤية الحاسوب والتعرف على الأنماط (CVPR'03). ماديسون، ويسكونسن.
  96. فاسيلسكو، ماو؛ تيرزوبولوس، د. (2002). التحليل متعدد الخطوط لمجموعات الصور: TensorFaces (ملف PDF) . سلسلة محاضرات في علوم الحاسوب 2350؛ (مقدمة في وقائع المؤتمر الأوروبي السابع لرؤية الحاسوب (ECCV'02)، كوبنهاغن، الدنمارك). سبرينغر، برلين، هايدلبرغ. doi : 10.1007/3-540-47969-4_30 . ISBN 978-3-540-43745-1.
  97. فاسيلسكو، ماو؛ تيرزوبولوس، د. (يونيو 2005). تحليل المكونات المستقلة متعددة الخطوط (ملف PDF) . وقائع مؤتمر IEEE حول رؤية الحاسوب والتعرف على الأنماط (CVPR'05). المجلد 1. سان دييغو، كاليفورنيا. الصفحات 547-553 .  
  98. كيريل سيمونوف، فيدور ف. فومين، بيتر أ. غولوفاتش، فهد بانولان (9-15 يونيو 2019). "تحسين تعقيد تحليل المكونات الرئيسية مع القيم الشاذة" . في: كماليكا تشودري، روسلان سالاخوتدينوف (محرران). وقائع المؤتمر الدولي السادس والثلاثين للتعلم الآلي (ICML 2019) . المجلد 97. لونغ بيتش، كاليفورنيا، الولايات المتحدة الأمريكية: PMLR. الصفحات 5818-5826 .  {{cite conference}}: صيانة CS1: أسماء متعددة: قائمة المؤلفين ( رابط )
  99. كريجل، إتش بي؛ كروجر، بي؛ شوبرت، إي؛ زيمك، إيه (2008). "إطار عام لزيادة متانة خوارزميات تجميع الارتباط القائمة على تحليل المكونات الرئيسية". إدارة قواعد البيانات العلمية والإحصائية . سلسلة محاضرات في علوم الحاسوب. المجلد 5069. الصفحات 418-435 . CiteSeerX 10.1.1.144.4864 . doi : 10.1007/978-3-540-69497-7_27 . ISBN    978-3-540-69476-2.
  100. إيمانويل ج. كانديس؛ شياودونغ لي؛ يي ما؛ جون رايت (2011). "تحليل المكونات الرئيسية القوي؟". مجلة ACM . 58 (3): 11. arXiv : 0912.3599 . doi : 10.1145/1970392.1970395 . S2CID 7128002 . 
  101. ت. بوومانز؛ إ. زحزة (2014). "تحليل المكونات الرئيسية القوي عبر البحث عن المكونات الرئيسية: مراجعة للتقييم المقارن في المراقبة بالفيديو". رؤية الحاسوب وفهم الصور . 122 : 22-34 . Bibcode : 2014CVIU..122...22B . doi : 10.1016/j.cviu.2013.11.009 .
  102. ت. بوومانز؛ أ. سوبرال؛ س. جافيد؛ س. جونغ؛ إ. زحزة (2015). "التحليل إلى مصفوفات منخفضة الرتبة بالإضافة إلى مصفوفات جمعية لفصل الخلفية/المقدمة: مراجعة لتقييم مقارن باستخدام مجموعة بيانات واسعة النطاق". مجلة علوم الحاسوب . 23 : 1-71 . arXiv : 1511.01245 . Bibcode : 2015arXiv151101245B . doi : 10.1016/j.cosrev.2016.11.001 . S2CID 10420698 . 
  103. لياو، جيه سي؛ بوسكولو، آر؛ يانغ، واي إل؛ تران، إل إم؛ ساباتي، سي ؛ رويتشودري، في بي (2003). " تحليل مكونات الشبكة: إعادة بناء الإشارات التنظيمية في الأنظمة البيولوجية" . وقائع الأكاديمية الوطنية للعلوم . 100 (26): 15522-15527 . Bibcode : 2003PNAS..10015522L . doi : 10.1073/pnas.2136632100 . PMC 307600. PMID 14673099 .  
  104. لياو، ت.؛ جومبارت، س.؛ ديفيلارد، ف.؛ بالو (2010). "التحليل التمييزي للمكونات الرئيسية: طريقة جديدة لتحليل التجمعات السكانية ذات البنية الجينية" . مجلة BMC Genetics . 11 : 11:94. doi : 10.1186/1471-2156-11-94 . PMC 2973851. PMID 20950446 .  
  105. جويسون، س. (2020). "بديل لتحليل المكونات الرئيسية لتقدير الأنماط السائدة لتقلبات المناخ والظواهر المتطرفة، مع تطبيق على هطول الأمطار الموسمية في الولايات المتحدة والصين" . الغلاف الجوي . 11 (4): 354. Bibcode : 2020Atmos..11..354J . doi : 10.3390/atmos11040354 .
  106. شير، س.؛ جويسون، س.؛ ميسوري، ج. (2021). "سيناريوهات أسوأ الحالات القوية من التنبؤات الجماعية" . الطقس والتنبؤ . 36 (4): 1357-1373 . Bibcode : 2021WtFor..36.1357S . doi : 10.1175/WAF-D-20-0219.1 . S2CID 236300040 . 
  107. جيوسون، س.؛ ميسوري، ج.؛ بارباتو، ج.؛ ميركوليانو، ب.؛ ميسياك، ج.؛ ساسي، م. (2022). "تطوير سيناريوهات تأثير تمثيلية من مجموعات إسقاطات المناخ، مع تطبيق على هطول الأمطار UKCP18 وEURO-CORDEX" . مجلة التقدم في نمذجة أنظمة الأرض . 15 (1) e2022MS003038. doi : 10.1029/2022MS003038 . S2CID 254965361 . 
  108. "تحليل المكونات الرئيسية" . معهد البحوث والتعليم الرقمي . جامعة كاليفورنيا في لوس أنجلوس . تم الاطلاع عليه بتاريخ 29 مايو 2018 .

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