Teskari-Wishart taqsimoti - Inverse-Wishart distribution

Teskari-istak

Notation	${displaystyle {mathcal {W}} ^ {- 1} ({mathbf {Psi}}, u)}$
Parametrlar	${displaystyle u> p-1}$ erkinlik darajasi (haqiqiy ) ${displaystyle mathbf {Psi}> 0}$, ${displaystyle p imes p}$ o'lchov matritsasi (pos. def. )
Qo'llab-quvvatlash	${displaystyle mathbf {X}}$ bu p × p ijobiy aniq
PDF	${displaystyle {frac {left | mathbf {Psi} ight | ^ {u / 2}} {2 ^ {up / 2} Gamma _ {p} ({frac {u} {2}})}} left | mathbf { x} ight | ^ {- (u + p + 1) / 2} e ^ {- {frac {1} {2}} operator nomi {tr} (mathbf {Psi} mathbf {x} ^ {- 1})} }$ ${displaystyle Gamma _ {p}}$ bo'ladi ko'p o'zgaruvchan gamma funktsiyasi ${displaystyle operator nomi {tr}}$ bo'ladi iz funktsiya
Anglatadi	${displaystyle {frac {mathbf {Psi}} {u -p-1}}}$Uchun ${displaystyle u> p + 1}$
Rejim	${displaystyle {frac {mathbf {Psi}} {u + p + 1}}}$[1]:406
Varians	pastga qarang

Yilda statistika, Wishart-ning teskari taqsimoti, shuningdek teskari Wishart tarqatish, a ehtimollik taqsimoti real qiymat bo'yicha aniqlangan ijobiy-aniq matritsalar. Yilda Bayes statistikasi u sifatida ishlatiladi oldingi konjugat kovaryans matritsasi uchun a ko'p o'zgaruvchan normal tarqatish.

Biz aytamiz ${displaystyle mathbf {X}}$ sifatida belgilangan teskari Wishart taqsimotiga amal qiladi ${displaystyle mathbf {X} sim {mathcal {W}} ^ {- 1} (mathbf {Psi}, u)}$, agar u bo'lsa teskari ${displaystyle mathbf {X} ^ {- 1}}$ bor Istaklarni tarqatish ${displaystyle {mathcal {W}} (mathbf {Psi} ^ {- 1}, u)}$. Teskari-Wishart taqsimoti uchun muhim identifikatorlar olingan.^[2]

Zichlik

The ehtimollik zichligi funktsiyasi teskari Wishart:^[3]

${displaystyle f_ {mathbf {x}} ({mathbf {x}}; {mathbf {Psi}}, u) = {frac {left | {mathbf {Psi}} ight | ^ {u / 2}} {2 ^ {up / 2} Gamma _ {p} ({frac {u} {2}})}} chap | mathbf {x} ight | ^ {- (u + p + 1) / 2} e ^ {- {frac {1} {2}} operator nomi {tr} (mathbf {Psi} mathbf {x} ^ {- 1})}}$

qayerda ${displaystyle mathbf {x}}$ va ${displaystyle {mathbf {Psi}}}$ bor ${displaystyle p imes p}$ ijobiy aniq matritsalar va Γ_p(·) bo'ladi ko'p o'zgaruvchan gamma funktsiyasi.

Teoremalar

Wishart tomonidan taqsimlangan matritsaning teskari taqsimoti

Agar ${displaystyle {mathbf {X}} sim {mathcal {W}} ({mathbf {Sigma}}, u)}$ va ${displaystyle {mathbf {Sigma}}}$ hajmi ${displaystyle p imes p}$, keyin ${displaystyle mathbf {A} = {mathbf {X}} ^ {- 1}}$ teskari Wishart taqsimotiga ega ${displaystyle mathbf {A} sim {mathcal {W}} ^ {- 1} ({mathbf {Sigma}} ^ {- 1}, u)}$ .^[4]

Teskari Wishart-taqsimlangan matritsadan marginal va shartli taqsimotlar

Aytaylik ${displaystyle {mathbf {A}} sim {mathcal {W}} ^ {- 1} ({mathbf {Psi}}, u)}$ teskari Wishart taqsimotiga ega. Matritsalarni taqsimlash ${displaystyle {mathbf {A}}}$ va ${displaystyle {mathbf {Psi}}}$ mos ravishda bir-birlari bilan

