توزيع ديريشليه متعدد الحدود السالب

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

إذا كانت معلمات توزيع ديريشليه هيα{\displaystyle {\boldsymbol {\alpha }}}وإذا

X|صنيو مكسيكو(x0،ص)،{\displaystyle X\mid p\sim \operatorname {NM} (x_{0},\mathbf {p} ),}

أين

صمدير(α0،α)،{\displaystyle \mathbf {p} \sim \operatorname {Dir} (\alpha _{0},{\boldsymbol {\alpha }}),}

إذن، التوزيع الهامشي لـ X هو توزيع متعدد الحدود سالب من نوع ديريشليه:

XDNM(x0،α0،α).{\displaystyle X\sim \operatorname {DNM} (x_{0},\alpha _{0},{\boldsymbol {\alpha }}).}

في ما سبق،نيو مكسيكو(x0،ص){\displaystyle \operatorname {NM} (x_{0},\mathbf {p} )}هو التوزيع متعدد الحدود السالب ومدير(α0،α){\displaystyle \operatorname {Dir} (\alpha _{0},{\boldsymbol {\alpha }})}هو توزيع ديريشليه .

تحفيز

توزيع ديريشليه متعدد الحدود السالب كتوزيع مركب

يُعد توزيع ديريشليه توزيعًا مترافقًا مع التوزيع متعدد الحدود السالب. وتؤدي هذه الحقيقة إلى توزيع مركب قابل للمعالجة التحليلية . بالنسبة لمتجه عشوائي لعدد الفئاتx=(x1،...،xم){\displaystyle \mathbf {x} =(x_{1},\dots ,x_{m})}، موزعة وفقًا لتوزيع متعدد الحدود السالب ، يتم الحصول على التوزيع المركب عن طريق التكامل على التوزيع لـ p والذي يمكن اعتباره متجهًا عشوائيًا يتبع توزيع ديريشليه:

برو(x|x0،α0،α)=صشمالهـزمuلت(x|x0،ص)دأنار(ص|α0،α)دص{\displaystyle \Pr(\mathbf {x} \mid x_{0},\alpha _{0},{\boldsymbol {\alpha }})=\int _{\mathbf {p} }\mathrm {NegMult} (\mathbf {x} \mid x_{0},\mathbf {p} )\mathrm {Dir} (\mathbf {p} \mid \alpha _{0},{\boldsymbol {\alpha }}){\textrm {d}}\mathbf {p} }
برو(x|x0،α0،α)=Γ(أنا=0مxأنا)Γ(x0)أنا=1مxأنا!1ب(α+)صأنا=0مصأناxأنا+αأنا-1دص{\displaystyle \Pr(\mathbf {x} \mid x_{0},\alpha _{0},{\boldsymbol {\alpha }})={\frac {\Gamma \left(\sum _{i=0}^{m}{x_{i}}\right)}{\Gamma (x_{0})\prod _{i=1}^{m}x_{i}!}}{\frac {1}{\mathrm {B} ({\boldsymbol {\alpha }}_{+})}}\int _{\mathbf {p} }\prod _{i=0}^{m}p_{i}^{x_{i}+\alpha _{i}-1}{\textrm {d}}\mathbf {p} }

مما ينتج عنه الصيغة التالية:

برو(x|x0،α0،α)=Γ(أنا=0مxأنا)Γ(x0)أنا=1مxأنا!ب(x++α+)ب(α+){\displaystyle \Pr(\mathbf {x} \mid x_{0},\alpha _{0},{\boldsymbol {\alpha }})={\frac {\Gamma \left(\sum _{i=0}^{m}{x_{i}}\right)}{\Gamma (x_{0})\prod _{i=1}^{m}x_{i}!}}{\frac {{\mathrm {B} }(\mathbf {x_{+}} +{\boldsymbol {\alpha }}_{+})}{\mathrm {B} ({\boldsymbol {\alpha }}_{+})}}}

