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

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

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

الصياغة الرياضية

لنفترض وجود مصفوفة بيانات ،X{\displaystyle X}، حيث كل منص{\displaystyle p}تمثل الأعمدة متغيرًا مُدخلًا، وكل منهان{\displaystyle n}تمثل الصفوف عينة مستقلة من مجموعة البيانات. يفترض المرء أن كل عمود منX{\displaystyle X}يكون متوسطها صفرًا، وإلا يمكن طرح متوسط ​​كل عمود من كل عنصر من عناصرX{\displaystyle X}. يتركΣ=1ن-1XX{\displaystyle \Sigma ={\frac {1}{n-1}}X^{\top }X}لتكن مصفوفة التغاير التجريبية لـX{\displaystyle X}، والتي لها أبعادص×ص{\displaystyle p\times p}.

بفرض عدد صحيحك{\displaystyle k}مع1كص{\displaystyle 1\leq k\leq p}يمكن صياغة مشكلة تحليل المكونات الرئيسية المتفرقة على أنها تعظيم التباين على طول اتجاه ممثل بواسطة متجهvRص{\displaystyle v\in \mathbb {R} ^{p}}مع تقييد عدد عناصرها:

الأعلىvتيΣvرهناً بـv2=1v0ك.{\displaystyle {\begin{aligned}\max \quad &v^{T}\Sigma v\\{\text{subject to}}\quad &\left\Vert v\right\Vert _{2}=1\\&\left\Vert v\right\Vert _{0}\leq k.\end{aligned}}}المعادلة 1

ينص القيد الأول على أن v متجه وحدة . أما في القيد الثاني،v0{\displaystyle \left\Vert v\right\Vert _{0}}يمثل0{\displaystyle \ell _{0}}المعيار الزائف للمصفوفة v ، والذي يُعرَّف بأنه عدد مكوناتها غير الصفرية. لذا، ينص القيد الثاني على أن عدد المكونات غير الصفرية في v أقل من أو يساوي k ، وهو عادةً عدد صحيح أصغر بكثير من البُعد p . تُعرف القيمة المثلى للمعادلة 1 باسم أكبر قيمة ذاتية متفرقة من الرتبة k .

إذا أخذنا k=p ، فإن المشكلة تختزل إلى PCA العادي ، وتصبح القيمة المثلى هي أكبر قيمة ذاتية لمصفوفة التغاير Σ .

بعد إيجاد الحل الأمثل v ، يتم تقليص Σ للحصول على مصفوفة جديدة

Σ1=Σ-(vتيΣv)vvتي،{\displaystyle \Sigma _{1}=\Sigma -(v^{T}\Sigma v)vv^{T},}

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

التعريف المكافئ التالي مكتوب بصيغة المصفوفة. ليكنV{\displaystyle V}إذا كانت مصفوفة متناظرة من الرتبة p×p ، فيمكن إعادة كتابة مسألة تحليل المكونات الرئيسية المتفرقة على النحو التالي:

الأعلىتير(ΣV)رهناً بـتير(V)=1V0ك2Rأنك(V)=1،V0.{\displaystyle {\begin{aligned}\max \quad &Tr(\Sigma V)\\{\text{subject to}}\quad &Tr(V)=1\\&\Vert V\Vert _{0}\leq k^{2}\\&Rank(V)=1,V\succeq 0.\end{aligned}}}المعادلة 2

Tr هو أثر المصفوفة ، وV0{\displaystyle \Vert V\Vert _{0}}يمثل العناصر غير الصفرية في المصفوفة V. يشير السطر الأخير إلى أن رتبة المصفوفة V تساوي واحدًا وأنها شبه موجبة . يعني السطر الأخير أن واحدًا لديهV=vvتي{\displaystyle V=vv^{T}}وبالتالي فإن المعادلة 2 تعادل المعادلة 1 .

علاوة على ذلك، فإن قيد الرتبة في هذه الصيغة هو في الواقع زائد عن الحاجة، وبالتالي يمكن صياغة PCA المتفرقة على أنها برنامج شبه محدد مختلط الأعداد الصحيحة التالي [ 1 ].

الأعلىتير(ΣV)رهناً بـتير(V)=1|Vأنا،أنا|zأنا،أنا{1،...،ص}،|Vأنا،ج|12zأنا،أنا،ج{1،...،ص}:أناج،V0،z{0،1}ص،أناzأناك{\displaystyle {\begin{aligned}\max \quad &Tr(\Sigma V)\\{\text{subject to}}\quad &Tr(V)=1\\&\vert V_{i,i}\vert \leq z_{i},\forall i\in \{1,...,p\},\vert V_{i,j}\vert \leq {\frac {1}{2}}z_{i},\forall i,j\in \{1,...,p\}:i\neq j,\\&V\succeq 0,z\in \{0,1\}^{p},\sum _{i}z_{i}\leq k\end{aligned}}}المعادلة 3