${displaystyle {mathbf {A}} = {egin {bmatrix} mathbf {A} _ {11} & mathbf {A} _ {12} mathbf {A} _ {21} & mathbf {A} _ {22} end {bmatrix }},; {mathbf {Psi}} = {egin {bmatrix} mathbf {Psi} _ {11} & mathbf {Psi} _ {12} mathbf {Psi} _ {21} & mathbf {Psi} _ {22} end {bmatrix}}}$

qayerda ${displaystyle {mathbf {A} _ {ij}}}$ va ${displaystyle {mathbf {Psi} _ {ij}}}$ bor ${displaystyle p_ {i} imes p_ {j}}$ matritsalar, keyin bizda bor

i) ${displaystyle mathbf {A} _ {11}}$ dan mustaqildir ${displaystyle mathbf {A} _ {11} ^ {- 1} mathbf {A} _ {12}}$ va ${displaystyle {mathbf {A}} _ {22cdot 1}}$, qayerda ${displaystyle {mathbf {A} _ {22cdot 1}} = {mathbf {A}} _ {22} - {mathbf {A}} _ {21} {mathbf {A}} _ {11} ^ {- 1} {mathbf {A}} _ {12}}$ bo'ladi Schur to'ldiruvchisi ning ${displaystyle {mathbf {A} _ {11}}}$ yilda ${displaystyle {mathbf {A}}}$;

ii) ${displaystyle {mathbf {A} _ {11}} sim {mathcal {W}} ^ {- 1} ({mathbf {Psi} _ {11}}, u -p_ {2})}$;

iii) ${displaystyle {mathbf {A}} _ {11} ^ {- 1} {mathbf {A}} _ {12} mid {mathbf {A}} _ {22cdot 1} sim MN_ {p_ {1} imes p_ {2 }} ({mathbf {Psi}} _ {11} ^ {- 1} {mathbf {Psi}} _ {12}, {mathbf {A}} _ {22cdot 1} otimes {mathbf {Psi}} _ {11 } ^ {- 1})}$, qayerda ${displaystyle MN_ {p imes q} (cdot, cdot)}$ a matritsaning normal taqsimlanishi;

iv) ${displaystyle {mathbf {A}} _ {22cdot 1} sim {mathcal {W}} ^ {- 1} ({mathbf {Psi}} _ {22cdot 1}, u)}$, qayerda ${displaystyle {mathbf {Psi} _ {22cdot 1}} = {mathbf {Psi}} _ {22} - {mathbf {Psi}} _ {21} {mathbf {Psi}} _ {11} ^ {- 1} {mathbf {Psi}} _ {12}}$;

Konjugatning tarqalishi

Kovaryans matritsasi haqida xulosa chiqarishni xohlaymiz ${displaystyle {mathbf {Sigma}}}$ kimning oldin ${displaystyle {p (mathbf {Sigma})}}$ bor ${displaystyle {mathcal {W}} ^ {- 1} ({mathbf {Psi}}, u)}$ tarqatish. Agar kuzatishlar bo'lsa ${displaystyle mathbf {X} = [mathbf {x} _ {1}, ldots, mathbf {x} _ {n}]}$ a dan olingan mustaqil p-variatsion Gauss o'zgaruvchilari ${displaystyle N (mathbf {0}, {mathbf {Sigma}})}$ taqsimot, keyin shartli taqsimot ${displaystyle {p (mathbf {Sigma} mid mathbf {X})}}$ bor ${displaystyle {mathcal {W}} ^ {- 1} ({mathbf {A}} + {mathbf {Psi}}, n + u)}$ tarqatish, qaerda ${displaystyle {mathbf {A}} = mathbf {X} mathbf {X} ^ {T}}$.

Oldingi va orqa taqsimotlari bir xil oila bo'lgani uchun, biz teskari Wishart taqsimoti deymiz birlashtirmoq ko'p o'zgaruvchan Gaussga.

Ko'p o'zgaruvchan Gauss bilan konjugatsiyasi tufayli, buni amalga oshirish mumkin marginallashtirish (integratsiya) Gauss parametrini ${displaystyle mathbf {Sigma}}$.

${displaystyle f_ {mathbf {X}, mid, Psi, u} (mathbf {x}) = int f_ {mathbf {X}, mid, mathbf {Sigma}, =, sigma} (mathbf {x}) f_ {mathbf {Sigma}, mid, mathbf {Psi}, u} (sigma), dsigma = {frac {| mathbf {Psi} | ^ {u / 2} Gamma _ {p} chap ({frac {u + n} {2) }} kech)} {pi ^ {np / 2} | mathbf {Psi} + mathbf {A} | ^ {(u + n) / 2} Gamma _ {p} ({frac {u} {2}}) }}}$