أينx+{\displaystyle \mathbf {x_{+}} }وα+{\displaystyle {\boldsymbol {\alpha }}_{+}}همم+1{\displaystyle m+1}المتجهات ذات الأبعاد التي تم إنشاؤها عن طريق إلحاق القيم العدديةx0{\displaystyle x_{0}}وα0{\displaystyle \alpha _{0}}إلىم{\displaystyle m}المتجهات متعددة الأبعادx{\displaystyle \mathbf {x} }وα{\displaystyle {\boldsymbol {\alpha }}} على التوالي وب{\displaystyle \mathrm {B} }هي الصيغة متعددة المتغيرات لدالة بيتا . يمكننا كتابة هذه المعادلة صراحةً على النحو التالي:

برو(x|x0،α0،α)=x0Γ(أنا=0مxأنا)Γ(أنا=0مαأنا)Γ(أنا=0م(xأنا+αأنا))أنا=0مΓ(xأنا+αأنا)Γ(xأنا+1)Γ(αأنا).{\displaystyle \Pr(\mathbf {x} \mid x_{0},\alpha _{0},{\boldsymbol {\alpha }})=x_{0}{\frac {\Gamma (\sum _{i=0}^{m}x_{i})\Gamma (\sum _{i=0}^{m}\alpha _{i})}{\Gamma (\sum _{i=0}^{m}(x_{i}+\alpha _{i}))}}\prod _{i=0}^{m}{\frac {\Gamma (x_{i}+\alpha _{i})}{\Gamma (x_{i}+1)\Gamma (\alpha _{i})}}.}

توجد صيغ بديلة. أحد التمثيلات الملائمة [ 1 ] هو

برو(x|x0،α0،α)=Γ(x)Γ(x0)أنا=1مΓ(xأنا+1)×Γ(α)أنا=0مΓ(αأنا)×أنا=0مΓ(xأنا+αأنا)Γ(x+α){\displaystyle \Pr(\mathbf {x} \mid x_{0},\alpha _{0},{\boldsymbol {\alpha }})={\frac {\Gamma (x_{\bullet })}{\Gamma (x_{0})\prod _{i=1}^{m}\Gamma (x_{i}+1)}}\times {\frac {\Gamma (\alpha _{\bullet })}{\prod _{i=0}^{m}\Gamma (\alpha _{i})}}\times {\frac {\prod _{i=0}^{m}\Gamma (x_{i}+\alpha _{i})}{\Gamma (x_{\bullet }+\alpha _{\bullet })}}}

أينx=x0+x1++xم{\displaystyle x_{\bullet }=x_{0}+x_{1}+\cdots +x_{m}}وα=α0+α1++αم{\displaystyle \alpha _{\bullet }=\alpha _{0}+\alpha _{1}+\cdots +\alpha _{m}}.

ويمكن كتابة هذا أيضًا

برو(x|x0،α0،α)=ب(x،α)ب(x0،α0)أنا=1مΓ(xأنا+αأنا)xأنا!Γ(αأنا).{\displaystyle \Pr(\mathbf {x} \mid x_{0},\alpha _{0},{\boldsymbol {\alpha }})={\frac {\mathrm {B} (x_{\bullet },\alpha _{\bullet })}{\mathrm {B} (x_{0},\alpha _{0})}}\prod _{i=1}^{m}{\frac {\Gamma (x_{i}+\alpha _{i})}{x_{i}!\Gamma (\alpha _{i})}}.}

ملكيات

التوزيعات الهامشية

للحصول على التوزيع الهامشي لمجموعة فرعية من المتغيرات العشوائية متعددة الحدود السالبة من نوع ديريشليه، يكفي حذف المتغيرات غير ذات الصلة.αأنا{\displaystyle \alpha _{i}}المتغيرات (التي يرغب المرء في استبعادها) منα{\displaystyle {\boldsymbol {\alpha }}}المتجه. التوزيع المشترك للمتغيرات العشوائية المتبقية هودشمالم(x0،α0،α(-)){\displaystyle \mathrm {DNM} (x_{0},\alpha _{0},{\boldsymbol {\alpha _{(-)}}})}أينα(-){\displaystyle {\boldsymbol {\alpha _{(-)}}}}هو المتجه بعد إزالةαأنا{\displaystyle \alpha _{i}}يقال إن التوزيعات الهامشية أحادية المتغير تتبع التوزيع الثنائي السالب بيتا .

التوزيعات الشرطية

إذا تم تقسيم x ذي الأبعاد m على النحو التالي

