確率分布
確率論 において 、 対数正規 分布 (または 対数正規分布 )は、 対数 が 正規分布する ランダム 変数 の 連続 確率分布 です。したがって、ランダム変数 X が対数正規分布する場合、 Y = ln( X ) は正規分布します。 [2] [3] 同様に、 Y が正規分布する場合、 Y の 指数関数 X = exp( Y ) は 対数正規分布します。対数正規分布するランダム変数は、正の実数値のみを取ります。これは、精密科学や工学だけでなく、医学、経済学、その他のトピック(エネルギー、濃度、長さ、金融商品の価格、その他の測定基準など)の測定にも 便利 で 有用 な モデル です 。
この分布はフランシス・ゴルトン にちなんで ゴルトン分布 または ゴルトンの分布 と呼ばれることもあります 。 [4] 対数正規分布は、 マカリスター分布 、 ジブラット分布 、 コブ・ダグラス 分布など、他の名前でも呼ばれています。 [4]
対数正規過程は、それぞれが正である 多数の 独立した ランダム変数 の乗法積の統計的実現である 。 これは、対数領域の 中心極限定理( ジブラの法則 と呼ばれることもある)を考慮することによって正当化される。対数正規分布は、 ランダム変数 Xの 最大エントロピー確率分布であり、 ln( X ) の平均と分散 が指定されている。 [5]
定義
生成とパラメータ
を標準正規変数 と し 、 と を 2つの実数( )とします 。このとき、確率変数の分布は
Z
{\displaystyle \ Z\ }
μ
{\displaystyle \mu }
σ
{\displaystyle \sigma }
σ
>
0
{\displaystyle \sigma >0}
X
=
e
μ
+
σ
Z
{\displaystyle X=e^{\mu +\sigma Z}}
は、パラメータおよび を持つ対数正規分布と呼ばれます 。これらは、変数の自然対数 の期待値 (または平均 ) と標準偏差であり 、 自体 の 期待 値 と 標準 偏差 で は ありません 。
μ
{\displaystyle \mu }
σ
{\displaystyle \sigma }
ln
(
X
)
{\displaystyle \ln(X)}
X
{\displaystyle \ X\ }
正規分布と対数正規分布の関係。 が正規分布している場合、 は 対数正規分布します。
Y
=
μ
+
σ
Z
{\displaystyle \ Y=\mu +\sigma Z\ }
X
∼
e
Y
{\displaystyle \ X\sim e^{Y}\ }
この関係は、対数関数または指数関数の底に関係なく成り立ちます。 が正規分布している場合、任意の 2 つの正の数に対して も が正規分布します。 同様に、 が対数正規分布している場合、 も正規分布し ます 。
log
a
(
X
)
{\displaystyle \ \log _{a}(X)\ }
log
b
(
X
)
{\displaystyle \ \log _{b}(X)\ }
a
,
b
≠
1
.
{\displaystyle \ a,b\neq 1~.}
e
Y
{\displaystyle \ e^{Y}\ }
a
Y
,
{\displaystyle \ a^{Y}\ ,}
0
<
a
≠
1
{\displaystyle 0<a\neq 1}
望ましい平均と分散を持つ 分布
を生成するには 、
μ
X
{\displaystyle \mu _{X}}
σ
X
2
,
{\displaystyle \ \sigma _{X}^{2}\ ,}
μ
=
ln
(
μ
X
2
μ
X
2
+
σ
X
2
)
{\displaystyle \ \mu =\ln \left({\frac {\mu _{X}^{2}}{\ {\sqrt {\mu _{X}^{2}+\sigma _{X}^{2}\ }}\ }}\right)\ }
σ
2
=
ln
(
1
+
σ
X
2
μ
X
2
)
.
{\displaystyle \ \sigma ^{2}=\ln \left(1+{\frac {\ \sigma _{X}^{2}\ }{\mu _{X}^{2}}}\right)~.}
あるいは、「乗法」または「幾何」パラメータ とを 使用することもできます。これらはより直接的な解釈が可能です。 は 分布の 中央値 であり、 「散布」間隔を決定するのに役立ちます (以下を参照)。
μ
∗
=
e
μ
{\displaystyle \ \mu ^{*}=e^{\mu }\ }
σ
∗
=
e
σ
{\displaystyle \ \sigma ^{*}=e^{\sigma }\ }
μ
∗
{\displaystyle \ \mu ^{*}\ }
σ
∗
{\displaystyle \ \sigma ^{*}\ }
確率密度関数
正の確率変数は 、 の自然対数が 平均 と分散で正規分布している場合 、対数正規分布(つまり )します。
X
{\displaystyle \ X\ }
X
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle \ X\sim \operatorname {Lognormal} \left(\ \mu ,\sigma ^{2}\ \right)\ }
X
{\displaystyle \ X\ }
μ
{\displaystyle \mu }
σ
2
:
{\displaystyle \ \sigma ^{2}\ :}
ln
(
X
)
∼
N
(
μ
,
σ
2
)
{\displaystyle \ln(X)\sim {\mathcal {N}}(\mu ,\sigma ^{2})}
をそれぞれ累積確率分布関数と標準正規分布の確率密度関数とする と 、 [ 2] [4] 対数正規分布の確率密度関数は次のように表さ
れる 。
Φ
{\displaystyle \ \Phi \ }
φ
{\displaystyle \ \varphi \ }
N
(
0
,
1
)
{\displaystyle \ {\mathcal {N}}(\ 0,1\ )\ }
f
X
(
x
)
=
d
d
x
P
X
[
X
≤
x
]
=
d
d
x
P
X
[
ln
X
≤
ln
x
]
=
d
d
x
Φ
(
ln
x
−
μ
σ
)
=
φ
(
ln
x
−
μ
σ
)
d
d
x
(
ln
x
−
μ
σ
)
=
φ
(
ln
x
−
μ
σ
)
1
σ
x
=
1
x
σ
2
π
exp
(
−
(
ln
x
−
μ
)
2
2
σ
2
)
.
{\displaystyle {\begin{aligned}f_{X}(x)&={\frac {\rm {d}}{{\rm {d}}x}}\ \operatorname {\mathbb {P} _{\mathit {X}}} \,\!{\bigl [}\ X\leq x\ {\bigr ]}\\[6pt]&={\frac {\rm {d}}{{\rm {d}}x}}\ \operatorname {\mathbb {P} _{\mathit {X}}} \,\!{\bigl [}\ \ln X\leq \ln x\ {\bigr ]}\\[6pt]&={\frac {\rm {d}}{{\rm {d}}x}}\operatorname {\Phi } \!\!\left({\frac {\ \ln x-\mu \ }{\sigma }}\right)\\[6pt]&=\operatorname {\varphi } \!\left({\frac {\ln x-\mu }{\sigma }}\right){\frac {\rm {d}}{{\rm {d}}x}}\left({\frac {\ \ln x-\mu \ }{\sigma }}\right)\\[6pt]&=\operatorname {\varphi } \!\left({\frac {\ \ln x-\mu \ }{\sigma }}\right){\frac {1}{\ \sigma \ x\ }}\\[6pt]&={\frac {1}{\ x\ \sigma {\sqrt {2\ \pi \ }}\ }}\exp \left(-{\frac {\ (\ln x-\mu )^{2}\ }{2\ \sigma ^{2}}}\right)~.\end{aligned}}}
累積分布関数
累積分布 関数 は
F
X
(
x
)
=
Φ
(
(
ln
x
)
−
μ
σ
)
{\displaystyle F_{X}(x)=\Phi \left({\frac {(\ln x)-\mu }{\sigma }}\right)}
ここで、は 標準正規分布の累積分布関数です(つまり、 )。
Φ
{\displaystyle \ \Phi \ }
N
(
0
,
1
)
{\displaystyle \ \operatorname {\mathcal {N}} (\ 0,\ 1)\ }
これは次のように表現することもできる。 [2]
1
2
[
1
+
erf
(
ln
x
−
μ
σ
2
)
]
=
1
2
erfc
(
−
ln
x
−
μ
σ
2
)
{\displaystyle {\frac {1}{2}}\left[1+\operatorname {erf} \left({\frac {\ln x-\mu }{\sigma {\sqrt {2}}}}\right)\right]={\frac {1}{2}}\operatorname {erfc} \left(-{\frac {\ln x-\mu }{\sigma {\sqrt {2}}}}\right)}
ここで、 erfcは 相補誤差関数 です 。
多変量対数正規分布
が多変量正規分布 である 場合 、 は 多変量対数正規分布に従う。 [6] [7] 指数関数はランダムベクトルに要素ごとに適用されます 。 の平均 は
X
∼
N
(
μ
,
Σ
)
{\displaystyle {\boldsymbol {X}}\sim {\mathcal {N}}({\boldsymbol {\mu }},\,{\boldsymbol {\Sigma }})}
Y
i
=
exp
(
X
i
)
{\displaystyle Y_{i}=\exp(X_{i})}
X
{\displaystyle {\boldsymbol {X}}}
Y
{\displaystyle {\boldsymbol {Y}}}
E
[
Y
]
i
=
e
μ
i
+
1
2
Σ
i
i
,
{\displaystyle \operatorname {E} [{\boldsymbol {Y}}]_{i}=e^{\mu _{i}+{\frac {1}{2}}\Sigma _{ii}},}
そしてその 共分散行列 は
Var
[
Y
]
i
j
=
e
μ
i
+
μ
j
+
1
2
(
Σ
i
i
+
Σ
j
j
)
(
e
Σ
i
j
−
1
)
.
{\displaystyle \operatorname {Var} [{\boldsymbol {Y}}]_{ij}=e^{\mu _{i}+\mu _{j}+{\frac {1}{2}}(\Sigma _{ii}+\Sigma _{jj})}(e^{\Sigma _{ij}}-1).}
多変量対数正規分布は広く使用されていないため、このエントリの残りの部分では 単変量分布 のみを扱います。
特性関数とモーメント生成関数
対数正規分布のすべてのモーメントが存在し、
E
[
X
n
]
=
e
n
μ
+
n
2
σ
2
/
2
{\displaystyle \operatorname {E} [X^{n}]=e^{n\mu +n^{2}\sigma ^{2}/2}}
これは、積分内に を入れることによって導くことができます 。しかし、対数正規分布はモーメントによって決まるわけではありません。 [8] これは、ゼロの近傍に定義されたモーメント生成関数を持つことができないことを意味します。 [9] 実際、定義積分が発散するため、
引数の任意の正の値に対して 期待値 は定義されません。
z
=
ln
(
x
)
−
(
μ
+
n
σ
2
)
σ
{\displaystyle z={\tfrac {\ln(x)-(\mu +n\sigma ^{2})}{\sigma }}}
E
[
e
t
X
]
{\displaystyle \operatorname {E} [e^{tX}]}
t
{\displaystyle t}
特性 関数は t の実数値に対して定義されるが、 負の虚数部を持つ t の複素数値に対しては定義されないため、特性関数は原点において 解析的で はない。したがって、対数正規分布の特性関数は無限収束級数として表すことはできない。 [10] 特に、そのテイラー 形式級数 は発散する。
E
[
e
i
t
X
]
{\displaystyle \operatorname {E} [e^{itX}]}
∑
n
=
0
∞
(
i
t
)
n
n
!
e
n
μ
+
n
2
σ
2
/
2
{\displaystyle \sum _{n=0}^{\infty }{\frac {(it)^{n}}{n!}}e^{n\mu +n^{2}\sigma ^{2}/2}}
しかしながら、いくつかの代替的な 発散級数 表現が得られている。 [10] [11] [12] [13]
収束領域内の特性関数 の閉形式の公式は 知られていない。比較的単純な近似式が閉形式で利用可能であり、 [14]で与えられる。
φ
(
t
)
{\displaystyle \varphi (t)}
t
{\displaystyle t}
φ
(
t
)
≈
exp
(
−
W
2
(
−
i
t
σ
2
e
μ
)
+
2
W
(
−
i
t
σ
2
e
μ
)
2
σ
2
)
1
+
W
(
−
i
t
σ
2
e
μ
)
{\displaystyle \varphi (t)\approx {\frac {\exp \left(-{\frac {W^{2}(-it\sigma ^{2}e^{\mu })+2W(-it\sigma ^{2}e^{\mu })}{2\sigma ^{2}}}\right)}{\sqrt {1+W(-it\sigma ^{2}e^{\mu })}}}}
ここで、 は Lambert W 関数 です 。この近似は漸近法によって導出されますが、 の収束領域全体にわたって鋭いままです 。
W
{\displaystyle W}
φ
{\displaystyle \varphi }
プロパティ
a. は対数正規変数で、 です 。は、 正規変数 に変換し 、その密度を (青い領域) で定義される領域で積分することで計算されます 。これには、レイトレーシングの数値法が使用されます。 [15] b & c. 対数正規変数の関数の pdf と cdf も この方法で計算できます。
y
{\displaystyle y}
μ
=
1
,
σ
=
0.5
{\displaystyle \mu =1,\sigma =0.5}
p
(
sin
y
>
0
)
{\displaystyle p(\sin y>0)}
x
=
ln
y
{\displaystyle x=\ln y}
sin
e
x
>
0
{\displaystyle \sin e^{x}>0}
sin
y
{\displaystyle \sin y}
さまざまな領域における確率
任意の領域における対数正規分布の確率内容は、まず変数を正規分布に変換し、次にレイトレース法を使用して数値積分することで、必要な精度で計算できます。 [15] (Matlabコード)
対数正規変数の関数の確率
対数正規分布の確率は任意の領域で計算できるため、対数正規分布の任意の関数の累積分布関数(および結果として累積分布関数と逆累積分布関数)も計算できることを意味します。 [15] (Matlabコード)
幾何モーメントまたは乗法モーメント
対数正規分布の幾何平均または乗法平均は です 。 これ は中央値に等しくなります。 幾何標準偏差または乗法標準偏差 は です 。 [16] [17]
GM
[
X
]
=
e
μ
=
μ
∗
{\displaystyle \operatorname {GM} [X]=e^{\mu }=\mu ^{*}}
GSD
[
X
]
=
e
σ
=
σ
∗
{\displaystyle \operatorname {GSD} [X]=e^{\sigma }=\sigma ^{*}}
算術統計学との類推により、幾何分散を定義することができ 、 幾何変動係数 [16] が提案されている。この用語は、対数正規データにおける乗法変動を記述するための変動係数に 類似したもの として意図されていたが、 このGCVの定義には、それ自体の推定値としての理論的根拠がない ( 変動係数 も参照)。
GVar
[
X
]
=
e
σ
2
{\displaystyle \operatorname {GVar} [X]=e^{\sigma ^{2}}}
GCV
[
X
]
=
e
σ
−
1
{\displaystyle \operatorname {GCV} [X]=e^{\sigma }-1}
CV
{\displaystyle \operatorname {CV} }
幾何平均は算術平均よりも小さいことに注意してください。これは AM-GM不等式によるもので、対数が 凹関数 であることの結果です 。実際、
E
[
X
]
=
e
μ
+
1
2
σ
2
=
e
μ
⋅
e
σ
2
=
GM
[
X
]
⋅
GVar
[
X
]
.
{\displaystyle \operatorname {E} [X]=e^{\mu +{\frac {1}{2}}\sigma ^{2}}=e^{\mu }\cdot {\sqrt {e^{\sigma ^{2}}}}=\operatorname {GM} [X]\cdot {\sqrt {\operatorname {GVar} [X]}}.}
[18]
金融においては、この用語は 凸性補正 として解釈されることがある。 確率論 の観点から見ると、これは 幾何ブラウン運動の伊藤の補題 における補正項と同じである 。
e
−
1
2
σ
2
{\displaystyle e^{-{\frac {1}{2}}\sigma ^{2}}}
算術モーメント
任意の実数または複素数 n に対して、 対数正規分布に従う変数 Xの n 次 モーメントは [4] で与えられる。
E
[
X
n
]
=
e
n
μ
+
1
2
n
2
σ
2
.
{\displaystyle \operatorname {E} [X^{n}]=e^{n\mu +{\frac {1}{2}}n^{2}\sigma ^{2}}.}
具体的には、対数正規分布変数 X の算術平均、期待二乗、算術分散、算術標準偏差はそれぞれ次のように表される: [2]
E
[
X
]
=
e
μ
+
1
2
σ
2
,
E
[
X
2
]
=
e
2
μ
+
2
σ
2
,
Var
[
X
]
=
E
[
X
2
]
−
E
[
X
]
2
=
(
E
[
X
]
)
2
(
e
σ
2
−
1
)
=
e
2
μ
+
σ
2
(
e
σ
2
−
1
)
,
SD
[
X
]
=
Var
[
X
]
=
E
[
X
]
e
σ
2
−
1
=
e
μ
+
1
2
σ
2
e
σ
2
−
1
,
{\displaystyle {\begin{aligned}\operatorname {E} [X]&=e^{\mu +{\tfrac {1}{2}}\sigma ^{2}},\\[4pt]\operatorname {E} [X^{2}]&=e^{2\mu +2\sigma ^{2}},\\[4pt]\operatorname {Var} [X]&=\operatorname {E} [X^{2}]-\operatorname {E} [X]^{2}=(\operatorname {E} [X])^{2}(e^{\sigma ^{2}}-1)=e^{2\mu +\sigma ^{2}}(e^{\sigma ^{2}}-1),\\[4pt]\operatorname {SD} [X]&={\sqrt {\operatorname {Var} [X]}}=\operatorname {E} [X]{\sqrt {e^{\sigma ^{2}}-1}}=e^{\mu +{\tfrac {1}{2}}\sigma ^{2}}{\sqrt {e^{\sigma ^{2}}-1}},\end{aligned}}}
算術 変動係数は 比である 。対数正規分布の場合、それは [3]に等しい。
CV
[
X
]
{\displaystyle \operatorname {CV} [X]}
SD
[
X
]
E
[
X
]
{\displaystyle {\tfrac {\operatorname {SD} [X]}{\operatorname {E} [X]}}}
CV
[
X
]
=
e
σ
2
−
1
.
{\displaystyle \operatorname {CV} [X]={\sqrt {e^{\sigma ^{2}}-1}}.}
この推定値は、幾何分散を使用しているため、「幾何CV」(GCV)と呼ばれることもあります [19] [20] 。算術標準偏差とは異なり、算術変動係数は算術平均とは無関係です。
算術平均と算術分散がわかれば、
パラメータ μ と σ を取得できます。
μ
=
ln
(
E
[
X
]
2
E
[
X
2
]
)
=
ln
(
E
[
X
]
2
Var
[
X
]
+
E
[
X
]
2
)
,
σ
2
=
ln
(
E
[
X
2
]
E
[
X
]
2
)
=
ln
(
1
+
Var
[
X
]
E
[
X
]
2
)
.
{\displaystyle {\begin{aligned}\mu &=\ln \left({\frac {\operatorname {E} [X]^{2}}{\sqrt {\operatorname {E} [X^{2}]}}}\right)=\ln \left({\frac {\operatorname {E} [X]^{2}}{\sqrt {\operatorname {Var} [X]+\operatorname {E} [X]^{2}}}}\right),\\[4pt]\sigma ^{2}&=\ln \left({\frac {\operatorname {E} [X^{2}]}{\operatorname {E} [X]^{2}}}\right)=\ln \left(1+{\frac {\operatorname {Var} [X]}{\operatorname {E} [X]^{2}}}\right).\end{aligned}}}
確率分布はモーメント E[ X n ] = e nμ + によって一意に決定されるわけではない。 1 / 2 n 2 σ 2 ( n ≥ 1 。つまり、同じモーメントを持つ他の分布も存在する。 [4] 実際、対数正規分布と同じモーメントを持つ分布の族全体が存在する。 [ 要出典 ]
歪度 が異なる 2 つの対数正規分布の 平均 、 中央値 、 最頻値 の比較 。
モード は 確率密度関数のグローバル最大値の点です。特に、方程式を解くと 次のようになります。
(
ln
f
)
′
=
0
{\displaystyle (\ln f)'=0}
Mode
[
X
]
=
e
μ
−
σ
2
.
{\displaystyle \operatorname {Mode} [X]=e^{\mu -\sigma ^{2}}.}
対数変換された 変数は 正規分布を持ち、単調変換の下では分位数が保存される ため、の分位数 は
Y
=
ln
X
{\displaystyle Y=\ln X}
X
{\displaystyle X}
q
X
(
α
)
=
e
μ
+
σ
q
Φ
(
α
)
=
μ
∗
(
σ
∗
)
q
Φ
(
α
)
,
{\displaystyle q_{X}(\alpha )=e^{\mu +\sigma q_{\Phi }(\alpha )}=\mu ^{*}(\sigma ^{*})^{q_{\Phi }(\alpha )},}
ここで、は 標準正規分布の分位数です。
q
Φ
(
α
)
{\displaystyle q_{\Phi }(\alpha )}
具体的には、対数正規分布の中央値はその乗法平均に等しい。 [21]
Med
[
X
]
=
e
μ
=
μ
∗
.
{\displaystyle \operatorname {Med} [X]=e^{\mu }=\mu ^{*}~.}
部分的な期待
閾値に関する 確率変数の部分期待値は 次のように定義される。
X
{\displaystyle X}
k
{\displaystyle k}
g
(
k
)
=
∫
k
∞
x
f
X
(
x
∣
X
>
k
)
d
x
.
{\displaystyle g(k)=\int _{k}^{\infty }xf_{X}(x\mid X>k)\,dx.}
あるいは、条件付き期待値 の定義を使用すると 、 と書くこともできます 。対数正規分布のランダム変数の場合、部分期待値は次のように表されます。
g
(
k
)
=
E
[
X
∣
X
>
k
]
P
(
X
>
k
)
{\displaystyle g(k)=\operatorname {E} [X\mid X>k]P(X>k)}
g
(
k
)
=
∫
k
∞
x
f
X
(
x
∣
X
>
k
)
d
x
=
e
μ
+
1
2
σ
2
Φ
(
μ
+
σ
2
−
ln
k
σ
)
{\displaystyle g(k)=\int _{k}^{\infty }xf_{X}(x\mid X>k)\,dx=e^{\mu +{\tfrac {1}{2}}\sigma ^{2}}\,\Phi \!\left({\frac {\mu +\sigma ^{2}-\ln k}{\sigma }}\right)}
ここで、 は 正規累積分布関数 です 。式の導出については、 トークページ で説明しています。偏期待値式は 保険 や 経済学に応用されており、 ブラック・ショールズ式 につながる偏微分方程式を解く際に使用されます 。
Φ
{\displaystyle \Phi }
条件付き期待値
対数正規分布のランダム変数の条件付き期待値は 、閾値に関して 、その部分期待値をその範囲内にある累積確率で割った値です。
X
{\displaystyle X}
k
{\displaystyle k}
E
[
X
∣
X
<
k
]
=
e
μ
+
σ
2
2
⋅
Φ
[
ln
(
k
)
−
μ
−
σ
2
σ
]
Φ
[
ln
(
k
)
−
μ
σ
]
E
[
X
∣
X
⩾
k
]
=
e
μ
+
σ
2
2
⋅
Φ
[
μ
+
σ
2
−
ln
(
k
)
σ
]
1
−
Φ
[
ln
(
k
)
−
μ
σ
]
E
[
X
∣
X
∈
[
k
1
,
k
2
]
]
=
e
μ
+
σ
2
2
⋅
Φ
[
ln
(
k
2
)
−
μ
−
σ
2
σ
]
−
Φ
[
ln
(
k
1
)
−
μ
−
σ
2
σ
]
Φ
[
ln
(
k
2
)
−
μ
σ
]
−
Φ
[
ln
(
k
1
)
−
μ
σ
]
{\displaystyle {\begin{aligned}E[X\mid X<k]&=e^{\mu +{\frac {\sigma ^{2}}{2}}}\cdot {\frac {\Phi \left[{\frac {\ln(k)-\mu -\sigma ^{2}}{\sigma }}\right]}{\Phi \left[{\frac {\ln(k)-\mu }{\sigma }}\right]}}\\[8pt]E[X\mid X\geqslant k]&=e^{\mu +{\frac {\sigma ^{2}}{2}}}\cdot {\frac {\Phi \left[{\frac {\mu +\sigma ^{2}-\ln(k)}{\sigma }}\right]}{1-\Phi \left[{\frac {\ln(k)-\mu }{\sigma }}\right]}}\\[8pt]E[X\mid X\in [k_{1},k_{2}]]&=e^{\mu +{\frac {\sigma ^{2}}{2}}}\cdot {\frac {\Phi \left[{\frac {\ln(k_{2})-\mu -\sigma ^{2}}{\sigma }}\right]-\Phi \left[{\frac {\ln(k_{1})-\mu -\sigma ^{2}}{\sigma }}\right]}{\Phi \left[{\frac {\ln(k_{2})-\mu }{\sigma }}\right]-\Phi \left[{\frac {\ln(k_{1})-\mu }{\sigma }}\right]}}\end{aligned}}}
代替パラメータ化 または
による特徴付けに加えて 、対数正規分布をパラメータ化する方法は複数あります。 確率分布 の知識ベースとオントロジーである ProbOnto [22] [23] では、次のような7つの形式がリストされています。
μ
,
σ
{\displaystyle \mu ,\sigma }
μ
∗
,
σ
∗
{\displaystyle \mu ^{*},\sigma ^{*}}
対数正規分布のパラメータ化の概要。
LogNormal1(μ,σ) は、 平均 μ と 標準偏差 σ で、両方とも対数スケールです [24]
P
(
x
;
μ
,
σ
)
=
1
x
σ
2
π
exp
[
−
(
ln
x
−
μ
)
2
2
σ
2
]
{\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {\sigma }})={\frac {1}{x\sigma {\sqrt {2\pi }}}}\exp \left[-{\frac {(\ln x-\mu )^{2}}{2\sigma ^{2}}}\right]}
LogNormal2(μ,υ) (平均 μ と分散 υ は両方とも対数スケール)
P
(
x
;
μ
,
v
)
=
1
x
v
2
π
exp
[
−
(
ln
x
−
μ
)
2
2
v
]
{\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {v}})={\frac {1}{x{\sqrt {v}}{\sqrt {2\pi }}}}\exp \left[-{\frac {(\ln x-\mu )^{2}}{2v}}\right]}
LogNormal3(m,σ)、 中央値 mは自然スケール、標準偏差σは対数スケール [24]
P
(
x
;
m
,
σ
)
=
1
x
σ
2
π
exp
[
−
ln
2
(
x
/
m
)
2
σ
2
]
{\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {\sigma }})={\frac {1}{x\sigma {\sqrt {2\pi }}}}\exp \left[-{\frac {\ln ^{2}(x/m)}{2\sigma ^{2}}}\right]}
LogNormal4(m,cv)、中央値m、 変動係数 cv、ともに自然スケール
P
(
x
;
m
,
c
v
)
=
1
x
ln
(
c
v
2
+
1
)
2
π
exp
[
−
ln
2
(
x
/
m
)
2
ln
(
c
v
2
+
1
)
]
{\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {cv}})={\frac {1}{x{\sqrt {\ln(cv^{2}+1)}}{\sqrt {2\pi }}}}\exp \left[-{\frac {\ln ^{2}(x/m)}{2\ln(cv^{2}+1)}}\right]}
LogNormal5(μ,τ) (平均 μ と精度 τ、両方とも対数スケール) [25]
P
(
x
;
μ
,
τ
)
=
τ
2
π
1
x
exp
[
−
τ
2
(
ln
x
−
μ
)
2
]
{\displaystyle P(x;{\boldsymbol {\mu }},{\boldsymbol {\tau }})={\sqrt {\frac {\tau }{2\pi }}}{\frac {1}{x}}\exp \left[-{\frac {\tau }{2}}(\ln x-\mu )^{2}\right]}
LogNormal6(m,σ g )、中央値m、幾何 標準偏差 σ g 、ともに自然尺度 [26]
P
(
x
;
m
,
σ
g
)
=
1
x
ln
(
σ
g
)
2
π
exp
[
−
ln
2
(
x
/
m
)
2
ln
2
(
σ
g
)
]
{\displaystyle P(x;{\boldsymbol {m}},{\boldsymbol {\sigma _{g}}})={\frac {1}{x\ln(\sigma _{g}){\sqrt {2\pi }}}}\exp \left[-{\frac {\ln ^{2}(x/m)}{2\ln ^{2}(\sigma _{g})}}\right]}
LogNormal7(μ N ,σ N )、平均 μ N および標準偏差 σ N 、どちらも自然スケール [27]
P
(
x
;
μ
N
,
σ
N
)
=
1
x
2
π
ln
(
1
+
σ
N
2
/
μ
N
2
)
exp
(
−
[
ln
x
−
ln
μ
N
1
+
σ
N
2
/
μ
N
2
]
2
2
ln
(
1
+
σ
N
2
/
μ
N
2
)
)
{\displaystyle P(x;{\boldsymbol {\mu _{N}}},{\boldsymbol {\sigma _{N}}})={\frac {1}{x{\sqrt {2\pi \ln \left(1+\sigma _{N}^{2}/\mu _{N}^{2}\right)}}}}\exp \left(-{\frac {{\Big [}\ln x-\ln {\frac {\mu _{N}}{\sqrt {1+\sigma _{N}^{2}/\mu _{N}^{2}}}}{\Big ]}^{2}}{2\ln(1+\sigma _{N}^{2}/\mu _{N}^{2})}}\right)}
再パラメータ化の例
2つの異なる最適設計ツール、例えばPFIM [28] とPopED [29] を使用してモデルを実行したい状況を考えてみましょう。 前者はそれぞれLN2パラメータ化をサポートし、後者はLN7パラメータ化をサポートしています。したがって、再パラメータ化が必要であり、そうしないと2つのツールが異なる結果を生成します。
遷移については 次の式 が成り立ちます 。
LN2
(
μ
,
v
)
→
LN7
(
μ
N
,
σ
N
)
{\displaystyle \operatorname {LN2} (\mu ,v)\to \operatorname {LN7} (\mu _{N},\sigma _{N})}
μ
N
=
exp
(
μ
+
v
/
2
)
{\textstyle \mu _{N}=\exp(\mu +v/2)}
σ
N
=
exp
(
μ
+
v
/
2
)
exp
(
v
)
−
1
{\textstyle \sigma _{N}=\exp(\mu +v/2){\sqrt {\exp(v)-1}}}
遷移については 次の式 が成り立ちます 。
LN7
(
μ
N
,
σ
N
)
→
LN2
(
μ
,
v
)
{\displaystyle \operatorname {LN7} (\mu _{N},\sigma _{N})\to \operatorname {LN2} (\mu ,v)}
μ
=
ln
(
μ
N
/
1
+
σ
N
2
/
μ
N
2
)
{\textstyle \mu =\ln \left(\mu _{N}/{\sqrt {1+\sigma _{N}^{2}/\mu _{N}^{2}}}\right)}
v
=
ln
(
1
+
σ
N
2
/
μ
N
2
)
{\textstyle v=\ln(1+\sigma _{N}^{2}/\mu _{N}^{2})}
残りの再パラメータ化の式はすべて、プロジェクトのウェブサイトの仕様書に記載されています。 [30]
多重、相互、力
定数による乗算: の 場合、
X
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}
a
X
∼
Lognormal
(
μ
+
ln
a
,
σ
2
)
{\displaystyle aX\sim \operatorname {Lognormal} (\mu +\ln a,\ \sigma ^{2})}
a
>
0.
{\displaystyle a>0.}
逆数: ならば
X
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}
1
X
∼
Lognormal
(
−
μ
,
σ
2
)
.
{\displaystyle {\tfrac {1}{X}}\sim \operatorname {Lognormal} (-\mu ,\ \sigma ^{2}).}
パワー:もし 、ならば 、
X
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}
X
a
∼
Lognormal
(
a
μ
,
a
2
σ
2
)
{\displaystyle X^{a}\sim \operatorname {Lognormal} (a\mu ,\ a^{2}\sigma ^{2})}
a
≠
0.
{\displaystyle a\neq 0.}
独立した対数正規分布のランダム変数の乗算と除算
2 つの独立した 対数正規変数 とを 乗算 [除算] すると 、積 [比] は再び対数正規となり、パラメータは [ ] および となります (ただし ) 。これは、このような変数の積に簡単に一般化できます 。
X
1
{\displaystyle X_{1}}
X
2
{\displaystyle X_{2}}
μ
=
μ
1
+
μ
2
{\displaystyle \mu =\mu _{1}+\mu _{2}}
μ
=
μ
1
−
μ
2
{\displaystyle \mu =\mu _{1}-\mu _{2}}
σ
{\displaystyle \sigma }
σ
2
=
σ
1
2
+
σ
2
2
{\displaystyle \sigma ^{2}=\sigma _{1}^{2}+\sigma _{2}^{2}}
n
{\displaystyle n}
より一般的には、 が 独立した対数正規分布に従う変数である場合、
X
j
∼
Lognormal
(
μ
j
,
σ
j
2
)
{\displaystyle X_{j}\sim \operatorname {Lognormal} (\mu _{j},\sigma _{j}^{2})}
n
{\displaystyle n}
Y
=
∏
j
=
1
n
X
j
∼
Lognormal
(
∑
j
=
1
n
μ
j
,
∑
j
=
1
n
σ
j
2
)
.
{\displaystyle Y=\textstyle \prod _{j=1}^{n}X_{j}\sim \operatorname {Lognormal} {\Big (}\textstyle \sum _{j=1}^{n}\mu _{j},\ \sum _{j=1}^{n}\sigma _{j}^{2}{\Big )}.}
乗法中心極限定理
独立で同一に分布する正のランダム変数 の幾何平均または乗法平均は 、 が有限であると仮定すると、 に対して、近似的にパラメータ および を持つ対数正規分布を示します 。
n
{\displaystyle n}
X
i
{\displaystyle X_{i}}
n
→
∞
{\displaystyle n\to \infty }
μ
=
E
[
ln
(
X
i
)
]
{\displaystyle \mu =E[\ln(X_{i})]}
σ
2
=
var
[
ln
(
X
i
)
]
/
n
{\displaystyle \sigma ^{2}={\mbox{var}}[\ln(X_{i})]/n}
σ
2
{\displaystyle \sigma ^{2}}
実際、ランダム変数は同一に分布している必要はありません。 の分布がすべて有限分散を持ち、 中心極限定理 のさまざまな変形のいずれかの他の条件を満たしていれば十分です 。
ln
(
X
i
)
{\displaystyle \ln(X_{i})}
これは一般にジブラの法則 として知られています 。
対数正規分布の裾の重なり
対数正規分布が真のヘビーテール分布であるかどうかは、まだ議論が続いている。主な理由は、その分散が常に有限であり、たとえば特定のパレート分布で起こることと異なるためである。しかし、最近の研究では、ロビンソン非標準分析を使用して、無限分散を持つ対数正規分布を作成することが可能であることが示されています。 [31]
他の
対数正規分布から生じるデータセットは対称的な ローレンツ曲線 を持ちます( ローレンツ非対称係数 も参照)。 [32]
この分布の 調和平均 、幾何平均 、算術平均は関連しており、 [33] その関係は次のように表される。
H
{\displaystyle H}
G
{\displaystyle G}
A
{\displaystyle A}
H
=
G
2
A
.
{\displaystyle H={\frac {G^{2}}{A}}.}
対数正規分布は 無限に分割可能 であるが [34] 、簡単に分布を作れる 安定した分布 ではない。 [35]
が正規分布 で ある場合 、
X
∼
N
(
μ
,
σ
2
)
{\displaystyle X\sim {\mathcal {N}}(\mu ,\sigma ^{2})}
exp
(
X
)
∼
Lognormal
(
μ
,
σ
2
)
.
{\displaystyle \exp(X)\sim \operatorname {Lognormal} (\mu ,\sigma ^{2}).}
が対数正規分布する 場合、 は 正規確率変数です。
X
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}
ln
(
X
)
∼
N
(
μ
,
σ
2
)
{\displaystyle \ln(X)\sim {\mathcal {N}}(\mu ,\sigma ^{2})}
を、変化する可能性のあるパラメータと パラメータ を持つ独立した対数正規分布変数、およびと します 。 の分布は閉じた形式の表現を持ちませんが、 右裾で 別の対数正規分布によって合理的に近似できます。 [36] 0 の近傍での確率密度関数は特徴付けられており [35] 、どの対数正規分布にも似ていません。LF Fenton による一般的な近似値 (ただし、以前は RI Wilkinson によって述べられ、Marlow [37] によって数学的に正当化されています) は、別の対数正規分布の平均と分散を一致させることによって得られます。すべてが 同じ分散パラメータ を持つ 場合 、これらの式は次のように簡略化されます。
X
j
∼
Lognormal
(
μ
j
,
σ
j
2
)
{\displaystyle X_{j}\sim \operatorname {Lognormal} (\mu _{j},\sigma _{j}^{2})}
σ
{\displaystyle \sigma }
μ
{\displaystyle \mu }
Y
=
∑
j
=
1
n
X
j
{\textstyle Y=\sum _{j=1}^{n}X_{j}}
Y
{\displaystyle Y}
Z
{\displaystyle Z}
σ
Z
2
=
ln
[
∑
e
2
μ
j
+
σ
j
2
(
e
σ
j
2
−
1
)
(
∑
e
μ
j
+
σ
j
2
/
2
)
2
+
1
]
,
μ
Z
=
ln
[
∑
e
μ
j
+
σ
j
2
/
2
]
−
σ
Z
2
2
.
{\displaystyle {\begin{aligned}\sigma _{Z}^{2}&=\ln \!\left[{\frac {\sum e^{2\mu _{j}+\sigma _{j}^{2}}(e^{\sigma _{j}^{2}}-1)}{(\sum e^{\mu _{j}+\sigma _{j}^{2}/2})^{2}}}+1\right],\\\mu _{Z}&=\ln \!\left[\sum e^{\mu _{j}+\sigma _{j}^{2}/2}\right]-{\frac {\sigma _{Z}^{2}}{2}}.\end{aligned}}}
X
j
{\displaystyle X_{j}}
σ
j
=
σ
{\displaystyle \sigma _{j}=\sigma }
σ
Z
2
=
ln
[
(
e
σ
2
−
1
)
∑
e
2
μ
j
(
∑
e
μ
j
)
2
+
1
]
,
μ
Z
=
ln
[
∑
e
μ
j
]
+
σ
2
2
−
σ
Z
2
2
.
{\displaystyle {\begin{aligned}\sigma _{Z}^{2}&=\ln \!\left[(e^{\sigma ^{2}}-1){\frac {\sum e^{2\mu _{j}}}{(\sum e^{\mu _{j}})^{2}}}+1\right],\\\mu _{Z}&=\ln \!\left[\sum e^{\mu _{j}}\right]+{\frac {\sigma ^{2}}{2}}-{\frac {\sigma _{Z}^{2}}{2}}.\end{aligned}}}
より正確な近似値を得るためには、 モンテカルロ法を 用いて累積分布関数、pdf、右裾を推定すればよい。 [38] [39]
相関対数正規分布するランダム変数の合計は、対数正規分布によって近似することもできる [ 要出典 ]
S
+
=
E
[
∑
i
X
i
]
=
∑
i
E
[
X
i
]
=
∑
i
e
μ
i
+
σ
i
2
/
2
σ
Z
2
=
1
/
S
+
2
∑
i
,
j
cor
i
j
σ
i
σ
j
E
[
X
i
]
E
[
X
j
]
=
1
/
S
+
2
∑
i
,
j
cor
i
j
σ
i
σ
j
e
μ
i
+
σ
i
2
/
2
e
μ
j
+
σ
j
2
/
2
μ
Z
=
ln
(
S
+
)
−
σ
Z
2
/
2
{\displaystyle {\begin{aligned}S_{+}&=\operatorname {E} \left[\sum _{i}X_{i}\right]=\sum _{i}\operatorname {E} [X_{i}]=\sum _{i}e^{\mu _{i}+\sigma _{i}^{2}/2}\\\sigma _{Z}^{2}&=1/S_{+}^{2}\,\sum _{i,j}\operatorname {cor} _{ij}\sigma _{i}\sigma _{j}\operatorname {E} [X_{i}]\operatorname {E} [X_{j}]=1/S_{+}^{2}\,\sum _{i,j}\operatorname {cor} _{ij}\sigma _{i}\sigma _{j}e^{\mu _{i}+\sigma _{i}^{2}/2}e^{\mu _{j}+\sigma _{j}^{2}/2}\\\mu _{Z}&=\ln \left(S_{+}\right)-\sigma _{Z}^{2}/2\end{aligned}}}
の 場合には、サポートを持つ 3パラメータ対数正規 分布 に従うと言われます 。 [40] 、 。
X
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle X\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}
X
+
c
{\displaystyle X+c}
x
∈
(
c
,
+
∞
)
{\displaystyle x\in (c,+\infty )}
E
[
X
+
c
]
=
E
[
X
]
+
c
{\displaystyle \operatorname {E} [X+c]=\operatorname {E} [X]+c}
Var
[
X
+
c
]
=
Var
[
X
]
{\displaystyle \operatorname {Var} [X+c]=\operatorname {Var} [X]}
対数正規分布は、半有界 ジョンソンSU分布 の特殊なケースである。 [41]
の 場合には 、 (鈴木分布)。
X
∣
Y
∼
Rayleigh
(
Y
)
{\displaystyle X\mid Y\sim \operatorname {Rayleigh} (Y)}
Y
∼
Lognormal
(
μ
,
σ
2
)
{\displaystyle Y\sim \operatorname {Lognormal} (\mu ,\sigma ^{2})}
X
∼
Suzuki
(
μ
,
σ
)
{\displaystyle X\sim \operatorname {Suzuki} (\mu ,\sigma )}
対数正規分布の代わりとなる、より基本的な関数[42] で表すことのできる 分布は、ロジスティック分布 に基づいて得られ、 CDF の近似値を得ることができます。 これは 対数ロジスティック分布 です。
F
(
x
;
μ
,
σ
)
=
[
(
e
μ
x
)
π
/
(
σ
3
)
+
1
]
−
1
.
{\displaystyle F(x;\mu ,\sigma )=\left[\left({\frac {e^{\mu }}{x}}\right)^{\pi /(\sigma {\sqrt {3}})}+1\right]^{-1}.}
統計的推論
パラメータの推定
対数正規分布パラメータ μ と σの 最大尤度 推定値を決定するには、 正規分布 の場合と 同じ手順を 使用できます 。
ここで、 は 正規分布の密度関数であることに注意してください 。したがって、対数尤度関数は
L
(
μ
,
σ
)
=
∏
i
=
1
n
1
x
i
φ
μ
,
σ
(
ln
x
i
)
,
{\displaystyle L(\mu ,\sigma )=\prod _{i=1}^{n}{\frac {1}{x_{i}}}\varphi _{\mu ,\sigma }(\ln x_{i}),}
φ
{\displaystyle \varphi }
N
(
μ
,
σ
2
)
{\displaystyle {\mathcal {N}}(\mu ,\sigma ^{2})}
ℓ
(
μ
,
σ
∣
x
1
,
x
2
,
…
,
x
n
)
=
−
∑
i
ln
x
i
+
ℓ
N
(
μ
,
σ
∣
ln
x
1
,
ln
x
2
,
…
,
ln
x
n
)
.
{\displaystyle \ell (\mu ,\sigma \mid x_{1},x_{2},\ldots ,x_{n})=-\sum _{i}\ln x_{i}+\ell _{N}(\mu ,\sigma \mid \ln x_{1},\ln x_{2},\dots ,\ln x_{n}).}
最初の項はμ と σ に関して一定であるため 、両方の対数尤度関数 とは 同じおよびで最大値に達します 。 したがって、最大尤度推定量は、観測値 に対する正規分布の推定量と同一です 。
ℓ
{\displaystyle \ell }
ℓ
N
{\displaystyle \ell _{N}}
μ
{\displaystyle \mu }
σ
{\displaystyle \sigma }
ln
x
1
,
ln
x
2
,
…
,
ln
x
n
)
{\displaystyle \ln x_{1},\ln x_{2},\dots ,\ln x_{n})}
μ
^
=
∑
i
ln
x
i
n
,
σ
^
2
=
∑
i
(
ln
x
i
−
μ
^
)
2
n
.
{\displaystyle {\widehat {\mu }}={\frac {\sum _{i}\ln x_{i}}{n}},\qquad {\widehat {\sigma }}^{2}={\frac {\sum _{i}\left(\ln x_{i}-{\widehat {\mu }}\right)^{2}}{n}}.}
有限の n の場合、 の推定値は 不偏ですが、 の推定値は 偏りがあります。正規分布の場合、 の不偏推定値は、 の式の分母 n を n −1に 置き換えることで得られます 。
μ
{\displaystyle \mu }
σ
{\displaystyle \sigma }
σ
{\displaystyle \sigma }
σ
^
2
{\displaystyle {\widehat {\sigma }}^{2}}
個々の値 が利用できないが、サンプルの平均 と 標準偏差 s が利用できる場合は、 モーメント法を 使用できます。対応するパラメータは、および の期待値 と分散の 方程式を解くことで得られる次の式によって決定されます 。
x
1
,
x
2
,
…
,
x
n
{\displaystyle x_{1},x_{2},\ldots ,x_{n}}
x
¯
{\displaystyle {\bar {x}}}
E
[
X
]
{\displaystyle \operatorname {E} [X]}
Var
[
X
]
{\displaystyle \operatorname {Var} [X]}
μ
{\displaystyle \mu }
σ
{\displaystyle \sigma }
μ
=
ln
(
x
¯
1
+
σ
^
2
/
x
¯
2
)
,
σ
2
=
ln
(
1
+
σ
^
2
/
x
¯
2
)
.
{\displaystyle \mu =\ln \left({\frac {\bar {x}}{\sqrt {1+{\widehat {\sigma }}^{2}/{\bar {x}}^{2}}}}\right),\qquad \sigma ^{2}=\ln \left(1+{{\widehat {\sigma }}^{2}}/{\bar {x}}^{2}\right).}
区間推定
対数正規分布データを分析する際に区間推定値 を取得する最も効率的な方法 は、正規分布に基づくよく知られた方法を対数変換されたデータに適用し、適切な場合は結果を逆変換することです。
予測区間
基本的な例として、 予測区間 が挙げられます。正規分布の場合、区間には 確率の約 3 分の 2 (68%) が含まれ、 95% が含まれます (または大きなサンプルの場合)。したがって、対数正規分布の場合、区間には
確率の 2/3 が含まれ、
95% が含まれます。推定パラメータを使用すると、これらの区間にはほぼ同じ割合のデータが含まれるはずです。
[
μ
−
σ
,
μ
+
σ
]
{\displaystyle [\mu -\sigma ,\mu +\sigma ]}
[
μ
−
2
σ
,
μ
+
2
σ
]
{\displaystyle [\mu -2\sigma ,\mu +2\sigma ]}
[
μ
∗
/
σ
∗
,
μ
∗
⋅
σ
∗
]
=
[
μ
∗
×
/
σ
∗
]
{\displaystyle [\mu ^{*}/\sigma ^{*},\mu ^{*}\cdot \sigma ^{*}]=[\mu ^{*}{}^{\times }\!\!/\sigma ^{*}]}
[
μ
∗
/
(
σ
∗
)
2
,
μ
∗
⋅
(
σ
∗
)
2
]
=
[
μ
∗
×
/
(
σ
∗
)
2
]
{\displaystyle [\mu ^{*}/(\sigma ^{*})^{2},\mu ^{*}\cdot (\sigma ^{*})^{2}]=[\mu ^{*}{}^{\times }\!\!/(\sigma ^{*})^{2}]}
信頼区間 eμ
この原理を用いると、 の 信頼区間 は となり 、 は 標準誤差、 q は自由度 n-1 の t 分布 の 97.5% 分位数となります 。逆変換により、 (中央値) の信頼区間は 次のようになります 。
μ
{\displaystyle \mu }
[
μ
^
±
q
⋅
s
e
^
]
{\displaystyle [{\widehat {\mu }}\pm q\cdot {\widehat {\mathop {se} }}]}
s
e
=
σ
^
/
n
{\displaystyle \mathop {se} ={\widehat {\sigma }}/{\sqrt {n}}}
μ
∗
=
e
μ
{\displaystyle \mu ^{*}=e^{\mu }}
[
μ
^
∗
×
/
(
sem
∗
)
q
]
{\displaystyle [{\widehat {\mu }}^{*}{}^{\times }\!\!/(\operatorname {sem} ^{*})^{q}]}
sem
∗
=
(
σ
^
∗
)
1
/
n
{\displaystyle \operatorname {sem} ^{*}=({\widehat {\sigma }}^{*})^{1/{\sqrt {n}}}}
E(X)の信頼区間
文献では、 (対数正規分布の平均) の 信頼区間を計算するためのいくつかのオプションが議論されています。これらには、 ブートストラップ法 やその他のさまざまな方法が含まれます。 [43] [44]
μ
{\displaystyle \mu }
コックス法 [a] は推定値をプラグインすることを提案している
μ
^
=
∑
i
ln
x
i
n
,
S
2
=
∑
i
(
ln
x
i
−
μ
^
)
2
n
−
1
{\displaystyle {\widehat {\mu }}={\frac {\sum _{i}\ln x_{i}}{n}},\qquad S^{2}={\frac {\sum _{i}\left(\ln x_{i}-{\widehat {\mu }}\right)^{2}}{n-1}}}
これらを使用して、 次のように
おおよその信頼区間を構築します。
C
I
(
E
(
X
)
)
:
e
(
μ
^
+
S
2
2
±
z
1
−
α
2
S
2
n
+
S
4
2
(
n
−
1
)
)
{\displaystyle CI(E(X)):e^{\left({\hat {\mu }}+{\frac {S^{2}}{2}}\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S^{2}}{n}}+{\frac {S^{4}}{2(n-1)}}}}\right)}}
Olsson 2005は、を に
置き換えた「修正Cox法」を提案した 。これは、小さなサンプルサイズでより良いカバレッジ結果をもたらすと思われる。 [43] :セクション3.4
z
1
−
α
2
{\displaystyle z_{1-{\frac {\alpha }{2}}}}
t
n
−
1
,
1
−
α
2
{\displaystyle t_{n-1,1-{\frac {\alpha }{2}}}}
2つの対数正規分布を比較するための信頼区間
2 つの対数正規分布を比較することは、多くの場合、たとえば治療群と対照群(たとえば、 A/B テスト )から行う必要があります。パラメータが と で、 サンプル サイズ がそれぞれ と である 2 つの独立した対数正規分布からのサンプルがあります 。
(
μ
1
,
σ
1
2
)
{\displaystyle (\mu _{1},\sigma _{1}^{2})}
(
μ
2
,
σ
2
2
)
{\displaystyle (\mu _{2},\sigma _{2}^{2})}
n
1
{\displaystyle n_{1}}
n
2
{\displaystyle n_{2}}
2 つの中央値を比較するには、それぞれのログを取得し、単純な信頼区間を作成して、それを指数スケールに戻すだけで簡単に行うことができます。
C
I
(
e
μ
1
−
μ
2
)
:
e
(
μ
^
1
−
μ
^
2
±
z
1
−
α
2
S
1
2
n
+
S
2
2
n
)
{\displaystyle CI(e^{\mu _{1}-\mu _{2}}):e^{\left({\hat {\mu }}_{1}-{\hat {\mu }}_{2}\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S_{1}^{2}}{n}}+{\frac {S_{2}^{2}}{n}}}}\right)}}
これらのCIは、疫学において相対リスク と オッズ比 のCIを計算するためによく使用されます 。 [47] そこでの計算方法は、2つの近似正規分布(例えば、RRの場合はp1とp2)があり、それらの比率を計算するというものです。 [b]
ただし、2 つのサンプルの期待値 (平均) の比率も興味深いかもしれませんが、その開発にはさらに作業が必要です。平均の比率は次のとおりです。
E
(
X
1
)
E
(
X
2
)
=
e
μ
1
+
σ
1
2
2
e
μ
2
+
σ
2
2
2
=
e
(
μ
1
−
μ
2
)
+
1
2
(
σ
1
2
−
σ
2
2
)
{\displaystyle {\frac {E(X_{1})}{E(X_{2})}}={\frac {e^{\mu _{1}+{\frac {\sigma _{1}^{2}}{2}}}}{e^{\mu _{2}+{\frac {\sigma _{2}^{2}}{2}}}}}=e^{(\mu _{1}-\mu _{2})+{\frac {1}{2}}(\sigma _{1}^{2}-\sigma _{2}^{2})}}
これらの各パラメータに推定値をプラグインすると、対数正規分布も生成されます。つまり、上で説明した Cox 法をこのユースケースに同様に使用できます。
C
I
(
E
(
X
1
)
E
(
X
2
)
=
e
μ
1
+
σ
1
2
2
e
μ
2
+
σ
2
2
2
)
:
e
(
(
μ
^
1
−
μ
^
2
+
1
2
S
1
2
−
1
2
S
2
2
)
±
z
1
−
α
2
S
1
2
n
1
+
S
2
2
n
2
+
S
1
4
2
(
n
1
−
1
)
+
S
2
4
2
(
n
2
−
1
)
)
{\displaystyle CI\left({\frac {E(X_{1})}{E(X_{2})}}={\frac {e^{\mu _{1}+{\frac {\sigma _{1}^{2}}{2}}}}{e^{\mu _{2}+{\frac {\sigma _{2}^{2}}{2}}}}}\right):e^{\left(({\hat {\mu }}_{1}-{\hat {\mu }}_{2}+{\frac {1}{2}}S_{1}^{2}-{\frac {1}{2}}S_{2}^{2})\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S_{1}^{2}}{n_{1}}}+{\frac {S_{2}^{2}}{n_{2}}}+{\frac {S_{1}^{4}}{2(n_{1}-1)}}+{\frac {S_{2}^{4}}{2(n_{2}-1)}}}}\right)}}
[証拠]
この比率の信頼区間を構築するには、まず が 正規分布に従い、 と は両方 とも カイ二乗分布 に従い 、これが( 関連 パラメータを持つ CLTにより) ほぼ 正規分布していることに注目します 。
μ
^
1
−
μ
^
2
{\displaystyle {\hat {\mu }}_{1}-{\hat {\mu }}_{2}}
S
1
2
{\displaystyle S_{1}^{2}}
S
2
2
{\displaystyle S_{2}^{2}}
これはつまり
(
μ
^
1
−
μ
^
2
+
1
2
S
1
2
−
1
2
S
2
2
)
∼
N
(
(
μ
1
−
μ
2
)
+
1
2
(
σ
1
2
−
σ
2
2
)
,
σ
1
2
n
1
+
σ
2
2
n
2
+
σ
1
4
2
(
n
1
−
1
)
+
σ
2
4
2
(
n
2
−
1
)
)
{\displaystyle ({\hat {\mu }}_{1}-{\hat {\mu }}_{2}+{\frac {1}{2}}S_{1}^{2}-{\frac {1}{2}}S_{2}^{2})\sim N\left((\mu _{1}-\mu _{2})+{\frac {1}{2}}(\sigma _{1}^{2}-\sigma _{2}^{2}),{\frac {\sigma _{1}^{2}}{n_{1}}}+{\frac {\sigma _{2}^{2}}{n_{2}}}+{\frac {\sigma _{1}^{4}}{2(n_{1}-1)}}+{\frac {\sigma _{2}^{4}}{2(n_{2}-1)}}\right)}
上記に基づいて、標準的な 信頼区間は( ピボタル量を 使用して)次のように 構築できます。
また、信頼区間は単調変換に対して保存されるため、次のようになります。
(
μ
^
1
−
μ
^
2
+
1
2
S
1
2
−
1
2
S
2
2
)
±
z
1
−
α
2
S
1
2
n
1
+
S
2
2
n
2
+
S
1
4
2
(
n
1
−
1
)
+
S
2
4
2
(
n
2
−
1
)
{\displaystyle ({\hat {\mu }}_{1}-{\hat {\mu }}_{2}+{\frac {1}{2}}S_{1}^{2}-{\frac {1}{2}}S_{2}^{2})\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S_{1}^{2}}{n_{1}}}+{\frac {S_{2}^{2}}{n_{2}}}+{\frac {S_{1}^{4}}{2(n_{1}-1)}}+{\frac {S_{2}^{4}}{2(n_{2}-1)}}}}}
C
I
(
E
(
X
1
)
E
(
X
2
)
=
e
μ
1
+
σ
1
2
2
e
μ
2
+
σ
2
2
2
)
:
e
(
(
μ
^
1
−
μ
^
2
+
1
2
S
1
2
−
1
2
S
2
2
)
±
z
1
−
α
2
S
1
2
n
1
+
S
2
2
n
2
+
S
1
4
2
(
n
1
−
1
)
+
S
2
4
2
(
n
2
−
1
)
)
{\displaystyle CI\left({\frac {E(X_{1})}{E(X_{2})}}={\frac {e^{\mu _{1}+{\frac {\sigma _{1}^{2}}{2}}}}{e^{\mu _{2}+{\frac {\sigma _{2}^{2}}{2}}}}}\right):e^{\left(({\hat {\mu }}_{1}-{\hat {\mu }}_{2}+{\frac {1}{2}}S_{1}^{2}-{\frac {1}{2}}S_{2}^{2})\pm z_{1-{\frac {\alpha }{2}}}{\sqrt {{\frac {S_{1}^{2}}{n_{1}}}+{\frac {S_{2}^{2}}{n_{2}}}+{\frac {S_{1}^{4}}{2(n_{1}-1)}}+{\frac {S_{2}^{4}}{2(n_{2}-1)}}}}\right)}}
ご希望通りです。
2つの期待値の比率に MLEを単純に使用して 比率推定値を作成すると、 一貫性 はあるものの偏りのある点推定値につながること に注意してください (比率の推定値が対数正規分布であるという事実を使用します) [c] :
E
[
E
^
(
X
1
)
E
^
(
X
2
)
]
=
E
[
e
(
μ
^
1
−
μ
^
2
)
+
1
2
(
S
1
2
−
S
2
2
)
]
=
e
(
μ
1
−
μ
2
)
+
1
2
(
σ
1
2
−
σ
2
2
)
+
1
2
(
σ
1
2
n
1
+
σ
2
2
n
2
+
σ
1
4
2
(
n
1
−
1
)
+
σ
2
4
2
(
n
2
−
1
)
)
{\displaystyle E\left[{\frac {{\widehat {E}}(X_{1})}{{\widehat {E}}(X_{2})}}\right]=E\left[e^{({\widehat {\mu }}_{1}-{\widehat {\mu }}_{2})+{\frac {1}{2}}(S_{1}^{2}-S_{2}^{2})}\right]=e^{(\mu _{1}-\mu _{2})+{\frac {1}{2}}(\sigma _{1}^{2}-\sigma _{2}^{2})+{\frac {1}{2}}\left({\frac {\sigma _{1}^{2}}{n_{1}}}+{\frac {\sigma _{2}^{2}}{n_{2}}}+{\frac {\sigma _{1}^{4}}{2(n_{1}-1)}}+{\frac {\sigma _{2}^{4}}{2(n_{2}-1)}}\right)}}
自由パラメータを固定するエントロピーの極限原理 σ
応用においては、 決定すべきパラメータである。生産と散逸がバランスした成長過程の場合、シャノンエントロピーの極限原理を用いると、 [48]
σ
{\displaystyle \sigma }
σ
=
1
6
{\displaystyle \sigma ={\frac {1}{\sqrt {6}}}}
この値は、対数正規分布の変曲点と最大点の間のスケーリング関係を与えるために使用することができます。 [48] この関係は自然対数の底によって決定され 、最小表面エネルギー原理と幾何学的に類似しています。これらのスケーリング関係は、多くの成長プロセス(伝染病の拡大、液滴の飛散、人口増加、浴槽渦の旋回速度、言語特性の分布、乱流の速度プロファイルなど)を予測するのに役立ちます。たとえば、このような対数正規関数は、 液滴衝突中に二次的に生成される液滴のサイズ [49] や伝染病の拡大によく適合します。 [50]
e
=
2.718
…
{\displaystyle e=2.718\ldots }
σ
{\displaystyle \sigma }
この値 はドレイク方程式の確率的解を与えるために使用される。 [51]
σ
=
1
/
6
{\textstyle \sigma =1{\big /}{\sqrt {6}}}
発生と応用
対数正規分布は自然現象の記述において重要である。多くの自然成長過程は、対数スケールで加算される多数の小さなパーセンテージ変化の蓄積によって駆動される。適切な規則性条件下では、結果として生じる蓄積された変化の分布は、上記の「乗法中心極限定理」のセクションで述べたように、対数正規分布によってますますよく近似される。これは、 企業のために定式化したロバート・ジブラット(1904-1980)にちなんで、 ジブラットの法則としても知られている。 [52] これらの小さな変化の蓄積率が時間の経過とともに変化しない場合、成長はサイズとは無関係になる。この仮定が真実でなくても、時間の経過とともに成長するもののどの年齢でもサイズ分布は対数正規分布になる傾向がある。 [ 引用が必要 ] したがって、健康な個人の測定値の 基準範囲は、 平均を中心に対称分布を仮定するよりも、対数正規分布を仮定することによってより正確に推定される。 [ 引用が必要 ]
2 つ目の正当化は、基本的な自然法則が正の変数の乗算と除算を意味するという観察に基づいています。例としては、質量と距離を結果として生じる力と結び付ける単純な万有引力の法則、または抽出物と生成物の濃度を結び付ける溶液中の化学物質の平衡濃度の式などがあります。これらのケースでは、関係する変数の対数正規分布を仮定すると、一貫したモデルが得られます。
具体的な例は以下のサブセクションで示します。 [53] には、地質学、生物学、医学、食品、生態学などの分野における対数正規分布のレビューと表が含まれています。 [54]は 、神経科学における対数正規分布に関するレビュー記事で、注釈付きの参考文献が付いています。
人間の行動
インターネットのディスカッションフォーラムに投稿されたコメントの長さは対数正規分布に従う。 [55]
ユーザーのオンライン記事(ジョーク、ニュースなど)への滞在時間は対数正規分布に従います。 [56]
チェスの ゲームの長さは 対数正規分布に従う傾向がある。 [57]
標準刺激と一致する音響比較刺激の開始持続時間は対数正規分布に従う。 [18]
生物学と医学
生体組織の大きさの測定(長さ、皮膚面積、重量)。 [58]
病気の潜伏期間。 [59]
バナナの葉の斑点の直径、大麦のうどんこ病。 [53]
2003年のSARSのような感染力の高い伝染病の場合、公的介入による制御政策が関与すると、エントロピーを仮定し、標準偏差をエントロピー生成の最大速度の原理によって決定すると、入院患者数は自由パラメータのない対数正規分布を満たすことが示されている 。 [ 60]
生物標本の不活性な付属肢(毛、爪、爪、歯)の成長方向の長さ。 [ 引用が必要 ]
任意のゲノム領域の正規化された RNA-Seq リードカウントは、対数正規分布によって十分に近似できます。
PacBio シーケン シングのリード長は対数正規分布に従う。 [61]
特定の生理学的測定、例えば成人の血圧(男性と女性の集団に分けた後)。 [62]
Cmax 、 消失半減期 、 消失速度定数 など の いくつかの 薬物動態 変数。 [63]
神経科学では、ニューロン集団全体の発火率の分布は、多くの場合、ほぼ対数正規分布になります。これは、最初に皮質と線条体で観察され [64] 、その後、海馬と嗅内皮質で [65] 、脳の他の部分でも観察されました。 [54] [66] また、固有のゲイン分布とシナプス重み分布も対数正規分布であるように見えます [67] 。
大脳皮質のニューロン密度は、神経発達中のノイズの多い細胞分裂プロセスによるものである。 [68]
手術室管理では、手術時間 の分布が重要です 。
生きた細胞の細胞骨格の破壊の雪崩の大きさは対数正規分布を示し、癌細胞では健康な細胞よりも著しく大きい。 [69]
化学
粒子サイズ分布 と モル質量分布 。
鉱物中の希少元素の濃度。 [70]
アイスクリームの結晶の直径、マヨネーズの油滴、ココアプレスケーキの孔。 [53]
年間最大日降水量に当てはめた累積対数正規分布。 分布の当てはめを参照。
水文学
水文学 では 、対数正規分布は、日降雨量や河川流量の月間および年間の最大値などの変数の極端な値を分析するために使用されます。 [71]
右の画像は CumFreq で作成されたもので、対数正規分布を年間最大日降水量の順位付けに当てはめた例を示しており、 二項分布 に基づく90% 信頼区間 も表示しています。 [72]
降雨量データは、 累積頻度分析 の一部として 位置をプロットする ことによって表されます。
社会科学と人口統計
経済学 では、 人口の97%~99%の 所得 が対数正規分布しているという証拠がある。 [73] (高所得者の所得分布は パレート分布 に従う)。 [74]
所得分布が標準偏差 の対数正規分布に従う場合 、所得の不平等を評価するために一般的に使用される ジニ係数 は 、次の ように計算できます。ここで は 誤差関数 です。 なぜなら 、 は 標準正規分布の累積分布関数だからです。
σ
{\displaystyle \sigma }
G
=
erf
(
σ
2
)
{\displaystyle G=\operatorname {erf} \left({\frac {\sigma }{2}}\right)}
erf
{\displaystyle \operatorname {erf} }
G
=
2
Φ
(
σ
2
)
−
1
{\displaystyle G=2\Phi \left({\frac {\sigma }{\sqrt {2}}}\right)-1}
Φ
(
x
)
{\displaystyle \Phi (x)}
金融 、特に ブラック・ショールズ・モデル では、 為替レート、価格指数、株価指数の対数の変化は正規分布であると仮定されている [ 75 ] ( これらの変数は単利ではなく複利のように振舞うため、乗法的である)。しかし、 ブノワ・マンデルブロなどの数学者の中には、 重い裾を持つ 対数レヴィ分布 の方が、特に 株式市場の暴落 の分析にはより適切なモデルであると 主張している [76] 。確かに、株価分布は典型的には 太い裾 を示している。 [77]株式市場の暴落時の変化の太い裾の分布は、 中心極限定理 の仮定を無効にする 。
科学計量学 では 、ジャーナル論文や特許の引用数は離散対数正規分布に従う。 [78] [79]
都市の規模 (人口)はジブラの法則を満たす。 [80] 都市規模の成長過程は、規模に比例し不変である。したがって、 中心極限定理 から、都市規模の対数は正規分布する。
性的パートナーの数は対数正規分布によって最もよく説明されると思われる。 [81]
テクノロジー
信頼性 分析では 、対数正規分布は保守可能なシステムの修復時間をモデル化するためによく使用されます。 [82]
無線通信 では 、「dBやネーパーなどの対数値で表される局所平均電力は、正規分布(ガウス分布)を示す。」 [83]また、大きな建物や丘陵による無線信号のランダムな遮蔽は、 シャドウイング と呼ばれ 、対数正規分布としてモデル化されることが多い。
ボールミル粉砕 などのランダムな衝撃による粉砕によって生成される粒子サイズ分布 。 [84]
公開されている音声・動画データファイル( MIMEタイプ )のファイル サイズ分布は、5 桁の対数 正規分布に従います 。 [85]
1999年に収集されたWindows OS搭載のパソコン上の1億4000万ファイルのファイルサイズ。 [86] [55]
テキストベースの電子メール(1990年代)とマルチメディアベースの電子メール(2000年代)のサイズ。 [55]
コンピュータネットワークと インターネットトラフィック 分析では、対数正規分布は単位時間あたりのトラフィック量を表すのに適した統計モデルとして示されています。これは、実際のインターネットトレースの大規模なグループに堅牢な統計的アプローチを適用することによって示されています。この文脈では、対数正規分布は、2つの主な使用例で優れたパフォーマンスを示しています。(1) トラフィックが特定のレベルを超える時間の割合を予測する (サービスレベル契約またはリンク容量の推定用)、つまり帯域幅プロビジョニングに基づくリンクの寸法設定、(2) 95パーセンタイル価格の予測。 [87]
物理試験 において、 特定の条件下でのアイテムの故障までの時間を測定する場合、データは対数正規分布を使用して分析するのが最も効果的であることが多い。 [88] [89]
参照
注記
^ コックス法は、ランド(1971) [45] で「個人的なコミュニケーション」として引用されており、また、周とガオ(1997) [46] とオルソン2005 [43]の引用でも引用されている :セクション3.3
^ 問題は、それを直接行う方法がわからないことです。そのため、それらの対数を取り、 デルタ法を 使用して、それらの対数自体が(近似的に)正規であるとします。このトリックにより、それらの exp が対数正規分布であると仮定し、その近似を使用して CI を構築できます。RR の場合、ベース分布(つまり、対数を取る前)の中央値と平均は実際には同一であることに注意してください(元々は正規分布であり、対数正規分布ではないため)。たとえば、および です。 したがって、対数に基づいて CI を構築し、それを逆変換すると、が得られます 。したがって、CI が中央値に対してであると予想されますが、この場合は、実際には元の分布の平均に対しても CI になります。つまり、元が対数 正規分布であれば、 であると予想されます 。しかし、実際には であることがわかっています 。したがって、近似値は 2 番目のステップ(デルタ法)にありますが、CI は実際には期待値(中央値だけでなく)に対してです。これは、正規分布の基本分布から始めて、対数分布の後に別の近似値を使用して正規分布に戻しているためです。つまり、CI の近似値の大部分はデルタ法によるものです。
p
^
1
∼
˙
N
(
p
1
,
p
1
(
1
−
p
1
)
/
n
)
{\displaystyle {\hat {p}}_{1}{\dot {\sim }}N(p_{1},p_{1}(1-p1)/n)}
l
o
g
(
p
^
1
)
∼
˙
N
(
l
o
g
(
p
1
)
,
(
1
−
p
1
)
/
(
p
1
∗
n
)
)
{\displaystyle log({\hat {p}}_{1}){\dot {\sim }}N(log(p_{1}),(1-p1)/(p_{1}*n))}
C
I
(
p
1
)
:
e
l
o
g
(
p
^
1
)
±
(
1
−
p
^
1
)
/
(
p
^
1
∗
n
)
)
{\displaystyle CI(p_{1}):e^{log({\hat {p}}_{1})\pm (1-{\hat {p}}_{1})/({\hat {p}}_{1}*n))}}
p
^
1
{\displaystyle {\hat {p}}_{1}}
E
[
p
^
1
]
=
e
l
o
g
(
p
1
)
+
1
/
2
∗
(
1
−
p
1
)
/
(
p
1
∗
n
)
{\displaystyle E[{\hat {p}}_{1}]=e^{log(p_{1})+1/2*(1-p1)/(p_{1}*n)}}
E
[
p
^
1
]
=
e
l
o
g
(
p
1
)
=
p
1
{\displaystyle E[{\hat {p}}_{1}]=e^{log(p_{1})}=p_{1}}
^ バイアスは以下を使用することで部分的に最小限に抑えることができます。
[
E
(
X
1
)
E
(
X
2
)
]
^
=
[
E
^
(
X
1
)
E
^
(
X
2
)
]
2
(
σ
1
2
n
1
+
σ
2
2
n
2
+
σ
1
4
2
(
n
1
−
1
)
+
σ
2
4
2
(
n
2
−
1
)
)
^
=
[
e
(
μ
^
1
−
μ
^
2
)
+
1
2
(
S
1
2
−
S
2
2
)
]
2
S
1
2
n
1
+
S
2
2
n
2
+
S
1
4
2
(
n
1
−
1
)
+
S
2
4
2
(
n
2
−
1
)
{\displaystyle {\widehat {\left[{\frac {E(X_{1})}{E(X_{2})}}\right]}}=\left[{\frac {{\widehat {E}}(X_{1})}{{\widehat {E}}(X_{2})}}\right]{\frac {2}{\widehat {\left({\frac {\sigma _{1}^{2}}{n_{1}}}+{\frac {\sigma _{2}^{2}}{n_{2}}}+{\frac {\sigma _{1}^{4}}{2(n_{1}-1)}}+{\frac {\sigma _{2}^{4}}{2(n_{2}-1)}}\right)}}}=\left[e^{({\widehat {\mu }}_{1}-{\widehat {\mu }}_{2})+{\frac {1}{2}}(S_{1}^{2}-S_{2}^{2})}\right]{\frac {2}{{\frac {S_{1}^{2}}{n_{1}}}+{\frac {S_{2}^{2}}{n_{2}}}+{\frac {S_{1}^{4}}{2(n_{1}-1)}}+{\frac {S_{2}^{4}}{2(n_{2}-1)}}}}}
参考文献
^ Norton, Matthew; Khokhlov, Valentyn; Uryasev, Stan (2019). 「ポートフォリオ最適化と密度推定への応用による共通確率分布のCVaRとbPOEの計算」 (PDF) . Annals of Operations Research . 299 (1–2). Springer: 1281–1315. arXiv : 1811.11301 . doi :10.1007/s10479-019-03373-1. S2CID 254231768. 2021-04-18にオリジナルから アーカイブ (PDF) . 2023-02-27 に取得– stonybrook.edu経由。
^ abcd Weisstein, Eric W. 「対数正規分布」。mathworld.wolfram.com 。 2020年9月13日 閲覧 。
^ ab "1.3.6.6.9. 対数正規分布". www.itl.nist.gov . 米国 国立標準技術研究所 (NIST) . 2020年9月13日 閲覧 。
^ abcde ジョンソン、ノーマン L.; コッツ、サミュエル; バラクリシュナン、N. (1994)、「14: 対数正規分布」、 連続一変量分布。第 1 巻 、ワイリー確率および数理統計シリーズ: 応用確率および統計 (第 2 版)、ニューヨーク: ジョン ワイリー & サンズ 、 ISBN 978-0-471-58495-7 、 MR 1299979
^ Park, Sung Y.; Bera, Anil K. (2009). 「最大エントロピー自己回帰条件付き異分散モデル」 (PDF) . Journal of Econometrics . 150 (2): 219–230, 特に表1、p. 221. CiteSeerX 10.1.1.511.9750 . doi :10.1016/j.jeconom.2008.12.014. 2016-03-07に オリジナル (PDF)からアーカイブ 。2011-06-02 に取得 。
^ Tarmast, Ghasem (2001). 多変量対数正規分布 (PDF) . ISI Proceedings: 53rd Session. ソウル. 2013-07-19 のオリジナルからアーカイブ (PDF) 。
^ Halliwell, Leigh (2015). The Lognormal Random Multivariate (PDF) . Casualty Actuarial Society E-Forum、2015年春。バージニア州アーリントン。 2015年9月30日時点のオリジナルよりアーカイブ (PDF) 。
^ Heyde, CC. (2010)、「対数正規分布の特性について」、 Journal of the Royal Statistical Society、シリーズ B 、第 25 巻、第 2 号、pp. 392–393、 doi : 10.1007/978-1-4419-5823-5_6 、 ISBN 978-1-4419-5822-8
^ ビリングスリー、パトリック(2012年)。確率と測定(記念版)。ホーボーケン、ニュージャージー:ワイリー 。p.415。ISBN 978-1-118-12237-2 . OCLC 780289503.
^ ab Holgate, P. (1989). 「対数正規特性関数、第18巻、pp. 4539–4548、1989」。 統計 における 通信 - 理論と方法 。18 (12): 4539–4548。doi :10.1080/03610928908830173。
^ Barakat, R. (1976). 「独立した対数正規分布ランダム変数の合計」. Journal of the Optical Society of America . 66 (3): 211–216. Bibcode :1976JOSA...66..211B. doi :10.1364/JOSA.66.000211.
^ Barouch, E.; Kaufman, GM.; Glasser, ML. (1986). 「対数正規分布のランダム変数の合計について」 (PDF) . 応用数学研究 . 75 (1): 37–55. doi :10.1002/sapm198675137. hdl : 1721.1/48703 .
^ Leipnik, Roy B. (1991年1月). 「対数正規分布のランダム変数について: I – 特性関数」 (PDF) . オーストラリア数学会誌、シリーズB. 32 ( 3): 327–347. doi : 10.1017/S0334270000006901 .
^ S. Asmussen、JL Jensen、L. Rojas-Nandayapa (2016)。「対数正規分布のラプラス変換について」、応用確率論および計算18(2)、441-458。Thieleレポート6(13)。
^ abc Das, Abhranil (2021). 「正規分布を積分および分類する方法」. Journal of Vision . 21 (10): 1. arXiv : 2012.14331 . doi : 10.1167/jov.21.10.1. PMC 8419883. PMID 34468706.
^ abカークウッド、トーマスBL(1979 年12 月 )。「幾何平均と分散の測定」 バイオメトリクス 。35 (4):908-9。JSTOR 2530139。
^ Limpert, E; Stahel, W; Abbt, M (2001). 「科学分野における対数正規分布:鍵と手がかり」. BioScience . 51 (5): 341–352. doi : 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2 .
^ ab Heil P, Friedrich B (2017). 「音響刺激のオンセット-持続時間のマッチングの再考:従来の算術と提案された精度と精密度の幾何学的尺度」。Frontiers in Psychology . 7 : 2013. doi : 10.3389/fpsyg.2016.02013 . PMC 5216879 . PMID 28111557.
^ Sawant, S.; Mohan, N. (2011)「FAQ: SAS を使用した臨床試験データの有効性分析に関する問題」 Wayback Machine に 2011 年 8 月 24 日にアーカイブ 、 PharmaSUG2011 、論文 PO08
^ Schiff, MH; et al. (2014). 「関節リウマチ患者におけるメトトレキサートの経口投与と皮下投与の直接比較ランダム化クロスオーバー試験:15 mg以上の用量での経口メトトレキサートの薬剤曝露限界は皮下投与で克服できる可能性がある」 Ann Rheum Dis . 73 (8): 1–3. doi :10.1136/annrheumdis-2014-205228. PMC 4112421. PMID 24728329 .
^ Daly, Leslie E.; Bourke, Geoffrey Joseph (2000). 医療統計の解釈と利用 . 第46巻(第5版). オックスフォード、イギリス:Wiley-Blackwell. p. 89. doi :10.1002/9780470696750. ISBN 978-0-632-04763-5 . PMC 1059583 ; 印刷版。オンライン電子書籍 ISBN 9780470696750
^ 「ProbOnto」 。 2017年 7月1日 閲覧 。
^ Swat, MJ; Grenon, P; Wimalaratne, S (2016). 「ProbOnto: 確率分布のオントロジーと知識ベース」. バイオインフォマティクス . 32 (17): 2719–21. doi :10.1093/bioinformatics/btw170. PMC 5013898. PMID 27153608 .
^ ab Forbes et al. Probability Distributions (2011)、John Wiley & Sons、Inc.
^ Lunn, D. (2012). BUGSブック:ベイズ分析の実践的入門。統計科学テキスト。CRCプレス。
^ Limpert, E.; Stahel, WA; Abbt, M. (2001). 「科学分野における対数正規分布:鍵と手がかり」. BioScience . 51 (5): 341–352. doi : 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2 .
^ Nyberg, J.; et al. (2012). 「PopED - 拡張された並列化された集団最適設計ツール」 Comput Methods Programs Biomed . 108 (2): 789–805. doi :10.1016/j.cmpb.2012.05.005. PMID 22640817.
^ Retout, S; Duffull, S; Mentré, F (2001). 「集団薬物動態設計の評価のための集団フィッシャー情報マトリックスの開発と実装」 Comp Meth Pro Biomed . 65 (2): 141–151. doi :10.1016/S0169-2607(00)00117-6. PMID 11275334.
^ PopED 開発チーム (2016)。PopED マニュアル、リリース バージョン 2.13。技術レポート、ウプサラ大学。
^ ProbOnto Webサイト、URL: http://probonto.org
^ Cococcioni, Marco; Fiorini, Francesco; Pagano, Michele (2023-04-06). 「数値的に検証可能な無限分散を持つ対数正規分布を使用したヘビーテール現象のモデル化」. 数学 . 11 (7): 1758. doi : 10.3390/math11071758 . hdl : 11568/1216554 . ISSN 2227-7390.
^ Damgaard, Christian; Weiner, Jacob (2000). 「植物の大きさや繁殖力の不平等の記述」. 生態学 . 81 (4): 1139–1142. doi :10.1890/0012-9658(2000)081[1139:DIIPSO]2.0.CO;2.
^ Rossman, Lewis A (1990 年 7 月). 「調和平均に基づく設計流量」. Journal of Hydraulic Engineering . 116 (7): 946–950. doi :10.1061/(ASCE)0733-9429(1990)116:7(946).
^ Thorin, Olof (1977). 「対数正規分布の無限分割可能性について」. Scandinavian Actuarial Journal . 1977 (3): 121–148. doi :10.1080/03461238.1977.10405635. ISSN 0346-1238.
^ ab Gao, Xin (2009). 「相関対数正規変数の合計に対する裾密度の漸近挙動」. 国際数学・数理科学ジャーナル . 2009 :1–28. doi : 10.1155/2009/630857 .
^ Asmussen, S.; Rojas-Nandayapa, L. (2008). 「ガウスコピュラによる対数正規確率変数の合計の漸近解析」 (PDF) . 統計および確率論文集 . 78 (16): 2709–2714. doi :10.1016/j.spl.2008.03.035.
^ Marlow, NA. (1967 年 11 月). 「独立した正規ランダム変数のべき乗和に対する正規極限定理」. Bell System Technical Journal . 46 (9): 2081–2089. doi :10.1002/j.1538-7305.1967.tb04244.x.
^ Botev, ZI; L'Ecuyer, P. (2017). 「従属対数正規変量の右裾の正確な計算」 2017 Winter Simulation Conference (WSC)、2017年12月3日~6日 。米国ネバダ州ラスベガス:IEEE。pp. 1880~1890。arXiv : 1705.03196。doi : 10.1109 / WSC.2017.8247924。ISBN 978-1-5386-3428-8 。
^ Asmussen, A.; Goffard, P.-O.; Laub, PJ (2016). 「正規直交多項式展開と対数正規和密度」. arXiv : 1601.01763v1 [math.PR].
^ Sangal, B.; Biswas, A. ( 1970). 「水文学における3パラメータ対数正規分布の応用」。 水資源研究 。6 (2): 505–515. doi :10.1029/WR006i002p00505。
^ Johnson, NL (1949). 「翻訳方法によって生成される頻度曲線のシステム」. Biometrika . 36 (1/2): 149–176. doi :10.2307/2332539. JSTOR 2332539. PMID 18132090.
^ Swamee, PK (2002). 「対数正規分布に近い」. 水文工学ジャーナル . 7 (6): 441–444. doi :10.1061/(ASCE)1084-0699(2002)7:6(441).
^ abc Olsson, Ulf. 「対数正規分布の平均の信頼区間」Journal of Statistics Education 13.1 (2005).pdf html
^ user10525、対数正規分布データセットの平均の信頼区間を計算するにはどうすればよいですか?、URL (バージョン: 2022-12-18): https://stats.stackexchange.com/q/33395
^ Land, CE (1971)、「正規平均と分散の線形関数の信頼区間」、Annals of Mathematical Statistics、42、1187–1205。
^ Zhou, XH.、Gao, S. (1997)、「対数正規平均の信頼区間」、Statistics in Medicine、16、783-790。
^ リスク比とオッズ比の信頼区間
^ ab Wu, Ziniu; Li, Juan; Bai, Chenyuan (2017). 「エントロピーの極限原理による対数正規型成長プロセスのスケーリング関係」. エントロピー . 19 (56): 1–14. Bibcode :2017Entrp..19...56W. doi : 10.3390/e19020056 .
^ Wu, Zi-Niu (2003). 「 固体壁に衝突する液滴のクラウンスプラッシュによる二次放出液滴のサイズ分布の予測」。 確率 論的工学力学 。18 (3): 241–249。Bibcode :2003PEngM..18..241W。doi : 10.1016/S0266-8920(03)00028-6。
^ Wang, WenBin; Wu, ZiNiu; Wang, ChunFeng; Hu, RuiFeng (2013). 「エントロピーベースの熱力学モデルによる制御された伝染病の拡散率のモデル化」. Science China 物理、力学、天文学 . 56 (11): 2143–2150. arXiv : 1304.5603 . Bibcode :2013SCPMA..56.2143W. doi :10.1007/s11433-013-5321-0. ISSN 1674-7348. PMC 7111546. PMID 32288765 .
^ Bloetscher, Frederick (2019). 「予測ベイジアンモンテカルロ-マルコフ連鎖法を用いたドレイク方程式の確率的解の提供」 Acta Astronautica . 155 :118–130. Bibcode :2019AcAau.155..118B. doi :10.1016/j.actaastro.2018.11.033. S2CID 117598888.
^サットン、ジョン(1997 年3月) 。 「ジブラットの遺産」。 経済文献ジャーナル 。32 (1):40-59。JSTOR 2729692。
^ abc Limpert, Eckhard; Stahel, Werner A.; Abbt, Markus (2001). 「科学分野における対数正規分布:鍵と手がかり」. BioScience . 51 (5): 341. doi : 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2 . ISSN 0006-3568.
^ ab Buzsáki, György; Mizuseki, Kenji (2017-01-06). 「対数動的脳:歪んだ分布がネットワーク操作に与える影響」. Nature Reviews. Neuroscience . 15 (4): 264–278. doi :10.1038/nrn3687. ISSN 1471-003X. PMC 4051294. PMID 24569488 .
^ abc Pawel, Sobkowicz; et al. (2013). 「インターネットディスカッションにおけるユーザー投稿の長さの対数正規分布 - ウェーバー・フェヒナーの法則の結果か?」 EPJ Data Science 。
^ Yin, Peifeng; Luo, Ping; Lee, Wang-Chien; Wang, Min (2013)。沈黙もまた証拠となる: 心理学的観点からの推薦のための滞在時間の解釈。ACM International Conference on KDD。
^ 「チェスのゲームの平均時間はどのくらいですか?」 chess.stackexchange.com 。 2018年 4月14日 閲覧 。
^ ハクスリー、ジュリアン・S.(1932年)。 相対的成長の問題 。ロンドン 。ISBN 978-0-486-61114-3 . OCLC 476909537.
^ サートウェル、フィリップ E.「感染症の潜伏期間の分布」 アメリカ衛生学誌 51 (1950): 310-318。
^ SK Chan, Jennifer; Yu, Philip LH (2006). 「閾値幾何プロセスを用いたSARSデータのモデリング」。 医学統計 . 25 (11): 1826–1839. doi :10.1002/sim.2376. PMID 16345017. S2CID 46599163.
^ 小野 幸輝; 浅井 潔; 浜田 道明 (2013-01-01). 「PBSIM: PacBio 読み取りシミュレーター - 正確なゲノムアセンブリに向けて」. バイオ インフォマティクス . 29 (1): 119–121. doi : 10.1093/bioinformatics/bts649 . ISSN 1367-4803. PMID 23129296.
^ Makuch, Robert W.; DH Freeman; MF Johnson (1979). 「血圧モデルとしての対数正規分布の正当性」. Journal of Chronic Diseases . 32 (3): 245–250. doi :10.1016/0021-9681(79)90070-5. PMID 429469.
^ Lacey, LF; Keene, ON; Pritchard, JF; Bye, A. (1997-01-01). 「一般的な非コンパートメント薬物動態変数:それらは正規分布か対数正規分布か?」 Journal of Biopharmaceutical Statistics . 7 (1): 171–178. doi :10.1080/10543409708835177. ISSN 1054-3406. PMID 9056596.
^ Scheler, Gabriele; Schumann, Johann (2006-10-08). 神経出力率の多様性と安定性 。第36回神経科学会議、アトランタ。
^ 水関 健二; Buzsáki György (2013-09-12). 「海馬と嗅内皮質における発火率の事前設定された偏った分布」. Cell Reports . 4 (5): 1010–1021. doi :10.1016/j.celrep.2013.07.039. ISSN 2211-1247. PMC 3804159. PMID 23994479 .
^ Wohrer, Adrien; Humphries, Mark D.; Machens, Christian K. (2013-04-01). 「知覚的意思決定中の神経活動の集団全体分布」。Progress in Neurobiology . 103 : 156–193. doi :10.1016/j.pneurobio.2012.09.004. ISSN 1873-5118. PMC 5985929. PMID 23123501 .
^ Scheler, Gabriele (2017-07-28). 「対数分布は内在的学習がヘブ分布であることを証明する」 F1000Research . 6 : 1222. doi : 10.12688/f1000research.12130.2 . PMC 5639933 . PMID 29071065.
^ Morales-Gregorio, Aitor; van Meegen, Alexander; van Albada, Sacha (2023). 「哺乳類大脳皮質におけるニューロン密度の普遍的対数正規分布」. 大脳皮質 . 33 (16): 9439–9449. doi :10.1093/cercor/bhad160. PMC 10438924. PMID 37409647 .
^ ポリッツィ、ステファノ;ラペルーザ、バスティアン。ペレス・レチェ、フランシスコ・J;ニコリーニ、フランク E;サッタ、ヴェロニク・マグエル。アルネオド、アラン。アルグール、フランソワーズ (2018-05-29)。 「生細胞可塑性のための最小破壊カスケードモデル」。 新しい物理学ジャーナル 。 20 (5): 053057。 ビブコード :2018NJPh...20e3057P。 土井 :10.1088/1367-2630/aac3c7。 hdl : 2164/10561 。 ISSN 1367-2630。
^ Ahrens, LH (1954-02-01). 「元素の対数正規分布(地球化学の基本法則とその補助法則)」. Geochimica et Cosmochimica Acta . 5 (2): 49–73. Bibcode :1954GeCoA...5...49A. doi :10.1016/0016-7037(54)90040-X. ISSN 0016-7037.
^ Oosterbaan, RJ (1994). 「6: 頻度と回帰分析」 (PDF) 。 Ritzema, HP (編) 排水原理と応用、出版物 16。 オランダ、ワーゲニンゲン: 国際土地再生・改善研究所 (ILRI)。 pp. 175–224。 ISBN 978-90-70754-33-4 。
^ CumFreq、分布フィッティング用無料ソフトウェア
^ クレメンティ、ファビオ、 ガレガティ、マウロ (2005)「パレートの所得分配法則: ドイツ、イギリス、アメリカの証拠」、EconWPA
^ 相馬渉 (2002-02-22). 「個人所得の物理学」. 高安秀樹 (編) 『金融変動の経験科学:経済物理学の到来』 . Springer. arXiv : cond-mat/0202388 . doi :10.1007/978-4-431-66993-7.
^ブラック、F.; ショールズ、M. ( 1973)。「オプションと企業 負債 の価格設定」。 政治経済ジャーナル 。81 (3): 637。doi :10.1086/260062。S2CID 154552078。
^ マンデルブロ、ブノワ (2004)。市場の(誤った)行動。ベーシックブックス 。ISBN 9780465043552 。
^ Bunchen, P., Advanced Option Pricing 、シドニー大学コースブック、2007年
^ Thelwall, Mike; Wilson, Paul (2014). 「引用データの回帰分析: さまざまな手法の評価」 Journal of Informetrics . 8 (4): 963–971. arXiv : 1510.08877 . doi :10.1016/j.joi.2014.09.011. S2CID 8338485.
^ シェリダン、ポール; 小野寺、拓 (2020). 「優先的愛着のパラドックス: 優先的愛着と成長が組み合わさって、対数正規入次数分布を持つネットワークが生成される仕組み」. Scientific Reports . 8 (1): 2811. arXiv : 1703.06645 . doi :10.1038/s41598-018-21133-2. PMC 5809396. PMID 29434232.
^ Eeckhout, Jan (2004). 「(すべての)都市に対するジブラの法則」 アメリカ経済評論 94 (5): 1429–1451. doi : 10.1257/0002828043052303. JSTOR 3592829 – JSTOR経由。
^ Kault, David (1996). 「性的パートナー数の分布の形状」。 医学統計 。15 (2): 221–230。doi : 10.1002 /(SICI)1097-0258(19960130)15:2<221::AID-SIM148>3.0.CO;2- Q。PMID 8614756 。
^ O'Connor, Patrick; Kleyner, Andre (2011). 実践的信頼性工学 . John Wiley & Sons. p. 35. ISBN 978-0-470-97982-2 。
^ 「Shadowing」 www.WirelessCommunication.NL . 2012年1月13日時点のオリジナルよりアーカイブ。
^ Dexter, AR; Tanner, DW (1972 年 7 月). 「対数正規分布を持つ球体の混合物の充填密度」. Nature Physical Science . 238 (80): 31–32. Bibcode :1972NPhS..238...31D. doi :10.1038/physci238031a0. ISSN 2058-1106.
^ Gros, C; Kaczor, G.; Markovic, D (2012). 「地球規模での人間のデータ生成に対する神経心理学的制約」. The European Physical Journal B . 85 (28): 28. arXiv : 1111.6849 . Bibcode :2012EPJB...85...28G. doi :10.1140/epjb/e2011-20581-3. S2CID 17404692.
^ Douceur, John R.; Bolosky, William J. (1999-05-01). 「ファイルシステムの内容に関する大規模な調査」. ACM SIGMETRICS パフォーマンス評価レビュー . 27 (1): 59–70. doi : 10.1145/301464.301480 . ISSN 0163-5999.
^ Alamsar, Mohammed; Parisis, George; Clegg, Richard; Zakhleniuk, Nickolay (2019). 「インターネットのトラフィック量の分布とその影響について」 arXiv : 1902.03853 [cs.NI].
^ ASTM D3654、感圧テープのせん断接着力の標準試験方法w
^ ASTM D4577、一定荷重下での容器の圧縮抵抗の標準試験方法>\
さらに読む
クロウ、エドウィン L.; 清水 邦夫編 (1988)、 対数正規分布、理論と応用 、統計学:教科書とモノグラフ、第 88 巻、ニューヨーク: マルセル デッカー社、pp. xvi+387、 ISBN 978-0-8247-7803-3 、 MR 0939191、 Zbl 0644.62014
Aitchison, J. および Brown, JAC (1957) 『対数正規分布』 、ケンブリッジ大学出版局。
Limpert, E; Stahel, W; Abbt, M (2001). 「科学分野における対数正規分布:鍵と手がかり」. BioScience . 51 (5): 341–352. doi : 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2 .
Holgate, P. (1989). 「対数正規特性関数」. 統計におけるコミュニケーション - 理論と方法 . 18 (12): 4539–4548. doi :10.1080/03610928908830173.
ブルックス、ロバート、コーソン、 ドナル、ウェールズ ( 1994)。「原資産がすべて対数正規分布に従う場合のインデックス オプションの価格設定」。 先物 および オプション研究の進歩 。7。SSRN 5735。
外部リンク
ウィキメディア・コモンズには、対数正規分布 に関連するメディアがあります 。