(bu foydali, chunki dispersiya matritsasi ${displaystyle mathbf {Sigma}}$ amalda ma'lum emas, lekin chunki ${displaystyle {mathbf {Psi}}}$ ma'lum aprioriva ${displaystyle {mathbf {A}}}$ ma'lumotlardan olish mumkin, o'ng tomonni to'g'ridan-to'g'ri baholash mumkin). Oldindan teskari-Wishart taqsimoti mavjud o'tkazilgan orqali tuzilishi mumkin oldingi bilim.^[5]

Lahzalar

Quyidagilar Press, S. J. (1982) "Amaliy ko'p o'zgaruvchan tahlil", 2-nashrga asoslangan. (Dover Publications, Nyu-York), erkinlik darajasini p.d.f.ga mos keladigan darajada qayta o'rnatgandan so'ng. yuqoridagi ta'rif.

O'rtacha:^[4]:85

${displaystyle operator nomi {E} (mathbf {X}) = {frac {mathbf {Psi}} {u -p-1}}.}$

Ning har bir elementining dispersiyasi ${displaystyle mathbf {X}}$:

${displaystyle operator nomi {Var} (x_ {ij}) = {frac {(u -p + 1) psi _ {ij} ^ {2} + (u -p-1) psi _ {ii} psi _ {jj} } {(u -p) (u -p-1) ^ {2} (u -p-3)}}}$

Diagonalning dispersiyasi yuqoridagi formuladan foydalanadi ${displaystyle i = j}$, bu quyidagilarni soddalashtiradi:

${displaystyle operator nomi {Var} (x_ {ii}) = {frac {2psi _ {ii} ^ {2}} {(u -p-1) ^ {2} (u -p-3)}}.}$

Elementlarining kovaryansiyasi ${displaystyle mathbf {X}}$ quyidagilar tomonidan beriladi:

${displaystyle operator nomi {Cov} (x_ {ij}, x_ {kell}) = {frac {2psi _ {ij} psi _ {kell} + (u -p-1) (psi _ {ik} psi _ {jell} + psi _ {iell} psi _ {kj})} {(u -p) (u -p-1) ^ {2} (u -p-3)}}}$

Natijalar fon Rozen tomonidan aniqroq Kronecker mahsuloti shaklida ifodalangan^[6] quyidagicha.

${displaystyle mathbf {E} chap (W ^ {- 1} otimes W ^ {- 1} ight) = c_ {1} Psi otimes Psi + c_ {2} Vec (Psi) Vec (Psi) ^ {T} + c_ {2} K_ {pp} Psi otimes Psi}$

${displaystyle mathbf {Cov} chap (W ^ {- 1} otimes W ^ {- 1} ight) = (c_ {1} -c_ {3}) Psi otimes Psi + c_ {2} Vec (Psi) Vec (Psi) ) ^ {T} + c_ {2} K_ {pp} Psi otimes Psi}$

qayerda ${displaystyle c_ {2} = chap [(u -p) (u -p-1) (u -p-3) ight] ^ {- 1}, ;; c_ {1} = (u -p-2) c_ {2} ,; c_ {3} = (u -p-1) ^ {- 2}}$va ${displaystyle K_ {pp} {ext {- bu a}} p ^ {2} imes p ^ {2}}$ kommutatsiya matritsasi. Qog'ozda xatoning xatosi mavjud bo'lib, uning koeffitsienti ${displaystyle K_ {pp} Psi otimes Psi}$ sifatida berilgan ${displaystyle c_ {1}}$ dan ko'ra ${displaystyle c_ {2}}$. Shuningdek, 3.1-xulosa o'rtacha Wishart kvadratiga teskari ifodani o'qish kerak ${displaystyle mathbf {E} chap [W ^ {- 1} W ^ {- 1} ight] = (c_ {1} + c_ {2}) Sigma ^ {- 1} Sigma ^ {- 1} + c_ {2 } Sigma ^ {- 1} mathbf {tr} (Sigma ^ {- 1})}$

Kovaryans diagonali bo'lganda o'zaro ta'sir qiluvchi atamalar qanday qilib siyraklashishini ko'rsatish uchun, ruxsat bering ${displaystyle Psi = mathbf {I} _ {3 imes 3}}$ va ba'zi bir o'zboshimchalik parametrlarini kiritish ${displaystyle u, v, w}$:

${displaystyle mathbf {E} chap (W ^ {- 1} otimes W ^ {- 1} ight) = uPsi otimes Psi + vVec (Psi) Vec (Psi) ^ {T} + wK_ {pp} Psi otimes Psi}$

