مدل بردلی-تری

**مدل برادلی-تری** یک مدل احتمال برای نتیجه مقایسه‌های دوتایی بین آیتم‌ها، تیم‌ها یا اشیا است. با توجه به یک جفت آیتم  و  که از یک توزیع خاص گرفته شده است، این احتمال را تخمین می‌زند که مقایسه زوجی  درست باشد، این فرض به احتمال زیر صحیح است:بطوریکه  یک امتیاز عدد حقیقی مثبت است که به هر یک از ها اختصاص داده شده است. مقایسه  می‌توان به صورت "  به  ترجیح داده شده"، "رتبه  بالاتر از  است"، یا "  , را می‌زند " بسته به کاربرد، خواند و با استفاده از "قضاوت" امکان ارزیابی ذهنی را فراهم می‌کند.

به عنوان مثال،  ممکن است نشان دهنده مهارت یک تیم در یک تورنمنت ورزشی باشد و P(i>j) به این معنای احتمال برد  در بازی مقابل  است. یا  ممکن است نشان دهنده کیفیت یا مطلوبیت یک محصول تجاری باشد و P(i>j) احتمال این است که مصرف‌کننده محصول  را بر محصول  ترجیح می‌دهد.

مدل بردلی-تری می‌تواند در جهت رو به جلو برای پیش‌بینی نتایج، همان‌طور که توضیح داده شد، استفاده شود، اما معمولاً به صورت معکوس برای استنتاج امتیازات  با توجه به مجموعه‌ای از نتایج مشاهده شده استفاده می‌شود. در این نوع کاربرد  نشان دهنده مقداری از قدرت یا کیفیت i است و این مدل به ما امکان می‌دهد نقاط قوت را از یک سری مقایسه‌های زوجی تخمین بزنیم. به عنوان مثال، در بررسی ترجیحات شراب، ممکن است برای پاسخ دهندگان دشوار باشد که رتبه‌بندی کاملی از مجموعه بزرگی از شراب‌ها را ارائه دهند، اما برای آنها نسبتاً آسان است که جفت‌های نمونه شراب را با هم مقایسه کنند و بگویند که کدام یک از شراب‌ها بهتر است. بر اساس مجموعه‌ای از این مقایسه‌های زوجی، می‌توان از مدل بردلی-تری برای به دست آوردن رتبه‌بندی کامل شراب‌ها استفاده کرد.

هنگامی که مقادیر امتیاز  محاسبه شد، می‌توان از مدل در جهت رو به جلو نیز استفاده کرد، به عنوان مثال برای پیش‌بینی نتیجه احتمالی مقایسه‌هایی که هنوز واقعاً رخ نداده‌اند. به عنوان مثال، در مثال بررسی شراب، می‌توان احتمال ترجیح شراب i  بر شراب j را محاسبه کرد، حتی اگر هیچ‌کس در نظرسنجی مستقیماً آن جفت خاص را مقایسه نکرده باشد.

## تاریخچه و کاربردها
این مدل به افتخار رالف آ. بردلی و میلتون ای. تری، که آن را در سال ۱۹۵۲ ارائه کردند، نامگذاری شده است، اگرچه قبلاً توسط ارنست زرملو در دهه ۱۹۲۰ مورد مطالعه قرار گرفته بود. کاربردهای این مدل شامل رتبه‌بندی رقبا در مسابقات ورزشی، شطرنج و سایر مسابقات، رتبه‌بندی محصولات در نظرسنجی‌های مقایسه زوجی انتخاب مصرف‌کننده، تجزیه و تحلیل سلسله مراتب سلطه در جوامع حیوانی و انسانی، رتبه‌بندی مجلات، رتبه‌بندی مدل‌های هوش مصنوعی، و برآورد ارتباط اسناد در موتورهای جستجوی ماشینی و غیره است.

## تعریف
مدل بردلی-تری را می‌توان به روش‌های مختلفی پارامتری کرد. معادله () شاید رایج‌ترین روش باشد، اما تعدادی دیگر نیز وجود دارد. بردلی و تری خود توابع امتیاز نمایی، p_i = e^{\beta_i} ، را تعریف کرده‌اند به طوری که

: P(i > j) = \frac{e^{\beta_i}}{e^{\beta_i} + e^{\beta_j}}.

به‌طور معادل می‌توان از یک logit استفاده کرد، به طوری که

