توزيع ديريشلي المعمم

في الإحصاء ، يُعد توزيع ديريشلي المعمم ( GD ) تعميمًا لتوزيع ديريشلي ، يتميز ببنية تباين أكثر عمومية وعدد معلمات يقارب ضعف عدد معلمات ديريشلي. وتكون المتجهات العشوائية التي تخضع لتوزيع ديريشلي المعمم محايدة تمامًا . [ 1 ]

دالة الكثافة لـص1،...،صك-1{\displaystyle p_{1},\ldots ,p_{k-1}}يكون

[أنا=1ك-1ب(أأنا،بأنا)]-1صكبك-1-1أنا=1ك-1[صأناأأنا-1(ج=أناكصج)بأنا-1-(أأنا+بأنا)]{\displaystyle \left[\prod _{i=1}^{k-1}B(a_{i},b_{i})\right]^{-1}p_{k}^{b_{k-1}-1}\prod _{i=1}^{k-1}\left[p_{i}^{a_{i}-1}\left(\sum _{j=i}^{k}p_{j}\right)^{b_{i-1}-(a_{i}+b_{i})}\right]}

حيث نحددصك=1-أنا=1ك-1صأنا{\textstyle p_{k}=1-\sum _{i=1}^{k-1}p_{i}}. هناب(x،y){\displaystyle B(x,y)}يرمز إلى دالة بيتا . ويختزل هذا إلى توزيع ديريشليه القياسي إذابأنا-1=أأنا+بأنا{\displaystyle b_{i-1}=a_{i}+b_{i}}ل2أناك-1{\displaystyle 2\leqslant i\leqslant k-1}(ب0{\displaystyle b_{0}}(اعتباطي).

على سبيل المثال، إذا كانت k=4 ، فإن دالة الكثافة لـص1،ص2،ص3{\displaystyle p_{1},p_{2},p_{3}}يكون

[أنا=13ب(أأنا،بأنا)]-1ص1أ1-1ص2أ2-1ص3أ3-1ص4ب3-1(ص2+ص3+ص4)ب1-(أ2+ب2)(ص3+ص4)ب2-(أ3+ب3){\displaystyle \left[\prod _{i=1}^{3}B(a_{i},b_{i})\right]^{-1}p_{1}^{a_{1}-1}p_{2}^{a_{2}-1}p_{3}^{a_{3}-1}p_{4}^{b_{3}-1}\left(p_{2}+p_{3}+p_{4}\right)^{b_{1}-\left(a_{2}+b_{2}\right)}\left(p_{3}+p_{4}\right)^{b_{2}-\left(a_{3}+b_{3}\right)}}

أينص1+ص2+ص3<1{\displaystyle p_{1}+p_{2}+p_{3}<1}وص4=1-ص1-ص2-ص3{\displaystyle p_{4}=1-p_{1}-p_{2}-p_{3}}.

عرّف كونور وموسيمان دالة كثافة الاحتمال (PDF) على النحو الذي فعلاه للسبب التالي: تعريف المتغيرات العشوائيةz1،...،zك-1{\displaystyle z_{1},\ldots ,z_{k-1}}معz1=ص1،z2=ص2/(1-ص1)،z3=ص3/(1-(ص1+ص2))،...،zأنا=صأنا/(1-(ص1++صأنا-1)){\displaystyle z_{1}=p_{1},z_{2}=p_{2}/\left(1-p_{1}\right),z_{3}=p_{3}/\left(1-(p_{1}+p_{2})\right),\ldots ,z_{i}=p_{i}/\left(1-\left(p_{1}+\cdots +p_{i-1}\right)\right)}. ثمص1،...،صك{\displaystyle p_{1},\ldots ,p_{k}}إذا كان لدينا توزيع ديريشلي المعمم كما هو موضح أعلاه،zأنا{\displaystyle z_{i}}هي بيتا مستقلة ذات معلماتأأنا،بأنا{\displaystyle a_{i},b_{i}}،أنا=1،...،ك-1{\displaystyle i=1,\ldots ,k-1}.

صيغة بديلة قدمها وونغ

يقدم وونغ [ 2 ] الصيغة الأكثر إيجازًا لـx1++xك1{\displaystyle x_{1}+\cdots +x_{k}\leq 1}

أنا=1كxأناαأنا-1(1-x1--xأنا)γأناب(αأنا،βأنا){\displaystyle \prod _{i=1}^{k}{\frac {x_{i}^{\alpha _{i}-1}\left(1-x_{1}-\cdots -x_{i}\right)^{\gamma _{i}}}{B(\alpha _{i},\beta _{i})}}}

أينγج=βج-αج+1-βج+1{\displaystyle \gamma _{j}=\beta _{j}-\alpha _{j+1}-\beta _{j+1}}ل1جك-1{\displaystyle 1\leq j\leq k-1}وγك=βك-1{\displaystyle \gamma _{k}=\beta _{k}-1}لاحظ أن وونغ يُعرّف توزيعًا علىك{\displaystyle k}الفضاء ذو ​​الأبعاد (الذي يُعرّف ضمنيًاxك+1=1-أنا=1كxأنا{\textstyle x_{k+1}=1-\sum _{i=1}^{k}x_{i}}) بينما يستخدم كونور وموسيمانك-1{\displaystyle k-1}فضاء متعدد الأبعاد معxك=1-أنا=1ك-1xأنا{\textstyle x_{k}=1-\sum _{i=1}^{k-1}x_{i}}.

دالة العزم العامة

لوX=(X1،...،Xك)جيدك(α1،...،αك;β1،...،βك){\displaystyle X=\left(X_{1},\ldots ,X_{k}\right)\sim GD_{k}\left(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k}\right)}، ثم