keyin ikkinchi moment matritsasi bo'ladi

${displaystyle mathbf {E} chap (W ^ {- 1} otimes W ^ {- 1} ight) = {egin {bmatrix} u + v + w & cdot & cdot & cdot & v & cdot & cdot & cdot & v cdot & u & cdot & w & cdot & cdot & cd & cd & cd & cdot & cdot & U & cdot & cdot & cdot & w & cdot & cdot cdot & w & cdot & U & cdot & cdot & cdot & cdot & cdot v & cdot & cdot & cdot & U + v + w & cdot & cdot & cdot & v cdot & cdot & cdot & cdot & cdot & U & cdot & w & cdot cdot & cdot & w & cdot & cdot & cdot & U & cdot & cdot cdot & cdot & cdot & cdot & cdot & w & cdot & U & cdot v & cdot & cdot & cdot & v & cdot & cdot & cdot & u + v + w end {bmatrix}}}$

Wishart mahsulotining farqlari, shuningdek, Kuk va boshq. al.^[7] alohida holatda va kengaytirilgan holda, to'liq martabali holatga qadar. Murakkab holatda "oq" teskari kompleks Wishart ${displaystyle {mathcal {W}} ^ {- 1} (mathbf {I}, u, p)}$ Shaman tomonidan namoyish etilgan^[8] etakchi diagonal elementlar o'zaro bog'liq bo'lgan diagonali statistik tuzilishga ega bo'lish, qolgan barcha elementlar esa o'zaro bog'liq emas. Bundan tashqari, uni Brennan va Rid ko'rsatgan^[9] matritsani ajratish protsedurasidan foydalangan holda, murakkab o'zgaruvchan domenda bo'lsa ham, ushbu matritsaning [1,1] diagonal elementining marginal pdf Teskari chi-kvadrat taqsimot. Bu barcha diagonal elementlarga osonlikcha tarqaladi ${displaystyle {mathcal {W}} ^ {- 1} (mathbf {I}, u, p)}$ diagonali elementlarning almashinuvini o'z ichiga olgan ortogonal transformatsiyalar ostida statistik o'zgarmasdir.

Teskari Chi kvadratik taqsimoti uchun ${displaystyle u _ {c}}$ erkinlik darajasi, pdf esa

${displaystyle {ext {Inv -}} chi ^ {2} (x; u _ {c}) = {frac {2 ^ {- u _ {c} / 2}} {Gamma (u _ {c} / 2) )}} x ^ {- u _ {c} / 2-1} e ^ {- 1 / (2x)}.}$

o'rtacha va o'zgaruvchanligi ${displaystyle {frac {1} {u _ {c} -2}} {ext {and}} {frac {2} {(u _ {c} -2) ^ {2} (u _ {c} -4) )}}}$ navbati bilan. Ushbu ikkita parametr mos keladigan teskari Wishart diagonal momentlariga mos keladi ${displaystyle u _ {c} = u -p + 1}$ va shuning uchun marginal pdf ning diagonal elementi ${displaystyle {mathcal {W}} ^ {- 1} (mathbf {I}, u, p)}$ bo'ladi:

${displaystyle f_ {x_ {11}} (x_ {11}; Psi, u, p) = {frac {2 ^ {- (u -p + 1) / 2}} {Gamma chap ({frac {u -p) +1} {2}} tun)}}, x_ {11} ^ {- (u -p + 1) / 2-1} e ^ {- 1 / (2x_ {11})}}$

quyida, barcha diagonal elementlar uchun umumlashtirilgan. E'tibor bering, murakkab teskari Wishart o'rtacha qiymati ${displaystyle {frac {mathbf {I}} {u -p}}}$ va haqiqiy bo'lgan Wishart ishidan farq qiladi ${displaystyle {frac {mathbf {I}} {u -p-1}}}$.

Tegishli tarqatishlar

A bir o'zgaruvchan teskari-Wishart taqsimotining ixtisoslashuvi teskari-gamma taqsimoti. Bilan ${displaystyle p = 1}$ (ya'ni bitta o'zgaruvchan) va ${displaystyle alfa = u / 2}$, ${displaystyle eta = mathbf {Psi} / 2}$ va ${displaystyle x = mathbf {X}}$ The ehtimollik zichligi funktsiyasi teskari-Wishart taqsimotiga aylanadi

${displaystyle p (xmid alfa, eta) = {frac {eta ^ {alfa}, x ^ {- alfa -1} exp (- eta / x)} {Gamma _ {1} (alfa)}}.}$

