|  |
| --- |
| بخشی از مجموعه مباحث دربارهٔ آمار |
| تحلیل رگرسیون |
| 
مدل‌ها |
| * رگرسیون خطی
* رگرسیون خطی ساده
* رگرسیون چندجمله‌ای
* رگرسیون چندمتغیره
 |
| * مدل خطی تعمیم‌یافته
* انتخاب گسسته
* رگرسیون لجستیک
* لوجیت چندجمله‌ای
* لوجیت آمیخته
* مدل پروبیت
* پروبیت چندجمله‌ای
* لوجیت مرتب
* پروبیت مرتب
* رگرسیون پواسون
 |
| * مدل سلسله‌مراتبی
* مدل اثرات ثابت
* مدل اثرهای تصادفی
* مدل آمیخته
 |
| * رگرسیون غیرخطی
* رگرسیون غیرپارامتریک
* رگرسیون نیمه‌پارامتریک
* رگرسیون باثبات
* رگرسیون چندک
* رگرسیون ایزوتونیک
* رگرسیون مولفه اصلی
* رگرسیون کمترین زاویه
* رگرسیون موضعی
* رگرسیون مقطع
 |
| * مدل خطا در متغیرها
 |
| 
تخمین |
| * کمترین مربعات
* کمترین مربعات خطی
* کمترین مربعات غیرخطی
 |
| * حداقل مربعات معمولی
* حداقل مربعات وزن‌دار
* روش تعمیم‌یافته کمترین مربعات
 |
| * رگرسیون پاره‌ای کمتری مربعات
* مجموع کمترین مربعات
* کمترین مربعات نامنفی
* تنظیم تیخونوف
* کمترین مربعات منظم
 |
| * کمترین انحرافات مطلق
* کمترین مربعات بازوزن‌داده مکرر
* رگرسیون خطی بیزی
* رگرسیون چندمتغیره خطی بیزی
 |
| 
پیش‌زمینه |
| * اعتبارسنجی مدل رگرسیون
* پاسخ میانگین و پیش‌بینی‌شده
* خطاها و باقی‌مانده‌ها در آمار
* نیکویی برازش
* باقی‌مانده استودنت‌شده
* قضیه گوس-مارکف
 |
| * نشان درگاه درگاه آمار
 |
| * ن
* ب
* و
 |

**حداقل مربعات جزئی**  (انگلیسی: **Partial least squares regression** (**PLS regression**)) یکی از روش‌های تحلیل رگرسیون است. در این روش، راه حل حداقل مربعات بر روی تعدادی مؤلفه متعامد که ترکیبی خطی از متغیرهای مستقل هستند و به صورت متناوب و با هدف بیشینه‌کردن کوواریانسِ تبدیل خطیِ متغیرهای مستقل و متغیرهای وابسته ایجاد شده‌اند، اعمال می‌شود.

حداقل مربعات جزئی توسط آماردان سوئدی، Herman O. A. Wold معرفی‌شد که سپس آن را با پسرش Svante Wold توسعه‌داد. یک اصطلاح جایگزین برای «حداقل مربعات جزئی»، **نگاشت به ساختارهای پنهان** است، اما هوز «حداقل مربعات جزئی» در بسیاری از حوزه ها اصطلاح رایج است. اگرچه کاربردهای اصلی آن در علوم اجتماعی بود، اما امروزه PLS بیشترین کاربرد را در شیمیِ آماری و حوزه های مرتبط دارد. همچنین در بیوانفورماتیک، حس سنجی، علوم اعصاب و انسان‌شناسی نیز استفاده می‌شود.

رابطه ریاضی
-----------

فرض می‌کنیم {\displaystyle n} داده داریم که هر کدام از یک متغیر وابسته و {\displaystyle p} متغیر مستقل تشکیل شده‌است؛ به کمک متغیرهای مستقل متغیر وابسته را پیش‌بینی می‌کنیم. مقادیر متغیرهای مستقل (به همراه یک بردار ثابت 1) و مقادیر متغیر وابسته را به ترتیب در ماتریس‌های {\displaystyle \mathbf {X} \_{n\times (p+1)}} و {\displaystyle \mathbf {Y} \_{n\times 1}} به شکل پایین نمایش می‌دهیم؛ در اینجا {\displaystyle \mathbf {X} \_{i}} یک ماتریس {\displaystyle n\times 1} از مقادیر {\displaystyle i} امین متغیر مستقل است:

{\displaystyle \mathbf {X} \_{n\times (p+1)}=\left(\mathbf {1} ,\mathbf {X} \_{1},\ldots ,\mathbf {X} \_{n}\right)^{\top }}
{\displaystyle \mathbf {Y} \_{n\times 1}=\left(y\_{1},\ldots ,y\_{n}\right)^{\top }}
هدف از رگرسیون خطی بدست آوردن پارامتر {\displaystyle {\boldsymbol {\beta }}\in \mathbb {R} ^{p+1}} است به شکلی که {\displaystyle \mathbf {Y} =\mathbf {X} {\boldsymbol {\beta }}+{\boldsymbol {\varepsilon }}} در اینجا {\displaystyle \varepsilon } یک متغیر تصادفی است که خطای مدل را نشان می‌دهد. این خطا از توزیع طبیعی با میانگین صفر و واریانس ثابت برای تمامی ابعاد پیروی می‌کند به این معنی که {\displaystyle \operatorname {E} \left({\boldsymbol {\varepsilon }}\right)=\mathbf {0} \;} و {\displaystyle \;\operatorname {Var} \left({\boldsymbol {\varepsilon }}\right)=\sigma ^{2}I\_{n\times n}}. از طریق روش حداقل مربعات می‌توان {\displaystyle ||\mathbf {Y} -\mathbf {X} {\boldsymbol {\beta }}||^{2}} یعنی مربع میزان خطاها را کمینه کرد و به پارامتر بهینه رسید. این پارامتر با {\displaystyle {\widehat {\boldsymbol {\beta }}}\_{\mathrm {ols} }=(\mathbf {X} ^{\top }\mathbf {X} )^{-1}\mathbf {X} ^{\top }\mathbf {Y} } برابراست. یکی از مشکلات اصلی این روش عدم وارونپذیری {\displaystyle (\mathbf {X} ^{\top }\mathbf {X} )^{-1}} است. برای حل این مشکل، حداقل مربعات جزئی، متناوباً مولفه‌هایی متعامد ({\displaystyle \mathbf {Z\_{m}} } در پایین) که ترکیبی خطی از متغیرهای مستقل هستند تولید می‌کند و در نهایت راه حل حداقل مربعات را بر روی این مولفه‌ها اعمال می‌کند. ضرایب متغیرهای مستقل در مولفه‌ها با ضرب نقطه‌ایِ متغیرهای وابسته و مستقل برابر است. متغیرهای مستقل در ابتدای کار استاندارد شده‌اند یعنی میانگین صفر و واریانس یک دارند. در پایان هر مرحله متغیرهای مستقل نسبت به مولفه آن مرحله متعامد می‌شوند. این کار باعث می‌شود که در پایان تمام مولفه‌ها نسبت به هم متعامد باشند. الگوریتم تولید مولفه‌ها برای {\displaystyle k\leq p} به شکل پایین است ({\displaystyle k} از روش اعتبارسنجی متقابل محاسبه می‌شود):

{\displaystyle {\begin{cases}\,\,{\mbox{ for }}\,\,j\,\,=\,\,1,\cdots ,p:\\\,\,\,\,\,\,\,\,\,\,\mathbf {X} \_{j}^{(0)}=\mathbf {X} \_{j}\\\,\,{\mbox{ for }}\,\,m=\,\,1,\cdots ,k:\\\,\,\,\,\,\,\,\,\,\,\mathbf {Z} \_{m}=\sum \_{j=1}^{p}\left\langle \mathbf {X} \_{j}^{(m-1)},\mathbf {Y} \right\rangle \mathbf {X} \_{j}^{(m-1)}\\\,\,\,\,\,\,\,\,\,{\mbox{ for }}\,\,j\,\,=\,\,1,\cdots ,p:\\\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\mathbf {X} \_{j}^{(m)}=\mathbf {X} \_{j}^{(m)}-{\frac {\langle \mathbf {Z} \_{m}\,,\,\mathbf {Y} \rangle }{\langle \mathbf {Z} \_{m}\,,\,\mathbf {Z} \_{m}\rangle }}\mathbf {Z} \_{m}\\\,\,{\hat {\mathbf {Y} }}={\bar {\mathbf {Y} }}\mathbf {1} +\sum \_{m=1}^{k}{\frac {\langle \mathbf {Z} \_{m}\,,\,\mathbf {Y} \rangle }{\langle \mathbf {Z} \_{m}\,,\,\mathbf {Z} \_{m}\rangle }}\mathbf {Z} \_{m}\\\end{cases}}}

