威沙特分布

统计学家约翰·威沙特为名的威沙特分布统计学上的一种半正定矩阵随机分布。[1]这个分布在多变量分析协方差矩阵估计上相当重要。

威沙特
参数 自由度 (实数)
尺度矩阵 (正定)
值域 是正定的
概率密度函数
期望值
众数
特征函数


定义

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

 

则威沙特分布为 散异矩阵英语Scatter matrix

 

机率分布

 有该机率分布通常记为

 

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

常见应用

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

机率密度函数

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

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

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

 

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

 

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

特征函数

威沙特分布的特征函数

 

也就是说

 

其中 为期望值

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

理论架构

 为一自由度为m,共变异矩阵为 的威沙特分布,记为— —其中 为一 q秩矩阵,则[4]

 

推论1

 为一非负 常数向量,则[4]  .

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

推论2

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

 

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

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

多变量常态分布的估计

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

分布抽样

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

  1. 生成一随机 三角矩阵  使得:
    •  ,意即  为一 卡方分布随机样本的平方根。
    •  其中 ,为一 常态分布的随机样本。[8]
  2. 计算 Cholesky分解
  3. 计算 。此时,  为一 的随机样本。

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

参考条目

参考资料

  1. ^ Wishart, J. The generalised product moment distribution in samples from a normal multivariate population. Biometrika. 1928, 20A (1–2): 32–52. JFM 54.0565.02. JSTOR 2331939. doi:10.1093/biomet/20A.1-2.32. 
  2. ^ Uhlig, H. On Singular Wishart and Singular Multivariate Beta Distributions. The Annals of Statistics. 1994, 22: 395–405. doi:10.1214/aos/1176325375. 
  3. ^ Anderson, T. W. An Introduction to Multivariate Statistical Analysis 3rd. Hoboken, N. J.: Wiley Interscience. 2003: 259. ISBN 0-471-36091-0. 
  4. ^ 4.0 4.1 Rao, C. R. Linear Statistical Inference and its Applications. Wiley. 1965: 535. 
  5. ^ Seber, George A. F. Multivariate Observations. Wiley. 2004. ISBN 978-0471691211. 
  6. ^ Chatfield, C.; Collins, A. J. Introduction to Multivariate Analysis. London: Chapman and Hall. 1980: 103–108. ISBN 0-412-16030-7. 
  7. ^ Smith, W. B.; Hocking, R. R. Algorithm AS 53: Wishart Variate Generator. Journal of the Royal Statistical Society, Series C. 1972, 21 (3): 341–345. JSTOR 2346290. 
  8. ^ Anderson, T. W. An Introduction to Multivariate Statistical Analysis 3rd. Hoboken, N. J.: Wiley Interscience. 2003: 257. ISBN 0-471-36091-0.