逆威沙特分布

✍ dations ◷ 2025-02-24 08:02:39 #连续分布,多变量统计


逆威沙特分布,也叫反威沙特分布作是统计学中出现的一类概率分布函数,定义在实值的正定矩阵上。在贝叶斯统计中,逆威沙特分布会用作多变量正态分布协方差矩阵的共轭先验分布。如果一个正定矩阵 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)。

相关

  • 沙池沙池是一个装载沙的范围,可作为游戏或跳远之用。较小的又叫沙盒。《田径竞赛规则》第四节B部分第一百八十四条第六及第七款: 落地区宽度最小2.75米,最大3米。助跑道应对准落地
  • 维利·鲍迈斯特维利·鲍迈斯特(Willi Baumeister,1889年1月22日-1955年8月31日),德国画家、景观设计师、艺术教授、印刷工艺师。他从师于斯图加特大学赫尔策教授,早期的作品受法国绘画影响,色彩轻
  • 黄纬禄黄纬禄(1916年12月18日-2011年11月23日 ),安徽省芜湖人,自动控制和导弹技术专家,中国科学院资深院士,中国导弹与航天技术的主要开拓者之一,两弹一星功勋奖章获得者,曾参与开拓中国核
  • 陈伯谅陈伯谅(1482年-1526年),字执之,福建福州府福清县人。明朝政治人物、进士出身。弘治十七年(1504年)甲子乡试第六名,正德三年(1508年)登戊辰科三甲一百四十九名。初授浙江富阳知县,官至四
  • 芒齿小檗芒齿小檗(学名:)为小檗科小檗属的植物,为中国的特有植物。分布于中国大陆的贵州、四川、陕西、湖南、湖北等地,生长于海拔500米至2,100米的地区,一般生于山坡杂木林中,目前尚未由人
  • 奇异化学志《奇异化学志》(拉丁语:,直译“奇异化学图书馆”),又称《惊异化学记》或《灵妙化学丛书》,是一部炼金术文献合集,于1702年在日内瓦出版,由瑞士医师兼作家尚-贾克·曼吉(英语:Jean-Jacq
  • 腓特烈王子 (安哈尔特-德绍)弗里德里希(Frederick;1769年12月27日-1814年5月27日),出生于德绍。安哈尔特王储。是安哈尔特-德绍亲王利奥波德三世与勃兰登堡-施韦特的路易丝的独生子。
  • 杨廷枢沧浪亭五百名贤像之杨廷枢石刻像杨廷枢(1595年-1647年),字维斗,号复庵,直隶苏州府吴县(今属江苏省苏州市)人。明朝末年学者,复社领袖之一。杨廷枢出身簪缨世家。早年为诸生,以气节自任
  • 威廉·明岑贝尔格威廉·明岑贝尔格(德语:Wilhelm "Willi" Münzenberg 1889年8月14日-1940年6月)德国共产党政治活动家,1919年至1920年担任青年共产国际首任主席,1921年建立国际工人救济会。魏玛
  • 牛津大学基督圣体学院基督圣体学院(Corpus Christi College)创校于1517年,按成立年份排名,是英国牛津大学第12所成员学院。基督圣体学院座落于墨顿街,位于墨顿学院和奥里尔学院之间,只有大约230名本科