هـ[X1ر1X2ر2Xكرك]=ج=1كΓ(αج+βج)Γ(αج+رج)Γ(βج+دلتاج)Γ(αج)Γ(βج)Γ(αج+βج+رج+دلتاج){\displaystyle E\left[X_{1}^{r_{1}}X_{2}^{r_{2}}\cdots X_{k}^{r_{k}}\right]=\prod _{j=1}^{k}{\frac {\Gamma \left(\alpha _{j}+\beta _{j}\right)\Gamma \left(\alpha _{j}+r_{j}\right)\Gamma \left(\beta _{j}+\delta _{j}\right)}{\Gamma \left(\alpha _{j}\right)\Gamma \left(\beta _{j}\right)\Gamma \left(\alpha _{j}+\beta _{j}+r_{j}+\delta _{j}\right)}}}

أيندلتاج=رج+1+رج+2++رك{\displaystyle \delta _{j}=r_{j+1}+r_{j+2}+\cdots +r_{k}}لج=1،2،،ك-1{\displaystyle j=1,2,\cdots ,k-1}ودلتاك=0{\displaystyle \delta _{k}=0}. هكذا

هـ(Xج)=αجαج+βجم=1ج-1βمαم+βم.{\displaystyle E\left(X_{j}\right)={\frac {\alpha _{j}}{\alpha _{j}+\beta _{j}}}\prod _{m=1}^{j-1}{\frac {\beta _{m}}{\alpha _{m}+\beta _{m}}}.}

الاختزال إلى توزيع ديريشليه القياسي

كما ذكر أعلاه، إذابأنا-1=أأنا+بأنا{\displaystyle b_{i-1}=a_{i}+b_{i}}ل2أناك{\displaystyle 2\leq i\leq k}عندئذٍ، يختزل التوزيع إلى توزيع ديريشليه القياسي. يختلف هذا الشرط عن الحالة المعتادة، حيث يؤدي ضبط المعاملات الإضافية للتوزيع المعمم على الصفر إلى الحصول على التوزيع الأصلي. مع ذلك، في حالة توزيع ديريشليه المعمم، ينتج عن ذلك دالة كثافة معقدة للغاية.

التحليل البايزي

يفترضX=(X1،...،Xك)جيدك(α1،...،αك;β1،...،βك){\displaystyle X=\left(X_{1},\ldots ,X_{k}\right)\sim GD_{k}\left(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k}\right)}هو ديريشليه معمّم، وY|X{\displaystyle Y\mid X}متعددة الحدود معن{\displaystyle n}التجارب (هنا)Y=(Y1،...،Yك){\displaystyle Y=\left(Y_{1},\ldots ,Y_{k}\right)}). كتابةYج=yج{\displaystyle Y_{j}=y_{j}}ل1جك{\displaystyle 1\leq j\leq k}وyك+1=ن-أنا=1كyأنا{\textstyle y_{k+1}=n-\sum _{i=1}^{k}y_{i}}المفصل الخلفي لـX|Y{\displaystyle X|Y}هو توزيع ديريشليه معمّم مع

X|Yجيدك(α1،...،αك;β1،...،βك){\displaystyle X\mid Y\sim GD_{k}\left({\alpha '}_{1},\ldots ,{\alpha '}_{k};{\beta '}_{1},\ldots ,{\beta '}_{k}\right)}

أينαج=αج+yج{\displaystyle {\alpha '}_{j}=\alpha _{j}+y_{j}}وβج=βج+أنا=ج+1ك+1yأنا{\displaystyle {\beta '}_{j}=\beta _{j}+\sum _{i=j+1}^{k+1}y_{i}}ل1جك.{\displaystyle 1\leqslant j\leqslant k.}

