توزيع بيتا المعمم

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

تعريف

يُعرَّف المتغير العشوائي بيتا المعمم، Y ، بواسطة دالة كثافة الاحتمال التالية (pdf):

جيب(y؛أ،ب،ج،ص،q)=|أ|yأص-1(1-(1-ج)(y/ب)أ)q-1بأصب(ص،q)(1+ج(y/ب)أ)ص+q ل 0<yأ<بأ1-ج،{\displaystyle GB(y;a,b,c,p,q)={\frac {|a|y^{ap-1}(1-(1-c)(y/b)^{a})^{q-1}}{b^{ap}B(p,q)(1+c(y/b)^{a})^{p+q}}}\quad \quad {\text{ for }}0<y^{a}<{\frac {b^{a}}{1-c}},}

وصفر فيما عدا ذلك. هنا تحقق المعاملاتأ0{\displaystyle a\neq 0}،0ج1{\displaystyle 0\leq c\leq 1}وب{\displaystyle b}،ص{\displaystyle p}، وq{\displaystyle q}موجب. الدالة B ( p,q ) هي دالة بيتا . المعاملب{\displaystyle b}هو مُعامل المقياس ، وبالتالي يمكن ضبطه على1{\displaystyle 1}دون الإخلال بعمومية الحل ، ولكن عادةً ما يتم توضيح ذلك صراحةً كما في الدالة أعلاه. أما مُعامل الموقع (غير المُدرج في الصيغة أعلاه) فيُترك عادةً ضمنيًا ويُضبط على0{\displaystyle 0}.

شجرة توزيع GB

ملكيات

لحظات

يمكن إثبات أن العزم من الرتبة h يمكن التعبير عنه على النحو التالي:

هـجيب(Yح)=بحب(ص+ح/أ،q)ب(ص،q)2F1[ص+ح/أ،ح/أ؛جص+q+ح/أ؛]،{\displaystyle \operatorname {E} _{GB}(Y^{h})={\frac {b^{h}B(p+h/a,q)}{B(p,q)}}{}_{2}F_{1}{\begin{bmatrix}p+h/a,h/a;c\\p+q+h/a;\end{bmatrix}},}

أين2F1{\displaystyle {}_{2}F_{1}}يشير إلى سلسلة فرط هندسية (التي تتقارب لجميع قيم h إذا كانت c < 1، أو لجميع قيم h / a < q إذا كانت c = 1).

يشمل التوزيع بيتا المعمم العديد من التوزيعات كحالات حدية أو خاصة. هذه التوزيعات موضحة في شجرة توزيع بيتا المعمم أعلاه. فيما يلي فروعه الثلاثة المباشرة، أو عائلاته الفرعية.

بيتا المعممة من النوع الأول (GB1)

يتم تعريف بيتا المعممة من النوع الأول بواسطة دالة الكثافة الاحتمالية التالية:

جيب1(y؛أ،ب،ص،q)=|أ|yأص-1(1-(y/ب)أ)q-1بأصب(ص،q){\displaystyle GB1(y;a,b,p,q)={\frac {|a|y^{ap-1}(1-(y/b)^{a})^{q-1}}{b^{ap}B(p,q)}}}

ل0<yأ<بأ{\displaystyle 0<y^{a}<b^{a}}أينب{\displaystyle b}،ص{\displaystyle p}، وq{\displaystyle q}نتائج إيجابية. من السهل التحقق من ذلك.

جيب1(y؛أ،ب،ص،q)=جيب(y؛أ،ب،ج=0،ص،q).{\displaystyle GB1(y;a,b,p,q)=GB(y;a,b,c=0,p,q).}

تُعطى عزوم GB1 بواسطة

هـجيب1(Yح)=بحب(ص+ح/أ،q)ب(ص،q).{\displaystyle \operatorname {E} _{GB1}(Y^{h})={\frac {b^{h}B(p+h/a,q)}{B(p,q)}}.}

يشمل GB1 بيتا من النوع الأول (B1)، وجاما المعممة (GG)، وباريتو (PA) كحالات خاصة:

ب1(y؛ب،ص،q)=جيب1(y؛أ=1،ب،ص،q)،{\displaystyle B1(y;b,p,q)=GB1(y;a=1,b,p,q),}
جيجي(y؛أ،β،ص)=ليمqجيب1(y؛أ،ب=q1/أβ،ص،q)،{\displaystyle GG(y;a,\beta ,p)=\lim _{q\to \infty }GB1(y;a,b=q^{1/a}\beta ,p,q),}
Pأ(y؛ب،ص)=جيب1(y؛أ=-1،ب،ص،q=1).{\displaystyle PA(y;b,p)=GB1(y;a=-1,b,p,q=1).}