بسبب قيد العدد، يصعب حل مسألة التعظيم بدقة، خاصةً عندما يكون البُعد p كبيرًا. في الواقع، تُعدّ مسألة تحليل المكونات الرئيسية المتفرقة في المعادلة 1 مسألة صعبة الحل من نوع NP بالمعنى الدقيق. [ 2 ]

الاعتبارات الحسابية

كما هو الحال مع معظم المشاكل المتفرقة، فإن اختيار المتغيرات في SPCA هو مشكلة غير محدبة صعبة حسابيًا من نوع NP-hard، [ 3 ] لذلك غالبًا ما يتم استخدام الخوارزميات الجشعة شبه المثلى لإيجاد الحلول.

تجدر الإشارة أيضًا إلى أن SPCA تُدخل معلمات فائقة تُحدد مدى تأثير معاقبة القيم الكبيرة للمعلمات. [ 4 ] قد يتطلب الأمر ضبط هذه المعلمات لتحقيق أداء مُرضٍ، مما يزيد من التكلفة الحسابية الإجمالية.

خوارزميات لـ SPCA

تم اقتراح العديد من المناهج البديلة ( للمعادلة 1 )، بما في ذلك

  • إطار الانحدار، [ 5 ]
  • إطار عمل لتحليل المصفوفة المعاقبة ، [ 6 ]
  • إطار عمل للبرمجة شبه المحددة/الاسترخاء المحدب، [ 7 ]
  • إطار عمل عام لطريقة القوة [ 8 ]
  • إطار عمل التحسين المتناوب [ 9 ]
  • البحث الجشع الأمامي والخلفي والأساليب الدقيقة باستخدام تقنيات التفرع والتقييد، [ 10 ]
  • نهج التفرع والتقييد الأمثل المعتمد [ 11 ]
  • إطار الصياغة البايزية. [ 12 ]
  • نهج التفرع والقطع شبه المحدد للأعداد الصحيحة المختلطة الأمثل المعتمد [ 1 ]

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

ملاحظات حول تخفيف البرمجة شبه المحددة

لقد اقتُرح أن تحليل المكونات الرئيسية المتفرقة يمكن تقريبه باستخدام البرمجة شبه المحددة (SDP). [ 7 ] إذا تم حذف قيد الرتبة وتخفيف قيد العددية بقيد محدب ذي معيار 1 ، فسيتم الحصول على استرخاء برمجة شبه محددة، والذي يمكن حله بكفاءة في وقت متعدد الحدود:

الأعلىتير(ΣV)رهناً بـتير(V)=11تي|V|1كV0.{\displaystyle {\begin{aligned}\max \quad &Tr(\Sigma V)\\{\text{subject to}}\quad &Tr(V)=1\\&\mathbf {1} ^{T}|V|\mathbf {1} \leq k\\&V\succeq 0.\end{aligned}}}المعادلة 3

في القيد الثاني،1{\displaystyle \mathbf {1} }هو متجه من الرتبة p×1 مكون من واحدات، و |V| هي المصفوفة التي عناصرها هي القيم المطلقة لعناصر V.

الحل الأمثلV{\displaystyle V}بالنسبة للمسألة المُخففة ، لا يُضمن أن تكون المعادلة 3 من الرتبة الأولى. في هذه الحالة،V{\displaystyle V}يمكن اقتطاعها للاحتفاظ فقط بالمتجه الذاتي المهيمن.

على الرغم من أن البرنامج شبه المحدد لا يتوسع إلى ما بعد n=300 متغير مشترك، فقد ثبت أن استرخاء المخروط من الدرجة الثانية للاسترخاء شبه المحدد يكون بنفس الدقة تقريبًا ويحل بنجاح المشكلات التي تحتوي على n=1000 من المتغيرات المشتركة [ 14 ] .

التطبيقات

تحليل البيانات المالية

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

علم الأحياء

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

اختبار الفرضيات متعددة الأبعاد

غالباً ما تحتوي مجموعات البيانات المعاصرة على عدد من متغيرات الإدخال (ص{\displaystyle p}) مماثلة أو حتى أكبر بكثير من عدد العينات (ن{\displaystyle n}). لقد ثبت أنه إذاص/ن{\displaystyle p/n}إذا لم تتقارب إلى الصفر، فإن تحليل المكونات الرئيسية الكلاسيكي غير متسق . بعبارة أخرى، إذا افترضناك=ص{\displaystyle k=p}في المعادلة 1 ، فإن القيمة المثلى لا تتقارب مع أكبر قيمة ذاتية لمجموعة البيانات عندما يكون حجم العينةن{\displaystyle n\rightarrow \infty }ولا يتقارب الحل الأمثل في اتجاه التباين الأقصى. لكن يمكن لتحليل المكونات الرئيسية المتفرقة أن يحافظ على اتساقه حتى لوصن.{\displaystyle p\gg n.}

