توزیع مزدوج پیشین

توزیع مزدوج پیشین (انگلیسی: Conjugate prior)
در تئوری احتمال بیزی، اگر توزیعات احتمال پسین با توزیع احتمال پیشین p () خانواده‌ای یکسان داشته باشد آن گاه به توزیعات پیشین و پسین توزیعات مزدوج گفته شده و توزیع پیشین را توزیع مزدوج پیشین برای تابع شباهت می‌نامند. به عنوان مثال، خانواده گاوسی مزدوج خودش (یا خود مزدوج) نسبت به یک تابع شباهت گاوسی است: اگر تابع شباهت گوسی باشد، انتخاب توزیع پیشین گاوسی برای میانگین تضمین می‌کند که توزیع پسین نیز گاوسی خواهد بود. این به این معناست که توزیع گاوسی، تابع مزدوج پیشین تابع شباهتی گوسی است. این مفهوم و عبارت «مزدوج پیشین» توسط هاوارد رایفا و رابرت اشلایفر در کارشان روی تئوری انتخاب بیزی معرفی شدند. مفهوم مشابهی نیز به‌طور مستقل توسط جورج آلفرد بارنارد کشف شده بود.
مسئله کلی یافتن توزیع برای یک پارامتر  به شرط یک داده یا مجموعه داده‌های x را در نظر بگیرید. از قضیه بیز می‌دانیم که توزیع پسین برابر با حاصلضرب تابع شباهت و توزیع پیشین نرمالیزه شده با (تقسیم بر) احتمال داده است:
فرض می‌کنیم تابع شباهت ثابت است. تابع شباهت معمولا به خوبی از بیانی از فرایند تولید داده مشخص می‌گردد. واضح است که انتخاب‌های متفاوت برای توزیع پیشین p() محاسبه انتگرال را سخت‌تر یا آسان‌تر می‌کند و حاصل ضرب p (x|) × p () فرم جبری متفاوتی بگیرد. برای انتخاب‌های مشخصی از توزیع پیشین، توزیع پسین فرم جبری مشابه توزیع پیشین خواهد داشت (عموما با مقادیر پارامترهای متفاوت). چنین انتخابی یک توزیع مزدوج پیشین است.
یک توزیع مزدوج پیشین یک آسودگی جبری است که توصیف بسته‌ای برای توزیع پسین ارائه می‌دهد. در غیر این صورت ممکن است انتگرال عددی لازم باشد. علاوه بر این، توزیعات پیشین مزدوج با بیشتر و شفاف نشان دادن اینکه چگونه یک تابع شباهت یک توزیع پیشین را به روز می‌کند درکی مستقیم ارائه می‌دهند.
تمام اعضای خانواده نمایی توزیعات پیشین مزدوج دارند.
مثال
فرم توزیع پیشین مزدوج عموما می‌تواند با بررسی تابع چگالی احتمال یا تابع جرم احتمال آن توزیع مشخص گردد. به عنوان مثال، یک متغیر تصادفی را در نظر بگیرید که شامل تعداد موفقیت‌های در آزمایش برنولی با احتمال موفقیت نامعلوم در بازه ۰،۱ است. این متغیر تصادفی با تابع جرم احتمال زیر از توزیع دوجمله‌ای تبعیت می‌کند:
توزیع مزدوج پیشین معمول توزیع بتا با پارامترهای ( , ):
که در آن و برای نمایش هرگونه باور یا اطلاعات به کار می‌رود (۱ = و ۲ = توزیع یکنواخت را به دست می‌دهند) و B(, ) تابع بتا است که به عنوان ثابت نرمالیزه کننده استفاده می‌شود.
در اینجا و هایپرپارامتر خوانده می‌شوند (پارامترهای توزیع پیشین) تا از پارامترهای مدل اساسی (در اینجا q) متمایز گردند. یکی از ویژگی‌های معمول توزیعات پیشین مزدوج این است که ابعاد هایپرپارامترها یک واحد بزرگتر از ابعاد پارامترهای توزیع اصلی‌اند. اگر همه پارامترها مقادیر اسکالر باشند آنگاه به این معناست که یک هایپرپارامتر بیشتر از پارامتر خواهیم داشت. البته این نکته درباره پارامترهای با مقادیر برداری و ماتریسی نیز صدق می‌کند. (به مقاله کلی درباره خانواده نمایی رجوع کنید و توزیع ویشارت را که مزدوج پیشین ماتریس کوواریانس توزیع نرمال چند متغیره است به عنوان مثالی برای ابعاد بزرگ در نظر بگیرید)
آن گاه اگر از این متغیر تصادفی نمونه برداری کنیم و s موفقیت و p شکست به دست آوریم داریم:
\begin{align
P(s, f \mid q=x) -= {s+f \choose s x^s(1-x)^f,\\
P(x) -= {x^{\alpha-1(1-x)^{\beta-1 \over \Beta(\alpha,\beta),\\
P(q=x \mid s,f) -= \frac{P(s, f \mid x)P(x){\int P(s, f \mid x)P(x)dx\\
- = \over \int_{y=0^1 \left({s+f \choose s y^{s+\alpha-1(1-y)^{f+\beta-1 / \Beta(\alpha,\beta)\right) dy \\
- = {x^{s+\alpha-1(1-x)^{f+\beta-1 \over \Beta(s+\alpha,f+\beta), \\
\end{align
که توزیع بتای دیگری با پارامترهای (s + , + f) است. این توزیع پسین نیز می‌تواند به عنوان توزیع پیشین نمونه‌های بیشتر با اضافه کردن هر چه بیشتر اطلاعات توسط هایپرپارامترها مورد استفاده قرار گیرد.
شبه مشاهدات
اغلب سودمند است که هایپرپارامترهای توزیع مزدوج پیشین را متناظر با انجام تعدادی شبه مشاهده با ویژگی‌هایی که توسط پارامترها مشخص می‌شود در نظر بگیریم. به عنوان مثال، مقادیر و برای یک توزیع بتا را می‌توان متناظر با ۱ - موفقیت و ۱ - شکست در نظر گرفت اگر مد توزیع پسین برای انتخاب یک تنظیمات پارامتر بهینه انتخاب شود یا موفقیت و شکست اگر میانگین توزیع پسین برای انتخاب یک تنظیمات پارامتر بهینه انتخاب شود. در حالت کلی، برای تقریبا تمام توزیعات مزدوج پیشین هایپر پارامترها می‌توانند بر حسب شبه مشاهدات تفسیر شوند. این عمل می‌تواند هم به فراهم آوردن درکی برای معادلات به‌روزرسانی اغلب درهم و شلوغ و هم برای انتخاب هایپرپارامترهای معقول برای یک توزیع پیشین کمک می‌کند.
تفسیرها
شباهت با توابع ویژه
توزیعات مزدوج پیشین مشابه تابع ویژه در نظریه عمل‌گرها هستند از این جهت که توزیعاتی هستند که در آن‌ها «عمل‌گر شرطی سازی»، در صورت در نظر گرفتن فرایند تغییر از توزیع پیشین به پسین به عنوان یک عمل‌گر، به روشی کاملا قابل درک عمل می‌کند.
در هر دو توابع ویژه و توزیعات مزدوج پیشین فضایی با بعد متناهی وجود دارد که توسط عمل‌گر حفظ و نگهداری می‌شود: خروجی فرمی مشابه (فضایی مشابه) ورودی دارد. این عمل به شدت تجزیه و تحلیل را ساده می‌کند چون در غیر این صورت فضایی با بعد نامتناهی (فضای تمام توابع، فضای تمام توزیع‌ها) را در نظر می‌گیرد.
با این وجود، فرایندها تنها مشابه و نه یکسان هستند: شرطی سازی خطی نیست زیرا فضای توزیعات تحت ترکیب خطی بسته نمی‌شود بلکه تحت ترکیب محدب بسته می‌شود و توزیع پسین تنها فرمی مشابه توزیع پیشین دارد و ضریبی اسکالر از آن نیست.
همان‌طور که به راحتی می‌توان تجزیه و تحلیل کرد که چگونه ترکیبی خطی از توابع ویژه تحت استفاده از یک عمل‌گر تکامل می‌یابد (چون عمل‌گر نسبت به این توابع، قطری شده است)، به راحتی می‌توان تجزیه و تحلیل کرد که چگونه یک ترکیب محدب از توزیعات مزدوج پیشین تحت شرط سازی کامل می‌یابد؛ این عمل به وسیله یک توزیع هایپرپیشین صورت می‌گیرد و متناظر با استفاده از چگالی مخلوط توزیعات مزدوج پیشین به جای تک توزیع مزدوج پیشین است.
سیستم دینامیکی
می‌توان عمل شرطی سازی روی توزیعات مزدوج پیشین را تعریف کردن نوعی سیستم دینامیکی (زمان گسسته) در نظر گرفت: به ازای مجموعه‌ای داده شده‌ای از هایپرپارامترها، داده ورودی این هایپرپارامترها را به روز می‌کند. پس می‌توان تغییر در هایپرپارامترها را به نوعی «تکامل زمانی» تعبیر کرد که با «یادگیری» متناظر است. شروع در نقاط متفاوت جریان‌هایی متفاوت در طول زمان را سبب می‌گردد. این نیز مشابه سیستم دینامیکی است که با یک عمل‌گر خطی تعریف می‌شود ولی توجه کنید چون نمونه‌های متفاوت به استنتاج‌های متفاوت ختم می‌شود، این عمل به تنهایی وابسته به زمان نیست بلکه بیشتر وابسته به داده‌است. برای رویکردهای مرتبط، تخمین بیزی بازگشتی و داده گواری را ببینید.
جدول توزیعات مزدوج
فرض کنید n نشان دهنده تعداد مشاهدات باشد. در تمام موارد زیر، فرض می‌شود که داده از n نقطه (که در موارد چند متغیره بردارهای تصادفی خواهند بود) تشکیل شود.
اگر تابع شباهت متعلق به خانواده نمایی باشد آن گاه توزیع مزدوج پیشین وجود دارد که اغلب اوقات آن نیز عضو خانواده نمایی است. به رجوع کنید.
توزیعات گسسته
شباهت
پارامترهای مدل
توزیع مزدوج پیشین
هایپرپارامترهای پیشین
هایپرپارامترهای پسین
تفسیر هایپر پارامترها
توزیع پسین پیشگویانه
برنولی
p (احتمال)
بتا
موفقیت، شکست
دوجمله‌ای
p (احتمال)
بتا
موفقیت، شکست
(بتا-دوجمله‌ای)
دوجمله‌ای منفیبا تعداد شکست‌های معلوم، r
p (احتمال)
بتا
total موفقیت، شکست (یعنی، آزمایش، با فرض اینکه ثابت بماند)
پواسون
(نرخ)
گاما
(دوجمله‌ای منفی)
تعداد رخدادها در بازه
(دوجمله‌ای منفی)
مطلق
p (بردار احتمال), k (تعداد دسته‌بندی‌ها؛ یعنی، اندازه p )
دیریکله
که تعداد مشاهدات در دسته‌بندی i است
رخداد در دسته‌بندی
\begin{align
p(\tilde{x=i) -= \frac{\sum_i {\alpha_i \\
-= \frac{\alpha_i + c_i{\sum_i \alpha_i + n
\end{align
چندجمله‌ای
p (بردار احتمال), k (تعداد دسته‌بندی‌های؛ یعنی، اندازه p )
دیریکله
رخدادهای دسته‌بندی
(دیریکله—چندجمله‌ای)
فوق‌هندسیبا اندازه جمعیت کل معلوم، N
M (تعداد اعضای هدف)
بتا-دوجمله‌ای
موفقیت، شکست
هندسی
p 0 (احتمال)
بتا
آزمایش، کل شکست‌ها
توزیعات پیوسته
شباهت
پارامترهای مدل
توزیع مزدوج پیشین
هایپرپارامترهای پیشین
هایپرپارامترهای پسین
تفسیر هایپرپارامترها
توزیع پسین پیشگویانه
نرمالبا واریانس معلوم  2
(میانگین)
نرمال
\frac{1{\frac{1{\sigma_0^2 + \frac{n{\sigma^2\left(\frac{\mu_0{\sigma_0^2 + \frac{\sum_{i=1^n x_i{\sigma^2\right),
\left(\frac{1{\sigma_0^2 + \frac{n{\sigma^2\right)^{-1
میانگین از مشاهدات با دقت کل (مجموع تمام دقت‌های منفرد) و میانگین نمونه تحمین زده شده‌است
نرمالبا دقت معلوم 
(mean)
نرمال
میانگین از مشاهدات با دقت کل (مجموع تمام دقت‌های منفرد) و میانگین نمونه تخمین زده شده‌است
نرمالبا میانگین معلوم 
2 (واریانس)
گامای وارونه
واریانس از مشاهدات با واریانس نمونه (یعنی جمع مجذور انحرافات , که انحرافات از میانگین معلوم هستند) تخمین زده شده‌است
نرمالبا میانگین معلوم 
2 (واریانس)
کی دو وارونه مقیاس شده
واریانس از مشاهدات با واریانس نمونه تخمین زده شده‌است
نرمالبا میانگین معلوم 
(دقت)
گاما
دقت از مشاهدات با واریانس نمونه (یعنی جمع مجذور انحرافات , که انحرافات از میانگین معلوم هستند) تخمین زده شده‌است
نرمال
و  2 با این فرض که تعویض پذیری دارند
گامای نرمال وارونه
\beta + \tfrac{1{2 \sum_{i=1^n (x_i - \bar{x)^2 + \frac{n\nu{\nu+n\frac{(\bar{x-\mu_0)^2{2
میانگین نمونه است
میانگین از مشاهدات با میانگین نمونه تخمین زده شده‌است؛ واریانس از مشاهدات با میانگین نمونه و مجموع مجذور انحرافات تخمین زده شده‌است
نرمال
و با فرض اینکه تعویض پذیری دارند
نرمال-گاما
\beta + \tfrac{1{2 \sum_{i=1^n (x_i - \bar{x)^2 + \frac{n\nu{\nu+n\frac{(\bar{x-\mu_0)^2{2
میانگین نمونه است
میانگین از مشاهدات با میانگین نمونه , و دقت از مشاهدات با میانگین نمونه و مجموع مجذور انحرافات تخمین زده شده‌است
نرمال چندمتغیره با ماتریس کوواریانس معلوم 
(بردار میانگین)
نرمال چند متغیره
میانگین نمونه است
میانگین از مشاهداتی با دقت کل (مجموع تمام دقت‌های منفرد) و میانگین نمونه تخمین زده شده‌است
نرمال چندمتغیره با ماتریس دقت معلوم 
(بردار میانگین)
نرمال چندمتغیره
میانگین نمونه است
میانگین از مشاهداتی با دقت کل (مجموع تمام دقت‌های منفرد) و میانگین نمونه تخمین زده شده‌است
نرمال چندمتغیره با میانگین معلوم 
(ماتریس کوواریانس)
ویشارت وارون
ماتریس کوواریانس از مشاهدات با مجموع حاصل‌ضرب‌های انحراف دو به دو تخمین زده شده‌است
نرمال چندمتغیره با میانگین معلوم 
(ماتریس دقت)
ویشارت
ماتریس کوواریانس از مشاهدات با مجموع حاصل‌ضرب‌های انحراف دو به دو تخمین زده شده‌است
نرمال چندمتغیره
(بردار میانگین) و  (ماتریس کوواریانس)
وارون نرمال-ویشارت
میانگین نمونه است
میانگین از مشاهدات با میانگین نمونه تخمین زده شده‌است؛ ماتریس کوواریانس از مشاهدات با میانگین نمونه و با مجموع حاصل‌ضرب‌های انحراف دو به دو تخمین زده شده‌است
نرمال چندمتغیره
(بردار میانگین) and  (ماتریس دقت)
نرمال-ویشارت
میانگین نمونه است
میانگین از مشاهدات با میانگین نمونه تخمین زده شده‌است؛ ماتریس کوواریانس از مشاهدات با میانگین نمونه و با مجموع حاصل‌ضرب‌های انحراف دو به دو تخمین زده شده‌است
یکنواخت
پارتو
مشاهدات با مقدار بیشینه
پارتوبا کمینه معلوم x m
k (شکل)
گاما
مشاهدات با مجموع مرتبه بزرگی هر مشاهده (یعنی لگاریتم نسبت هر مشاهده به کمینه )
توزیع وایبولبا شکل معلوم 
(مقیاس)
توزیع گامای وارونه
مشاهدات با مجموع  امین توان هر مشاهده
توزیع لگاریتمی نرمالبا دقت معلوم 
(میانگین)
نرمال
«میانگین» از مشاهداتی با دقت کل (مجموع تمام دقت‌های منفرد) و با میانگین نمونه تخمین زده شده‌است
توزیع لگاریتمی نرمالبا میانگین معلوم 
(دقت)
گاما
دقت از مشاهدات با واریانس نمونه (یعنی مجموع مجذور انحرافات لگاریتمی — یعنی انحرافات لگاریتم نقاط داده از "میانگین")
نمایی
(نرخ)
گاما
مشاهدات با مجموع
(توزیع لوماکس)
گامابا شکل معلوم 
(نرخ)
گاما
مشاهدات با مجموع
توزیع گامای وارونهبا شکل معلوم 
(نرخ معکوس)
گاما
مشاهدات با مجموع
گامابا نرخ معلوم 
(شکل)
or مشاهدات ( برای تخمین , برای تخمین ) با حاصل‌ضرب
گاما
(شکل),  (نرخ معکوس)
از مشاهدات با حاصل‌ضرب تخمین زده شده‌است؛ از مشاهدات با مجموع تخمین زده شده‌است
توزیع بتا-دوجمله‌ای
نکات