: \operatorname{logit} P(i > j) = \log \frac{P(i > j)}{1 - P(i > j)} = \log \frac{P(i > j)}{P(j > i)} = \beta_i - \beta_j,

یعنی  \operatorname{logit} p = \log\frac p{1-p}  برای  0

این فرمول شباهت بین مدل بردلی-تری و رگرسیون لجستیک را برجسته می‌کند. هر دو اساساً از یک مدل واحد، اما به روش‌های مختلف استفاده می‌کنند. در رگرسیون لجستیک فرد معمولاً پارامترهای \beta_i را می‌داند و تلاش می‌کند تا شکل عملکردی P(i>j) را استنباط کند؛ در رتبه‌بندی تحت مدل بردلی-تری، فرد شکل عملکردی را می‌شناسد و سعی می‌کند پارامترها را استنتاج کند.

### تخمین پارامترها
رایج‌ترین کاربرد مدل برادلی-تری، استنتاج مقادیر پارامترهای p_i با توجه به مجموعه ای از نتایج مشاهده شده i>j است مانند برد و باخت در یک مسابقه. ساده‌ترین راه برای تخمین پارامترها ، تخمین حداکثر درستنمایی است، یعنی با به حداکثر رساندن احتمال نتایج مشاهده شده با توجه به مقادیر مدل و پارامتر.

فرض کنید نتایج مجموعه ای از رقابت‌های دوتایی بین گروه خاصی از افراد را می‌دانیم و اجازه دهید  تعداد دفعاتی باشد که  فردی  را شکست می‌دهد. سپس احتمال این مجموعه از نتایج در مدل بردلی-تری\prod_{ij} [P(i>j)]^{w_{ij}} است و درستنمایی پارامتر  برابر است با

: 
\begin{align}
L(\mathbf{p}) & = \ln \prod_{ij} [P(i>j)]^{w_{ij}} = \sum_{i=1}^n \sum_{j=1}^n \ln \left(\frac{p_i}{p_i+p_j}\right)^{w_{ij}} \\[6pt]
& = \sum_{ij} w_{ij} \ln \frac{p_i}{p_i+p_j} = \sum_{ij} [w_{ij} \ln p_i - w_{ij} \ln(p_i + p_j)].
\end{align}


Zermelo نشان داد که این عبارت تنها دارای یک ماکزیمم واحد است که با اگر نسبت به p_i مشتق گرفته و نتیجه را مساوی صفر قرار دهیم خواهیم داشت:{\sum_{j} (w_{ij}+w_{ji})/ (p_i+p_j)}.|}}این معادله هیچ راه حل بسته شناخته شده‌ای ندارد، اما زرملو حل آن را با تکرار ساده پیشنهاد کرد. با شروع از هر مجموعه مناسبی از مقادیر اولیه (مثبت) برای p_i ، یکی به‌طور مکرر به روز رسانی را انجام می‌دهد{\sum_{j} (w_{ij}+w_{ji})/ (p_i+p_j)}|}}برای همه ‌ها به نوبت این مقادیر محاسبه می‌شود. پارامترهای به دست آمده تا یک ضریب ثابت، دلخواه هستند، بنابراین پس از محاسبه همه مقادیر جدید باید با تقسیم بر میانگین هندسی آنها نرمال سازی شوند:.|}}این فرایند تخمین، احتمال ورود به سیستم را در هر تکرار بهبود می‌بخشد و تضمین می‌شود که در نهایت به حداکثر منحصر به فرد برسد. با این حال، همگرایی کند است. اخیراً اشاره شده است که معادله () را نیز می‌توان به صورت بازآرایی کرد.

: p_i = \frac{\sum_{j} w_{ij} p_j/(p_i+p_j)}{\sum_{j} w_{ji}/(p_i+p_j)},

که با تکرار قابل حل استپس از هر دور به روز رسانی با استفاده از معادله () دوباره عادی می‌شود. این تکرار نتایج یکسانی را با خروجی معادله () می‌دهد، اما خیلی سریعتر همگرا می‌شود و از این رو معمولاً بر () ترجیح داده می‌شود.

#### نمونه کار شده از روش حل
یک رقابت ورزشی بین چهار تیم را در نظر بگیرید که در مجموع ۲۲ بازی را بین خود انجام می‌دهند. بردهای هر تیم در ردیف‌های جدول زیر و حریفان به صورت ستون آورده شده‌اند:

به عنوان مثال، تیم A دو بار تیم B را شکست داده و سه بار به تیم B باخته است و با تیم C اصلاً بازی نکرده است. یک برد و چهار باخت مقابل تیم D دارد.

