版权归原作者所有,如有侵权,请联系我们

[科普中国]-威沙特分布

科学百科
原创
科学百科为用户提供权威科普内容,打造知识科普阵地
收藏

以统计学家约翰·威沙特为名的威沙特分布是统计学上的一种半正定矩阵随机分布。这个分布在多变量分析的共变异矩阵估计上相当重要。

定义假设X为一n×p矩阵,其各行(row)来自同一均值向量为 维多变量正态分布且彼此独立。

则威沙特分布为 散异矩阵

的机率分布。

有该机率分布通常记

其中正整数 为自由度。有时亦记号为 。若 则该分布退化为一自由度为 的单变量卡方分布。

常见应用威沙特分布常用于多变量的概似比检定,亦用于随机矩阵的频谱理论中。

机率密度函数威沙特分布具有下述的机率密度函数:

令' 为一 正定对称随机变数矩阵。令 为一特定正定 矩阵。

如此,若 ,则 服从于一具自由度n的威沙特分布且有机率度函数

其中 为多变量Gamma分布,其定义为

上述定义可推广至任一实数

特征函数威沙特分布的特征函数为

也就是说

其中 为期望值

(这里的 皆为与 维度相同的矩阵。 为单位矩阵,而i为-1的平方根)

理论架构若 为一自由度为m,共变异矩阵为 的威沙特分布,记为— —其中 为一 的q秩矩阵,则

推论1若 为一非负 常数向量,则

则在此情形下, 为一卡方分布且 (因 为正定,所以 为一正常数)。

推论2在 的情形下(亦即第j个元素为1其他为0),推论1可导出

为矩阵的每一个对对角元素的边际分布。

统计学家George Seber曾论证威沙特分布并非多变量卡方分布,这是因为非对角元素的边际分布并非卡方分布,Seber倾向于将某某多变量分布此一遣词用于所有元素的边际分布皆相同的情形1。

多变量正态分布的估计由于威沙特分布可视为一多变量正态分布其共变异矩阵的最大概似估计量(MLE)的的分布,其衍自MLE的计算可为令人惊喜地简约而优雅。基于频谱理论,可将一标量视为一 矩阵的迹(trace)。请参考共变异矩阵的估计。

分布抽样以下的算法取材自 Smith & Hocking (1972)2。一个来自自由度为n及共变异矩阵为 的威沙特分布的 (其中 )随机样本可以如下方式抽样而得:

生成一随机 下三角矩阵 使得:

,意即 为一 卡方分布随机样本的平方根。

其中 ,为一 正态分布的随机样本。

计算 的Cholesky分解。

计算 。此时, 为一 的随机样本。

,则因 ,可以直接以 进行抽样。

参考条目共变异矩阵的估计

Hotelling的T平方分布

逆威沙特分布

本词条内容贡献者为:

李嘉骞 - 博士 - 同济大学