方差

描述一個隨機變數的離散程度

方差(英語:variance)又稱變異數[1]變方[2],在概率論統計學中,描述的是一個隨機變量的離散程度,即一組數字與其平均值之間的距離的度量,是隨機變量與其總體均值或樣本均值的離差的平方的期望值。方差在統計中有非常核心的地位,其應用領域包括敘述統計學推論統計學假設檢定、度量適合度,以及蒙地卡羅採樣。由於科學分析經常涉及統計,方差也是重要的科研工具。方差是標準差的平方、分佈的二階動差,以及隨機變量與其自身的協方差,其常用的符號表示有,以及[3]

「Variance」的各地常用名稱
中國大陸方差
臺灣變異數
港澳方差
日本、韓國分散
越南分散(phương sai
兩個來自不同總體的樣本的案例,均值相同但方差不同。兩者總體均值都是100,紅色總體的方差為100(標準差SD=10),而藍色總體的方差為2500(SD=50)。

方差作為離散度量的優點是,它比其他離散度量(如平均差)更易於代數運算;例如,一組不相關的隨機變量和的方差等於它們方差的和。在實際應用中,方差的一個缺點是它與隨機變量的單位不同,而標準差則單位相同,這就是計算完成後通常採用標準差來衡量離散程度的原因。

有兩個不同的概念都被稱為「方差」。一種如上所述,是理論概率分佈的方差。而另一種方差是一組觀測值的特徵。觀測值通常是從真實世界的系統中測量的。如果給出系統的所有可能的觀測,則它們算出的方差稱為總體方差;然而,一般情況下我們只使用總體的一個子集(樣本),由此計算出的方差稱為樣本方差。用樣本計算出的方差可認為是對整個總體的方差的估計量。

方差的正平方根稱為該隨機變量的標準差;方差除以期望值歸一化的值叫分散指數;標準差除以平均值歸一化的值叫變異系數

定義

編輯

X為服從分佈F隨機變量,如果E[X]是隨機變量X期望值(均值μ=E[X]),則隨機變量X或者分佈F的方差為X離差平方的期望值:

 

這個定義涵蓋了連續、離散,或兩者皆非的隨機變量。方差亦可視作隨機變量與自身的協方差

 

方差也等價於生成X的概率分佈的二階累積量。方差的常用的表達有 ,有時作  ,也可寫作符號  (讀作「sigma方」)。方差的表達式可展開如下:

 

也就是說,X的方差等於X平方的均值減去X均值的平方。該等式不應該用於浮點運算,因為如果等式的兩個成分大小相似,將會造成災難性抵消

離散隨機變量

編輯

如果隨機變量X是具有概率質量函數的離散隨機分佈x1 ↦ p1, ..., xn ↦ pn,則:

 

此處 是其期望值,即:

 

  表示實現值(realized value)

X為有n個相等概率值的離散型均勻分佈時:

 
 

n個相等概率值的方差亦可以點對點間的方變量表示為:

 

連續型隨機變量

編輯

如果隨機變量X是連續分佈,概率密度函數f(x),相應的累積分佈函數F(x),則其方差為:

 

或等價地:

 

其中  的期望值,其計算方法如下:

 

這些公式中,  的積分分別為勒貝格積分勒貝格-斯蒂爾吉斯積分英語Lebesgue–Stieltjes integration

若函數 在每個有限區間 都是黎曼可積的,則:

 

該積分為非正常黎曼積分

常見概率分佈

編輯

下表列出了一些常用概率分佈的方差。

概率分佈類型 概率分佈函數 均值 方差
二項式分佈      
幾何分佈      
正態分佈      
連續型均勻分佈      
指數分佈      
泊松分佈      

特性

編輯

方差不會是負的,因為平方運算結果為非負數:

 

一個常數隨機變量的方差為零。反之,若有限個數組成的資料集方差為零,則其內所有數皆相等。對於一般隨機變量,也有類似結論,即方差為零推出該變量幾乎總是取同一個值:

 

方差不變於定位參數的變動。也就是說,如果一個常數被加至一個數列中的所有變量值,此數列的方差不會改變:

 

如果所有數值被放大一個常數倍,方差會放大此常數的平方倍:

 

兩個隨機變量合的方差為:

 
 

此處Cov(X, Y)代表協方差

對於 個隨機變量 的總和:

 

在樣本空間Ω上存在有限期望值和方差的隨機變量構成一個希爾伯特空間L2 (Ω, dP),不過這裏的內積和長度跟協方差,標準差還是不大一樣。所以,我們得把這個空間「除」常變量構成的子空間,也就是說把相差一個常數的所有原來那個空間的隨機變量做成一個等價類。這還是一個新的無窮維線性空間,並且有一個從舊空間內積誘導出來的新內積,而這個內積就是協方差。

總體方差和樣本方差

編輯

總體方差

編輯

一般而言,一個有限的容量為N、元素的值為xi總體總體方差為:

 

其中總體均值為:

 

總體方差也可用下式計算:

 

該式成立,是因為:

 

總體方差與生成該總體的概率分佈的方差相匹配。因此,「總體」的概念可推廣到具有無限總體的連續隨機變量。

樣本方差

編輯

有偏樣本方差

編輯

在許多實際情況下,總體的真實方差無法事先知道,必須以某種方式計算出來。在面對非常大的總體時,不可能計算總體中的每一個元素,因此必須從總體中抽取樣本進行計算。[4]樣本方差還可以應用於用連續分佈的樣本來估計該分佈的方差。

下面我們從總體中有放回抽取n個數值Y1, ..., Yn,其中n < N,並用該樣本來估計總體的方差。[5]直接使用樣本數據的方差,得到的是離差平方英語squared deviations的均值:

 

此處, 表示樣本均值

 

由於Yi是隨機選取的,  都是隨機變量。它們的期望值可以用從總體中抽取的所有可能的容量為n{Yi}的樣本集合來估計。對於 即為:

 

因此, 給出的是總體方差的有偏估計量,偏差為 。因此, 稱為有偏樣本方差

無偏樣本方差

編輯

將偏差糾正後,可得到無偏樣本方差,記為 

 

當語境明確時,兩個估計量都可以簡稱為「樣本方差」。同樣的證明也適用於取自連續概率分佈的樣本。

其中,對n − 1的使用稱為貝塞爾校正英語Bessel's correction,它也用於樣本協方差英語sample covariance樣本標準差(方差的平方根)。平方根是一個凹函數,因此會引入負偏差(根據簡森不等式),具體取決於分佈,因此校正的樣本標準差(使用貝塞爾校正)是有偏的。標準差的無偏估計英語unbiased estimation of standard deviation是一個技術上複雜的問題,不過對於正態分佈,使用n - 1.5能得到幾乎無偏的估計值。

無偏樣本方差是函數ƒ(y1y2) = (y1 − y2)2/2U-統計量

一般化

編輯

如果X是一個向量其取值範圍在實數空間Rn,並且其每個元素都是一個一維隨機變量,我們就把X稱為隨機向量。隨機向量的方差是一維隨機變量方差的自然推廣,其定義為E[(Xμ)(Xμ)T],其中μ = E(X)XTX的轉置。這個方差是一個非負定方陣,通常稱為協方差矩陣

如果X是一個複數隨機變量的向量(向量中每個元素均為複數的隨機變量),那麼其方差定義則為E[(X − μ)(X − μ)*],其中X*X共軛轉置向量或稱為埃爾米特向量。根據這個定義,方差為實數。

歷史

編輯

方差」(variance)這個名詞率先由羅納德·費雪(英語:Ronald Fisher)在論文《The Correlation between Relatives on the Supposition of Mendelian Inheritance[6]中提出。

後來方差逐漸衍生出了「半方差」(semivariance)、「亞方差」(hypo variance)、「超方差」(super variance)、「圓方差英語circular variance」(circular variance)與「倒方差」(inverse variance)等概念。

半方差

編輯

半方差的計算方式與方差類似,但是只包括了低於均值的觀測值:

 

半方差在不同應用領域也被用作特殊的量度。對於偏態分佈,半方差能提供方差所不能提供的額外資訊。[7]

參見

編輯

方差類型

編輯

參考文獻

編輯
  1. ^ 存档副本. [2023-07-25]. (原始內容存檔於2023-07-25). 
  2. ^ 存档副本. [2023-07-25]. (原始內容存檔於2023-07-25). 
  3. ^ Wasserman, Larry. All of Statistics: a concise course in statistical inference. Springer texts in statistics. 2005: 51. ISBN 9781441923226. 
  4. ^ Navidi, William (2006) Statistics for Engineers and Scientists, McGraw-Hill, pg 14.
  5. ^ Montgomery, D. C. and Runger, G. C. (1994) Applied statistics and probability for engineers, page 201. John Wiley & Sons New York
  6. ^ Ronald Fisher(1918)The correlation between relatives on the supposition of Mendelian Inheritance頁面存檔備份,存於互聯網檔案館
  7. ^ Fama, Eugene F.; French, Kenneth R. Q&A: Semi-Variance: A Better Risk Measure?. Fama/French Forum. 2010-04-21 [2022-06-10]. (原始內容存檔於2021-07-25).