区間内の一様分布
確率論 と 統計学 において 、 連続一様分布 または 矩形分布は 対称 確率分布 の族である 。このような分布は、特定の境界の間に任意の結果がある実験を説明する。 [1] 境界はパラメータによって定義され、 最小 値と最大値である。区間は 閉じた (つまり )か 開いた (つまり )かのいずれかである。 [2] したがって、分布はしばしば と省略され 、 は 一様分布を表す。 [1] 境界間の差は区間の長さを定義する。 分布の サポート上の同じ長さの 区間 はすべて、確率が等しい。これは、分布のサポートに含まれるという制約以外の制約がない ランダム変数 の 最大エントロピー確率分布 である。 [3]
a
{\displaystyle a}
b
,
{\displaystyle b,}
[
a
,
b
]
{\displaystyle [a,b]}
(
a
,
b
)
{\displaystyle (a,b)}
U
(
a
,
b
)
,
{\displaystyle U(a,b),}
U
{\displaystyle U}
X
{\displaystyle X}
定義
確率密度関数
連続一様分布の
確率密度関数 は
f
(
x
)
=
{
1
b
−
a
for
a
≤
x
≤
b
,
0
for
x
<
a
or
x
>
b
.
{\displaystyle f(x)={\begin{cases}{\frac {1}{b-a}}&{\text{for }}a\leq x\leq b,\\[8pt]0&{\text{for }}x<a\ {\text{ or }}\ x>b.\end{cases}}}
2 つの境界 とにおける の値は、 通常重要ではありません。なぜなら、それらの値は、 任意の区間におけるの値も 、 の値も 、任意の高次のモーメントの値も変更しないからです。 ときには、それらの値はゼロに選ばれ、ときには に選ばれます。 後者は、 最大尤度 法による推定のコンテキストでは適切です。 フーリエ解析 のコンテキストでは、 または の値を とすることができます 。その理由は、その場合、この一様関数の多くの 積分変換 の逆変換によって、 関数自体が返されるからです。関数は、「 ほとんどすべての場所」、つまり 測度 がゼロの点の集合を除いて等しい関数ではありません。 また、それは 、そのような曖昧さのない
符号関数 と一致しています。
f
(
x
)
{\displaystyle f(x)}
a
{\displaystyle a}
b
{\displaystyle b}
∫
c
d
f
(
x
)
d
x
{\textstyle \int _{c}^{d}f(x)dx}
[
c
,
d
]
,
{\displaystyle [c,d],}
∫
a
b
x
f
(
x
)
d
x
,
{\textstyle \int _{a}^{b}xf(x)dx,}
1
b
−
a
.
{\displaystyle {\tfrac {1}{b-a}}.}
f
(
a
)
{\displaystyle f(a)}
f
(
b
)
{\displaystyle f(b)}
1
2
(
b
−
a
)
,
{\displaystyle {\tfrac {1}{2(b-a)}},}
任意の確率密度関数は積分すると連続一様分布の確率密度関数は 、底辺の長さと 高 さ を それぞれ 表す長方形としてグラフ化されます 。底辺の長さが長くなると、高さ(分布境界内の特定の値における密度)は減少します。 [4]
1
,
{\displaystyle 1,}
b
−
a
{\displaystyle b-a}
1
b
−
a
{\displaystyle {\tfrac {1}{b-a}}}
平均と分散の 観点から見ると 、連続一様分布の確率密度関数は
μ
{\displaystyle \mu }
σ
2
,
{\displaystyle \sigma ^{2},}
f
(
x
)
=
{
1
2
σ
3
for
−
σ
3
≤
x
−
μ
≤
σ
3
,
0
otherwise
.
{\displaystyle f(x)={\begin{cases}{\frac {1}{2\sigma {\sqrt {3}}}}&{\text{for }}-\sigma {\sqrt {3}}\leq x-\mu \leq \sigma {\sqrt {3}},\\0&{\text{otherwise}}.\end{cases}}}
累積分布関数
連続一様分布の
累積分布関数 は次のようになります。
F
(
x
)
=
{
0
for
x
<
a
,
x
−
a
b
−
a
for
a
≤
x
≤
b
,
1
for
x
>
b
.
{\displaystyle F(x)={\begin{cases}0&{\text{for }}x<a,\\[8pt]{\frac {x-a}{b-a}}&{\text{for }}a\leq x\leq b,\\[8pt]1&{\text{for }}x>b.\end{cases}}}
その逆は次のようになります。
F
−
1
(
p
)
=
a
+
p
(
b
−
a
)
for
0
<
p
<
1.
{\displaystyle F^{-1}(p)=a+p(b-a)\quad {\text{ for }}0<p<1.}
平均と分散 の観点から見ると、 連続一様分布の累積分布関数は次のようになります。
μ
{\displaystyle \mu }
σ
2
,
{\displaystyle \sigma ^{2},}
F
(
x
)
=
{
0
for
x
−
μ
<
−
σ
3
,
1
2
(
x
−
μ
σ
3
+
1
)
for
−
σ
3
≤
x
−
μ
<
σ
3
,
1
for
x
−
μ
≥
σ
3
;
{\displaystyle F(x)={\begin{cases}0&{\text{for }}x-\mu <-\sigma {\sqrt {3}},\\{\frac {1}{2}}\left({\frac {x-\mu }{\sigma {\sqrt {3}}}}+1\right)&{\text{for }}-\sigma {\sqrt {3}}\leq x-\mu <\sigma {\sqrt {3}},\\1&{\text{for }}x-\mu \geq \sigma {\sqrt {3}};\end{cases}}}
その逆は次のようになります。
F
−
1
(
p
)
=
σ
3
(
2
p
−
1
)
+
μ
for
0
≤
p
≤
1.
{\displaystyle F^{-1}(p)=\sigma {\sqrt {3}}(2p-1)+\mu \quad {\text{ for }}0\leq p\leq 1.}
ランダム変数を 見つける には
X
∼
U
(
0
,
23
)
,
{\displaystyle X\sim U(0,23),}
P
(
2
<
X
<
18
)
:
{\displaystyle P(2<X<18):}
P
(
2
<
X
<
18
)
=
(
18
−
2
)
⋅
1
23
−
0
=
16
23
.
{\displaystyle P(2<X<18)=(18-2)\cdot {\frac {1}{23-0}}={\frac {16}{23}}.}
連続一様分布関数のグラフ表現では 、確率を示す指定された境界内の曲線の下の領域は長方形になります。上記の特定の例では、底辺は 、高さは となります [5]
[
f
(
x
)
vs
x
]
,
{\displaystyle [f(x){\text{ vs }}x],}
16
,
{\displaystyle 16,}
1
23
.
{\displaystyle {\tfrac {1}{23}}.}
ランダム変数を 見つけるには
X
∼
U
(
0
,
23
)
,
{\displaystyle X\sim U(0,23),}
P
(
X
>
12
|
X
>
8
)
:
{\displaystyle P(X>12\ |\ X>8):}
P
(
X
>
12
|
X
>
8
)
=
(
23
−
12
)
⋅
1
23
−
8
=
11
15
.
{\displaystyle P(X>12\ |\ X>8)=(23-12)\cdot {\frac {1}{23-8}}={\frac {11}{15}}.}
上記の例は、 連続一様分布の 条件付き確率の例です。 が
X
>
8
{\displaystyle X>8}
真であるとすると、 と
X
>
12
?
{\displaystyle X>12?}
なる確率はどれくらいでしょうか。条件付き確率は標本空間を変えるので、新しい区間の長さ を
b
−
a
′
{\displaystyle b-a'}
計算する必要があります。ここで 、および [5] です。グラフィカルな表現は例 1 に従い、指定された境界内の曲線の下の領域が確率を表示します。長方形の底辺は で、高さは [5] です。
b
=
23
{\displaystyle b=23}
a
′
=
8.
{\displaystyle a'=8.}
11
,
{\displaystyle 11,}
1
15
.
{\displaystyle {\tfrac {1}{15}}.}
生成関数
モーメント生成関数
連続一様分布の モーメント生成関数は次 の通りである: [6]
M
X
=
E
(
e
t
X
)
=
∫
a
b
e
t
x
d
x
b
−
a
=
e
t
b
−
e
t
a
t
(
b
−
a
)
=
B
t
−
A
t
t
(
b
−
a
)
,
{\displaystyle M_{X}=\mathrm {E} (\mathrm {e} ^{tX})=\int _{a}^{b}\mathrm {e} ^{tx}{\frac {dx}{b-a}}={\frac {\mathrm {e} ^{tb}-\mathrm {e} ^{ta}}{t(b-a)}}={\frac {B^{t}-A^{t}}{t(b-a)}},}
そこから 生のモーメントを計算することができる
m
k
:
{\displaystyle m_{k}:}
m
1
=
a
+
b
2
,
{\displaystyle m_{1}={\frac {a+b}{2}},}
m
2
=
a
2
+
a
b
+
b
2
3
,
{\displaystyle m_{2}={\frac {a^{2}+ab+b^{2}}{3}},}
m
k
=
∑
i
=
0
k
a
i
b
k
−
i
k
+
1
.
{\displaystyle m_{k}={\frac {\sum _{i=0}^{k}a^{i}b^{k-i}}{k+1}}.}
連続一様分布に従うランダム変数の場合、 期待値 は 、 分散 は
m
1
=
a
+
b
2
,
{\displaystyle m_{1}={\tfrac {a+b}{2}},}
m
2
−
m
1
2
=
(
b
−
a
)
2
12
.
{\displaystyle m_{2}-m_{1}^{2}={\tfrac {(b-a)^{2}}{12}}.}
特殊なケースでは、 連続一様分布の確率密度関数は次のようになります。
a
=
−
b
,
{\displaystyle a=-b,}
f
(
x
)
=
{
1
2
b
for
−
b
≤
x
≤
b
,
0
otherwise
;
{\displaystyle f(x)={\begin{cases}{\frac {1}{2b}}&{\text{for }}-b\leq x\leq b,\\[8pt]0&{\text{otherwise}};\end{cases}}}
モーメント生成関数は次のような単純な形に簡約される。
M
X
=
sinh
b
t
b
t
.
{\displaystyle M_{X}={\frac {\sinh bt}{bt}}.}
キュムラント生成関数
n
≥
2
,
{\displaystyle n\geq 2,}
に対して、 区間 上の連続一様分布の- 番目 の キュムラントは 、 - 番目の ベルヌーイ数 です 。 [ 7]
n
{\displaystyle n}
[
−
1
2
,
1
2
]
{\displaystyle [-{\tfrac {1}{2}},{\tfrac {1}{2}}]}
B
n
n
,
{\displaystyle {\tfrac {B_{n}}{n}},}
B
n
{\displaystyle B_{n}}
n
{\displaystyle n}
パラメータおよび を持つ連続一様分布は、 標準 一様分布 と呼ばれます 。
a
=
0
{\displaystyle a=0}
b
=
1
,
{\displaystyle b=1,}
U
(
0
,
1
)
,
{\displaystyle U(0,1),}
標準一様分布の興味深い性質の1つは、 が 標準一様分布を持つ場合、 も を持つということです。この性質は、 とりわけ反対の変量 を 生成するために使用できます 。言い換えると、この性質は 反転法 として知られており、連続標準一様分布を使用して、他の任意の連続分布の 乱数 を生成することができます。 [4] が標準一様分布を持つ一様乱数、つまり を持つ場合、 は 指定 された 累積分布関数 を持つ任意の連続分布から 乱数を生成します [4]
u
1
{\displaystyle u_{1}}
1
−
u
1
.
{\displaystyle 1-u_{1}.}
u
1
{\displaystyle u_{1}}
U
(
0
,
1
)
,
{\displaystyle U(0,1),}
x
=
F
−
1
(
u
1
)
{\displaystyle x=F^{-1}(u_{1})}
x
{\displaystyle x}
F
.
{\displaystyle F.}
他の機能との関係
遷移点で同じ規則に従う限り、連続一様分布の確率密度関数は、 ヘヴィサイドのステップ関数 で次のように表現することもできます。
f
(
x
)
=
H
(
x
−
a
)
−
H
(
x
−
b
)
b
−
a
,
{\displaystyle f(x)={\frac {\operatorname {H} (x-a)-\operatorname {H} (x-b)}{b-a}},}
または長方形関数 では 次のようになります。
f
(
x
)
=
1
b
−
a
rect
(
x
−
a
+
b
2
b
−
a
)
.
{\displaystyle f(x)={\frac {1}{b-a}}\ \operatorname {rect} \left({\frac {x-{\frac {a+b}{2}}}{b-a}}\right).}
符号関数 の遷移点には曖昧さはありません 。遷移点における半最大値規則を使用すると、連続一様分布は符号関数で次のように表すことができます。
f
(
x
)
=
sgn
(
x
−
a
)
−
sgn
(
x
−
b
)
2
(
b
−
a
)
.
{\displaystyle f(x)={\frac {\operatorname {sgn} {(x-a)}-\operatorname {sgn} {(x-b)}}{2(b-a)}}.}
プロパティ
瞬間
連続一様分布の
平均(第1原 モーメント)は次のようになります。
E
(
X
)
=
∫
a
b
x
d
x
b
−
a
=
b
2
−
a
2
2
(
b
−
a
)
=
b
+
a
2
.
{\displaystyle E(X)=\int _{a}^{b}x{\frac {dx}{b-a}}={\frac {b^{2}-a^{2}}{2(b-a)}}={\frac {b+a}{2}}.}
この分布の 2 番目の生のモーメントは次のとおりです。
E
(
X
2
)
=
∫
a
b
x
2
d
x
b
−
a
=
b
3
−
a
3
3
(
b
−
a
)
.
{\displaystyle E(X^{2})=\int _{a}^{b}x^{2}{\frac {dx}{b-a}}={\frac {b^{3}-a^{3}}{3(b-a)}}.}
一般に、 この分布の - 次生モーメントは次のようになります。
n
{\displaystyle n}
E
(
X
n
)
=
∫
a
b
x
n
d
x
b
−
a
=
b
n
+
1
−
a
n
+
1
(
n
+
1
)
(
b
−
a
)
.
{\displaystyle E(X^{n})=\int _{a}^{b}x^{n}{\frac {dx}{b-a}}={\frac {b^{n+1}-a^{n+1}}{(n+1)(b-a)}}.}
この分布の
分散(第2 中心モーメント)は次のようになります。
V
(
X
)
=
E
(
(
X
−
E
(
X
)
)
2
)
=
∫
a
b
(
x
−
a
+
b
2
)
2
d
x
b
−
a
=
(
b
−
a
)
2
12
.
{\displaystyle V(X)=E\left({\big (}X-E(X){\big )}^{2}\right)=\int _{a}^{b}\left(x-{\frac {a+b}{2}}\right)^{2}{\frac {dx}{b-a}}={\frac {(b-a)^{2}}{12}}.}
注文統計
を からの iid サンプルと し 、をこのサンプルからの 番目の 順序 統計量 とします 。
X
1
,
.
.
.
,
X
n
{\displaystyle X_{1},...,X_{n}}
U
(
0
,
1
)
,
{\displaystyle U(0,1),}
X
(
k
)
{\displaystyle X_{(k)}}
k
{\displaystyle k}
X
(
k
)
{\displaystyle X_{(k)}}
はベータ分布 、パラメータ および を 持つ 。
k
{\displaystyle k}
n
−
k
+
1.
{\displaystyle n-k+1.}
期待値は次のとおりです。
E
(
X
(
k
)
)
=
k
n
+
1
.
{\displaystyle \operatorname {E} \left(X_{(k)}\right)={k \over n+1}.}
この事実はQ-Qプロットを 作成するときに役立ちます 。
差異は次のとおりです。
V
(
X
(
k
)
)
=
k
(
n
−
k
+
1
)
(
n
+
1
)
2
(
n
+
2
)
.
{\displaystyle \operatorname {V} \left(X_{(k)}\right)={k(n-k+1) \over (n+1)^{2}(n+2)}.}
連続的に一様分布するランダム変数が固定長の任意の区間内に入る確率は、区間が分布のサポート内に含まれている限り、区間自体の位置とは無関係です(ただし、区間のサイズには依存します )。
(
ℓ
)
{\displaystyle (\ell )}
実際、 が の部分区間で が固定されている場合、 次のようになります。
X
∼
U
(
a
,
b
)
{\displaystyle X\sim U(a,b)}
[
x
,
x
+
ℓ
]
{\displaystyle [x,x+\ell ]}
[
a
,
b
]
{\displaystyle [a,b]}
ℓ
>
0
,
{\displaystyle \ell >0,}
P
(
X
∈
[
x
,
x
+
ℓ
]
)
=
∫
x
x
+
ℓ
d
y
b
−
a
=
ℓ
b
−
a
,
{\displaystyle P{\big (}X\in [x,x+\ell ]{\big )}=\int _{x}^{x+\ell }{\frac {dy}{b-a}}={\frac {\ell }{b-a}},}
これは独立しており、 この事実がこの分布の名前の由来となっています。
x
.
{\displaystyle x.}
一様分布は、区間よりも一般的な集合に一般化できます。
正式には、 を 正の有限 ルベーグ測度の ボレル集合 とします。すなわち、 上の一様分布は、 確率密度関数を外側ではゼロ 、外側では常に と等しいと定義することで指定できます 。
S
{\displaystyle S}
λ
(
S
)
,
{\displaystyle \lambda (S),}
0
<
λ
(
S
)
<
+
∞
.
{\displaystyle 0<\lambda (S)<+\infty .}
S
{\displaystyle S}
S
{\displaystyle S}
1
λ
(
S
)
{\displaystyle {\tfrac {1}{\lambda (S)}}}
S
.
{\displaystyle S.}
興味深い特殊なケースとして、集合 S が 単体 である場合が挙げられます。標準的な n 頂点単体上で一様分布を次の方法で得ることができます。 [8] : Thm.4.1 同じ 指数分布に従う n 個の独立したランダム変数を取り 、それらを X 1 ,...,X n で表し、Y i := X i / (sum i X i ) とします。すると、ベクトル Y 1 ,...,Y n は単体上で一様分布します。
Xが 標準一様分布に従う 場合、 逆変換サンプリング 法によれば、 Y = − λ −1 ln( X )は (速度)パラメータ λを持つ 指数分布 に従います。
Xが 標準一様分布に従う 場合、 Y = X nは パラメータ(1/ n ,1)を持つ ベータ分布 に従います 。したがって、
アーウィン ・ホール分布は、 n個の i.id U (0,1)分布の合計です 。
ベイツ分布は 、 n i.id U (0,1)分布の平均です 。
標準一様分布は、パラメータが(1,1)である ベータ分布 の特殊なケースです。
2 つの独立した一様分布 U 1 (a,b)+ U 2 (c,d)の和は、サポート [a+c,b+d] 上の平均に関して対称な 台形分布になります。台形分布の幅は、 U 1 と U 2 の幅の絶対差に等しくなります 。傾斜部分の幅は、最も狭い一様分布の幅に対応します。
一様分布の幅 w が同じ場合、結果は サポート [a+c,a+c+2w] 上の、平均を中心に対称な 三角分布になります。
2つの独立した均等に分布する一様分布 U 1 (a,b)+ U 2 (a,b)の合計は、サポート[2a,2b]上の 対称 三角分布を生成します。
2つのiid 一様ランダム変数 間の距離| U 1 (a,b)- U 2 (a,b)|も、対称ではありませんが、サポート[0,ba]上で 三角分布を 持ちます。
統計的推論
パラメータの推定
最大値の推定
最小分散不偏推定量
上の一様分布が 未知数である場合、最大値の 最小 分散不偏推定量 (UMVUE) は 次のようになります。
[
0
,
b
]
{\displaystyle [0,b]}
b
,
{\displaystyle b,}
b
^
UMVU
=
k
+
1
k
m
=
m
+
m
k
,
{\displaystyle {\hat {b}}_{\text{UMVU}}={\frac {k+1}{k}}m=m+{\frac {m}{k}},}
ここで、 は サンプルの最大値 、 は サンプル サイズ で、非復元サンプリングです (ただし、連続分布の場合は、この区別はほとんど確実に違いを生じません)。これは、 離散分布 の推定 と同じ理由で得られ、 最大間隔推定の非常に単純なケースと見なすことができます。この問題は、最大値推定を 第二次世界大戦 中のドイツの戦車生産の推定に適用したことから、 一般に ドイツの戦車問題 として知られています。
m
{\displaystyle m}
k
{\displaystyle k}
モーメント推定法
モーメント推定法は次 のとおり です。
b
^
M
M
=
2
X
¯
,
{\displaystyle {\hat {b}}_{MM}=2{\bar {X}},}
サンプル平均は
どこですか。
X
¯
{\displaystyle {\bar {X}}}
最大尤度推定量
最大 尤度 推定量は次のとおりです。
b
^
M
L
=
m
,
{\displaystyle {\hat {b}}_{ML}=m,}
ここで、は サンプルの最大値 であり、 サンプルの
最大 順序統計量 とも呼ばれます。
m
{\displaystyle m}
m
=
X
(
n
)
,
{\displaystyle m=X_{(n)},}
最小値の推定
a が 未知で 上の一様分布が与えられた場合、 a の最大尤度推定値は 次のようになります。
[
a
,
b
]
{\displaystyle [a,b]}
a
^
M
L
=
min
{
X
1
,
…
,
X
n
}
{\displaystyle {\hat {a}}_{ML}=\min\{X_{1},\dots ,X_{n}\}}
、
サンプル 最小値 [ 9]
中間点の推定
分布の中点は、 一様分布の平均と中央値の両方です。標本平均と標本中央値はどちらも中点の 不偏推定値 ですが、どちらも 標本 中範囲、つまり標本最大値と標本最小値の算術平均ほど 効率的で はありません。これは中点の UMVU 推定値 (および 最大尤度推定値 ) です。
a
+
b
2
,
{\displaystyle {\tfrac {a+b}{2}},}
信頼区間
最大限に
が母集団の最大値で ある サンプルであると する。すると、 ルベーグ・ボレル密度 [10]は
X
1
,
X
2
,
X
3
,
.
.
.
,
X
n
{\displaystyle X_{1},X_{2},X_{3},...,X_{n}}
U
[
0
,
L
]
,
{\displaystyle U_{[0,L]},}
L
{\displaystyle L}
X
(
n
)
=
max
(
X
1
,
X
2
,
X
3
,
.
.
.
,
X
n
)
{\displaystyle X_{(n)}=\max(X_{1},X_{2},X_{3},...,X_{n})}
f
=
d
Pr
X
(
n
)
d
λ
:
{\displaystyle f={\frac {d\Pr _{X_{(n)}}}{d\lambda }}:}
f
(
t
)
=
n
1
L
(
t
L
)
n
−
1
=
n
t
n
−
1
L
n
1
1
[
0
,
L
]
(
t
)
,
{\displaystyle f(t)=n{\frac {1}{L}}\left({\frac {t}{L}}\right)^{n-1}\!=n{\frac {t^{n-1}}{L^{n}}}1\!\!1_{[0,L]}(t),}
の 指標関数は どこに あるか
1
1
[
0
,
L
]
{\displaystyle 1\!\!1_{[0,L]}}
[
0
,
L
]
.
{\displaystyle [0,L].}
前述の信頼区間は数学的に正しくありません。
Pr
(
[
θ
^
,
θ
^
+
ε
]
∋
θ
)
≥
1
−
α
{\displaystyle \Pr {\big (}[{\hat {\theta }},{\hat {\theta }}+\varepsilon ]\ni \theta {\big )}\geq 1-\alpha }
を知らない と解くことはできない 。しかし、
ε
{\displaystyle \varepsilon }
θ
{\displaystyle \theta }
Pr
(
[
θ
^
,
θ
^
(
1
+
ε
)
]
∋
θ
)
≥
1
−
α
{\displaystyle \Pr {\big (}[{\hat {\theta }},{\hat {\theta }}(1+\varepsilon )]\ni \theta {\big )}\geq 1-\alpha }
不明だが有効 な
ε
≥
(
1
−
α
)
−
1
/
n
−
1
{\displaystyle \varepsilon \geq (1-\alpha )^{-1/n}-1}
θ
;
{\displaystyle \theta ;}
次に、上記の条件を満たす 最小のものを選択する。区間の長さはランダム変数に依存することに注意する。
ε
{\displaystyle \varepsilon }
θ
^
.
{\displaystyle {\hat {\theta }}.}
発生と応用
一様分布関数の確率は、関数形式が単純なため、計算が簡単です。 [2] そのため、この分布は、仮説検定、ランダムサンプリング、金融など、以下に示すようにさまざまな用途に使用できます。さらに、一般に、物理的な起源の実験は一様分布に従います(例:放射性 粒子 の放出)。 [1] ただし、どの用途でも、固定長の間隔に入る確率は一定であるという不変の仮定があることに注意することが重要です。 [2]
経済学の分野では、通常、 需要 と補充は予想される正規分布に従わない場合があります。その結果、 ベルヌーイ過程 などの他の分布モデルが、確率と傾向をより正確に予測するために使用されます。 [11] しかし、Wanke(2008)によると、完全に新しい製品を分析している ライフサイクル の開始時に在庫管理の リードタイム を調査するという特定のケースでは、一様分布の方が有用であることが証明されています。 [11] この状況では、新製品に関する既存のデータがないか、需要履歴が利用できないため、実際には適切または既知の分布がないため、他の分布は実行できない可能性があります。 [11] この状況では、リードタイム(需要に関連)のランダム変数は新製品については不明ですが、結果は2つの値の妥当な範囲内になる可能性が高いため、一様分布が理想的です。 [11] したがって、リード タイムはランダム変数を表します。一様分布モデルから、サイクルサービスレベルやサイクルあたりの不足など、 リードタイム に関連する他の要因を 計算することができました。また、計算の単純さから均一分布も使用されたことも指摘されている。 [11]
任意の分布からのサンプリング
一様分布は、任意の分布からサンプリングする場合に便利です。一般的な方法は、ターゲットランダム変数の 累積分布関数 (CDF) を使用する逆変換サンプリング法です。この方法は理論的な作業で非常に便利です。この方法を使用したシミュレーションでは、ターゲット変数の CDF を反転する必要があるため、CDF が閉じた形式でわからない場合のために、代替方法が考案されています。そのような方法の 1 つが 拒否サンプリング です。
正規 分布は 、逆変換法が効率的でない重要な例です。ただし、 逆変換を使用して 2 つの独立した一様 ランダム変数を2 つの独立した 正規分布 ランダム変数
に変換する正確な方法である Box–Muller 変換があります。
量子化誤差
アナログからデジタルへの変換では、量子化誤差が発生します。この誤差は、丸めまたは切り捨てによって発生します。元の信号が 1 つの 最下位ビット (LSB) よりもはるかに大きい場合、量子化誤差は信号と有意な相関関係になく、ほぼ均一に分布します。したがって、 RMS 誤差は この分布の分散から生じます。
ランダム変量生成
シミュレーション実験を実行すると便利なアプリケーションは数多くあります。多くの プログラミング言語 には、標準の一様分布に従って効果的に分布する
疑似乱数 を生成する実装が付属しています。
一方、均一に分布した数は、 非均一な乱数変量生成 の基礎としてよく使用されます。
が標準一様分布からサンプリングされた値である 場合、その値は 上記のように
および によってパラメータ化された一様分布に従います。
u
{\displaystyle u}
a
+
(
b
−
a
)
u
{\displaystyle a+(b-a)u}
a
{\displaystyle a}
b
,
{\displaystyle b,}
歴史
一様分布の概念の歴史的起源は確定的ではないが、「一様」という用語はサイコロゲームにおける 等確率 の概念から生じたと推測されている(サイコロゲームは 連続一様標本空間ではなく 離散的であることに注意)。 等確率は ジェロラモ・カルダノ の 『Liber de Ludo Aleae』 で言及されている。これは 16世紀に書かれたマニュアルで、サイコロに関する高度な確率計算について詳述されている。 [12]
参照
参考文献
^ abc デッキング、ミシェル(2005)。 確率と統計の現代的入門:なぜ、どのように理解するか 。ロンドン、英国:シュプリンガー。pp. 60–61。ISBN 978-1-85233-896-1 。
^ abc ウォルポール、ロナルド、他 (2012)。 エンジニアと科学者のための確率と統計 。米国ボストン:プレンティスホール。pp. 171–172。ISBN 978-0-321-62911-1 。
^ Park, Sung Y.; Bera, Anil K. (2009). 「最大エントロピー自己回帰条件付き異分散モデル」. Journal of Econometrics . 150 (2): 219–230. CiteSeerX 10.1.1.511.9750 . doi :10.1016/j.jeconom.2008.12.014.
^ abc 「一様分布(連続)」。 MathWorks 。2019年。 2019年 11月22日 閲覧 。
^ abc Illowsky, Barbara; et al. (2013). 統計学入門. ライス大学、テキサス州ヒューストン、米国: OpenStax College. pp. 296–304. ISBN 978-1-938168-20-8 。
^ カセラ&バーガー 2001、626 ページ
^ ウィチュラ、マイケル J. (2001 年 1 月 11 日)。 「キュムラント」 (PDF) 。 Stat 304 配布資料 。シカゴ大学。
^ 非均一ランダム変量生成。doi :10.1007 / 978-1-4613-8643-8。
^
。 因数 は 、 によって 制限される最大の a によって最大化されます 。したがって、 は の最大値です 。
L
n
(
a
,
b
)
=
∏
i
=
1
n
f
(
X
i
)
=
1
(
b
−
a
)
n
1
[
a
,
b
]
(
X
1
,
…
,
X
n
)
{\displaystyle L_{n}(a,b)=\prod _{i=1}^{n}f(X_{i})={\frac {1}{(b-a)^{n}}}\mathbf {1} _{[a,b]}(X_{1},\dots ,X_{n})}
=
1
(
b
−
a
)
n
1
{
a
≤
min
{
X
1
,
…
,
X
n
}
}
1
{
max
{
X
1
,
…
,
X
n
}
≤
b
}
{\displaystyle ={\frac {1}{(b-a)^{n}}}\mathbf {1} _{\{a\leq \min\{X_{1},\dots ,X_{n}\}\}}\mathbf {1} _{\{\max\{X_{1},\dots ,X_{n}\}\leq b\}}}
n
≥
1
{\displaystyle n\geq 1}
1
(
b
−
a
)
n
{\displaystyle {\frac {1}{(b-a)^{n}}}}
L
n
(
a
,
b
)
{\displaystyle L_{n}(a,b)}
min
{
X
1
,
…
,
X
n
}
{\displaystyle \min\{X_{1},\dots ,X_{n}\}}
a
=
min
{
X
1
,
…
,
X
n
}
{\displaystyle a=\min\{X_{1},\dots ,X_{n}\}}
L
n
(
a
,
b
)
{\displaystyle L_{n}(a,b)}
^ Nechval KN、Nechval NA、Vasamanis EK、Makeev VY (2002) 最短の長さの信頼区間の構築。運輸および電気通信 3 (1) 95-103
^ abcde Wanke, Peter (2008). 「新製品在庫管理への最初の実践的アプローチとしての均一分布」 International Journal of Production Economics . 114 (2): 811–819. doi :10.1016/j.ijpe.2008.04.004 – Research Gate経由。
^ デヴィッド・ベルハウス (2005 年 5 月)。 「カルダノのLiber de Ludoを解読する」。 ヒストリア マセマティカ 。 32 :180-202。 土井 : 10.1016/j.hm.2004.04.001 。
さらに読む
カセラ、ジョージ、ロジャー・L・バーガー(2001)、 統計的推論 (第2版)、トムソンラーニング、 ISBN 978-0-534-24312-8 、 LCCN 2001025794
外部リンク