تجربة أخذ العينات

يقدم وونغ النظام التالي كمثال على كيفية اختلاف توزيعات ديريشليه وتوزيعات ديريشليه المعممة. يفترض أن جرة كبيرة تحتوي على كرات منك+1{\displaystyle k+1}ألوان مختلفة. نسبة كل لون غير معروفة. اكتبX=(X1،...،Xك){\displaystyle X=(X_{1},\ldots ,X_{k})}نسبة الكرات الملونةج{\displaystyle j}في الجرة.

التجربة 1. يعتقد المحلل 1 أنXد(α1،...،αك،αك+1){\displaystyle X\sim D(\alpha _{1},\ldots ,\alpha _{k},\alpha _{k+1})}(أي،X{\displaystyle X}دالة ديريشلي ذات المعاملاتαأنا{\displaystyle \alpha _{i}}ثم يقوم المحلل بـك+1{\displaystyle k+1}صناديق زجاجية ويضعαأنا{\displaystyle \alpha _{i}}كرات ملونةأنا{\displaystyle i}في الصندوقأنا{\displaystyle i}(يفترض أنαأنا{\displaystyle \alpha _{i}}هي أعداد صحيحة1{\displaystyle \geq 1}ثم يقوم المحلل 1 بسحب كرة من الجرة، ويلاحظ لونها (قل اللون).ج{\displaystyle j}) ويضعها في الصندوقج{\displaystyle j}يستطيع تحديد الصندوق الصحيح لأنه شفاف، وبالتالي تظهر ألوان الكرات الزجاجية بداخله. وتستمر العملية حتىن{\displaystyle n}تم سحب الكرات. وبالتالي فإن التوزيع الاحتمالي اللاحق هو توزيع ديريشليه، حيث تمثل المعلمات عدد الكرات في كل صندوق.

التجربة الثانية . يعتقد المحلل الثاني أنX{\displaystyle X}يتبع توزيع ديريشلي المعمم:Xجيد(α1،...،αك;β1،...،βك){\displaystyle X\sim GD(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k})}يُفترض مجدداً أن جميع المعاملات أعداد صحيحة موجبة. يقوم المحلل بـك+1{\displaystyle k+1}صناديق خشبية. تحتوي الصناديق على قسمين: قسم للكرات وقسم للكرات الزجاجية. الكرات ملونة، أما الكرات الزجاجية فليست ملونة. ثم لـج=1،...،ك{\displaystyle j=1,\ldots ,k}، يضعαج{\displaystyle \alpha _{j}}كرات ملونةج{\displaystyle j}، وβج{\displaystyle \beta _{j}}كرات زجاجية، في الصندوقج{\displaystyle j}ثم يضع كرة من اللونك+1{\displaystyle k+1}في الصندوقك+1{\displaystyle k+1}ثم يسحب المحلل كرة من الجرة. ولأن الصناديق مصنوعة من الخشب، لا يستطيع المحلل تحديد الصندوق الذي يضع فيه الكرة (كما فعل في التجربة الأولى أعلاه)؛ كما أن ذاكرته ضعيفة ولا يتذكر لون الكرات في كل صندوق. عليه أن يكتشف الصندوق الصحيح. يفعل ذلك بفتح الصندوق الأول ومقارنة الكرات الموجودة فيه بالكرة المسحوبة. إذا اختلفت الألوان، فالصندوق خاطئ. يضع المحلل كرة زجاجية في الصندوق الأول ثم ينتقل إلى الصندوق الثاني. يكرر العملية حتى تتطابق الكرات في الصندوق مع الكرة المسحوبة، وعندها يضع الكرة في الصندوق مع الكرات الأخرى من نفس اللون. ثم يسحب المحلل كرة أخرى من الجرة ويكرر العملية حتىن{\displaystyle n}يتم سحب الكرات. ثم يتم تعميم التوزيع الاحتمالي اللاحق باستخدام دالة ديريشليه ذات المعاملات.α{\displaystyle \alpha }عدد الكرات، وβ{\displaystyle \beta }عدد الكرات الزجاجية في كل صندوق.

لاحظ أنه في التجربة الثانية، كان لتغيير ترتيب الصناديق تأثير غير تافه، على عكس التجربة الأولى.

انظر أيضاً

مراجع

  1. آر جيه كونور وجيه إي موسيمان 1969. مفاهيم الاستقلال للنسب مع تعميم لتوزيع ديريشليه . مجلة الجمعية الإحصائية الأمريكية، المجلد 64، الصفحات 194-206
  2. تي. تي. وونغ 1998. توزيع ديريشلي المعمم في التحليل البايزي . الرياضيات التطبيقية والحساب، المجلد 97، الصفحات 165-181