بيتا المعممة من النوع الثاني (GB2)

يتم تعريف GB2 بواسطة ملف PDF التالي:

جيب2(y؛أ،ب،ص،q)=|أ|yأص-1بأصب(ص،q)(1+(y/ب)أ)ص+q{\displaystyle GB2(y;a,b,p,q)={\frac {|a|y^{ap-1}}{b^{ap}B(p,q)(1+(y/b)^{a})^{p+q}}}}

ل0<y<{\displaystyle 0<y<\infty }وصفر فيما عدا ذلك. يمكن للمرء أن يتحقق من ذلك.

جيب2(y؛أ،ب،ص،q)=جيب(y؛أ،ب،ج=1،ص،q).{\displaystyle GB2(y;a,b,p,q)=GB(y;a,b,c=1,p,q).}

تُعطى عزوم GB2 بواسطة

هـجيب2(Yح)=بحب(ص+ح/أ،q-ح/أ)ب(ص،q).{\displaystyle \operatorname {E} _{GB2}(Y^{h})={\frac {b^{h}B(p+h/a,qh/a)}{B(p,q)}}.}

يُعرف GB2 أيضًا باسم بيتا برايم المعمم (باتيل، بوسويل، راتناباريخي (1984))، [ 2 ] بيتا المحول (فينتر، 1983)، [ 3 ] F المعمم (كالفلايش وبرينتيس، 1980)، [ 4 ] وهو حالة خاصة (μ≡0) من توزيع فيلر-باريتو (أرنولد، 1983) [ 5 ] . يضم GB2 توزيعات شائعة مثل جاما المعمم (GG)، وبور من النوع 3، وبور من النوع 12 ، وداغوم ، واللوغاريتمي الطبيعي ، وويبول ، وجاما ، ولوماكس ، وإحصائية F ، وفيسك أو رايلي ، ومربع كاي ، ونصف الطبيعي ، ونصف t للطالب ، والأسي ، ولوغاريتم لابلاس غير المتماثل، ولوغاريتم لابلاس ، ودالة القوة ، واللوغاريتم اللوجستي . [ 6 ]

بيتا

يتم تعريف عائلة التوزيعات بيتا (B) على النحو التالي: [ 1 ]

ب(y؛ب،ج،ص،q)=yص-1(1-(1-ج)(y/ب))q-1بصب(ص،q)(1+ج(y/ب))ص+q{\displaystyle B(y;b,c,p,q)={\frac {y^{p-1}(1-(1-c)(y/b))^{q-1}}{b^{p}B(p,q)(1+c(y/b))^{p+q}}}}

ل0<y<ب/(1-ج){\displaystyle 0<y<b/(1-c)}وصفر فيما عدا ذلك. وتتضح علاقتها بالبريطانية العظمى أدناه:

ب(y؛ب،ج،ص،q)=جيب(y؛أ=1،ب،ج،ص،q).{\displaystyle B(y;b,c,p,q)=GB(y;a=1,b,c,p,q).}

تشمل عائلة بيتا بيتا من النوع الأول والثاني [ 7 ] (B1 وB2، حيث يُشار إلى B2 أيضًا باسم بيتا الأولية )، والتي تتوافق مع c = 0 و c = 1 على التوالي.ج=0{\displaystyle c=0}،ب=1{\displaystyle b=1}ينتج عنه توزيع بيتا القياسي ذو المعلمتين .

جاما المعممة

يُعد توزيع جاما المعمم ( GG) حالة حدية لتوزيع جاما المعمم من النوع الثاني (GB2). ويتم تعريف دالة كثافة الاحتمال الخاصة به على النحو التالي: [ 8 ]

جيجي(y؛أ،β،ص)=ليمqجيب2(y،أ،ب=q1/أβ،ص،q)=|أ|yأص-1هـ-(y/β)أβأصΓ(ص){\displaystyle GG(y;a,\beta ,p)=\lim _{q\rightarrow \infty }GB2(y,a,b=q^{1/a}\beta ,p,q)={\frac {|a|y^{ap-1}e^{-(y/\beta )^{a}}}{\beta ^{ap}\Gamma (p)}}}

معح{\displaystyle h}اللحظات المعطاة بواسطة

هـ(Yجيجيح)=βحΓ(ص+ح/أ)Γ(ص).{\displaystyle \operatorname {E} (Y_{GG}^{h})={\frac {\beta ^{h}\Gamma (p+h/a)}{\Gamma (p)}}.}

كما ذكرنا سابقًا، فإن شجرة عائلة التوزيع GB تصور بشكل مرئي الحالات الخاصة والحدودية (انظر McDonald و Xu (1995)).