ya'ni teskari-gamma taqsimoti, bu erda ${displaystyle Gamma _ {1} (cdot)}$ oddiy Gamma funktsiyasi.

Teskari Wishart taqsimoti - bu alohida holat teskari matritsa gamma taqsimoti qachon shakl parametri ${displaystyle alfa = {frac {u} {2}}}$ va o'lchov parametri ${displaystyle eta = 2}$.

Boshqa bir umumlashma umumlashtirilgan teskari Wishart taqsimoti deb nomlandi, ${displaystyle {mathcal {GW}} ^ {- 1}}$. A ${displaystyle p imes p}$ ijobiy aniq matritsa ${displaystyle mathbf {X}}$ sifatida tarqatilishi aytilmoqda ${displaystyle {mathcal {GW}} ^ {- 1} (mathbf {Psi}, u, mathbf {S})}$ agar ${displaystyle mathbf {Y} = mathbf {X} ^ {1/2} mathbf {S} ^ {- 1} mathbf {X} ^ {1/2}}$ sifatida taqsimlanadi ${displaystyle {mathcal {W}} ^ {- 1} (mathbf {Psi}, u)}$. Bu yerda ${displaystyle mathbf {X} ^ {1/2}}$ ning nosimmetrik matritsasi kvadrat ildizini bildiradi ${displaystyle mathbf {X}}$, parametrlari ${displaystyle mathbf {Psi}, mathbf {S}}$ bor ${displaystyle p imes p}$ ijobiy aniq matritsalar va parametr ${displaystyle u}$ dan kattaroq ijobiy skalar hisoblanadi ${displaystyle 2p}$. Qachon ekanligini unutmang ${displaystyle mathbf {S}}$ identifikatsiya matritsasiga teng, ${displaystyle {mathcal {GW}} ^ {- 1} (mathbf {Psi}, u, mathbf {S}) = {mathcal {W}} ^ {- 1} (mathbf {Psi}, u)}$. Ushbu umumlashtirilgan teskari Wishart taqsimoti ko'p o'zgaruvchan avtoregressiv jarayonlarning taqsimlanishini baholash uchun qo'llanilgan.^[10]

Umumlashtirishning boshqa turi bu normal-teskari-Wishart taqsimoti, asosan a mahsuloti ko'p o'zgaruvchan normal taqsimot teskari Wishart taqsimoti bilan.

Agar o'lchov matritsasi identifikatsiya matritsasi bo'lsa, ${displaystyle {mathcal {Psi}} = mathbf {I}, {ext {and}} {mathcal {Phi}}}$ o'zboshimchalik bilan ortogonal matritsa, o'rnini bosish ${displaystyle mathbf {X}}$ tomonidan ${displaystyle {Phi} mathbf {X} {mathcal {Phi}} ^ {T}}$ ning pdf-ni o'zgartirmaydi ${displaystyle mathbf {X}}$ shunday ${displaystyle {mathcal {W}} ^ {- 1} (mathbf {I}, u, p)}$ qaysidir ma'noda sferik o'zgarmas tasodifiy jarayonlar (SIRPs) oilasiga tegishli.
Shunday qilib, o'zboshimchalik bilan p-vektor ${displaystyle V}$ bilan ${displaystyle l_ {2} {ext {length}} V ^ {T} V = 1}$ vektorga aylantirilishi mumkin ${displaystyle mathbf {Phi} V = [1; 0; 0cdots] ^ {T}}$ pdf-ni o'zgartirmasdan ${displaystyle V ^ {T} mathbf {X} V}$, bundan tashqari ${displaystyle mathbf {Phi}}$ diagonal elementlarni almashtiradigan almashtirish matritsasi bo'lishi mumkin. Bundan kelib chiqadiki, ning diagonal elementlari ${displaystyle mathbf {X}}$ pdf bilan bir xil teskari chi kvadratga taqsimlanadi ${displaystyle f_ {x_ {11}}}$ oldingi bo'limda bo'lsa ham, ular o'zaro mustaqil emas. Natija maqbul portfel statistikasida ma'lum, chunki Bodnar va boshqalarning 2-teoremasi, 1-xulosa.^[11] bu erda teskari shaklda ifodalangan ${displaystyle {frac {V ^ {T} mathbf {Psi} V} {V ^ {T} mathbf {X} V}} sim chi _ {u -p + 1} ^ {2}}$.

Shuningdek qarang

• Teskari matritsali gamma-taqsimot
• Matritsaning normal taqsimlanishi
• Istaklarni tarqatish
• Wishart-ning teskari taqsimoti

Adabiyotlar