يمكن استخدام أكبر قيمة ذاتية متفرقة من الرتبة k (القيمة المثلى للمعادلة 1 ) للتمييز بين نموذج متساوي القياس، حيث يكون لكل اتجاه نفس التباين، ونموذج التغاير ذي التباين المتغير في بيئة عالية الأبعاد. [ 15 ] لنفترض اختبار فرضية حيث تنص الفرضية الصفرية على أن البياناتX{\displaystyle X}يتم توليدها من توزيع طبيعي متعدد المتغيرات بمتوسط ​​0 وتغاير يساوي مصفوفة الوحدة ، وتنص الفرضية البديلة على أن البياناتX{\displaystyle X}يتم توليدها من نموذج ذي قوة إشارة مرتفعةθ{\displaystyle \theta }:

ح0:Xشمال(0،أناص)،ح1:Xشمال(0،أناص+θvvتي)،{\displaystyle H_{0}:X\sim N(0,I_{p}),\quad H_{1}:X\sim N(0,I_{p}+\theta vv^{T}),}

أينvRص{\displaystyle v\in \mathbb {R} ^{p}}يحتوي على k إحداثيات غير صفرية فقط . يمكن لأكبر قيمة ذاتية متفرقة k أن تميز بين الفرضيتين إذا وفقط إذاθ>Θ(كسجل(ص)/ن){\displaystyle \theta >\Theta ({\sqrt {k\log(p)/n}})}.

بما أن حساب القيمة الذاتية للمصفوفة المتفرقة من الرتبة k هو مسألة صعبة حسابيًا (NP-hard)، يمكن تقريبها بالقيمة المثلى لتقريب البرمجة شبه المحددة ( المعادلة 3 ). في هذه الحالة، يمكننا التمييز بين الفرضيتين إذاθ>Θ(ك2سجل(ص)/ن){\displaystyle \theta >\Theta ({\sqrt {k^{2}\log(p)/n}})}. الإضافيك{\displaystyle {\sqrt {k}}}لا يمكن تحسين هذا المصطلح بواسطة أي خوارزمية أخرى ذات وقت متعدد الحدود إذا كانت فرضية الزمرة المزروعة صحيحة.

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

  • amanpg - حزمة R لتحليل المكونات الرئيسية المتفرقة باستخدام طريقة التدرج التقريبي للمشعب المتناوب [ 16 ]
  • elasticnet – حزمة R للتقدير المتفرق وتحليل المكونات الرئيسية المتفرق باستخدام الشبكات المرنة [ 17 ]
  • epca – حزمة R لتحليل المكونات الرئيسية الاستكشافي لمجموعات البيانات واسعة النطاق، بما في ذلك تحليل المكونات الرئيسية المتفرقة وتقريب المصفوفة المتفرقة. [ 18 ]
  • nsprcomp - حزمة R لتحليل المكونات الرئيسية المتفرقة و/أو غير السالبة بناءً على تكرارات القوة المحددة [ 19 ]
  • scikit-learn – مكتبة بايثون للتعلم الآلي تحتوي على Sparse PCA وتقنيات أخرى في وحدة التفكيك. [ 20 ]

ملحوظات

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