ما می‌خواهیم نقاط قوت نسبی تیم‌ها را تخمین بزنیم که این کار را با محاسبه پارامترهای p_i انجام می‌دهیم، بطوریکه پارامترهای بالاتر نشان دهنده مهارت بیشتر است. برای انجام این کار، چهار ورودی را در بردار پارامتر  به‌طور دلخواه مقداردهی می‌کنیم، به عنوان مثال مقدار ۱ را به هر تیم اختصاص می‌دهیم: . سپس معادله () را برای به روز رسانی p_1 اعمال می‌کنیم، که نتیجه زیر حاصل می‌شود:p_1 = \frac{\sum_{j(\ne 1)} w_{1j} p_j/(p_1+p_j)}{\sum_{j(\ne 1)} w_{j1}/(p_1+p_j)} =
\frac{2\frac{1}{1+1} + 0\frac{1}{1+1} + 1\frac{1}{1+1}}{3\frac{1}{1+1}+0\frac{1}{1+1} + 4\frac{1}{1+1}} = 0.429.اکنون () را دوباره برای به روز رسانی p_2 اعمال می‌کنیم، مطمئن شوید که از مقدار جدید p_1 که محاسبه شد استفاده می‌کنید:p_2 = \frac{\sum_{j(\ne 2)} w_{2j} p_j/(p_2+p_j)}{\sum_{j(\ne 2)} w_{j2}/(p_2+p_j)} =
\frac{3\frac{0.429}{1+0.429} + 5\frac{1}{1+1} + 0\frac{1}{1+1}}{2\frac{1}{1+0.429}+3\frac{1}{1+1} + 0\frac{1}{1+1}} = 1.172به‌طور مشابه برای p_3 و p_4 داریم:p_3 = \frac{\sum_{j(\ne 3)} w_{3j} p_j/(p_3+p_j)}{\sum_{j(\ne 3)} w_{j3}/(p_3+p_j)} =
\frac{0\frac{0.429}{1+0.429} + 3\frac{1.172}{1+1.172} + 1\frac{1}{1+1}}{0\frac{1}{1+0.429}+5\frac{1}{1+1.172} + 3\frac{1}{1+1}} = 0.557p_4 = \frac{\sum_{j(\ne 4)} w_{4j} p_j/(p_4+p_j)}{\sum_{j(\ne 4)} w_{j4}/(p_4+p_j)} =
\frac{4\frac{0.429}{1+0.429} + 0\frac{1.172}{1+1.172} + 3\frac{0.557}{1+0.557}}{1\frac{1}{1+0.429}+0\frac{1}{1+1.172} + 1\frac{1}{1+0.557}} = 1.694سپس تمام پارامترها را با تقسیم بر میانگین هندسی آنها نرمال می‌کنیم (0.429\times1.172\times0.557\times1.694)^{1/4} = 0.830 برای بدست آوردن پارامترهای تخمینی .

برای بهبود بیشتر تخمین‌ها، با استفاده از مقادیر  جدید، فرایند را تکرار می‌کنیم؛ مثلاً، p_1 =
\frac{2\cdot\frac{1.413}{0.516+1.413} + 0\cdot\frac{0.672}{0.516+0.672} + 1\cdot\frac{2.041}{0.516+2.041}}{3\cdot\frac{1}{0.516+1.413}+0\cdot\frac{1}{0.516+0.672} + 4\cdot\frac{1}{0.516+2.041}} = 0.725.با تکرار این فرایند برای پارامترهای باقیمانده و عادی سازی،  را دریافت می‌کنیم. تکرار ۱۰ بار دیگر همگرایی سریع به سمت حل نهایی  می‌دهد. این نشان می‌دهد که تیم D قوی‌ترین و تیم B دومین قوی‌تر است، در حالی که تیم‌های A و C تقریباً از نظر قدرت برابر هستند اما کمتر از تیم‌های B و D هستند. حتی اگر همه تیم‌ها با هم بازی نکرده باشند.

## جستارهای وابسته
* رگرسیون ترتیبی
* مدل راش
* مقیاس (علوم اجتماعی)
* سیستم رتبه‌بندی الو
* مدل تورستونی

## منابع

رده:رگرسیون لجستیک
رده:صفحه‌های با ترجمه بازبینی‌نشده
رده:مدل‌های آماری
رده:مدل‌های رگرسیون
رده:یادگیری ماشین