باريتو

يُعد توزيع باريتو ( PA) الحالة الحدية التالية لتوزيع غاما المعمم:

Pأ(y؛β،θ)=ليمأ-جيجي(y؛أ،β،ص=-θ/أ)=ليمأ-(θy-θ-1هـ-(y/β)أβ-θ(-θ/أ)Γ(-θ/أ))={\displaystyle PA(y;\beta ,\theta )=\lim _{a\rightarrow -\infty }GG(y;a,\beta ,p=-\theta /a)=\lim _{a\rightarrow -\infty }\left({\frac {\theta y^{-\theta -1}e^{-(y/\beta )^{a}}}{\beta ^{-\theta }(-\theta /a)\Gamma (-\theta /a)}}\right)=}
ليمأ-(θy-θ-1هـ-(y/β)أβ-θΓ(1-θ/أ))=θy-θ-1β-θ{\displaystyle \lim _{a\rightarrow -\infty }\left({\frac {\theta y^{-\theta -1}e^{-(y/\beta )^{a}}}{\beta ^{-\theta }\Gamma (1-\theta /a)}}\right)={\frac {\theta y^{-\theta -1}}{\beta ^{-\theta }}}}لβ<y{\displaystyle \beta <y}و0{\displaystyle 0}خلاف ذلك.

وظيفة الطاقة

يمثل توزيع دالة القوة (P) الحالة الحدية التالية لتوزيع جاما المعمم:

P(y؛β،θ)=ليمأجيجي(y؛أ=θ/ص،β،ص)=ليمأ|θص|yθ-1هـ-(y/β)أβθΓ(ص)=ليمأθyθ-1صΓ(ص)βθهـ-(y/β)أ={\displaystyle P(y;\beta ,\theta )=\lim _{a\rightarrow \infty }GG(y;a=\theta /p,\beta ,p)=\lim _{a\rightarrow \infty }{\frac {\mid {\frac {\theta }{p}}|y^{\theta -1}e^{-(y/\beta )^{a}}}{\beta ^{\theta }\Gamma (p)}}=\lim _{a\rightarrow \infty }{\frac {\theta y^{\theta -1}}{p\Gamma (p)\beta ^{\theta }}}e^{-(y/\beta )^{a}}=}
ليمأθyθ-1Γ(ص+1)βθهـ-(y/β)أ=ليمأθyθ-1Γ(θأ+1)βθهـ-(y/β)أ=θyθ-1βθ{\displaystyle \lim _{a\rightarrow \infty }{\frac {\theta y^{\theta -1}}{\Gamma (p+1)\beta ^{\theta }}}e^{-(y/\beta )^{a}}=\lim _{a\rightarrow \infty }{\frac {\theta y^{\theta -1}}{\Gamma ({\frac {\theta }{a}}+1)\beta ^{\theta }}}e^{-(y/\beta )^{a}}={\frac {\theta y^{\theta -1}}{\beta ^{\theta }}}}ل0<y<β{\displaystyle 0<y<\beta }وθ>0{\displaystyle \theta >0}.

لوغاريتم لابلاس غير المتماثل

يتم تعريف توزيع لوغاريتم لابلاس غير المتماثل (يشار إليه أيضًا باسم توزيع باريتو المزدوج [ 9 ] ) على النحو التالي: [ 10 ]

