逆威沙特分布

✍ dations ◷ 2025-07-06 18:55:43 #连续分布,多变量统计


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

相关

  • DNA 序列核酸序列(英语:Nucleic acid sequence,亦称为核酸的一级结构)使用一串字母表示的真实的或者假设的携带基因信息的DNA分子的一级结构。每个字母代表一种核碱基,两个碱基形成一个碱
  • 埃德温·萨瑟恩埃德温·麦勒·萨瑟恩爵士 (1938年6月7日-)是英国分子生物学家,英国皇家学会院士,爱丁堡皇家学会会员,拉斯克奖获得者,牛津大学生物化学名誉教授,牛津大学三一学院研究员。 他最广为
  • 行政院新闻局行政院新闻局(简称新闻局)为已裁撤之中华民国行政院附属机关,具有部会级地位,成立于1947年5月2日,后随着政府迁台而有多次的缩编与恢复。其主要负责行政院的公共关系、政策宣传、
  • 保罗·索伦蒂诺保罗·索伦蒂诺(意大利语:Paolo Sorrentino,1970年5月31日-)是一名意大利电影导演和编剧。2013年艺术作品《绝美之城》赢得欧洲电影奖最佳影片、导演、男演员和剪辑四个奖,以及第7
  • 贝尔县贝尔县(Bell County, Texas)位美国德克萨斯州中部偏东的一个县。面积2,818平方公里。根据美国2000年人口普查,共有人口237,974人。县治贝尔顿(Belton)。成立于1850年1月22日,县政
  • 迷失 (第四季)美国电视连续剧《迷失》的第四季在美国及加拿大地区于2008年1月31日首播,并于2008年5月29日结束。在故事时间线内,本季开始的超过90天前,一架飞机意外堕落在南太平洋一个岛屿上
  • 周在周在(?年-?年),字善卿,直隶太仓州人。明朝政治人物。正德九年(1514年)甲戌科进士。授宝坻知县,因忤中贵,逮捕,慷慨陈词,吏无法刁难。还职,县人薛凤鸣以御史附刘瑾,犯罪,周在绳之以法。刘瑾因
  • 白符站白符站(日语:白符駅/しらふえき  */?)是位于北海道松前郡福岛町字白符的北海道旅客铁道(JR北海道)松前线的铁路站。在1988年停用。车站名称源自阿伊努语的“cir-o-p”(鸟很多的地
  • 贡萨洛·桑切斯·德洛萨达贡萨洛·桑切斯·德洛萨达(西班牙语:Gonzalo Sánchez de Lozada y Sánchez de Bustamante,1930年7月1日-),生于玻利维亚拉巴斯市,政治人物,曾两度出任玻利维亚总统(1993-1997、2002
  • 刘业翔刘业翔(1930年9月1日-1930年9月1日),湖北武汉人,有色金属冶金专家,中国工程院院士,曾任中南工业大学(现中南大学)党委书记、书记兼校长。