逆威沙特分布

✍ dations ◷ 2025-10-22 23:35:35 #连续分布,多变量统计


逆威沙特分布,也叫反威沙特分布作是统计学中出现的一类概率分布函数,定义在实值的正定矩阵上。在贝叶斯统计中,逆威沙特分布会用作多变量正态分布协方差矩阵的共轭先验分布。如果一个正定矩阵 B {\displaystyle {\mathbf {B} }} (·) 则是多变量伽马分布(英语:Multivariate gamma function)。函数

指的是迹函数。

设矩阵 A W ( Σ , m ) {\displaystyle {\mathbf {A} }\sim W({\mathbf {\Sigma } },m)} 并且 Σ {\displaystyle {\mathbf {\Sigma } }} p × p {\displaystyle p\times p} 的矩阵,那么 B = A 1 {\displaystyle {\mathbf {B} }={\mathbf {A} }^{-1}} 遵从逆威沙特分布: B W 1 ( Σ 1 , m ) {\displaystyle {\mathbf {B} }\sim W^{-1}({\mathbf {\Sigma } }^{-1},m)} 。它的概率密度函数是:

其中 Ψ = Σ 1 {\displaystyle \mathbf {\Psi } =\mathbf {\Sigma } ^{-1}} ,而 Γ p ( ) {\displaystyle \Gamma _{p}(\cdot )} 是多变量伽马分布。

设矩阵 A W 1 ( Ψ , m ) {\displaystyle {\mathbf {A} }\sim W^{-1}({\mathbf {\Psi } },m)} 遵从逆威沙特分布。并且假设矩阵 A {\displaystyle {\mathbf {A} }} Ψ {\displaystyle {\mathbf {\Psi } }} 都有相适合的分块矩阵表示方式:

其中子矩阵 A i j {\displaystyle {\mathbf {A} _{ij}}} Ψ i j {\displaystyle {\mathbf {\Psi } _{ij}}} p i × p j {\displaystyle p_{i}\times p_{j}} 的矩阵,那么会有:

甲) A 11 {\displaystyle {\mathbf {A} _{11}}} A 11 1 A 12 {\displaystyle {\mathbf {A} }_{11}^{-1}{\mathbf {A} }_{12}} A 22 1 {\displaystyle {\mathbf {A} }_{22\cdot 1}} 相互独立,其中 A 22 1 = A 22 A 21 A 11 1 A 12 {\displaystyle {\mathbf {A} _{22\cdot 1}}={\mathbf {A} }_{22}-{\mathbf {A} }_{21}{\mathbf {A} }_{11}^{-1}{\mathbf {A} }_{12}} 是子矩阵 A 11 {\displaystyle {\mathbf {A} _{11}}} A {\displaystyle {\mathbf {A} }} 中的舒尔补。

乙) A 11 W 1 ( Ψ 11 , m p 2 ) {\displaystyle {\mathbf {A} _{11}}\sim W^{-1}({\mathbf {\Psi } _{11}},m-p_{2})} ;

丙) A 11 1 A 12 | A 22 1 M N p 1 × p 2 ( Ψ 11 1 Ψ 12 , A 22 1 Ψ 11 1 ) {\displaystyle {\mathbf {A} }_{11}^{-1}{\mathbf {A} }_{12}|{\mathbf {A} }_{22\cdot 1}\sim MN_{p_{1}\times p_{2}}({\mathbf {\Psi } }_{11}^{-1}{\mathbf {\Psi } }_{12},{\mathbf {A} }_{22\cdot 1}\otimes {\mathbf {\Psi } }_{11}^{-1})} ,其中 M N p × q ( , ) {\displaystyle MN_{p\times q}(\cdot ,\cdot )} 是矩阵正态分布。

丁) A 22 1 W 1 ( Ψ 22 1 , m ) {\displaystyle {\mathbf {A} }_{22\cdot 1}\sim W^{-1}({\mathbf {\Psi } }_{22\cdot 1},m)}