ألل(y؛ب،λ1،λ2)=ليمأجيب2(y؛أ،ب،ص=λ1/أ،q=λ2/أ)=λ1λ2y(λ1+λ2){(yب)λ1ل 0<y<ب(بy)λ2ل yب{\displaystyle ALL(y;b,\lambda _{1},\lambda _{2})=\lim _{a\rightarrow \infty }GB2(y;a,b,p=\lambda _{1}/a,q=\lambda _{2}/a)={\frac {\lambda _{1}\lambda _{2}}{y(\lambda _{1}+\lambda _{2})}}{\begin{cases}({\frac {y}{b}})^{\lambda _{1}}&{\mbox{for }}0<y<b\\({\frac {b}{y}})^{\lambda _{2}}&{\mbox{for }}y\geq b\end{cases}}}

حيثح{\displaystyle h}تُعطى اللحظات من خلال

هـ(Yأللح)=بحλ1λ2(λ1+ح)(λ2-ح).{\displaystyle \operatorname {E} (Y_{ALL}^{h})={\frac {b^{h}\lambda _{1}\lambda _{2}}{(\lambda _{1}+h)(\lambda _{2}-h)}}.}

متىλ1=λ2{\displaystyle \lambda _{1}=\lambda _{2}}وهذا يعادل توزيع لوغاريتم لابلاس .

توزيع بيتا المعمم الأسي

تأجيرYجيب(y؛أ،ب،ج،ص،q){\displaystyle Y\sim GB(y;a,b,c,p,q)}(بدون معلمة الموقع)، المتغير العشوائيZ=ln(Y){\displaystyle Z=\ln(Y)}، مع إعادة تحديد المعلماتدلتا=ln(ب){\displaystyle \delta =\ln(b)}وσ=1/أ{\displaystyle \sigma =1/a}، يتم توزيعها كتوزيع بيتا معمّم أُسّي (EGB)، مع دالة كثافة الاحتمال التالية:

هـجيب(z؛دلتا،σ،ج،ص،q)=هـص(z-دلتا)/σ(1-(1-ج)هـ(z-دلتا)/σ)q-1|σ|ب(ص،q)(1+جهـ(z-دلتا)/σ)ص+q{\displaystyle EGB(z;\delta ,\sigma ,c,p,q)={\frac {e^{p(z-\delta )/\sigma }(1-(1-c)e^{(z-\delta )/\sigma })^{q-1}}{|\sigma |B(p,q)(1+ce^{(z-\delta )/\sigma })^{p+q}}}}

ل-<z-دلتاσ<ln(11-ج){\displaystyle -\infty <{\frac {z-\delta }{\sigma }}<\ln({\frac {1}{1-c}})}وصفر فيما عدا ذلك. يشمل التوزيع الاحتمالي الموسع تعميمات لتوزيعات جومبيرتز ، وجامبل ، والقيمة القصوى من النوع الأول ، واللوجستية ، وبور-2، والأسية ، والطبيعية . المعلمةدلتا=ln(ب){\displaystyle \delta =\ln(b)}يمثل معلمة الموقع لـ EGB (بينماب{\displaystyle b}(وهو معامل المقياس لـ GB)، وσ=1/أ{\displaystyle \sigma =1/a}يمثل معامل المقياس لـ EGB (بينماأ{\displaystyle a}(وهو أحد معلمات شكل GB)؛ وبالتالي فإن EGB له ثلاثة معلمات شكل .

يتضمن الشكل المرفق العلاقة بين قاعدة البيانات الأوروبية وحالاتها الخاصة والحدودية. [ 11 ]

عائلة توزيعات EGB

دالة توليد العزوم

باستخدام رموز مشابهة لما سبق، يمكن التعبير عن دالة توليد العزوم لـ EGB على النحو التالي:

مهـجيب(Z)=هـدلتاتب(ص+تσ،q)ب(ص،q)2F1[ص+تσ،تσ؛جص+q+تσ؛].{\displaystyle M_{EGB}(Z)={\frac {e^{\delta t}B(p+t\sigma ,q)}{B(p,q)}}{}_{2}F_{1}{\begin{bmatrix}p+t\sigma ,t\sigma ;c\\p+q+t\sigma ;\end{bmatrix}}.}

توزيع بيتا المعمم متعدد المتغيرات

تُوسّع دالة كثافة الاحتمال بيتا المعممة متعددة المتغيرات التوزيعات أحادية المتغير المذكورة أعلاه.ن{\displaystyle n}المتغيراتy=(y1،...،yن){\displaystyle y=(y_{1},...,y_{n})}، يُعرِّف1xن{\displaystyle 1xn}متجهات المعلمات بواسطةأ=(أ1،...،أن){\displaystyle a=(a_{1},...,a_{n})}،ب=(ب1،...،بن){\displaystyle b=(b_{1},...,b_{n})}،ج=(ج1،...،جن){\displaystyle c=(c_{1},...,c_{n})}، وص=(ص1،...،صن){\displaystyle p=(p_{1},...,p_{n})}حيث كلبأنا{\displaystyle b_{i}}وصأنا{\displaystyle p_{i}}إيجابي، و0{\displaystyle 0}{\displaystyle \leq }جأنا{\displaystyle c_{i}}{\displaystyle \leq }1{\displaystyle 1}المعلمةq{\displaystyle q}يُفترض أن تكون موجبة، ونُعرّف الدالةب(ص1،...،صن،q){\displaystyle B(p_{1},...,p_{n},q)}=Γ(ص1)...Γ(صن)Γ(q)Γ(ص¯+q){\displaystyle {\frac {\Gamma (p_{1})...\Gamma (p_{n})\Gamma (q)}{\Gamma ({\bar {p}}+q)}}}لص¯{\displaystyle {\bar {p}}}=أنا=1نصأنا{\displaystyle \sum _{i=1}^{n}p_{i}}.

دالة كثافة الاحتمال لبيتا المعممة متعددة المتغيرات (مجيب{\displaystyle MGB}يمكن كتابة ) على النحو التالي:

مجيب(y؛أ،ب،ص،q،ج)=(أنا=1ن|أأنا|yأناأأناصأنا-1)(1-أنا=1ن(1-جأنا)(yأنابأنا)أأنا)q-1(أنا=1نبأناأأناصأنا)ب(ص1،...،صن،q)(1+أنا=1نجأنا(yأنابأنا)أأنا)ص¯+q{\displaystyle MGB(y;a,b,p,q,c)={\frac {(\prod _{i=1}^{n}|a_{i}|y_{i}^{a_{i}p_{i}-1})(1-\sum _{i=1}^{n}(1-c_{i})({\frac {y_{i}}{b_{i}}})^{a_{i}})^{q-1}}{(\prod _{i=1}^{n}b_{i}^{a_{i}p_{i}})B(p_{1},...,p_{n},q)(1+\sum _{i=1}^{n}c_{i}({\frac {y_{i}}{b_{i}}})^{a_{i}})^{{\bar {p}}+q}}}}

أين0{\displaystyle 0}<{\displaystyle <}أنا=1ن(1-جأنا)(yأنابأنا)أأنا{\displaystyle \sum _{i=1}^{n}(1-c_{i})({\frac {y_{i}}{b_{i}}})^{a_{i}}}<{\displaystyle <}1{\displaystyle 1}ل0{\displaystyle 0}{\displaystyle \leq }جأنا{\displaystyle c_{i}}<{\displaystyle <}1{\displaystyle 1}و0{\displaystyle 0}<{\displaystyle <}yأنا{\displaystyle y_{i}}متىجأنا{\displaystyle c_{i}}=1{\displaystyle 1}.

على غرار توزيع بيتا المعمم أحادي المتغير، يشمل توزيع بيتا المعمم متعدد المتغيرات عدة توزيعات في عائلته كحالات خاصة. ومن خلال فرض قيود معينة على متجهات المعلمات، يمكن اشتقاق التوزيعات التالية بسهولة. [ 12 ]

بيتا المعمم متعدد المتغيرات من النوع الأول (MGB1)

عندما كلجأنا{\displaystyle c_{i}}إذا كانت قيمة MGB تساوي صفرًا، فإن دالة MGB تتبسط إلى بيتا المعممة متعددة المتغيرات من النوع الأول (MGB1)، والتي تُعرَّف على النحو التالي:

مجيب1(y؛أ،ب،ص،q)=(أنا=1ن|أأنا|yأناأأناصأنا-1)(1-أنا=1ن(yأنابأنا)أأنا)q-1(أنا=1نبأناأأناصأنا)ب(ص1،...،صن،q){\displaystyle MGB1(y;a,b,p,q)={\frac {(\prod _{i=1}^{n}|a_{i}|y_{i}^{a_{i}p_{i}-1})(1-\sum _{i=1}^{n}({\frac {y_{i}}{b_{i}}})^{a_{i}})^{q-1}}{(\prod _{i=1}^{n}b_{i}^{a_{i}p_{i}})B(p_{1},...,p_{n},q)}}}

أين0{\displaystyle 0}<{\displaystyle <}أنا=1ن(yأنابأنا)أأنا{\displaystyle \sum _{i=1}^{n}({\frac {y_{i}}{b_{i}}})^{a_{i}}}<{\displaystyle <}1{\displaystyle 1}.

بيتا المعمم متعدد المتغيرات من النوع الثاني (MGB2)

في حالة كلجأنا{\displaystyle c_{i}}إذا كانت قيمة MGB تساوي 1، فإنها تتبسط إلى بيتا المعممة متعددة المتغيرات من النوع الثاني (MGB2)، مع تعريف دالة كثافة الاحتمال أدناه:

مجيب2(y؛أ،ب،ص،q)=(أنا=1ن|أأنا|yأناأأناصأنا-1)(أنا=1نبأناأأناصأنا)ب(ص1،...،صن،q)(1+أنا=1ن(yأنابأنا)أأنا)ص¯+q{\displaystyle MGB2(y;a,b,p,q)={\frac {(\prod _{i=1}^{n}|a_{i}|y_{i}^{a_{i}p_{i}-1})}{(\prod _{i=1}^{n}b_{i}^{a_{i}p_{i}})B(p_{1},...,p_{n},q)(1+\sum _{i=1}^{n}({\frac {y_{i}}{b_{i}}})^{a_{i}})^{{\bar {p}}+q}}}}

متى0{\displaystyle 0}<{\displaystyle <}yأنا{\displaystyle y_{i}}للجميعyأنا{\displaystyle y_{i}}.

جاما المعمم متعدد المتغيرات

يمكن اشتقاق دالة كثافة الاحتمال لتوزيع جاما المعمم متعدد المتغيرات (MGG) من دالة كثافة الاحتمال لتوزيع جاما المعمم متعدد المتغيرات (MGB) عن طريق الاستبدالبأنا{\displaystyle b_{i}}=βأناq1أأنا{\displaystyle \beta _{i}q^{\frac {1}{a_{i}}}}وباعتبار الحد كماq{\displaystyle q}{\displaystyle \to }{\displaystyle \infty }، مع تقريب ستيرلينغ لدالة غاما، مما ينتج عنه الدالة التالية:

مجيجي(y؛أ،β،ص)=((أنا=1ن|أأنا|yأناأأناصأنا-1)(أنا=1نβأناأأناصأنا)Γ(صأنا))هـ-أنا=1ن(yأناβأنا)أأنا=أنا=1نجيجي(yأنا؛أأنا،βأنا،صأنا){\displaystyle MGG(y;a,\beta ,p)=({\frac {(\prod _{i=1}^{n}|a_{i}|y_{i}^{a_{i}p_{i}-1})}{(\prod _{i=1}^{n}\beta _{i}^{a_{i}p_{i}})\Gamma (p_{i})}})e^{-\sum _{i=1}^{n}({\frac {y_{i}}{\beta _{i}}})^{a_{i}}}=\prod _{i=1}^{n}GG(y_{i};a_{i},\beta _{i},p_{i})}

وهو ناتج عن متغيرات عشوائية غاما معممة موزعة بشكل مستقل ولكن ليس بالضرورة بشكل متطابق.

توزيعات متعددة المتغيرات أخرى

يمكن إنشاء دوال كثافة احتمالية مماثلة لمتغيرات أخرى في شجرة العائلة الموضحة أعلاه، وذلك ببساطة عن طريق وضع الحرف M أمام اسم كل دالة كثافة احتمالية، ثم إيجاد الحالات الحدية والخاصة المناسبة لتوزيع MGB وفقًا للقيود والحدود المحددة للتوزيع أحادي المتغير. وتشمل دوال كثافة الاحتمالية متعددة المتغيرات الإضافية في المراجع توزيع ديريشليه (الصيغة القياسية) الموضح في:مجيب1(y؛أ=1،ب=1،ص،q){\displaystyle MGB1(y;a=1,b=1,p,q)}، التوزيع بيتا المعكوس متعدد المتغيرات وتوزيع ديريشليه المعكوس (نوع ديريشليه 2) المعطى بواسطةمجيب2(y؛أ=1،ب=1،ص،q){\displaystyle MGB2(y;a=1,b=1,p,q)}وتوزيع بور متعدد المتغيرات المعطى بواسطةمجيب2(y؛أ،ب،ص،q=1){\displaystyle MGB2(y;a,b,p,q=1)}.

دوال الكثافة الحدية

إن دوال الكثافة الحدية لـ MGB1 و MGB2، على التوالي، هي توزيعات بيتا المعممة من النوع الأول والثاني، ويتم تقديمها على النحو التالي:

جيب1(yأنا؛أأنا،بأنا،صأنا،ص¯-صأنا+q)=|أأنا|yأناأأناصأنا-1(1-(yأنابأنا)أأنا)ص¯-صأنا+q-1بأناأأناصأناب(صأنا،ص¯-صأنا+q){\displaystyle GB1(y_{i};a_{i},b_{i},p_{i},{\bar {p}}-p_{i}+q)={\frac {|a_{i}|y_{i}^{a_{i}p_{i}-1}(1-({\frac {y_{i}}{b_{i}}})^{a_{i}})^{{\bar {p}}-p_{i}+q-1}}{b_{i}^{a_{i}p_{i}}B(p_{i},{\bar {p}}-p_{i}+q)}}}
جيب2(yأنا؛أأنا،بأنا،صأنا،q)=|أأنا|yأناأأناصأنا-1بأناأأناصأناب(صأنا،q)(1+(yأنابأنا)أأنا)صأنا+q{\displaystyle GB2(y_{i};a_{i},b_{i},p_{i},q)={\frac {|a_{i}|y_{i}^{a_{i}p_{i}-1}}{b_{i}^{a_{i}p_{i}}B(p_{i},q)(1+({\frac {y_{i}}{b_{i}}})^{a_{i}})^{p_{i}+q}}}}

التطبيقات

تُستخدم المرونة التي توفرها عائلة GB في نمذجة توزيع:

  • توزيع الدخل
  • وظائف المخاطر
  • عوائد المخزون
  • خسائر التأمين

تشمل التطبيقات التي تتضمن أعضاء عائلة EGB ما يلي: [ 1 ] [ 6 ]

  • التقدير التكيفي الجزئي لنماذج الانحدار
  • نماذج السلاسل الزمنية
  • نماذج (G)ARCH

توزيع الدخل

استُخدم نموذج GB2 والعديد من حالاته الخاصة والمحدودة على نطاق واسع كنماذج لتوزيع الدخل. للاطلاع على بعض الأمثلة المبكرة، انظر: Thurow (1970)، [ 13 ] Dagum (1977)، [ 14 ] Singh and Maddala (1976)، [ 15 ] وMcDonald (1984). [ 6 ] ويمكن إجراء تقديرات الاحتمال الأقصى بسهولة باستخدام بيانات فردية أو مجمعة أو مُرمّزة أعلى مستوى مع هذه التوزيعات.

يمكن التعبير عن مقاييس عدم المساواة، مثل مؤشر جيني (G) ومؤشر بيترا (P) ومؤشر ثيل (T) من حيث معلمات التوزيع، كما ورد في ماكدونالد ورانسوم (2008): [ 16 ]

جي=(12μ)هـ(|Y-X|)=(P12μ)00|x-y|و(x)و(y)دxدy=1-0(1-F(y))2دy0(1-F(y))دyP=(12μ)هـ(|Y-μ|)=(12μ)0|y-μ|و(y)دyتي=هـ(ln(Y/μ)Y/μ)=0(y/μ)ln(y/μ)و(y)دy{\displaystyle {\begin{aligned}G=\left({\frac {1}{2\mu }}\right)\operatorname {E} (|Y-X|)=\left(P{\frac {1}{2\mu }}\right)\int _{0}^{\infty }\int _{0}^{\infty }|x-y|f(x)f(y)\,dxdy\\=1-{\frac {\int _{0}^{\infty }(1-F(y))^{2}\,dy}{\int _{0}^{\infty }(1-F(y))\,dy}}\\P=\left({\frac {1}{2\mu }}\right)\operatorname {E} (|Y-\mu |)=\left({\frac {1}{2\mu }}\right)\int _{0}^{\infty }|y-\mu |f(y)\,dy\\T=\operatorname {E} (\ln(Y/\mu )^{Y/\mu })=\int _{0}^{\infty }(y/\mu )\ln(y/\mu )f(y)\,dy\end{aligned}}}

وظائف المخاطر

دالة المخاطرة ، h(s)، حيث f(s) هي دالة كثافة الاحتمال وF(s) هي دالة التوزيع التراكمي المقابلة، تُعرَّف على النحو التالي:

ح(s)=و(s)1-F(s){\displaystyle h(s)={\frac {f(s)}{1-F(s)}}}

تُعدّ دوال المخاطر مفيدة في العديد من التطبيقات، مثل نمذجة مدة البطالة، ووقت فشل المنتجات، ومتوسط ​​العمر المتوقع. على سبيل المثال، إذا رمزنا لطول العمر بـ s، فإن h(s) هو معدل الوفاة عند العمر s، بافتراض أن الفرد قد عاش حتى العمر s. قد يكون شكل دالة المخاطر لبيانات وفيات البشر كما يلي: انخفاض معدل الوفيات في الأشهر الأولى من العمر، ثم فترة من معدل وفيات ثابت نسبيًا، وأخيرًا ازدياد احتمالية الوفاة مع التقدم في العمر.

تُتيح الحالات الخاصة لتوزيع بيتا المعمم مرونةً أكبر في نمذجة شكل دالة المخاطر، والتي قد تتطلب أشكالًا على هيئة "∪" أو "∩"، أو خطوطًا متزايدة تمامًا (يرمز لها بـ I) أو متناقصة تمامًا (يرمز لها بـ D). يكون توزيع غاما المعمم على هيئة "∪" عندما يكون a>1 و p<1/a، وعلى هيئة "∩" عندما يكون a<1 و p>1/a، وعلى هيئة I عندما يكون a>1 و p>1/a، وعلى هيئة D عندما يكون a<1 و p>1/a. [ 17 ] يُلخص ذلك في الشكل أدناه. [ 18 ] [ 19 ]

أشكال دالة المخاطر المحتملة باستخدام جاما المعممة

مراجع

  1. 1 2 3 ماكدونالد، جيمس ب. وشو، يكسياو ج. (1995) "تعميم لتوزيع بيتا مع تطبيقات"، مجلة الاقتصاد القياسي ، 66(1-2)، 133-152 doi : 10.1016/0304-4076(94)01612-4
  2. باتيل، جي بي، بوسويل، إم تي، وراتنابارخي، إم في، قاموس وفهرس مصنف للتوزيعات الإحصائية في سلسلة الأعمال العلمية، محرر جي بي باتيل، دار النشر التعاونية الداخلية، بيرتونزفيل، ماريلاند، 1984.
  3. فينتر، جي، توزيعات بيتا وجاما المحولة والخسائر الإجمالية، وقائع جمعية الخبراء الاكتواريين في مجال الحوادث، 1983.
  4. كالفلايش، جيه دي وبرينتيس، آر إل، التحليل الإحصائي لبيانات وقت الفشل، نيويورك: جيه وايلي، 1980
  5. أرنولد، بي سي، توزيعات باريتو، المجلد 5 في التوزيعات الإحصائية في سلسلة الأعمال العلمية، دار النشر التعاونية الدولية، بيرتونزفيل، ماريلاند 1983.
  6. 1 2 3 ماكدونالد، جيه بي (1984) "بعض الدوال المعممة لتوزيعات حجم الدخل"، Econometrica 52، 647-663.
  7. ستيوارت، أ. وأورد، جيه كيه (1987): نظرية كيندال المتقدمة في الإحصاء، نيويورك: مطبعة جامعة أكسفورد.
  8. ستايسي، إي دبليو (1962). "تعميم لتوزيع غاما". حوليات الإحصاء الرياضي 33(3): 1187-1192. JSTOR 2237889 
  9. ريد، دبليو جيه (2001). "قوانين باريتو، وزيبف، وغيرها من قوانين القوة". رسائل اقتصادية 74: 15-19. doi : 10.1016/S0165-1765(01)00524-9
  10. هيغبي، جيه دي، وجينسن، جيه إي، وماكدونالد، جيه بي (2019). "توزيع لابلاس اللوغاريتمي غير المتماثل كحالة حدية لتوزيع بيتا المعمم." رسائل الإحصاء والاحتمالات 151: 73-78. doi : 10.1016/j.spl.2019.03.018
  11. ماكدونالد، جيمس ب. وكيرمان، شون س. (2013) "حدود الالتواء والتفرطح لـ EGB1 وEGB2 والحالات الخاصة"، قيد النشر
  12. ويليام إم. كوكريل وجيمس بي. ماكدونالد (2017): عائلتان من بيتا المعممة متعددة المتغيرات، الاتصالات في الإحصاء - النظرية والأساليب، doi : 10.1080/03610926.2017.1400058
  13. ثورو، إل سي (1970) "تحليل توزيع الدخل الأمريكي"، أوراق ومحاضر، الجمعية الاقتصادية الأمريكية ، 60، 261-269
  14. داغوم، سي. (1977) "نموذج جديد لتوزيع الدخل الشخصي: التحديد والتقدير"، الاقتصاد التطبيقي ، 30، 413-437
  15. سينغ، إس كيه ومادالا، جي إس (1976) "دالة لتوزيع حجم الدخول"، Econometrica ، 44، 963-970
  16. ماكدونالد، جيه بي ورانسوم، إم. (2008) "توزيع بيتا المعمم كنموذج لتوزيع الدخل: تقدير مقاييس عدم المساواة ذات الصلة"، نمذجة التوزيعات ومنحنيات لورنز ، "دراسات اقتصادية في عدم المساواة: الإقصاء الاجتماعي والرفاه"، سبرينغر: نيويورك، المحرر جاك سيلبر، 5، 147-166
  17. جلاسر، رونالد إي. (1980) "أحواض الاستحمام وخصائص معدل الفشل ذات الصلة"، مجلة الجمعية الإحصائية الأمريكية ، 75(371)، 667-672 doi : 10.1080/01621459.1980.10477530
  18. ماكدونالد، جيمس ب. (1987) "منهجية عامة لتحديد أشكال التوزيع مع تطبيقات في الموثوقية"، مجلة التخطيط والاستدلال الإحصائي ، 16، 365-376، doi : 10.1016/0378-3758(87)90089-9
  19. ماكدونالد، جيه بي وريتشاردز، دي أو (1987) "دوال الخطر وتوزيعات بيتا المعممة"، معاملات IEEE في الموثوقية ، 36، 463-466

فهرس

  • سي. كلايبر وإس. كوتز (2003) التوزيعات الإحصائية للأحجام في الاقتصاد والعلوم الاكتوارية . نيويورك: وايلي
  • جونسون، إن إل، إس. كوتز، وإن. بالاكريشنان (1994) التوزيعات أحادية المتغير المستمرة . المجلد 2، هوبوكين، نيوجيرسي: وايلي-إنترساينس.