از آنجا که {\displaystyle {\hat {\mathbf {Y} }}} یا همان مقدار پیش‌بینی شده، ترکیبی خطی از {\displaystyle \mathbf {Z} \_{m}}‌ها است و خود {\displaystyle \mathbf {Z} \_{m}}‌ها هم ترکیبی خطی از متغیرهای مستقل هستند، در نهایت مدل رگرسیون ترکیبی خطی از متغیرهای مستقل خواهد بود.

الگوریتم PLS1
-------------

PLS1 یک الگوریتم پرکاربرد است که وقتی از آن استفاده می‌شود که Y یک بردار باشد. این الگوریتم در شبه‌کد زیر بیان می‌شود (حروف بزرگ ماتریس هستند، حروف کوچک اگر بالانویس داشته باشند بردار، و اگر اندیس داشته‌باشند اسکالر هستند):

```
 1 **function** PLS1(X, y, l)
 2     {\displaystyle X^{(0)}\gets X}
 3     {\displaystyle w^{(0)}\gets X^{\mathrm {T} }y/||X^{\mathrm {T} }y||}, an initial estimate of w.
 4     **for** {\displaystyle k=0} **to** {\displaystyle l-1}
 5         {\displaystyle t^{(k)}\gets X^{(k)}w^{(k)}}
 6         {\displaystyle t\_{k}\gets {t^{(k)}}^{\mathrm {T} }t^{(k)}} (note this is a scalar)
 7         {\displaystyle t^{(k)}\gets t^{(k)}/t\_{k}}
 8         {\displaystyle p^{(k)}\gets {X^{(k)}}^{\mathrm {T} }t^{(k)}}
 9         {\displaystyle q\_{k}\gets {y}^{\mathrm {T} }t^{(k)}} (note this is a scalar)
10         **if** {\displaystyle q\_{k}=0}
11             {\displaystyle l\gets k}, **break** the **for loop**
12         **if** {\displaystyle k<(l-1)}
13             {\displaystyle X^{(k+1)}\gets X^{(k)}-t\_{k}t^{(k)}{p^{(k)}}^{\mathrm {T} }}
14             {\displaystyle w^{(k+1)}\gets {X^{(k+1)}}^{\mathrm {T} }y}
15     **end** **for**
16     **define** W to be the matrix with columns {\displaystyle w^{(0)},w^{(1)},...w^{(l-1)}}.
       Do the same to form the P matrix and q vector.
17     {\displaystyle B\gets W{(P^{\mathrm {T} }W)}^{-1}q}
18     {\displaystyle B\_{0}\gets q\_{0}-{P^{(0)}}^{\mathrm {T} }B}
19     **return** {\displaystyle B,B\_{0}}
```

این شکل از الگوریتم نیازی به مرکزیت ورودی X و Y (کم کردن میانگین هر ستون از تک تک درایه‌های آن ستون در ماتریس) ندارد، زیرا این کار به طور ضمنی توسط الگوریتم انجام می‌شود.متغیر *l* (فراهم‌شده توسط کاربر) محدودیتی برای تعداد متغییرهای پنهان در رگرسیون است و اگر با رتبه ماتریس *X* برابر باشد، الگوریتم تخمین‌های رگرسیون کمترین مربعات را برای *B و* {\displaystyle B\_{0}} به دست می‌دهد.

جستارهای وابسته
---------------

تحلیل مولفه‌های اصلی

تحلیل واریانس

رگرسیون خطی

یادداشت‌ها
----------

1. ↑  partial least squares

منابع
-----

1. 1 2  Hastie, Trevor; Tibshirani, Robert; Friedman, Jerome (2009). "The Elements of Statistical Learning". *Springer Series in Statistics* (به انگلیسی): 80–82. doi:10.1007/978-0-387-84858-7. ISSN 0172-7397.
2. ↑  de Jong, S.; ter Braak, C.J.F. (1994). "Comments on the PLS kernel algorithm". *J. Chemometrics*. **8** (2): 169–174. doi:10.1002/cem.1180080208.