値が平均からどれだけ離れているかを示す統計的尺度
平均は同じだが分散が異なる 2 つの母集団からのサンプルの例。赤い母集団の平均は 100、分散は 100 (SD=10) で、青い母集団の平均は 100、分散は 2500 (SD=50) です。ここで、SD は標準偏差を表します。
確率論 と 統計学 において 、 分散とは ランダム変数 の 平均からの偏差 の二乗の 期待値 です 。 標準偏差 (SD) は分散の平方根として得られます。分散は 散布度 の尺度であり、つまり、一連の数値が平均値から どれだけ離れているかの尺度です。これは 分布 の2番目の 中心モーメント であり、ランダム変数とそれ自身の 共分散 であり、、、、、 または で表されます 。 [1 ]
σ
2
{\displaystyle \sigma ^{2}}
s
2
{\displaystyle s^{2}}
ヴァール
(
バツ
)
{\displaystyle \operatorname {Var} (X)}
五
(
バツ
)
{\displaystyle V(X)}
五
(
バツ
)
{\displaystyle \mathbb {V} (X)}
分散を分散の尺度として用いる利点は、期待絶対偏差 などの他の分散の尺度よりも代数的操作がしやすいことです 。たとえば、相関のないランダム変数の合計の分散は、それらの分散の合計に等しくなります。実際の用途における分散の欠点は、標準偏差とは異なり、ランダム変数とは単位が異なることです。そのため、計算が完了すると、分散の尺度として標準偏差が報告されるのが一般的です。もう 1 つの欠点は、分散が多くの分布に対して有限ではないことです。
どちらも「分散」と呼ばれる 2 つの異なる概念があります。1 つは、上で説明したように、理論的な 確率分布 の一部であり、方程式によって定義されます。もう 1 つの分散は、一連の観測の特性です。分散が観測から計算される場合、それらの観測は通常、現実世界のシステムから測定されます。システムのすべての可能な観測が存在する場合、計算された分散は母集団分散と呼ばれます。ただし、通常はサブセットのみが利用可能であり、これから計算された分散はサンプル分散と呼ばれます。サンプルから計算された分散は、完全な母集団分散の推定値と見なされます。母集団分散の推定値を計算する方法は複数あり、これについては以下のセクションで説明します。
2 種類の分散は密接に関連しています。どのように関連しているかを確認するには、理論上の確率分布を仮説的観測の生成元として使用できることを考えてみましょう。分布を使用して無限の数の観測が生成されると、その無限セットから計算されたサンプル分散は、分布の分散方程式を使用して計算された値と一致します。分散は統計において中心的な役割を果たしており、分散を使用するアイデアには、 記述統計 、 統計的推論 、 仮説検定 、 適合度 、 モンテカルロ サンプリング などがあります。
任意の分布(2、4、4、4、5、5、7、9)の分散の幾何学的視覚化: 頻度分布が構築されます。 分布の重心はその平均を与えます。 各値に対して、平均値からの差に等しい辺を持つ正方形が形成されます。 正方形を長方形に配置し、その一辺を値の数 n に等しくすると、もう一辺は分布の分散 σ 2 になります。
意味
ランダム変数の分散は、 平均
からの偏差の二乗 の 期待 値 です 。
バツ
{\displaystyle X}
バツ
{\displaystyle X}
μ
=
え
[
バツ
]
{\displaystyle \mu =\operatorname {E} [X]}
ヴァール
(
バツ
)
=
え
[
(
バツ
−
μ
)
2
]
。
{\displaystyle \operatorname {Var} (X)=\operatorname {E} \left[(X-\mu )^{2}\right].}
この定義には、離散的 、 連続的 、 どちらでもない 、または混合のプロセスによって生成されるランダム変数が含まれます。分散は 、ランダム変数とそれ自身の
共分散 と考えることもできます。
ヴァール
(
バツ
)
=
カバー
(
バツ
、
バツ
)
。
{\displaystyle \operatorname {Var} (X)=\operatorname {Cov} (X,X).}
分散は、 を生成する確率分布の2 番目の キュムラント と同等です 。分散は 、通常 、または 、あるいは 記号 的に 、または 単に (「 シグマ スクエアード」と発音)と表されます。分散の式は、次のように展開できます。
バツ
{\displaystyle X}
ヴァール
(
バツ
)
{\displaystyle \operatorname {Var} (X)}
五
(
バツ
)
{\displaystyle V(X)}
五
(
バツ
)
{\displaystyle \mathbb {V} (X)}
σ
バツ
2
{\displaystyle \sigma _{X}^{2}}
σ
2
{\displaystyle \sigma ^{2}}
ヴァール
(
バツ
)
=
え
[
(
バツ
−
え
[
バツ
]
)
2
]
=
え
[
バツ
2
−
2
バツ
え
[
バツ
]
+
え
[
バツ
]
2
]
=
え
[
バツ
2
]
−
2
え
[
バツ
]
え
[
バツ
]
+
え
[
バツ
]
2
=
え
[
バツ
2
]
−
2
え
[
バツ
]
2
+
え
[
バツ
]
2
=
え
[
バツ
2
]
−
え
[
バツ
]
2
{\displaystyle {\begin{aligned}\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]-2\operatorname {E} [X]^{2}+\operatorname {E} [X]^{2}\\[4pt]&=\operatorname {E} \left[X^{2}\right]-\operatorname {E} [X]^{2}\end{aligned}}}
言い換えると、 Xの分散は、 X の二乗の平均から X の平均の二乗を引いたものに等しくなります。この式は、 浮動小数点演算 を使用した計算には使用しないでください。式の 2 つの要素の大きさが同じである場合、 壊滅的な相殺 が発生するからです 。数値的に安定した他の代替方法については、「 分散を計算するためのアルゴリズム」 を参照してください。
離散確率変数
確率変数の生成器が 確率質量関数 を持つ 離散 的で ある場合、
バツ
{\displaystyle X}
x
1
↦
p
1
、
x
2
↦
p
2
、
…
、
x
ん
↦
p
ん
{\displaystyle x_{1}\mapsto p_{1},x_{2}\mapsto p_{2},\ldots ,x_{n}\mapsto p_{n}}
ヴァール
(
バツ
)
=
∑
私
=
1
ん
p
私
⋅
(
x
私
−
μ
)
2
、
{\displaystyle \operatorname {Var} (X)=\sum _{i=1}^{n}p_{i}\cdot (x_{i}-\mu )^{2},}
ここで 期待値です。つまり、
μ
{\displaystyle \mu}
μ
=
∑
私
=
1
ん
p
私
x
私
。
{\displaystyle \mu =\sum _{i=1}^{n}p_{i}x_{i}.}
(このような離散 加重分散が、 合計が 1 でない重みによって指定される場合、重みの合計で割ります。)
等確率の値の集合の分散 は次のように表される。
ん
{\displaystyle n}
ヴァール
(
バツ
)
=
1
ん
∑
私
=
1
ん
(
x
私
−
μ
)
2
{\displaystyle \operatorname {Var} (X)={\frac {1}{n}}\sum _{i=1}^{n}(x_{i}-\mu )^{2}}
ここで 平均値です。つまり、
μ
{\displaystyle \mu}
μ
=
1
ん
∑
私
=
1
ん
x
私
。
{\displaystyle \mu ={\frac {1}{n}}\sum _{i=1}^{n}x_{i}.}
等確率の値の集合の分散は 、平均値を直接参照することなく、点同士の距離の二乗偏差の二乗で表現することができる。 [2]
ん
{\displaystyle n}
ヴァール
(
バツ
)
=
1
ん
2
∑
私
=
1
ん
∑
じ
=
1
ん
1
2
(
x
私
−
x
じ
)
2
=
1
ん
2
∑
私
∑
じ
>
私
(
x
私
−
x
じ
)
2
。
{\displaystyle \operatorname {Var} (X)={\frac {1}{n^{2}}}\sum _{i=1}^{n}\sum _{j=1}^{n}{\frac {1}{2}}(x_{i}-x_{j})^{2}={\frac {1}{n^{2}}}\sum _{i}\sum _{j>i}(x_{i}-x_{j})^{2}.}
絶対連続ランダム変数
確率変数が 確率密度関数 を持ち 、 対応する 累積分布関数 である場合、
バツ
{\displaystyle X}
ふ
(
x
)
{\displaystyle f(x)}
ふ
(
x
)
{\displaystyle F(x)}
ヴァール
(
バツ
)
=
σ
2
=
∫
R
(
x
−
μ
)
2
ふ
(
x
)
d
x
=
∫
R
x
2
ふ
(
x
)
d
x
−
2
μ
∫
R
x
ふ
(
x
)
d
x
+
μ
2
∫
R
ふ
(
x
)
d
x
=
∫
R
x
2
d
ふ
(
x
)
−
2
μ
∫
R
x
d
ふ
(
x
)
+
μ
2
∫
R
d
ふ
(
x
)
=
∫
R
x
2
d
ふ
(
x
)
−
2
μ
⋅
μ
+
μ
2
⋅
1
=
∫
R
x
2
d
ふ
(
x
)
−
μ
2
、
{\displaystyle {\begin{aligned}\operatorname {Var} (X)=\sigma ^{2}&=\int _{\mathbb {R} }(x-\mu )^{2}f(x)\,dx\\[4pt]&=\int _{\mathbb {R} }x^{2}f(x)\,dx-2\mu \int _{\mathbb {R} }xf(x)\,dx+\mu ^{2}\int _{\mathbb {R} }f(x)\,dx\\[4pt]&=\int _{\mathbb {R} }x^{2}\,dF(x)-2\mu \int _{\mathbb {R} }x\,dF(x)+\mu ^{2}\int _{\mathbb {R} }\,dF(x)\\[4pt]&=\int _{\mathbb {R} }x^{2}\,dF(x)-2\mu \cdot \mu +\mu ^{2}\cdot 1\\[4pt]&=\int _{\mathbb {R} }x^{2}\,dF(x)-\mu ^{2},\end{aligned}}}
または同等に、
ヴァール
(
バツ
)
=
∫
R
x
2
ふ
(
x
)
d
x
−
μ
2
、
{\displaystyle \operatorname {Var} (X)=\int _{\mathbb {R} }x^{2}f(x)\,dx-\mu ^{2},}
ここで の期待値は 次のように与えられる。
μ
{\displaystyle \mu}
バツ
{\displaystyle X}
μ
=
∫
R
x
ふ
(
x
)
d
x
=
∫
R
x
d
ふ
(
x
)
。
{\displaystyle \mu =\int _{\mathbb {R} }xf(x)\,dx=\int _{\mathbb {R} }x\,dF(x).}
これらの式では、および に関する積分はそれぞれ ルベーグ積分 と ルベーグ・スティルチェス 積分
です。
d
x
{\displaystyle dx}
d
ふ
(
x
)
{\displaystyle dF(x)}
関数が すべての有限区間で リーマン積分可能 である場合 、
x
2
ふ
(
x
)
{\displaystyle x^{2}f(x)}
[
1つの
、
b
]
⊂
R
、
{\displaystyle [a,b]\subset \mathbb {R} ,}
ヴァール
(
バツ
)
=
∫
−
∞
+
∞
x
2
ふ
(
x
)
d
x
−
μ
2
、
{\displaystyle \operatorname {Var} (X)=\int _{-\infty }^{+\infty }x^{2}f(x)\,dx-\mu ^{2},}
ここで積分は 不定リーマン積分 である。
例
指数分布
パラメータ λ を持つ指数 分布は、 確率密度関数 が次のように与えられる
連続分布である。
ふ
(
x
)
=
λ
e
−
λ
x
{\displaystyle f(x)=\lambda e^{-\lambda x}}
区間 [0,∞) 上の平均は
え
[
バツ
]
=
∫
0
∞
x
λ
e
−
λ
x
d
x
=
1
λ
。
{\displaystyle \operatorname {E} [X]=\int _{0}^{\infty }x\lambda e^{-\lambda x}\,dx={\frac {1}{\lambda }}.}
部分積分を 使用し 、すでに計算された期待値を利用すると、次のようになります。
え
[
バツ
2
]
=
∫
0
∞
x
2
λ
e
−
λ
x
d
x
=
[
−
x
2
e
−
λ
x
]
0
∞
+
∫
0
∞
2
x
e
−
λ
x
d
x
=
0
+
2
λ
え
[
バツ
]
=
2
λ
2
。
{\displaystyle {\begin{aligned}\operatorname {E} \left[X^{2}\right]&=\int _{0}^{\infty }x^{2}\lambda e^{-\lambda x}\,dx\\&=\left[-x^{2}e^{-\lambda x}\right]_{0}^{\infty }+\int _{0}^{\infty }2xe^{-\lambda x}\,dx\\&=0+{\frac {2}{\lambda }}\operatorname {E} [X]\\&={\frac {2}{\lambda ^{2}}}.\end{aligned}}}
したがって、 X の分散は 次のように表される。
Var
(
X
)
=
E
[
X
2
]
−
E
[
X
]
2
=
2
λ
2
−
(
1
λ
)
2
=
1
λ
2
.
{\displaystyle \operatorname {Var} (X)=\operatorname {E} \left[X^{2}\right]-\operatorname {E} [X]^{2}={\frac {2}{\lambda ^{2}}}-\left({\frac {1}{\lambda }}\right)^{2}={\frac {1}{\lambda ^{2}}}.}
フェアダイ
公平な 6面サイコロは、 結果が1から6までで、それぞれが1/6の等確率である離散確率変数 Xとしてモデル化できます 。X の期待値は 、したがって、 X の分散は 、
(
1
+
2
+
3
+
4
+
5
+
6
)
/
6
=
7
/
2.
{\displaystyle (1+2+3+4+5+6)/6=7/2.}
Var
(
X
)
=
∑
i
=
1
6
1
6
(
i
−
7
2
)
2
=
1
6
(
(
−
5
/
2
)
2
+
(
−
3
/
2
)
2
+
(
−
1
/
2
)
2
+
(
1
/
2
)
2
+
(
3
/
2
)
2
+
(
5
/
2
)
2
)
=
35
12
≈
2.92.
{\displaystyle {\begin{aligned}\operatorname {Var} (X)&=\sum _{i=1}^{6}{\frac {1}{6}}\left(i-{\frac {7}{2}}\right)^{2}\\[5pt]&={\frac {1}{6}}\left((-5/2)^{2}+(-3/2)^{2}+(-1/2)^{2}+(1/2)^{2}+(3/2)^{2}+(5/2)^{2}\right)\\[5pt]&={\frac {35}{12}}\approx 2.92.\end{aligned}}}
n 面 サイコロの 結果 X の分散の一般的な式は
、
Var
(
X
)
=
E
(
X
2
)
−
(
E
(
X
)
)
2
=
1
n
∑
i
=
1
n
i
2
−
(
1
n
∑
i
=
1
n
i
)
2
=
(
n
+
1
)
(
2
n
+
1
)
6
−
(
n
+
1
2
)
2
=
n
2
−
1
12
.
{\displaystyle {\begin{aligned}\operatorname {Var} (X)&=\operatorname {E} \left(X^{2}\right)-(\operatorname {E} (X))^{2}\\[5pt]&={\frac {1}{n}}\sum _{i=1}^{n}i^{2}-\left({\frac {1}{n}}\sum _{i=1}^{n}i\right)^{2}\\[5pt]&={\frac {(n+1)(2n+1)}{6}}-\left({\frac {n+1}{2}}\right)^{2}\\[4pt]&={\frac {n^{2}-1}{12}}.\end{aligned}}}
よく使われる確率分布
次の表は、一般的に使用されるいくつかの確率分布の分散を示しています。
プロパティ
基本的なプロパティ
分散は、平方が正またはゼロであるため、負ではありません。
Var
(
X
)
≥
0.
{\displaystyle \operatorname {Var} (X)\geq 0.}
定数の分散はゼロです。
Var
(
a
)
=
0.
{\displaystyle \operatorname {Var} (a)=0.}
逆に、ランダム変数の分散が 0 の場合、それは ほぼ確実に 定数です。つまり、常に同じ値を持ちます。
Var
(
X
)
=
0
⟺
∃
a
:
P
(
X
=
a
)
=
1.
{\displaystyle \operatorname {Var} (X)=0\iff \exists a:P(X=a)=1.}
有限性の問題
分布が有限の期待値を持たない場合、 コーシー分布の場合のように、分散も有限にはならない。しかし、 分布 によっては、期待値が有限であるにもかかわらず、有限の分散を持たない 場合がある。例としては、 指数 が
k
{\displaystyle k}
1
<
k
≤
2.
{\displaystyle 1<k\leq 2.}
分解
分散分解または 全分散の法則 の一般的な公式は次の通りである。 と が 2つのランダム変数であり、 の分散が 存在する場合、
X
{\displaystyle X}
Y
{\displaystyle Y}
X
{\displaystyle X}
Var
[
X
]
=
E
(
Var
[
X
∣
Y
]
)
+
Var
(
E
[
X
∣
Y
]
)
.
{\displaystyle \operatorname {Var} [X]=\operatorname {E} (\operatorname {Var} [X\mid Y])+\operatorname {Var} (\operatorname {E} [X\mid Y]).}
の 条件 付き期待値 と 条件 付き分散 は、次のように理解できます。 ランダム変数 Yの特定の値 y が 与えられた場合、 イベント Y = y が与えられたときの条件付き期待値があります 。この量は特定の値 y に依存し、関数です。ランダム変数 Y で評価された同じ関数は 、条件付き期待値です。
E
(
X
∣
Y
)
{\displaystyle \operatorname {E} (X\mid Y)}
X
{\displaystyle X}
Y
{\displaystyle Y}
Var
(
X
∣
Y
)
{\displaystyle \operatorname {Var} (X\mid Y)}
E
(
X
∣
Y
=
y
)
{\displaystyle \operatorname {E} (X\mid Y=y)}
g
(
y
)
=
E
(
X
∣
Y
=
y
)
{\displaystyle g(y)=\operatorname {E} (X\mid Y=y)}
E
(
X
∣
Y
)
=
g
(
Y
)
.
{\displaystyle \operatorname {E} (X\mid Y)=g(Y).}
特に、が 離散確率変数であり、可能な値が 対応する確率で存在すると仮定すると 、全分散の式において、右辺第1項は次のようになる。
Y
{\displaystyle Y}
y
1
,
y
2
,
y
3
…
{\displaystyle y_{1},y_{2},y_{3}\ldots }
p
1
,
p
2
,
p
3
…
,
{\displaystyle p_{1},p_{2},p_{3}\ldots ,}
E
(
Var
[
X
∣
Y
]
)
=
∑
i
p
i
σ
i
2
,
{\displaystyle \operatorname {E} (\operatorname {Var} [X\mid Y])=\sum _{i}p_{i}\sigma _{i}^{2},}
ここで 、同様に、右辺の2番目の項は
σ
i
2
=
Var
[
X
∣
Y
=
y
i
]
{\displaystyle \sigma _{i}^{2}=\operatorname {Var} [X\mid Y=y_{i}]}
Var
(
E
[
X
∣
Y
]
)
=
∑
i
p
i
μ
i
2
−
(
∑
i
p
i
μ
i
)
2
=
∑
i
p
i
μ
i
2
−
μ
2
,
{\displaystyle \operatorname {Var} (\operatorname {E} [X\mid Y])=\sum _{i}p_{i}\mu _{i}^{2}-\left(\sum _{i}p_{i}\mu _{i}\right)^{2}=\sum _{i}p_{i}\mu _{i}^{2}-\mu ^{2},}
ここで 、およびである 。したがって、総分散は次のように与えられる。
μ
i
=
E
[
X
∣
Y
=
y
i
]
{\displaystyle \mu _{i}=\operatorname {E} [X\mid Y=y_{i}]}
μ
=
∑
i
p
i
μ
i
{\displaystyle \mu =\sum _{i}p_{i}\mu _{i}}
Var
[
X
]
=
∑
i
p
i
σ
i
2
+
(
∑
i
p
i
μ
i
2
−
μ
2
)
.
{\displaystyle \operatorname {Var} [X]=\sum _{i}p_{i}\sigma _{i}^{2}+\left(\sum _{i}p_{i}\mu _{i}^{2}-\mu ^{2}\right).}
同様の式が 分散分析 にも適用され、対応する式は次のようになる。
M
S
total
=
M
S
between
+
M
S
within
;
{\displaystyle {\mathit {MS}}_{\text{total}}={\mathit {MS}}_{\text{between}}+{\mathit {MS}}_{\text{within}};}
ここでの 二乗平均は、二乗平均を意味します。 線形回帰 分析では、対応する式は次のようになります。
M
S
{\displaystyle {\mathit {MS}}}
M
S
total
=
M
S
regression
+
M
S
residual
.
{\displaystyle {\mathit {MS}}_{\text{total}}={\mathit {MS}}_{\text{regression}}+{\mathit {MS}}_{\text{residual}}.}
これは分散の加法性からも導き出すことができます。合計(観測)スコアは予測スコアとエラー スコアの合計であり、後者 2 つは相関していないためです。
偏差の二乗和(平方和、 )についても同様の分解が可能です。
S
S
{\displaystyle {\mathit {SS}}}
S
S
total
=
S
S
between
+
S
S
within
,
{\displaystyle {\mathit {SS}}_{\text{total}}={\mathit {SS}}_{\text{between}}+{\mathit {SS}}_{\text{within}},}
S
S
total
=
S
S
regression
+
S
S
residual
.
{\displaystyle {\mathit {SS}}_{\text{total}}={\mathit {SS}}_{\text{regression}}+{\mathit {SS}}_{\text{residual}}.}
CDFからの計算
非負の確率変数の母分散は、 累積分布関数 F で次
のように表すことができます。
2
∫
0
∞
u
(
1
−
F
(
u
)
)
d
u
−
(
∫
0
∞
(
1
−
F
(
u
)
)
d
u
)
2
.
{\displaystyle 2\int _{0}^{\infty }u(1-F(u))\,du-\left(\int _{0}^{\infty }(1-F(u))\,du\right)^{2}.}
この式は、密度 ではなく CDF を便利に表現できる
状況で分散を計算するために使用できます。
特性プロパティ
ランダム変数の2次 モーメントは、 ランダム変数の1次モーメント(つまり平均)の周りでとったときに最小値、つまり になります 。逆に、連続関数が すべてのランダム変数 X に対して を満たす場合、その関数は必ず (a > 0 ) の 形になります 。これは多次元の場合にも当てはまります。 [3]
a
r
g
m
i
n
m
E
(
(
X
−
m
)
2
)
=
E
(
X
)
{\displaystyle \mathrm {argmin} _{m}\,\mathrm {E} \left(\left(X-m\right)^{2}\right)=\mathrm {E} (X)}
φ
{\displaystyle \varphi }
a
r
g
m
i
n
m
E
(
φ
(
X
−
m
)
)
=
E
(
X
)
{\displaystyle \mathrm {argmin} _{m}\,\mathrm {E} (\varphi (X-m))=\mathrm {E} (X)}
φ
(
x
)
=
a
x
2
+
b
{\displaystyle \varphi (x)=ax^{2}+b}
測定単位
期待される絶対偏差 とは異なり 、変数の分散の単位は、変数自体の単位の 2 乗です。たとえば、メートルで測定された変数の分散は、メートルの 2 乗で測定されます。このため、データ セットを分散よりも 標準偏差 または 二乗平均平方根偏差 で記述することが好まれることがよくあります。サイコロの例では、標準偏差は √ 2.9 ≈ 1.7 で、期待される絶対偏差の 1.5 よりわずかに大きくなります。
標準偏差と期待絶対偏差は、どちらも分布の「広がり」の指標として使用できます。標準偏差は期待絶対偏差よりも代数的な操作が容易で、分散とその一般化 共分散 とともに、理論統計で頻繁に使用されます。ただし、期待絶対偏差は、 測定の異常 や過度に 裾の重い分布 から生じる 外れ値 の影響を受けにくいため、より 堅牢な 傾向があります。
伝搬
定数による加算と乗算
分散は 位置パラメータ の変化に対して 不変 です。つまり、変数のすべての値に定数を追加しても、分散は変化しません。
Var
(
X
+
a
)
=
Var
(
X
)
.
{\displaystyle \operatorname {Var} (X+a)=\operatorname {Var} (X).}
すべての値が定数でスケーリングされる場合、分散はその定数の 2 乗でスケーリングされます。
Var
(
a
X
)
=
a
2
Var
(
X
)
.
{\displaystyle \operatorname {Var} (aX)=a^{2}\operatorname {Var} (X).}
2つのランダム変数の和の分散は次のように表される。
Var
(
a
X
+
b
Y
)
=
a
2
Var
(
X
)
+
b
2
Var
(
Y
)
+
2
a
b
Cov
(
X
,
Y
)
{\displaystyle \operatorname {Var} (aX+bY)=a^{2}\operatorname {Var} (X)+b^{2}\operatorname {Var} (Y)+2ab\,\operatorname {Cov} (X,Y)}
Var
(
a
X
−
b
Y
)
=
a
2
Var
(
X
)
+
b
2
Var
(
Y
)
−
2
a
b
Cov
(
X
,
Y
)
{\displaystyle \operatorname {Var} (aX-bY)=a^{2}\operatorname {Var} (X)+b^{2}\operatorname {Var} (Y)-2ab\,\operatorname {Cov} (X,Y)}
ここで 共分散 は です 。
Cov
(
X
,
Y
)
{\displaystyle \operatorname {Cov} (X,Y)}
線形結合
一般に、ランダム変数 の合計の場合 、分散は次のようになります。
N
{\displaystyle N}
{
X
1
,
…
,
X
N
}
{\displaystyle \{X_{1},\dots ,X_{N}\}}
Var
(
∑
i
=
1
N
X
i
)
=
∑
i
,
j
=
1
N
Cov
(
X
i
,
X
j
)
=
∑
i
=
1
N
Var
(
X
i
)
+
∑
i
≠
j
Cov
(
X
i
,
X
j
)
,
{\displaystyle \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\neq j}\operatorname {Cov} (X_{i},X_{j}),}
ビエネメ将軍の身元 も参照 。
これらの結果から、線形結合 の分散は 次のようになります。
Var
(
∑
i
=
1
N
a
i
X
i
)
=
∑
i
,
j
=
1
N
a
i
a
j
Cov
(
X
i
,
X
j
)
=
∑
i
=
1
N
a
i
2
Var
(
X
i
)
+
∑
i
≠
j
a
i
a
j
Cov
(
X
i
,
X
j
)
=
∑
i
=
1
N
a
i
2
Var
(
X
i
)
+
2
∑
1
≤
i
<
j
≤
N
a
i
a
j
Cov
(
X
i
,
X
j
)
.
{\displaystyle {\begin{aligned}\operatorname {Var} \left(\sum _{i=1}^{N}a_{i}X_{i}\right)&=\sum _{i,j=1}^{N}a_{i}a_{j}\operatorname {Cov} (X_{i},X_{j})\\&=\sum _{i=1}^{N}a_{i}^{2}\operatorname {Var} (X_{i})+\sum _{i\not =j}a_{i}a_{j}\operatorname {Cov} (X_{i},X_{j})\\&=\sum _{i=1}^{N}a_{i}^{2}\operatorname {Var} (X_{i})+2\sum _{1\leq i<j\leq N}a_{i}a_{j}\operatorname {Cov} (X_{i},X_{j}).\end{aligned}}}
ランダム変数 が
X
1
,
…
,
X
N
{\displaystyle X_{1},\dots ,X_{N}}
Cov
(
X
i
,
X
j
)
=
0
,
∀
(
i
≠
j
)
,
{\displaystyle \operatorname {Cov} (X_{i},X_{j})=0\ ,\ \forall \ (i\neq j),}
それらは 無相関で あると言われます。前述の式から、ランダム変数 が無相関である場合、それらの合計の分散はそれらの分散の合計に等しくなります。記号で表すと次のようになります。
X
1
,
…
,
X
N
{\displaystyle X_{1},\dots ,X_{N}}
Var
(
∑
i
=
1
N
X
i
)
=
∑
i
=
1
N
Var
(
X
i
)
.
{\displaystyle \operatorname {Var} \left(\sum _{i=1}^{N}X_{i}\right)=\sum _{i=1}^{N}\operatorname {Var} (X_{i}).}
独立したランダム変数は常に無相関であるため( 共分散 § 無相関と独立性 を参照)、上記の式は特にランダム変数 が独立している場合に当てはまります。したがって、合計の分散が分散の合計と等しくなるためには独立性は十分ですが、必要ではありません。
X
1
,
…
,
X
n
{\displaystyle X_{1},\dots ,X_{n}}
線形結合の分散の行列表記
をランダム変数 の列ベクトルとして 定義し 、を スカラー の列ベクトルとして定義します 。したがって、は これらのランダム変数の 線形結合 であり、 は の 転置を 表します 。また、 は の 共分散行列 とします 。 の分散は 次のように与えられます。 [4]
X
{\displaystyle X}
n
{\displaystyle n}
X
1
,
…
,
X
n
{\displaystyle X_{1},\ldots ,X_{n}}
c
{\displaystyle c}
n
{\displaystyle n}
c
1
,
…
,
c
n
{\displaystyle c_{1},\ldots ,c_{n}}
c
T
X
{\displaystyle c^{\mathsf {T}}X}
c
T
{\displaystyle c^{\mathsf {T}}}
c
{\displaystyle c}
Σ
{\displaystyle \Sigma }
X
{\displaystyle X}
c
T
X
{\displaystyle c^{\mathsf {T}}X}
Var
(
c
T
X
)
=
c
T
Σ
c
.
{\displaystyle \operatorname {Var} \left(c^{\mathsf {T}}X\right)=c^{\mathsf {T}}\Sigma c.}
これは、平均の分散が(1の列ベクトルで)次のように表されることを意味します。
Var
(
x
¯
)
=
Var
(
1
n
1
′
X
)
=
1
n
2
1
′
Σ
1.
{\displaystyle \operatorname {Var} \left({\bar {x}}\right)=\operatorname {Var} \left({\frac {1}{n}}1'X\right)={\frac {1}{n^{2}}}1'\Sigma 1.}
変数の合計
他の分散の尺度よりも分散が優先される理由の 1 つは、相関のない ランダム変数の合計 (または差) の分散 が、それらの分散の合計であるためです。
Var
(
∑
i
=
1
n
X
i
)
=
∑
i
=
1
n
Var
(
X
i
)
.
{\displaystyle \operatorname {Var} \left(\sum _{i=1}^{n}X_{i}\right)=\sum _{i=1}^{n}\operatorname {Var} (X_{i}).}
この式はビエネメの 公式 [5] と呼ばれ 、1853年に発見されました。 [6] [7]この式は、変数が 独立して いるというより強い条件で行われることが多いです が、無相関であれば十分です。したがって、すべての変数が同じ分散σ 2 を持つ場合、 n による除算は線形変換であるため 、この式は、それらの平均の分散が
Var
(
X
¯
)
=
Var
(
1
n
∑
i
=
1
n
X
i
)
=
1
n
2
∑
i
=
1
n
Var
(
X
i
)
=
1
n
2
n
σ
2
=
σ
2
n
.
{\displaystyle \operatorname {Var} \left({\overline {X}}\right)=\operatorname {Var} \left({\frac {1}{n}}\sum _{i=1}^{n}X_{i}\right)={\frac {1}{n^{2}}}\sum _{i=1}^{n}\operatorname {Var} \left(X_{i}\right)={\frac {1}{n^{2}}}n\sigma ^{2}={\frac {\sigma ^{2}}{n}}.}
つまり、 n が増加すると平均の分散は減少します。この平均の分散の式は、 中心極限定理 で使用される標本平均の 標準誤差 の定義で使用されます 。
最初の主張を証明するには、次のことを示せば十分である。
Var
(
X
+
Y
)
=
Var
(
X
)
+
Var
(
Y
)
.
{\displaystyle \operatorname {Var} (X+Y)=\operatorname {Var} (X)+\operatorname {Var} (Y).}
一般的な結果は帰納的に導かれる。定義から始めると、
Var
(
X
+
Y
)
=
E
[
(
X
+
Y
)
2
]
−
(
E
[
X
+
Y
]
)
2
=
E
[
X
2
+
2
X
Y
+
Y
2
]
−
(
E
[
X
]
+
E
[
Y
]
)
2
.
{\displaystyle {\begin{aligned}\operatorname {Var} (X+Y)&=\operatorname {E} \left[(X+Y)^{2}\right]-(\operatorname {E} [X+Y])^{2}\\[5pt]&=\operatorname {E} \left[X^{2}+2XY+Y^{2}\right]-(\operatorname {E} [X]+\operatorname {E} [Y])^{2}.\end{aligned}}}
期待演算子 の線形性と、 X と Y の独立性(または無相関性)の仮定を使用すると 、これは次のようにさらに単純化されます。
Var
(
X
+
Y
)
=
E
[
X
2
]
+
2
E
[
X
Y
]
+
E
[
Y
2
]
−
(
E
[
X
]
2
+
2
E
[
X
]
E
[
Y
]
+
E
[
Y
]
2
)
=
E
[
X
2
]
+
E
[
Y
2
]
−
E
[
X
]
2
−
E
[
Y
]
2
=
Var
(
X
)
+
Var
(
Y
)
.
{\displaystyle {\begin{aligned}\operatorname {Var} (X+Y)&=\operatorname {E} \left[X^{2}\right]+2\operatorname {E} [XY]+\operatorname {E} \left[Y^{2}\right]-\left(\operatorname {E} [X]^{2}+2\operatorname {E} [X]\operatorname {E} [Y]+\operatorname {E} [Y]^{2}\right)\\[5pt]&=\operatorname {E} \left[X^{2}\right]+\operatorname {E} \left[Y^{2}\right]-\operatorname {E} [X]^{2}-\operatorname {E} [Y]^{2}\\[5pt]&=\operatorname {Var} (X)+\operatorname {Var} (Y).\end{aligned}}}
一般に、 n 個の変数の合計の分散は、それらの 共分散 の合計です 。
Var
(
∑
i
=
1
n
X
i
)
=
∑
i
=
1
n
∑
j
=
1
n
Cov
(
X
i
,
X
j
)
=
∑
i
=
1
n
Var
(
X
i
)
+
2
∑
1
≤
i
<
j
≤
n
Cov
(
X
i
,
X
j
)
.
{\displaystyle \operatorname {Var} \left(\sum _{i=1}^{n}X_{i}\right)=\sum _{i=1}^{n}\sum _{j=1}^{n}\operatorname {Cov} \left(X_{i},X_{j}\right)=\sum _{i=1}^{n}\operatorname {Var} \left(X_{i}\right)+2\sum _{1\leq i<j\leq n}\operatorname {Cov} \left(X_{i},X_{j}\right).}
(注: 2 番目の等式は、 Cov( X i , X i ) = Var( X i ) という事実から生じます 。)
ここで、は 共分散 であり 、独立したランダム変数(存在する場合)に対しては 0 になります。この式は、和の分散が、成分の共分散行列のすべての要素の合計に等しいことを示しています。次の式は、和の分散が、共分散行列の対角線の合計に、その上三角要素(または下三角要素)の合計の 2 倍を加えたものであることを示しています。これは、共分散行列が対称であることを強調しています。この式は、 古典的な検定理論 の Cronbach のアルファ 理論で使用されます。
Cov
(
⋅
,
⋅
)
{\displaystyle \operatorname {Cov} (\cdot ,\cdot )}
したがって、変数の分散が σ 2 で 、異なる変数の 平均 相関が ρ である場合、それらの平均の分散は
Var
(
X
¯
)
=
σ
2
n
+
n
−
1
n
ρ
σ
2
.
{\displaystyle \operatorname {Var} \left({\overline {X}}\right)={\frac {\sigma ^{2}}{n}}+{\frac {n-1}{n}}\rho \sigma ^{2}.}
これは、平均の分散が相関の平均とともに増加することを意味します。言い換えれば、追加の相関観測は、追加の独立した観測ほど 平均の不確実性を 減らすのに効果的ではありません。さらに、変数が単位分散を持っている場合、たとえば標準化されている場合、これは次のように単純化されます。
Var
(
X
¯
)
=
1
n
+
n
−
1
n
ρ
.
{\displaystyle \operatorname {Var} \left({\overline {X}}\right)={\frac {1}{n}}+{\frac {n-1}{n}}\rho .}
この式は、古典的な検定理論のスピアマン・ブラウン予測式 で使用されている 。平均相関が一定であるか収束する限り、 nが 無限大になると ρ に収束する。したがって、等しい相関または収束する平均相関を持つ標準化変数の平均の分散については、
lim
n
→
∞
Var
(
X
¯
)
=
ρ
.
{\displaystyle \lim _{n\to \infty }\operatorname {Var} \left({\overline {X}}\right)=\rho .}
したがって、多数の標準化変数の平均の分散は、それらの平均相関とほぼ等しくなります。これにより、 大数の法則 では標本平均は独立変数に対して収束すると規定されているにもかかわらず、相関変数の標本平均は一般に母平均に収束しないことが明らかになります。
事前に、ある基準に従って許容される観測数がわからないままサンプルが採取される場合がある。そのような場合、サンプルサイズ Nはランダム変数であり、その変動は X の変動に加算される 。
Var
(
∑
i
=
1
N
X
i
)
=
E
[
N
]
Var
(
X
)
+
Var
(
N
)
(
E
[
X
]
)
2
{\displaystyle \operatorname {Var} \left(\sum _{i=1}^{N}X_{i}\right)=\operatorname {E} \left[N\right]\operatorname {Var} (X)+\operatorname {Var} (N)(\operatorname {E} \left[X\right])^{2}}
[8]
これは全分散の法則 に従います 。
N が ポアソン分布 に従う 場合 、推定値 n = N となります 。したがって、 の推定値 は となり 、次のようになります
( 標本平均の標準誤差 を参照)。
E
[
N
]
=
Var
(
N
)
{\displaystyle \operatorname {E} [N]=\operatorname {Var} (N)}
Var
(
∑
i
=
1
n
X
i
)
{\displaystyle \operatorname {Var} \left(\sum _{i=1}^{n}X_{i}\right)}
n
S
x
2
+
n
X
¯
2
{\displaystyle n{S_{x}}^{2}+n{\bar {X}}^{2}}
SE
(
X
¯
)
=
S
x
2
+
X
¯
2
n
{\displaystyle \operatorname {SE} ({\bar {X}})={\sqrt {\frac {{S_{x}}^{2}+{\bar {X}}^{2}}{n}}}}
変数の加重合計
スケーリング特性とビエネイメの公式、そして共 分散 Cov( aX , bY )= abCov ( X , Y ) の特性を組み合わせると 、次のことが分かります。
Var
(
a
X
±
b
Y
)
=
a
2
Var
(
X
)
+
b
2
Var
(
Y
)
±
2
a
b
Cov
(
X
,
Y
)
.
{\displaystyle \operatorname {Var} (aX\pm bY)=a^{2}\operatorname {Var} (X)+b^{2}\operatorname {Var} (Y)\pm 2ab\,\operatorname {Cov} (X,Y).}
これは、変数の加重合計において、最大の重みを持つ変数が、合計の分散において不釣り合いに大きな重みを持つことを意味します。たとえば、 X と Yが無相関で、 X の重みが Y の重みの 2 倍である場合、 X の分散の重みは Y の分散の重みの 4 倍になります 。
上記の式は、複数の変数の加重合計に拡張できます。
Var
(
∑
i
n
a
i
X
i
)
=
∑
i
=
1
n
a
i
2
Var
(
X
i
)
+
2
∑
1
≤
i
∑
<
j
≤
n
a
i
a
j
Cov
(
X
i
,
X
j
)
{\displaystyle \operatorname {Var} \left(\sum _{i}^{n}a_{i}X_{i}\right)=\sum _{i=1}^{n}a_{i}^{2}\operatorname {Var} (X_{i})+2\sum _{1\leq i}\sum _{<j\leq n}a_{i}a_{j}\operatorname {Cov} (X_{i},X_{j})}
変数の積
独立変数の積
2つの変数XとYが 独立して いる場合、それらの積の分散は [9]で与えられる。
Var
(
X
Y
)
=
[
E
(
X
)
]
2
Var
(
Y
)
+
[
E
(
Y
)
]
2
Var
(
X
)
+
Var
(
X
)
Var
(
Y
)
.
{\displaystyle \operatorname {Var} (XY)=[\operatorname {E} (X)]^{2}\operatorname {Var} (Y)+[\operatorname {E} (Y)]^{2}\operatorname {Var} (X)+\operatorname {Var} (X)\operatorname {Var} (Y).}
同様に、期待値の基本的な性質を用いると、次のように表される。
Var
(
X
Y
)
=
E
(
X
2
)
E
(
Y
2
)
−
[
E
(
X
)
]
2
[
E
(
Y
)
]
2
.
{\displaystyle \operatorname {Var} (XY)=\operatorname {E} \left(X^{2}\right)\operatorname {E} \left(Y^{2}\right)-[\operatorname {E} (X)]^{2}[\operatorname {E} (Y)]^{2}.}
統計的従属変数の積
一般に、2 つの変数が統計的に依存している場合、それらの積の分散は次のように表されます。
Var
(
X
Y
)
=
E
[
X
2
Y
2
]
−
[
E
(
X
Y
)
]
2
=
Cov
(
X
2
,
Y
2
)
+
E
(
X
2
)
E
(
Y
2
)
−
[
E
(
X
Y
)
]
2
=
Cov
(
X
2
,
Y
2
)
+
(
Var
(
X
)
+
[
E
(
X
)
]
2
)
(
Var
(
Y
)
+
[
E
(
Y
)
]
2
)
−
[
Cov
(
X
,
Y
)
+
E
(
X
)
E
(
Y
)
]
2
{\displaystyle {\begin{aligned}\operatorname {Var} (XY)={}&\operatorname {E} \left[X^{2}Y^{2}\right]-[\operatorname {E} (XY)]^{2}\\[5pt]={}&\operatorname {Cov} \left(X^{2},Y^{2}\right)+\operatorname {E} (X^{2})\operatorname {E} \left(Y^{2}\right)-[\operatorname {E} (XY)]^{2}\\[5pt]={}&\operatorname {Cov} \left(X^{2},Y^{2}\right)+\left(\operatorname {Var} (X)+[\operatorname {E} (X)]^{2}\right)\left(\operatorname {Var} (Y)+[\operatorname {E} (Y)]^{2}\right)\\[5pt]&-[\operatorname {Cov} (X,Y)+\operatorname {E} (X)\operatorname {E} (Y)]^{2}\end{aligned}}}
任意の関数
デルタ 法は、 2次の テイラー展開 を使用して、1つ以上のランダム変数の関数の分散を近似します。 ランダム変数の関数のモーメントについてはテイラー展開を 参照してください。たとえば、1つの変数の関数の近似分散は次のように表されます。
Var
[
f
(
X
)
]
≈
(
f
′
(
E
[
X
]
)
)
2
Var
[
X
]
{\displaystyle \operatorname {Var} \left[f(X)\right]\approx \left(f'(\operatorname {E} \left[X\right])\right)^{2}\operatorname {Var} \left[X\right]}
ただし、 f は2 回微分可能であり、 X の平均と分散は 有限であるものとする。
母分散と標本分散
昨日の一日の雨量の測定などの現実世界の観測は、通常、実行可能なすべての観測の完全なセットにはなり得ません。そのため、有限セットから計算された分散は、通常、可能性のある観測の全母集団から計算された分散とは一致しません。つまり、 推定 方程式を使用して、限られた観測セットから平均と分散を 推定します 。推定値は、潜在的な観測の全 母集団から観測バイアスなしで抽出された n 個の 観測の サンプル の関数です 。この例では、サンプルは、対象地域内にある利用可能な雨量計から昨日の降雨量の実際の測定値のセットになります。
母平均と母分散の最も単純な推定量は、単に標本の平均と分散、 標本平均 と (補正されていない)標本分散 です。これらは 矛盾のない推定量 (標本数が増えるにつれて母集団全体の値に収束する)ですが、改善することができます。最も単純な場合、標本分散は、 (標本)平均 の偏差 の二乗和を標本数 nで割って計算されます 。ただし、 n 以外の値を使用すると、 推定量がさまざまな方法で改善されます。分母の一般的な値は 、n、 n − 1、 n + 1、および n − 1.5 の4 つです 。n は最も単純(標本の分散)で、 n − 1 はバイアスを排除し、 n + 1 は正規分布の 平均二乗誤差 を最小化し、 n − 1.5 は正規分布の
標準偏差の不偏推定 におけるバイアスをほぼ排除します。
まず、真の母平均が不明な場合、標本分散(真の平均の代わりに標本平均を使用する)は 偏りのある推定値 になります。つまり、分散を係数 ( n − 1) / n で過小評価します。この係数を修正して、標本平均を中心とした偏差の二乗和を n ではなく n -1 で割ることを ベッセル補正 と呼びます。結果として得られる推定値は偏りがなく、 (修正)標本分散 または 偏りのない標本分散 と呼ばれます 。平均が、分散の推定に使用されたのと同じ標本以外の方法で決定される場合、この偏りは発生せず、分散は(独立に既知の)平均を中心とした標本の分散として安全に推定できます。
第二に、標本分散は一般に 標本分散と母集団分散の間の 平均二乗誤差 を最小化しません。バイアスを修正すると、状況は悪化することがよくあります。修正された標本分散よりもパフォーマンスのよいスケール係数を常に選択できますが、最適なスケール係数は母集団の 過剰尖度( 平均二乗誤差: 分散を参照) に依存し、バイアスを導入します。これは常に、不偏推定量をスケールダウン ( n − 1より大きい数で割る) することから成り、 収縮推定量 の簡単な例です。つまり 、不偏推定量をゼロに向かって「収縮」します。正規分布の場合、 n + 1 ( n − 1 または n ではなく ) で割ると、平均二乗誤差が最小化されます。ただし、結果として得られる推定量にはバイアスがあり、 バイアス標本分散 と呼ばれます。
人口分散
一般に、 値 x i を持つサイズ Nの 有限 母集団 の 母分散は 、母平均が 、が 期待値 演算子で ある で与えられます 。
σ
2
=
1
N
∑
i
=
1
N
(
x
i
−
μ
)
2
=
1
N
∑
i
=
1
N
(
x
i
2
−
2
μ
x
i
+
μ
2
)
=
(
1
N
∑
i
=
1
N
x
i
2
)
−
2
μ
(
1
N
∑
i
=
1
N
x
i
)
+
μ
2
=
E
[
x
i
2
]
−
μ
2
{\displaystyle {\begin{aligned}\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 {1}{N}}\sum _{i=1}^{N}x_{i}^{2}\right)-2\mu \left({\frac {1}{N}}\sum _{i=1}^{N}x_{i}\right)+\mu ^{2}\\[5pt]&=\operatorname {E} [x_{i}^{2}]-\mu ^{2}\end{aligned}}}
μ
=
E
[
x
i
]
=
1
N
∑
i
=
1
N
x
i
{\textstyle \mu =\operatorname {E} [x_{i}]={\frac {1}{N}}\sum _{i=1}^{N}x_{i}}
E
[
x
i
2
]
=
(
1
N
∑
i
=
1
N
x
i
2
)
{\textstyle \operatorname {E} [x_{i}^{2}]=\left({\frac {1}{N}}\sum _{i=1}^{N}x_{i}^{2}\right)}
E
{\textstyle \operatorname {E} }
母分散は [10]を使って計算することもできる。
σ
2
=
1
N
2
∑
i
<
j
(
x
i
−
x
j
)
2
=
1
2
N
2
∑
i
,
j
=
1
N
(
x
i
−
x
j
)
2
.
{\displaystyle \sigma ^{2}={\frac {1}{N^{2}}}\sum _{i<j}\left(x_{i}-x_{j}\right)^{2}={\frac {1}{2N^{2}}}\sum _{i,j=1}^{N}\left(x_{i}-x_{j}\right)^{2}.}
(右側には合計に重複する項がありますが、中央の側には合計する一意の項のみがあります。) これは 、母分散が生成確率分布の分散と一致するため当てはまります。この意味で、母集団の概念は、無限の母集団を持つ連続ランダム変数に拡張できます。
1
2
N
2
∑
i
,
j
=
1
N
(
x
i
−
x
j
)
2
=
1
2
N
2
∑
i
,
j
=
1
N
(
x
i
2
−
2
x
i
x
j
+
x
j
2
)
=
1
2
N
∑
j
=
1
N
(
1
N
∑
i
=
1
N
x
i
2
)
−
(
1
N
∑
i
=
1
N
x
i
)
(
1
N
∑
j
=
1
N
x
j
)
+
1
2
N
∑
i
=
1
N
(
1
N
∑
j
=
1
N
x
j
2
)
=
1
2
(
σ
2
+
μ
2
)
−
μ
2
+
1
2
(
σ
2
+
μ
2
)
=
σ
2
.
{\displaystyle {\begin{aligned}&{\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{aligned}}}
サンプル分散
偏ったサンプル分散
多くの実際の状況では、母集団の真の分散は 事前には わからないため、何らかの方法で計算する必要があります。非常に大きな母集団を扱う場合、母集団内のすべてのオブジェクトを数えることは不可能であるため、計算は 母集団の サンプルに対して実行する必要があります。 [11]これは一般に サンプル分散 または 経験的分散 と呼ばれます 。サンプル分散は、連続分布の分散をその分布のサンプルから推定するためにも適用できます。
サイズの母集団から n個 の 値 Y 1 , ..., Y n を置換した標本( n < N ) を採取し 、この標本に基づいて分散を推定する。 [12]標本データの分散を直接取ると、 偏差の二乗 の平均が得られる 。
N
{\textstyle N}
S
~
Y
2
=
1
n
∑
i
=
1
n
(
Y
i
−
Y
¯
)
2
=
(
1
n
∑
i
=
1
n
Y
i
2
)
−
Y
¯
2
=
1
n
2
∑
i
,
j
:
i
<
j
(
Y
i
−
Y
j
)
2
.
{\displaystyle {\tilde {S}}_{Y}^{2}={\frac {1}{n}}\sum _{i=1}^{n}\left(Y_{i}-{\overline {Y}}\right)^{2}=\left({\frac {1}{n}}\sum _{i=1}^{n}Y_{i}^{2}\right)-{\overline {Y}}^{2}={\frac {1}{n^{2}}}\sum _{i,j\,:\,i<j}\left(Y_{i}-Y_{j}\right)^{2}.}
[13]
(この式の導出については、母分散のセクションを参照してください。)ここで、は 標本平均 を表します 。
Y
¯
{\displaystyle {\overline {Y}}}
Y
¯
=
1
n
∑
i
=
1
n
Y
i
.
{\displaystyle {\overline {Y}}={\frac {1}{n}}\sum _{i=1}^{n}Y_{i}.}
Y i は ランダムに選択される ため、 と は両方とも ランダム変数 です。それらの期待値は 、母集団からのサイズ nのすべての可能なサンプル { Y i } の集合を平均することで評価できます。 これにより、次の式が得られます。
Y
¯
{\displaystyle {\overline {Y}}}
S
~
Y
2
{\displaystyle {\tilde {S}}_{Y}^{2}}
S
~
Y
2
{\displaystyle {\tilde {S}}_{Y}^{2}}
E
[
S
~
Y
2
]
=
E
[
1
n
∑
i
=
1
n
(
Y
i
−
1
n
∑
j
=
1
n
Y
j
)
2
]
=
1
n
∑
i
=
1
n
E
[
Y
i
2
−
2
n
Y
i
∑
j
=
1
n
Y
j
+
1
n
2
∑
j
=
1
n
Y
j
∑
k
=
1
n
Y
k
]
=
1
n
∑
i
=
1
n
(
E
[
Y
i
2
]
−
2
n
(
∑
j
≠
i
E
[
Y
i
Y
j
]
+
E
[
Y
i
2
]
)
+
1
n
2
∑
j
=
1
n
∑
k
≠
j
n
E
[
Y
j
Y
k
]
+
1
n
2
∑
j
=
1
n
E
[
Y
j
2
]
)
=
1
n
∑
i
=
1
n
(
n
−
2
n
E
[
Y
i
2
]
−
2
n
∑
j
≠
i
E
[
Y
i
Y
j
]
+
1
n
2
∑
j
=
1
n
∑
k
≠
j
n
E
[
Y
j
Y
k
]
+
1
n
2
∑
j
=
1
n
E
[
Y
j
2
]
)
=
1
n
∑
i
=
1
n
[
n
−
2
n
(
σ
2
+
μ
2
)
−
2
n
(
n
−
1
)
μ
2
+
1
n
2
n
(
n
−
1
)
μ
2
+
1
n
(
σ
2
+
μ
2
)
]
=
n
−
1
n
σ
2
.
{\displaystyle {\begin{aligned}\operatorname {E} [{\tilde {S}}_{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 {1}{n}}\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 {1}{n}}\sum _{i=1}^{n}\left(\operatorname {E} \left[Y_{i}^{2}\right]-{\frac {2}{n}}\left(\sum _{j\neq i}\operatorname {E} \left[Y_{i}Y_{j}\right]+\operatorname {E} \left[Y_{i}^{2}\right]\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 {1}{n}}\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 {1}{n}}\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{aligned}}}
ここでは、 および の独立性により 、 セクションで導出された および が使用されます。
σ
2
=
E
[
Y
i
2
]
−
μ
2
{\textstyle \sigma ^{2}=\operatorname {E} [Y_{i}^{2}]-\mu ^{2}}
E
[
Y
i
Y
j
]
=
E
[
Y
i
]
E
[
Y
j
]
=
μ
2
{\textstyle \operatorname {E} [Y_{i}Y_{j}]=\operatorname {E} [Y_{i}]\operatorname {E} [Y_{j}]=\mu ^{2}}
Y
i
{\textstyle Y_{i}}
Y
j
{\textstyle Y_{j}}
したがって、 の期待値は母分散(真の分散)よりもその係数だけ小さいため、の係数 によって偏りのある母分散の推定値が得られます 。このため、は 偏りのある標本分散 と呼ばれます 。
S
~
Y
2
{\textstyle {\tilde {S}}_{Y}^{2}}
n
−
1
n
{\textstyle {\frac {n-1}{n}}}
S
~
Y
2
{\textstyle {\tilde {S}}_{Y}^{2}}
S
~
Y
2
{\textstyle {\tilde {S}}_{Y}^{2}}
偏りのないサンプル分散
この偏りを補正すると、 次のように表される 不偏標本分散 が得られます。
S
2
{\displaystyle S^{2}}
S
2
=
n
n
−
1
S
~
Y
2
=
n
n
−
1
[
1
n
∑
i
=
1
n
(
Y
i
−
Y
¯
)
2
]
=
1
n
−
1
∑
i
=
1
n
(
Y
i
−
Y
¯
)
2
{\displaystyle S^{2}={\frac {n}{n-1}}{\tilde {S}}_{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}}
どちらの推定値も、文脈によってバージョンを決定できる場合は、
単に サンプル分散 と呼ばれることがあります。同じ証明は、連続確率分布から取得されたサンプルにも適用できます。
n − 1という項の使用は ベッセル補正 と呼ばれ 、 標本共分散 や 標本標準偏差 (分散の平方根) でも使用されます。平方根は 凹関数 であるため、分布に依存する負のバイアス ( ジェンセンの不等式 による) が導入され、そのため、補正された標本標準偏差 (ベッセル補正を使用) はバイアスされます。 標準偏差の不偏推定は技術的に複雑な問題ですが、正規分布の場合、 n − 1.5という項を使用すると、 ほぼ不偏推定値が得られます。
不偏標本分散は 関数 ƒ ( y 1 , y 2 ) = ( y 1 − y 2 ) 2 /2の U統計量 であり、これは母集団の2要素サブセットにわたって2標本統計量を平均することによって得られることを意味します。
例
数値セット {10、15、30、45、57、52、63、72、81、93、102、105} の場合、このセットが何らかの測定のデータ母集団全体である場合、分散は、このセットの平均に関する偏差の二乗の合計をセット メンバーの数である 12 で割った母分散 932.743 になります。セットが母集団全体のサンプルである場合、不偏サンプル分散は、サンプルの平均に関する偏差の二乗の合計を 12 ではなく 11 で割った 1017.538 として計算できます。Microsoft Excel の関数 VAR.S は不偏サンプル分散を返し、VAR.P は母分散を表します。
標本分散の分布
ランダム変数 の関数であるため 、標本分散自体もランダム変数であり、その分布を調べるのは自然なことです。Y i が正規分布からの独立した観測値である場合 、 コクラン の 定理 は 、 不偏標本分散 S 2がスケールされた カイ2乗分布 に従うことを示しています ( 漸近的性質 と 基本的な証明 も参照)。 [14]
(
n
−
1
)
S
2
σ
2
∼
χ
n
−
1
2
{\displaystyle (n-1){\frac {S^{2}}{\sigma ^{2}}}\sim \chi _{n-1}^{2}}
ここで σ 2 は 母分散である。直接的な帰結として、
E
(
S
2
)
=
E
(
σ
2
n
−
1
χ
n
−
1
2
)
=
σ
2
,
{\displaystyle \operatorname {E} \left(S^{2}\right)=\operatorname {E} \left({\frac {\sigma ^{2}}{n-1}}\chi _{n-1}^{2}\right)=\sigma ^{2},}
そして [15]
Var
[
S
2
]
=
Var
(
σ
2
n
−
1
χ
n
−
1
2
)
=
σ
4
(
n
−
1
)
2
Var
(
χ
n
−
1
2
)
=
2
σ
4
n
−
1
.
{\displaystyle \operatorname {Var} \left[S^{2}\right]=\operatorname {Var} \left({\frac {\sigma ^{2}}{n-1}}\chi _{n-1}^{2}\right)={\frac {\sigma ^{4}}{(n-1)^{2}}}\operatorname {Var} \left(\chi _{n-1}^{2}\right)={\frac {2\sigma ^{4}}{n-1}}.}
Y i が独立かつ同一分布に従うが、必ずしも正規分布するとは限らない 場合、 [16]
E
[
S
2
]
=
σ
2
,
Var
[
S
2
]
=
σ
4
n
(
κ
−
1
+
2
n
−
1
)
=
1
n
(
μ
4
−
n
−
3
n
−
1
σ
4
)
,
{\displaystyle \operatorname {E} \left[S^{2}\right]=\sigma ^{2},\quad \operatorname {Var} \left[S^{2}\right]={\frac {\sigma ^{4}}{n}}\left(\kappa -1+{\frac {2}{n-1}}\right)={\frac {1}{n}}\left(\mu _{4}-{\frac {n-3}{n-1}}\sigma ^{4}\right),}
ここで、 κ は 分布の 尖度、 μ 4 は4次 中心モーメント です。
大数の法則 の条件が 観測値の二乗に当てはまる場合、 S 2 は σ 2 の 一貫した推定値 になります 。推定値の分散が漸近的にゼロに近づくことがわかります。漸近的に同等の式は、Kenney and Keeping (1951:164)、Rose and Smith (2002:264)、およびWeisstein (nd)で与えられています。 [17] [18] [19]
サミュエルソンの不等式
サミュエルソンの不等式 は、サンプル平均と(偏りのある)分散が計算されていると仮定して、サンプル内の個々の観測値が取り得る値の境界を示す結果である。 [20] 値は、
y
¯
±
σ
Y
(
n
−
1
)
1
/
2
.
{\displaystyle {\bar {y}}\pm \sigma _{Y}(n-1)^{1/2}.}
調和平均と算術平均との関係
[21] によれば、 正の実数の
サンプル{ y i }に対して、
σ
y
2
≤
2
y
max
(
A
−
H
)
,
{\displaystyle \sigma _{y}^{2}\leq 2y_{\max }(A-H),}
ここで、 y max はサンプルの最大値、 A は算術平均、 H は サンプルの 調和平均 、はサンプルの (偏りのある) 分散です。
σ
y
2
{\displaystyle \sigma _{y}^{2}}
この境界は改善されており、分散は次のように制限されることが分かっています。
σ
y
2
≤
y
max
(
A
−
H
)
(
y
max
−
A
)
y
max
−
H
,
{\displaystyle \sigma _{y}^{2}\leq {\frac {y_{\max }(A-H)(y_{\max }-A)}{y_{\max }-H}},}
σ
y
2
≥
y
min
(
A
−
H
)
(
A
−
y
min
)
H
−
y
min
,
{\displaystyle \sigma _{y}^{2}\geq {\frac {y_{\min }(A-H)(A-y_{\min })}{H-y_{\min }}},}
ここで ymin は サンプルの最小値である。 [22]
分散の等価性の検定
分散の等価性の F 検定 と カイ 二乗検定は、 サンプルが正規分布している場合に適しています。非正規性の場合、2 つ以上の分散の等価性を検定することはより困難になります。
ノンパラメトリック検定がいくつか提案されています。これには、Barton–David–Ansari–Freund–Siegel–Tukey 検定、Capon 検定、 Mood 検定 、Klotz 検定、Sukhatme 検定が含まれます。Sukhatme 検定は 2 つの分散に適用され、両方の 中央値 が既知でゼロに等しいことが必要です。Mood 検定、Klotz 検定、Capon 検定、Barton–David–Ansari–Freund–Siegel–Tukey 検定も 2 つの分散に適用されます。これらの検定では中央値が不明でもかまいませんが、2 つの中央値が等しいことが必要です。
Lehmann 検定は、2 つの分散のパラメトリック検定です。この検定にはいくつかのバリエーションが知られています。分散の等価性に関するその他の検定には、Box 検定、Box–Anderson 検定、Moses 検定などがあります。
ブートストラップ法 や ジャックナイフ法 などの再サンプリング法は 、分散の等価性を検定するために使用できます。
慣性モーメント
確率分布の分散は、 古典力学 における 、質量分布の線に沿った質量中心の周りの回転に関する 慣性モーメントに類似しています。 [ 引用が必要 ] この類似性のため、分散のようなものは 確率分布 の モーメント と呼ばれます。 [ 引用が必要 ]共分散行列は、多変量分布の 慣性モーメントテンソル に関連しています。 の共分散 行列を持つ n 点のクラウドの慣性モーメントは、 次のように与えられます。 [ 引用が必要 ]
Σ
{\displaystyle \Sigma }
I
=
n
(
1
3
×
3
tr
(
Σ
)
−
Σ
)
.
{\displaystyle I=n\left(\mathbf {1} _{3\times 3}\operatorname {tr} (\Sigma )-\Sigma \right).}
物理学と統計学における慣性モーメントの違いは、直線上に集まった点については明らかです。多くの点が x 軸の近くにあり、それに沿って分布しているとします。共分散行列は次のようになります
。
Σ
=
[
10
0
0
0
0.1
0
0
0
0.1
]
.
{\displaystyle \Sigma ={\begin{bmatrix}10&0&0\\0&0.1&0\\0&0&0.1\end{bmatrix}}.}
つまり、 x方向の分散が最も大きい。物理学者は これをx 軸 周りの モーメントが低いと考える ので、慣性モーメントテンソルは
I
=
n
[
0.2
0
0
0
10.1
0
0
0
10.1
]
.
{\displaystyle I=n{\begin{bmatrix}0.2&0&0\\0&10.1&0\\0&0&10.1\end{bmatrix}}.}
半分散
半 分散は 分散と同じ方法で計算されますが、平均を下回る観測値のみが計算に含まれます。 また、さまざまな応用分野では特定の尺度として説明されています。歪んだ分布の場合、半分散は分散では得られない追加情報を提供することができます。 [23]
Semivariance
=
1
n
∑
i
:
x
i
<
μ
(
x
i
−
μ
)
2
{\displaystyle {\text{Semivariance}}={1 \over {n}}\sum _{i:x_{i}<\mu }(x_{i}-\mu )^{2}}
半分散に関連する不等式については、 チェビシェフの不等式 § 半分散を 参照してください。
語源
分散 という用語は、 ロナルド・フィッシャー が1918年に発表した論文 「メンデル遺伝の仮定に基づく親族間の相関関係」 で初めて導入された 。 [24]
利用可能な統計の大部分は、 人間の測定値 の平均からの偏差が 誤差の正規法則 に非常に近いこと、したがって変動は 平均二乗誤差 の 平方根 に対応する 標準偏差 によって一様に測定できることを示しています。 標準偏差 および を持つ均一な母集団分布を生成する可能性のある 2 つの独立した変動原因がある場合 、 両方の原因が一緒に作用すると、分布の標準偏差は になります 。 したがって、変動原因を分析するには、変動の尺度として標準偏差の二乗を扱うことが望ましいです。 この量を分散と呼びます...
σ
1
{\displaystyle \sigma _{1}}
σ
2
{\displaystyle \sigma _{2}}
σ
1
2
+
σ
2
2
{\displaystyle {\sqrt {\sigma _{1}^{2}+\sigma _{2}^{2}}}}
一般化
複素変数の場合
が スカラー 複素数 値ランダム変数で、値が の 場合、その分散 は の 複素 共役 です 。この分散は実数スカラーです。
x
{\displaystyle x}
C
,
{\displaystyle \mathbb {C} ,}
E
[
(
x
−
μ
)
(
x
−
μ
)
∗
]
,
{\displaystyle \operatorname {E} \left[(x-\mu )(x-\mu )^{*}\right],}
x
∗
{\displaystyle x^{*}}
x
.
{\displaystyle x.}
ベクトル値ランダム変数の場合
マトリックスとして
がベクトル 値のランダム変数で 、値が にあり 、が列ベクトルと考えられる 場合、分散の自然な一般化は となり 、 は の転置であり 、したがって は行ベクトルとなります。結果は 半正定値正方行列となり、一般に 分散共分散行列 (または単に 共分散行列 )
と呼ばれます。
X
{\displaystyle X}
R
n
,
{\displaystyle \mathbb {R} ^{n},}
E
[
(
X
−
μ
)
(
X
−
μ
)
T
]
,
{\displaystyle \operatorname {E} \left[(X-\mu )(X-\mu )^{\operatorname {T} }\right],}
μ
=
E
(
X
)
{\displaystyle \mu =\operatorname {E} (X)}
X
T
{\displaystyle X^{\operatorname {T} }}
X
,
{\displaystyle X,}
が ベクトルおよび複素数値のランダム変数で、値が の場合 、 共分散行列は となり、 は の 共役転置 です [ 要出典 ]。 この行列は半正定値で正方行列でもあります。
X
{\displaystyle X}
C
n
,
{\displaystyle \mathbb {C} ^{n},}
E
[
(
X
−
μ
)
(
X
−
μ
)
†
]
,
{\displaystyle \operatorname {E} \left[(X-\mu )(X-\mu )^{\dagger }\right],}
X
†
{\displaystyle X^{\dagger }}
X
.
{\displaystyle X.}
スカラーとして
ベクトル値ランダム変数の分散のもう一つの一般化は 、行列ではなくスカラー値になる 一般化分散 であり、共分散行列の 行列式 である。一般化分散は、平均の周りの点の多次元散布と関連していることが示される。 [25]
X
{\displaystyle X}
det
(
C
)
{\displaystyle \det(C)}
スカラー分散の式 を考察し、ランダム変数とその平均の間の ユークリッド距離 の二乗 、または、単にベクトル とそれ自身のスカラー積として 再解釈すると、別の一般化が得られます。この結果は、共分散行列の トレース である となります 。
Var
(
X
)
=
E
[
(
X
−
μ
)
2
]
{\displaystyle \operatorname {Var} (X)=\operatorname {E} \left[(X-\mu )^{2}\right]}
(
X
−
μ
)
2
{\displaystyle (X-\mu )^{2}}
X
−
μ
{\displaystyle X-\mu }
E
[
(
X
−
μ
)
T
(
X
−
μ
)
]
=
tr
(
C
)
,
{\displaystyle \operatorname {E} \left[(X-\mu )^{\operatorname {T} }(X-\mu )\right]=\operatorname {tr} (C),}
参照
数学ポータル
無料辞書のウィクショナリーで 「variance」 を調べてください。
差異の種類
参考文献
^ Wasserman, Larry (2005). All of Statistics: a concise course in statistics in statistical inference . Springer texts in statistics. p. 51. ISBN 978-1-4419-2322-6 。
^ Yuli Zhang; Huaiyu Wu; Lei Cheng (2012 年 6 月). 分散と共分散に関するいくつかの新しい変形式 。モデリング、識別、制御に関する第 4 回国際会議 (ICMIC2012) の議事録。pp. 987–992。
^ Kagan, A.; Shepp, LA (1998). 「なぜ分散があるのか?」. 統計と確率の手紙 . 38 (4): 329–333. doi :10.1016/S0167-7152(98)00041-8.
^ ジョンソン、リチャード、ウィチャーン、ディーン (2001)。 応用多変量統計分析 。プレンティスホール。p. 76。ISBN 0-13-187715-1 。
^ Loève, M. (1977)「確率論」、 Graduate Texts in Mathematics 、第45巻、第4版、Springer-Verlag、p. 12。
^ ビエネイメ、I.-J. (1853) "Considérations à l'appui de la découverte de Laplace sur la loi de probabilité dans la methode des moindres carrés"、 Comptes rendus de l'Académie des Sciences Paris 、37、p. 309–317;デジタルコピーが利用可能 [1]
^ ビエネイメ、I.-J. (1867) "Considérations à l'appui de la découverte de Laplace sur la loi de probabilité dans la methode des moindres carrés"、 Journal de Mathématiques Pures et Appliquées、Série 2 、Tome 12、p. 158–167;デジタルコピーが利用可能 [2][3]
^ Cornell, JR、Benjamin, CA、 「土木技術者のための確率、統計、意思決定」、 McGraw-Hill、NY、1970年、178-9ページ。
^ Goodman, Leo A. (1960 年 12 月). 「製品の正確な分散について」. アメリカ統計学会誌 . 55 (292): 708–713. doi :10.2307/2281592. JSTOR 2281592.
^ Yuli Zhang; Huaiyu Wu; Lei Cheng (2012 年 6 月). 分散と共分散に関するいくつかの新しい変形式 。モデリング、識別、制御に関する第 4 回国際会議 (ICMIC2012) の議事録。pp. 987–992。
^ Navidi, William (2006) Statistics for Engineers and Scientists 、McGraw-Hill、p. 14。
^ Montgomery, DC および Runger, GC (1994) エンジニアのための応用統計と確率 、201 ページ。John Wiley & Sons ニューヨーク
^ Yuli Zhang; Huaiyu Wu; Lei Cheng (2012 年 6 月). 分散と共分散に関するいくつかの新しい変形式 。モデリング、識別、制御に関する第 4 回国際会議 (ICMIC2012) の議事録。pp. 987–992。
^ Knight K. (2000)、 数学統計学 、Chapman and Hall、ニューヨーク。(命題2.11)
^ Casella and Berger (2002) 統計的推論 、例7.3.3、p. 331 [ 全文引用が必要 ]
^ Mood, AM, Graybill, FA, Boes, DC (1974) 統計理論入門 、第 3 版、McGraw-Hill、ニューヨーク、p. 229
^ Kenney, John F.; Keeping, ES (1951). Mathematics of Statistics. Part Two (PDF) (第2版). Princeton, New Jersey: D. Van Nostrand Company, Inc. 2018年11月17日のオリジナル (PDF) からのアーカイブ– KrishiKosh経由。
^ Rose, Colin; Smith, Murray D. (2002). 「Mathematica による数理統計」 Springer-Verlag、ニューヨーク。
^ Weisstein, Eric W. 「サンプル分散分布」。MathWorld Wolfram。
^サミュエルソン、ポール (1968)。「 あなたはどれだけ 逸脱 できるか?」 アメリカ統計学会誌 。63 (324): 1522–1525。doi :10.1080/01621459.1968.10480944。JSTOR 2285901 。
^ Mercer, A. McD. (2000). 「一般的な方法を用いた A–G、A–H、G–H、および Ky Fan 型の不等式の族の境界」 J. Math. Anal. Appl . 243 (1): 163–173. doi : 10.1006/jmaa.1999.6688 .
^ Sharma, R. (2008). 「算術平均、調和平均、分散に関するその他の不等式」. Journal of Mathematical Inequalities . 2 (1): 109–114. CiteSeerX 10.1.1.551.9397 . doi :10.7153/jmi-02-11.
^ Fama, Eugene F.; French, Kenneth R. (2010-04-21). 「Q&A: セミバリアンス: より優れたリスク測定法か?」 Fama/French フォーラム 。
^ ロナルド・フィッシャー (1918) メンデル遺伝の仮定に基づく親族間の相関関係
^ Kocherlakota, S.; Kocherlakota, K. (2004). 「一般化分散」。 統計科学百科事典 。Wileyオンラインライブラリ。doi : 10.1002 /0471667196.ess0869。ISBN 0-471-66719-6 。