{{Other uses|subject=變異數,又稱-{zh-cn: 变异数; zh-hk: 變異數; zh-tw:方差}-、變方|other=變異係數(coefficient of variation)|变异系数}}
SD=10),而蓝色总体的方差为2500(SD=50)。]]
變異數()又稱-{zh-cn: 变异数; zh-hk: 變異數; zh-tw:方差}-、變方,在概率论及统计学中,描述的是一个随机变量的离散程度,即一组数字与其平均值之间的距离的度量,是随机变量与其总体均值或样本均值的离差的平方的期望值。方差在统计中有非常核心的地位,其应用领域包括描述统计学、推論統計學、假說檢定、度量拟合优度,以及蒙特卡洛采样。由于科学分析经常涉及统计,方差也是重要的科研工具。方差是標準差的平方、分布的二阶矩,以及随机变量与其自身的协方差,其常用的符号表示有\sigma^2、s^2、\operatorname{Var}(X)、V(X),以及\mathbb{V}(X)。
方差作为离散度量的优点是,它比其他离散度量(如平均差)更易于代数运算;例如,一组不相关的随机变量和的方差等于它们方差的和。在实际应用中,方差的一个缺点是它与随机变量的单位不同,而標準差则单位相同,这就是计算完成后通常采用标准差来衡量离散程度的原因。
有两个不同的概念都被称为“方差”。一种如上所述,是理论概率分布的方差。而另一种方差是一组观测值的特征。观测值通常是从真实世界的系统中测量的。如果给出系统的所有可能的观测,则它们算出的方差称为总体方差;然而,一般情况下我们只使用总体的一个子集(样本),由此计算出的方差称为样本方差。用样本计算出的方差可认为是对整个总体的方差的估计量。
方差的正平方根称为该随机变量的標準差;方差除以期望值归一化的值叫分散指数;标准差除以平均值归一化的值叫变异系数。
定义
设为服从分布的随机变量,如果是随机变量的期望值(均值),则随机变量或者分布的方差为的离差平方的期望值:
:\operatorname{Var}(X) = \operatorname{E}\left[(X - \mu)^2 \right]
这个定义涵盖了连续、离散,或两者皆非的随机变量。方差亦可視作随机变量与自身的协方差:
:\operatorname{Var}(X) = \operatorname{Cov}(X, X)
方差也等价于生成的概率分布的二阶累积量。方差的常用的表达有\operatorname{Var}(X),有时作V(X)或\mathbb{V}(X),也可写作符号\sigma^2_X或\sigma^2(读作“sigma方”)。方差的表达式可展开如下:
:\begin{align}
\operatorname{Var}(X) &= \operatorname{E}\left[(X - \operatorname{E}[X])^2\right] \\[4pt]
&= \operatorname{E}\left[X^2 - 2X\operatorname{E}[X] + \operatorname{E}[X]^2\right] \\[4pt]
&= \operatorname{E}\left[X^2\right] - 2\operatorname{E}[X]\operatorname{E}[X] + \operatorname{E}[X]^2 \\[4pt]
&= \operatorname{E}\left[X^2 \right] - \operatorname{E}[X]^2
\end{align}
也就是说,的方差等于平方的均值减去均值的平方。该等式不应该用于浮点运算,因为如果等式的两个成分大小相似,将会造成灾难性抵消。
离散随机变量
如果随机变量是具有概率质量函數的离散随机分布,則:
:\operatorname{Var}(X) = \sum_{i=1}^n p_i\cdot(x_i - \mu)^2 = \sum_{i=1}^n (p_i\cdot x_i^2) - \mu^2
此處\mu是其期望值,即:
:\mu = \sum_{i=1}^n p_i\cdot x_i.
x_i 表示實現值(realized value)
當為有個相等機率值的離散型均勻分佈時:
:\mu = \frac{1}{n}\sum_{i=1}^n x_i,
:\operatorname{Var}(X) = \sigma^{2} =\frac{1}{n} \sum_{i=1}^n (x_i - \mu)^2 = \frac{1}{n} \left(\sum_{i=1}^n x_i^2 - n\mu^2 \right) = \frac{ \sum_{ i = 1 }^n x_i^2 }{ n } - \mu^2.
個相等機率值的方差亦可以點對點間的方變量表示為:
: \operatorname{Var}(X) = \frac{1}{n^2} \sum_{i=1}^n \sum_{j=1}^n \frac{1}{2}(x_i - x_j)^2.
连续型随机变量
如果随机变量是連續分布,機率密度函數為,相應的累积分布函数為,則其方差為:
:\begin{align}
\operatorname{Var}(X) = \sigma^2 &= \int_{\R} (x-\mu)^2 f(x) \, dx \\[4pt]
&= \int_{\R} x^2f(x)\,dx -2\mu\int_{\R} xf(x)\,dx + \mu^2\int_{\R} f(x)\,dx \\[4pt]
&= \int_{\R} x^2 \,dF(x) - 2 \mu \int_{\R} x \,dF(x) + \mu^2 \int_{\R} \,dF(x) \\[4pt]
&= \int_{\R} x^2 \,dF(x) - 2 \mu \cdot \mu + \mu^2 \cdot 1 \\[4pt]
&= \int_{\R} x^2 \,dF(x) - \mu^2,
\end{align}
或等價地:
:\operatorname{Var}(X) = \int_{\R} x^2 f(x) \,dx - \mu^2 ,
其中\mu為X的期望值,其計算方法如下:
:\mu = \int_{\R} x f(x) \, dx = \int_{\R} x \, d F(x).
這些公式中,dx和dF(x)的積分分別為勒貝格積分和。
若函數x^2f(x)在每個有限區間[a,b]\subset\R都是黎曼可積的,則:
:\operatorname{Var}(X) = \int^{+\infty}_{-\infty} x^2 f(x) \, dx - \mu^2,
該積分為非正常黎曼积分。
常见概率分布
下表列出了一些常用概率分布的方差。
特性
方差不會是負的,因為平方運算結果為非負數:
:\operatorname{Var}(X)\ge 0
一個常數隨機變數的方差為零。反之,若有限個數組成的資料集方差為零,則其內所有數皆相等。對於一般隨機變數,也有類似結論,即方差為零推出該變數幾乎總是取同一個值:
:P(X=a) = 1\Leftrightarrow \operatorname{Var}(X)= 0
方差不變於定位參數的變動。也就是說,如果一個常數被加至一個數列中的所有變數值,此數列的方差不會改變:
:\operatorname{Var}(X+a)=\operatorname{Var}(X).
如果所有數值被放大一個常數倍,方差會放大此常數的平方倍:
:\operatorname{Var}(aX)=a^2\operatorname{Var}(X)
兩個隨機變數合的方差為:
:\operatorname{Var}(aX+bY)=a^2\operatorname{Var}(X)+b^2\operatorname{Var}(Y)+2ab\, \operatorname{Cov}(X,Y),
:\operatorname{Var}(X-Y)=\operatorname{Var}(X)+\operatorname{Var}(Y)-2\, \operatorname{Cov}(X,Y),
此處代表共變異數。
對於N個隨機變數\{X_1,\dots,X_N\}的總和:
:\operatorname{Var}\left(\sum_{i=1}^N X_i\right)=\sum_{i,j=1}^N\operatorname{Cov}(X_i,X_j)=\sum_{i=1}^N\operatorname{Var}(X_i)+\sum_{i\ne j}\operatorname{Cov}(X_i,X_j)
在样本空间Ω上存在有限期望和方差的随机变量构成一个希尔伯特空间:,不过这裡的内积和长度跟协方差,标准差还是不大一样。所以,我们得把这个空间“除”常变量构成的子空间,也就是说把相差一个常数的所有原来那个空间的随机变量做成一个等价类。这还是一个新的无穷维线性空间,并且有一个从旧空间内积诱导出来的新内积,而这个内积就是协方差。
总体方差和样本方差
总体方差
一般而言,一个有限的容量为、元素的值为的总体的总体方差为:
: \begin{align}
\sigma^2 &= \frac{1}{N} \sum_{i=1}^N \left(x_i - \mu\right)^2 = \frac{1}{N} \sum_{i=1}^N \left(x_i^2 - 2\mu x_i + \mu^2 \right) \\[5pt]
&= \left(\frac 1N \sum_{i=1}^N x_i^2\right) - 2\mu \left(\frac{1}{N} \sum_{i=1}^N x_i\right) + \mu^2 \\[5pt]
&= \left(\frac{1}{N} \sum_{i=1}^N x_i^2\right) - \mu^2
\end{align}
其中总体均值为:
: \mu = \frac 1N \sum_{i=1}^N x_i.
总体方差也可用下式计算:
:\sigma^2 = \frac {1} {N^2}\sum_{i
该式成立,是因为:
: \begin{align}
&\frac{1}{2N^2} \sum_{i, j=1}^N\left( x_i - x_j \right)^2 \\[5pt]
={} &\frac{1}{2N^2} \sum_{i, j=1}^N\left( x_i^2 - 2x_i x_j + x_j^2 \right) \\[5pt]
={} &\frac{1}{2N} \sum_{j=1}^N\left(\frac{1}{N} \sum_{i=1}^N x_i^2\right) - \left(\frac{1}{N} \sum_{i=1}^N x_i\right)\left(\frac{1}{N} \sum_{j=1}^N x_j\right) + \frac{1}{2N} \sum_{i=1}^N\left(\frac{1}{N} \sum_{j=1}^N x_j^2\right) \\[5pt]
={} &\frac{1}{2} \left( \sigma^2 + \mu^2 \right) - \mu^2 + \frac{1}{2} \left( \sigma^2 + \mu^2 \right) \\[5pt]
={} &\sigma^2
\end{align}
总体方差与生成该总体的概率分布的方差相匹配。因此,“总体”的概念可推广到具有无限总体的连续随机变量。
样本方差
有偏样本方差
在许多实际情况下,总体的真实方差无法事先知道,必须以某种方式计算出来。在面对非常大的总体时,不可能计算总体中的每一个元素,因此必须从总体中抽取样本进行计算。样本方差还可以应用于用连续分布的样本来估计该分布的方差。
下面我们从总体中有放回抽取个数值,其中,并用该样本来估计总体的方差。直接使用样本数据的方差,得到的是的均值:
:\sigma_Y^2 =
\frac{1}{n} \sum_{i=1}^n \left(Y_i - \overline{Y}\right)^2 = \left(\frac 1n \sum_{i=1}^n Y_i^2\right) - \overline{Y}^2 =
\frac{1}{n^2} \sum_{i,j\,:\,i
此处,\overline{Y}表示:
:\overline{Y} = \frac{1}{n} \sum_{i=1}^n Y_i .
由于是随机选取的,\overline{Y}和\sigma_Y^2都是随机变量。它们的期望值可以用从总体中抽取的所有可能的容量为的{{Math|{Yi}}}的样本集合来估计。对于\sigma_Y^2即为:
:\begin{align}
\operatorname{E}[\sigma_Y^2]
&= \operatorname{E}\left[ \frac{1}{n} \sum_{i=1}^n \left(Y_i - \frac{1}{n} \sum_{j=1}^n Y_j \right)^2 \right] \\[5pt]
&= \frac 1n \sum_{i=1}^n \operatorname{E}\left[ Y_i^2 - \frac{2}{n} Y_i \sum_{j=1}^n Y_j + \frac{1}{n^2} \sum_{j=1}^n Y_j \sum_{k=1}^n Y_k \right] \\[5pt]
&= \frac 1n \sum_{i=1}^n \left( \frac{n - 2}{n} \operatorname{E}\left[Y_i^2\right] - \frac{2}{n} \sum_{j \neq i} \operatorname{E}\left[Y_i Y_j\right] + \frac{1}{n^2} \sum_{j=1}^n \sum_{k \neq j}^n \operatorname{E}\left[Y_j Y_k\right] +\frac{1}{n^2} \sum_{j=1}^n \operatorname{E}\left[Y_j^2\right] \right) \\[5pt]
&= \frac 1n \sum_{i=1}^n \left[ \frac{n - 2}{n} \left(\sigma^2 + \mu^2\right) - \frac{2}{n} (n - 1)\mu^2 + \frac{1}{n^2} n(n - 1)\mu^2 + \frac{1}{n} \left(\sigma^2 + \mu^2\right) \right] \\[5pt]
&= \frac{n - 1}{n} \sigma^2.
\end{align}
因此,\sigma_Y^2给出的是总体方差的有偏估计量,偏差为\frac{n - 1}{n}。因此,\sigma_Y^2称为有偏样本方差。
无偏样本方差
将偏差纠正后,可得到无偏样本方差,记为s^2:
:s^2 = \frac{n}{n - 1} \sigma_Y^2 = \frac{n}{n - 1} \left[ \frac{1}{n} \sum_{i=1}^n \left(Y_i - \overline{Y}\right)^2 \right] = \frac{1}{n - 1} \sum_{i=1}^n \left(Y_i - \overline{Y} \right)^2
当语境明确时,两个估计量都可以简称为“样本方差”。同样的证明也适用于取自连续概率分布的样本。
其中,对的使用称为,它也用于和样本标准差(方差的平方根)。平方根是一个凹函数,因此会引入负偏差(根据簡森不等式),具体取决于分布,因此校正的样本标准差(使用贝塞尔校正)是有偏的。是一个技术上复杂的问题,不过对于正态分布,使用能得到几乎无偏的估计值。
无偏样本方差是函数的U-统计量。
一般化
如果是一个向量其取值范围在實數空间,并且其每个元素都是一个一维随机变量,我们就把称为随机向量。随机向量的方差是一维随机变量方差的自然推广,其定义为,其中,是的转置。这个方差是一个非负定的方阵,通常称为协方差矩阵。
如果是一个複數随机变量的向量(向量中每個元素均為複數的隨機變數),那么其方差定义则为,其中是的共轭转置向量或稱為埃尔米特向量。根据这个定义,變異數为实数。
历史
「方差」(variance)这个名词率先由羅納德·費雪()在论文《》中提出。
后来方差逐渐衍生出了「半方差」(semivariance)、「亚方差」(hypo variance)、「超方差」(super variance)、「」(circular variance)与「倒方差」(inverse variance)等概念。
半方差
半方差的計算方式與方差類似,但是只包括了低於均值的觀測值:
: \text{Semivariance} = {1\over{n}}\sum_{i:x_{i}
半方差在不同应用领域也被用作特殊的量度。对于偏态分布,半方差能提供方差所不能提供的额外信息。
参见
- 方差分析
- 标准差
- 标准离差率
- 变异系数
- 异方差
- 最小平方頻譜分析法
- 离散程度
*
方差类型
- 相关 (概率论)
- 可解释变异
- 合并方差
*
参考文献
评论 (0)