假设要求先验分布 p ( Σ ) {\displaystyle {p(\mathbf {\Sigma } )}} 为逆威沙特分布 W 1 ( Ψ , m ) {\displaystyle W^{-1}({\mathbf {\Psi } },m)} 的协方差矩阵 Σ {\displaystyle {\mathbf {\Sigma } }} 。如果观测值 X = {\displaystyle \mathbf {X} =} 是从互相独立的 p-变量正态分布 N ( 0 , Σ ) {\displaystyle N(\mathbf {0} ,{\mathbf {\Sigma } })} 的随机变量得到的,那么条件分布 p ( Σ | X ) {\displaystyle {p(\mathbf {\Sigma } |\mathbf {X} )}} 遵从的是逆威沙特分布: W 1 ( A + Ψ , n + m ) {\displaystyle W^{-1}({\mathbf {A} }+{\mathbf {\Psi } },n+m)} 。其中 A = X X T {\displaystyle {\mathbf {A} }=\mathbf {X} \mathbf {X} ^{T}} 是样本协方差矩阵的 n {\displaystyle n} 倍。

因此,逆威沙特矩阵是多变量正态分布的共轭先验分布。

期望值::85

矩阵 B {\displaystyle \mathbf {B} } 的每一个系数的方差:

对角系数的方差是在上式中令 i = j {\displaystyle i=j} 得到,化简后变成:

当变量数目减到一个的时候,逆威沙特分布会变成特例:逆伽马分布(英语:Inverse-gamma distribution)。也就是说,当 p = 1 {\displaystyle p=1} α = m / 2 {\displaystyle \alpha =m/2} β = Ψ / 2 {\displaystyle \beta =\mathbf {\Psi } /2} 以及 x = B {\displaystyle x=\mathbf {B} } 的时候,逆威沙特分布的概率密度函数是:


这正是逆伽马分布。其中 Γ 1 ( ) {\displaystyle \Gamma _{1}(\cdot )} 是通常的伽马函数。


而逆威沙特分布也有推广,其中一个是正态逆威沙特分布(英语:Normal-inverse-Wishart distribution)。

相关

  • 扁桃体炎扁桃体炎(Tonsillitis),或称扁桃腺炎,通常会快速发病。扁桃体炎属于咽炎的一种。其症状包括咽喉痛、发烧、扁桃腺肿大、吞咽困难、颈部的淋巴结肿大(英语:Lymphadenopathy)。并发症
  • 恋尿淫乐排尿(英语:urolagnia),又称恋尿癖(英语:urophilia),是一种恋物癖形式,指人对尿感到特别的偏好。在这项领域中,会以圣水(英语:golden shower)代指支配者或施虐者排出的尿液。
  • 中士中士是军人的职衔,士官的一种,在大多数国家的军衔制度中,中士之下的阶级为下士,中士之上的阶级为上士。军队编制上多以士官长、上士担任副排长,中士担任班长,下士担任副班长,下辖三
  • 祖籍祖籍,是指一个家族族群以父系共同认定的一位祖先的生长地。祖籍与籍贯、户籍的意义不太相同,祖籍通常是追认极遥远的祖先生长地,籍贯通常是指父亲或祖父的长居之地,户籍则是其本
  • 网络存取互联网存取(英语:Internet access),又译为互联网接入,俗称上网、连网、连线,将个人电脑终端机、个人电脑、移动设备与电脑网络,透过各种连结技术与管道,连结到互联网(Internet)上,让使
  • 阿卜杜勒·卡迪尔·汗阿卜杜勒·卡迪尔·汗(乌尔都语:‏عبد القدیر خان‎‎,1936年4月1日-),出生于英属印度博帕尔(今属中央邦)。他作为巴基斯坦核技术及冶金工程领域的专家,并掌握着铀浓缩离
  • 七冲门饮食从进入人体内到排出体外,共需经过七个主要部位,这七个部位在人体饮食的代谢过程中就如七个关卡,在中医中称作“七冲门”。《难经》中描述这“七冲门”为:口唇是饮食进入人体
  • 田厦路田厦路(英语:Tin Ha Road),位于元朗区的西部,连接屏厦路近厦村市及青山公路-洪水桥段近田心,是一条南北行双向道路。田厦路近青山公路-洪水桥段交界处一段是元朗区及屯门区的分界线,
  • 周芝沅周芝沅(1783年-1835年),表字荣坤,号逸坡。江苏吴江人。乾隆工部尚书周元理之后,周元德之孙,周一士之子。监生出身,入粟得候选江西按察使司知事。生平乐善好施。道光三年(1823年)夏雨水
  • BPEL4PeopleBPEL4People的全称是WS-BPEL Extension for People,是BPEL在人工活动方面的扩展。2005年7月,IBM和SAP在一个联合白皮书中提出BPEL4People。2007年6月,Active Endpoints, Adobe,