x=[x(1)x(2)] مع المقاسات [q×1(م-q)×1]{\displaystyle \mathbf {x} ={\begin{bmatrix}\mathbf {x} ^{(1)}\\\mathbf {x} ^{(2)}\end{bmatrix}}{\text{ with sizes }}{\begin{bmatrix}q\times 1\\(m-q)\times 1\end{bmatrix}}}

وبناءً على ذلك α{\displaystyle {\boldsymbol {\alpha }}}

α=[α(1)α(2)] مع المقاسات [q×1(م-q)×1]{\displaystyle {\boldsymbol {\alpha }}={\begin{bmatrix}{\boldsymbol {\alpha }}^{(1)}\\{\boldsymbol {\alpha }}^{(2)}\end{bmatrix}}{\text{ with sizes }}{\begin{bmatrix}q\times 1\\(m-q)\times 1\end{bmatrix}}}

ثم التوزيع الشرطي لـX(1){\displaystyle \mathbf {X} ^{(1)}}علىX(2)=x(2){\displaystyle \mathbf {X} ^{(2)}=\mathbf {x} ^{(2)}}يكوندشمالم(x0،α0،α(1)){\displaystyle \mathrm {DNM} (x_{0}^{\prime },\alpha _{0}^{\prime },{\boldsymbol {\alpha }}^{(1)})}أين

x0=x0+أنا=1م-qxأنا(2){\displaystyle x_{0}^{\prime }=x_{0}+\sum _{i=1}^{m-q}x_{i}^{(2)}}

و

α0=α0+أنا=1م-qαأنا(2){\displaystyle \alpha _{0}^{\prime }=\alpha _{0}+\sum _{i=1}^{m-q}\alpha _{i}^{(2)}}.

إنه،

برو(x(1)|x(2)،x0،α0،α)=ب(x،α)ب(x0،α0)أنا=1qΓ(xأنا(1)+αأنا(1))(xأنا(1)!)Γ(αأنا(1)){\displaystyle \Pr(\mathbf {x} ^{(1)}\mid \mathbf {x} ^{(2)},x_{0},\alpha _{0},{\boldsymbol {\alpha }})={\frac {\mathrm {B} (x_{\bullet },\alpha _{\bullet })}{\mathrm {B} (x_{0}^{\prime },\alpha _{0}^{\prime })}}\prod _{i=1}^{q}{\frac {\Gamma (x_{i}^{(1)}+\alpha _{i}^{(1)})}{(x_{i}^{(1)}!)\Gamma (\alpha _{i}^{(1)})}}}

بشرط المجموع

التوزيع الشرطي لتوزيع ديريشليه متعدد الحدود السالب علىأنا=1مxأنا=ن{\displaystyle \sum _{i=1}^{m}x_{i}=n}هو توزيع ديريشليه متعدد الحدود ذو معلماتن{\displaystyle n} وα{\displaystyle {\boldsymbol {\alpha }}}. إنه

برو(x|أنا=1مxأنا=ن،x0،α0،α)=ن!Γ(أنا=1مαأنا)Γ(ن+أنا=1مαأنا)أنا=1مΓ(xأنا+αأنا)xأنا!Γ(αأنا){\displaystyle \Pr(\mathbf {x} \mid \sum _{i=1}^{m}x_{i}=n,x_{0},\alpha _{0},{\boldsymbol {\alpha }})={\frac {n!\Gamma \left(\sum _{i=1}^{m}\alpha _{i}\right)}{\Gamma \left(n+\sum _{i=1}^{m}\alpha _{i}\right)}}\prod _{i=1}^{m}{\frac {\Gamma (x_{i}+\alpha _{i})}{x_{i}!\Gamma (\alpha _{i})}}}.

لاحظ أن التعبير لا يعتمد علىx0{\displaystyle x_{0}}أوα0{\displaystyle \alpha _{0}}.

تجميع

لو

X=(X1،...،Xم)DNM(x0،α0،α1،...،αم){\displaystyle X=(X_{1},\ldots ,X_{m})\sim \operatorname {DNM} (x_{0},\alpha _{0},\alpha _{1},\ldots ,\alpha _{m})}

