威沙特分布

$n > 0\!$ 自由度 (实数)$\mathbf{V} > 0\,$ 尺度矩阵 (正定)

$\mathbf{W}\!$是正定的

$$\frac{\left|\mathbf{W}\right|^\frac{n-p-1}{2}} <pre> {2^\frac{np}{2}\left|{\mathbf V}\right|^\frac{n}{2}\Gamma_p(\frac{n}{2})} \exp\left(-\frac{1}{2}{\rm Tr}({\mathbf V}^{-1}\mathbf{W})\right)$$

$$n \mathbf{V}$$

$$(n-p-1)\mathbf{V}\text{ for }n \geq p+1$$

$$\Theta \mapsto \left|{\mathbf I} - 2i\,{\mathbf\Theta}{\mathbf V}\right|^{-n/2}$$

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

## 定义

假设X为一n × p矩阵，其各行（row）来自同一均值向量为$\mathbf 0$的$p$维多变量常态分布且彼此独立。

$$X_{(i)}{=}(x_i^1,\dots,x_i^p)^T\sim N_p(0,V),$$

则威沙特分布为$p\times p$散异矩阵

$$S=X^T X = \sum_{i = 1}^{n} X_{(i)} X_{(i)}^T, \,\!$$

几率分布。

$\mathbf S$有该几率分布通常记为

$$\mathbf S\sim W_p(\mathbf V,n).$$

其中正整数$n$为自由度。有时亦记号为$W(\mathbf V, p, n)$。若$p = 1$且$\mathbf V=1$则该分布退化为一自由度为$n$的单变量卡方分布。

## 常见应用

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

## 几率密度函数

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

令'$\mathbf W$为一$p\times p$正定对称随机变数矩阵。令$\mathbf V$为一特定正定$p\times p$矩阵。

如此，若$n > p$，则$\mathbf W$服从于一具自由度n的威沙特分布且有几率度函数$f_W$

$$f_{\mathbf W}(w)= \frac{ <pre> \left|w\right|^{(n-p-1)/2} \exp\left[ - {\rm trace}({\mathbf V}^{-1}w/2 )\right] </pre> }{ 2^{np/2}\left|{\mathbf V}\right|^{n/2}\Gamma_p(n/2) }$$

其中$\Gamma_p(\cdot)$为多变量Gamma分布，其定义为

$$\Gamma_p(n/2)= \pi^{p(p-1)/4}\Pi_{j=1}^p \Gamma\left[ (n+1-j)/2\right].$$

上述定义可推广至任一实数$n> p -1$

## 特征函数

威沙特分布的特征函数为

$$\Theta \mapsto \left|{\mathbf I} - 2i\,{\mathbf\Theta}{\mathbf V}\right|^{-n/2}.$$

也就是说

$$\Theta \mapsto {\mathcal E}\left\{\mathrm{exp}\left[i\cdot\mathrm{trace}({\mathbf W}{\mathbf\Theta})\right]\right\} = \left|{\mathbf I} - 2i{\mathbf\Theta}{\mathbf V}\right|^{-n/2}$$

其中${\mathcal E}(\cdot)$为期望值

（这里的$\Theta$及${\mathbf I}$ 皆为与${\mathbf V}$维度相同的矩阵。（${\mathbf I}$ 为单位矩阵，而$i$为－1的平方根）.

## 理论架构

若$\scriptstyle {\mathbf W}$为一自由度为m，共变异矩阵为$\scriptstyle {\mathbf V}$的威沙特分布，记为—$\scriptstyle {\mathbf W}\sim{\mathbf W}_p({\mathbf V},m)$—其中$\scriptstyle{\mathbf C}$为一$q\times p$的q秩矩阵，则

$${\mathbf C}{\mathbf W}{\mathbf C'} \sim {\mathbf W}_q\left({\mathbf C}{\mathbf V}{\mathbf C'},m\right).$$

### 推论1

若${\mathbf z}$为一非负$p\times 1$常数向量，则${\mathbf z'}{\mathbf W}{\mathbf z}\sim\sigma_z^2\chi_m^2$.

则在此情形下，$\chi_m^2$为一卡方分布且$\sigma_z^2={\mathbf z'}{\mathbf V}{\mathbf z}$（因${\mathbf V}$为正定，所以$\sigma_z^2$为一正常数）。

### 推论2

在${\mathbf z'}=(0,\ldots,0,1,0,\ldots,0)$ 的情形下（亦即第j个元素为1其他为0），推论1可导出

$$w_{jj}\sim\sigma_{jj}\chi^2_m$$

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

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

## 多变量常态分布的估计

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

## 分布抽样

以下的算法取材自 Smith & Hocking (1972)。一个来自自由度为n及共变异矩阵为$\mathbf V$的威沙特分布的$p\times p$（其中$n \geq p$）随机样本可以如下方式抽样而得：

- 生成一随机$p\times p$下三角矩阵 ${\textbf A}$使得：
- $a_{ii}=(\chi^2_{n-i+1})^{1/2}$，意即 $a_{ii}$为一$\chi^2_{n-i+1}$卡方分布随机样本的平方根。
- $a_{ij}$其中$j<i$，为一$N_1(0,1)$常态分布的随机样本。
- 计算${\textbf V} = {\textbf L}{\textbf L}^T$的Cholesky分解。
- 计算${\textbf X} = {\textbf L}{\textbf A}{\textbf A}^T{\textbf L}^T$。此时，${\textbf X}$ 为一$W_p({\textbf V},n)$的随机样本。

若${\textbf V}={\textbf I}$，则因${\textbf V}={\textbf I}{\textbf I}^T$，可以直接以${\textbf X} = {\textbf A}{\textbf A}^T$进行抽样。

## 参考条目

- 共变异矩阵的估计
- Hotelling的T平方分布
- 逆威沙特分布