مراجع

  1. 1 2 ديميتريس بيرتسيماس؛ رايان كوري-رايت؛ جان بوفيليه (2020). "حل مشكلة تحليل المكونات الرئيسية المتفرقة واسعة النطاق للوصول إلى الأمثلية (شبه) القابلة للتحقق". arXiv : 2005.05195 [ math.OC ].
  2. أندرياس م. تيلمان؛ مارك إي. بفيتش (2013). "التعقيد الحسابي لخاصية التماثل المقيد، وخاصية الفضاء الصفري، والمفاهيم ذات الصلة في الاستشعار المضغوط". معاملات IEEE في نظرية المعلومات . 60 (2): 1248-1259 . arXiv : 1205.2081 . CiteSeerX 10.1.1.760.2559 . doi : 10.1109 /TIT.2013.2290112 . S2CID 2788088 .  
  3. مقدم، باباك؛ فايس، يائير؛ أفيدان، شاي (2007-09-07). شولكوف، برنارد؛ بلات، جون؛ هوفمان، توماس (محررون). التطورات في أنظمة معالجة المعلومات العصبية 19: وقائع مؤتمر 2006. مطبعة معهد ماساتشوستس للتكنولوجيا. الصفحات 915-922 . doi : 10.7551/mitpress/7503.001.0001 . ISBN  978-0-262-25691-9.
  4. زو، هوي؛ تريفور، هاستي؛ تيبشيراني، روبرت (1 يونيو 2006). "تحليل المكونات الرئيسية المتفرقة" . مجلة الإحصاءات الحاسوبية والرسومية . 15 (2): 265-286 . doi : 10.1198/106186006X113430 . ISSN 1061-8600 . 
  5. هوي زو؛ تريفور هاستي؛ روبرت تيبشيراني (2006). "تحليل المكونات الرئيسية المتفرقة" (ملف PDF) . مجلة الإحصاءات الحاسوبية والرسومية . 15 (2): 262-286 . CiteSeerX 10.1.1.62.580 . doi : 10.1198/106186006x113430 . S2CID 5730904 .  
  6. فان تشين؛ كارل روه (2021). "أساس جديد لتحليل المكونات الرئيسية المتفرقة" . مجلة الإحصاءات الحاسوبية والرسومية . 33 (2): 421-434 . arXiv : 2007.00596 . doi : 10.1080/10618600.2023.2256502 .
  7. 1 2 ألكسندر داسبريمون؛ لوران الغاوي؛ مايكل آي. جوردان؛ جيرت آر جي لانكريت (2007). "صياغة مباشرة لتحليل المكونات الرئيسية المتفرقة باستخدام البرمجة شبه المحددة" (ملف PDF) . مجلة SIAM Review . 49 (3): 434-448 . arXiv : cs/0406021 . doi : 10.1137/050645506 . S2CID 5490061 . 
  8. ميشيل جورني؛ يوري نيستيروف؛ بيتر ريشتاريك؛ رودولف سيبولكر (2010). "طريقة القوة المعممة لتحليل المكونات الرئيسية المتفرقة" (ملف PDF) . مجلة أبحاث تعلم الآلة . 11 : 517-553 . arXiv : 0811.4724 . Bibcode : 2008arXiv0811.4724J . ورقة نقاش CORE 2008/70.
  9. بيتر ريشتاريك؛ ماجد جهاني؛ إس. داملا أهيباساوغلو؛ مارتن تاكاك (2021). "التعظيم المتناوب: إطار موحد لثمانية صيغ تحليل المكونات الرئيسية المتفرقة ورموز متوازية فعالة" . الهندسة والتحسين . 22 (3): 1493-1519 . arXiv : 1212.4137 . doi : 10.1007/s11081-020-09562-3 . S2CID 2549610 . 
  10. باباك مقدم؛ يائير فايس؛ شاي أفيدان (2005). "الحدود الطيفية لتحليل المكونات الرئيسية المتفرقة: خوارزميات دقيقة وجشعة" (ملف PDF) . التطورات في أنظمة معالجة المعلومات العصبية . المجلد 18. مطبعة معهد ماساتشوستس للتكنولوجيا. 
  11. لورين بيرك؛ ديميتريس بيرتسيماس (2019). "تحليل المكونات الرئيسية المتفرقة الأمثل القابل للتحقق". الحوسبة في البرمجة الرياضية . 11 (3). سبرينغر: 381-420 . doi : 10.1007/s12532-018-0153-6 . hdl : 1721.1/131566 . S2CID 126998398 . 
  12. يو غوان؛ جينيفر داي (2009). "تحليل المكونات الرئيسية الاحتمالية المتفرقة" (ملف PDF) . وقائع ورشة عمل ومؤتمر مجلة أبحاث تعلم الآلة . 5 : 185.
  13. هوي زو؛ لينغتشو شيو (2018). "نظرة عامة انتقائية على تحليل المكونات الرئيسية المتفرقة" . وقائع معهد مهندسي الكهرباء والإلكترونيات . 106 (8): 1311-1320 . doi : 10.1109/jproc.2018.2846588 .
  14. ديميتريس بيرتسيماس؛ رايان كوري-رايت (2020). "حول تحليلات المخروط متعدد السطوح والتحليلات المخروطية من الدرجة الثانية لمسائل التحسين شبه المحددة" . رسائل بحوث العمليات . 48 (1). إلسيفير: 78-85 . arXiv : 1910.03143 . doi : 10.1016/j.orl.2019.12.003 .
  15. كوينتين بيرتيه؛ فيليب ريغوليه (2013). "الكشف الأمثل عن المكونات الرئيسية المتفرقة في الأبعاد العالية". حوليات الإحصاء . 41 (1): 1780-1815 . arXiv : 1202.5070 . doi : 10.1214/13-aos1127 . S2CID 7162068 . 
  16. https://cran.r-project.org/web/packages/amanpg/index.html
  17. https://cran.r-project.org/web/packages/elasticnet/index.html
  18. https://cran.r-project.org/web/packages/epca/index.html
  19. https://cran.r-project.org/web/packages/nsprcomp/index.html
  20. http://scikit-learn.org/stable/modules/generated/sklearn.decomposition.SparsePCA.html