ثم، إذا تم حذف المتغيرات العشوائية ذات المؤشرات الموجبة i و j من المتجه واستبدالها بمجموعها،

X=(X1،...،Xأنا+Xج،...،Xم)DNM(x0،α0،α1،...،αأنا+αج،...،αم).{\displaystyle X'=(X_{1},\ldots ,X_{i}+X_{j},\ldots ,X_{m})\sim \operatorname {DNM} \left(x_{0},\alpha _{0},\alpha _{1},\ldots ,\alpha _{i}+\alpha _{j},\ldots ,\alpha _{m}\right).}

مصفوفة الارتباط

لα0>2{\displaystyle \alpha _{0}>2}عناصر مصفوفة الارتباط هي

ρ(Xأنا،Xأنا)=1.{\displaystyle \rho (X_{i},X_{i})=1.}
ρ(Xأنا،Xج)=كوف(Xأنا،Xج)متغير(Xأنا)متغير(Xج)=αأناαج(α0+αأنا-1)(α0+αج-1).{\displaystyle \rho (X_{i},X_{j})={\frac {\operatorname {cov} (X_{i},X_{j})}{\sqrt {\operatorname {var} (X_{i})\operatorname {var} (X_{j})}}}={\sqrt {\frac {\alpha _{i}\alpha _{j}}{(\alpha _{0}+\alpha _{i}-1)(\alpha _{0}+\alpha _{j}-1)}}}.}

ذي الذيل الثقيل

يُعد توزيع ديريشليه متعدد الحدود السالب توزيعًا ذو ذيل سميك . وليس له متوسط ​​محدود لـα01{\displaystyle \alpha _{0}\leq 1}ولها مصفوفة تباين لا نهائية لـα02{\displaystyle \alpha _{0}\leq 2}لذلك، فإن دالة توليد العزوم غير موجودة.

التطبيقات

Dirichlet متعدد الحدود السلبي كنموذج Pólya urn

في حالة عندمام+2{\displaystyle m+2}حدودx0،α0{\displaystyle x_{0},\alpha _{0}}وα{\displaystyle {\boldsymbol {\alpha }}}إذا كانت أعدادًا صحيحة موجبة، فيمكن أيضًا تبرير متعددة الحدود السالبة لـ Dirichlet بنموذج الجرة - أو بشكل أكثر تحديدًا نموذج جرة Pólya الأساسي . لنفترض جرة تحتوي في البداية علىأنا=0مαأنا{\displaystyle \sum _{i=0}^{m}{\alpha _{i}}}كرات منم+1{\displaystyle m+1}ألوان متنوعة تشملα0{\displaystyle \alpha _{0}}الكرات الحمراء (لون التوقف). المتجهα{\displaystyle {\boldsymbol {\alpha }}}يُعطي العدد النسبي للكرات الأخرى من مختلف الأنواعم{\displaystyle m}الألوان غير الحمراء. في كل خطوة من خطوات النموذج، تُسحب كرة عشوائيًا من الجرة وتُعاد إلى مكانها، مع كرة إضافية من نفس اللون. تُكرر العملية مرارًا وتكرارًا حتىx0{\displaystyle x_{0}}يتم رسم كرات حمراء اللون. المتجه العشوائيX{\displaystyle \mathbf {X} }من السحوبات الملحوظة للآخرم{\displaystyle m}يتم توزيع الألوان غير الحمراء وفقًا لـدشمالم(x0،α0،α){\displaystyle \mathrm {DNM} (x_{0},\alpha _{0},{\boldsymbol {\alpha }})}لاحظ أنه في نهاية التجربة، تحتوي الجرة دائمًا على العدد الثابتx0+α0{\displaystyle x_{0}+\alpha _{0}}من الكرات الحمراء التي تحتوي على الرقم العشوائيX+α{\displaystyle \mathbf {X} +{\boldsymbol {\alpha }}}من الآخرم{\displaystyle m}الألوان.

انظر أيضاً

مراجع

  1. فارويل، دانيال وفارويل، فيرنون. (2012). انحدار ديريشلي متعدد الحدود السالب لبيانات العد المترابطة ذات التشتت الزائد. الإحصاء الحيوي (أكسفورد، إنجلترا). 14. 10.1093/biostatistics/kxs050.