logo

ラプラス分布の平均と分散 📂確率分布論

ラプラス分布の平均と分散

公式

X∼X \sim Laplace⁡(μ,b)\operatorname{Laplace}(\mu, b)のとき、XXの平均と分散はそれぞれ次のようになる。

E(X)=μ E(X) = \mu Var⁡(X)=2b2 \Var(X) = 2b^{2}

証明

直接計算

期待値の定義と部分積分法によって、

E(X)=∫−∞∞x12be−∣x−μ∣/bdx=∫−∞μx2be(x−μ)/bdx+∫μ∞x2be−(x−μ)/bdx=([x2b(be(x−μ)/b)]−∞μ−∫−∞μ12e(x−μ)/bdx)+([x2b(−be−(x−μ)/b)]μ∞+∫μ∞12e−(x−μ)/bdx)=(μ2−[b2e(x−μ)/b]−∞μ)+(μ2+[−b2e−(x−μ)/b]μ∞)=(μ2−b2)+(μ2+b2)=μ \begin{align*} E(X) &= \int\limits_{-\infty}^{\infty} x \dfrac{1}{2b} e^{-|x - \mu|/b} dx \\ &= \int\limits_{-\infty}^{\mu} \dfrac{x}{2b} e^{(x - \mu)/b} dx + \int\limits_{\mu}^{\infty} \dfrac{x}{2b} e^{-(x - \mu)/b} dx \\ &= \left( \left[\dfrac{x}{2b}\left(b e^{(x - \mu)/b}\right)\right]_{-\infty}^{\mu} - \int\limits_{-\infty}^{\mu} \dfrac{1}{2} e^{(x - \mu)/b} dx \right) \\ &\qquad+ \left( \left[\dfrac{x}{2b}\left(-b e^{-(x - \mu)/b}\right)\right]_{\mu}^{\infty} + \int\limits_{\mu}^{\infty} \dfrac{1}{2} e^{-(x - \mu)/b} dx \right) \\ &= \left( \dfrac{\mu}{2} - \left[ \dfrac{b}{2} e^{(x - \mu)/b} \right]_{-\infty}^{\mu} \right) + \left( \dfrac{\mu}{2} + \left[ -\dfrac{b}{2} e^{-(x-\mu)/b} \right]_{\mu}^{\infty} \right) \\ &= \left( \dfrac{\mu}{2} - \dfrac{b}{2} \right) + \left( \dfrac{\mu}{2} + \dfrac{b}{2} \right) \\ &= \mu \end{align*}

分散を得るためにE(X2)E(X^{2})を計算しよう。

E(X2)=∫−∞∞x212be−∣x−μ∣/bdx=∫−∞μx22be(x−μ)/bdx+∫μ∞x22be−(x−μ)/bdx=([x22(e(x−μ)/b)]−∞μ−∫−∞μxe(x−μ)/bdx)+([−x22e−(x−μ)/b]μ∞+∫μ∞xe−(x−μ)/bdx)=[μ22−(μb−b2)]+[μ22+(μb−b2)]=μ2+2b2 \begin{align*} E(X^{2}) &= \int\limits_{-\infty}^{\infty} x^{2} \dfrac{1}{2b} e^{-|x - \mu|/b} dx \\ &= \int\limits_{-\infty}^{\mu} \dfrac{x^{2}}{2b} e^{(x - \mu)/b} dx + \int\limits_{\mu}^{\infty} \dfrac{x^{2}}{2b} e^{-(x - \mu)/b} dx \\ &= \left( \left[\dfrac{x^{2}}{2}\left( e^{(x - \mu)/b}\right)\right]_{-\infty}^{\mu} - \int\limits_{-\infty}^{\mu} x e^{(x - \mu)/b} dx \right) \\ &\qquad+ \left( \left[ - \dfrac{x^{2}}{2} e^{-(x - \mu)/b}\right]_{\mu}^{\infty} + \int\limits_{\mu}^{\infty} x e^{-(x - \mu)/b} dx \right) \\ &= \left[ \dfrac{\mu^{2}}{2} - (\mu b - b^{2}) \right] + \left[ \dfrac{\mu^{2}}{2} + (\mu b - b^{2}) \right] \\ &= \mu^{2} + 2b^{2} \end{align*}

三つ目の等号は部分積分法、四つ目の等号は期待値で計算した積分の結果を利用したものである。分散はVar⁡(X)=E(X2)−E(X)2\Var(X) = E(X^{2}) - E(X)^{2}なので、

Var⁡(X)=E(X2)−E(X)2=μ2+2b2−μ2=2b2 \Var(X) = E(X^{2}) - E(X)^{2} = \mu^{2} + 2b^{2} - \mu^{2} = 2b^{2}

■

積率生成函数より

ラプラス分布の積率生成函数は次のようになる。

m(t)=11−b2t2eμtfor ∣t∣<1b m(t) = \dfrac{1}{1 - b^{2}t^{2}} e^{\mu t} \qquad \text{for } |t| < \dfrac{1}{b}

期待値はm′(0)m^{\prime}(0)なので微分してみると、

m′(t)=ddt(11−b2t2)eμt+11−b2t2ddteμt=dd(1−b2t2)(11−b2t2)d(1−b2t2)dteμt+11−b2t2μeμt=2b2t(1−b2t2)2eμt+μ1−b2t2eμt \begin{align*} m^{\prime}(t) &= \dfrac{d}{dt} \left( \dfrac{1}{1 - b^{2}t^{2}} \right) e^{\mu t} + \dfrac{1}{1 - b^{2}t^{2}} \dfrac{d}{dt} e^{\mu t} \\ &= \dfrac{d}{d(1-b^{2}t^{2})} \left( \dfrac{1}{1 - b^{2}t^{2}} \right) \dfrac{d(1-b^{2}t^{2})}{dt} e^{\mu t} + \dfrac{1}{1 - b^{2}t^{2}} \mu e^{\mu t} \\ &= \dfrac{2b^{2}t}{(1 - b^{2}t^{2})^{2}} e^{\mu t} + \dfrac{\mu}{1 - b^{2}t^{2}} e^{\mu t} \end{align*}

したがって次を得る。

E(X)=m′(0)=μ E(X) = m^{\prime}(0) = \mu

分散を求めるためにもう一度微分してみると、

m′′(t)=ddt(2b2t(1−b2t2)2eμt+μ1−b2t2eμt) m^{\prime\prime}(t) = \dfrac{d}{dt} \left( \dfrac{2b^{2}t}{(1 - b^{2}t^{2})^{2}} e^{\mu t} + \dfrac{\mu}{1 - b^{2}t^{2}} e^{\mu t} \right)

二番目の項の微分は上の結果からすぐに導かれる。

m′′(t)=ddt(2b2t(1−b2t2)2eμt)+2μb2t(1−b2t2)2eμt+μ21−b2t2eμt m^{\prime\prime}(t) = \dfrac{d}{dt} \left( \dfrac{2b^{2}t}{(1 - b^{2}t^{2})^{2}} e^{\mu t} \right) + \dfrac{2\mu b^{2}t}{(1 - b^{2}t^{2})^{2}} e^{\mu t} + \dfrac{\mu^{2}}{1 - b^{2}t^{2}} e^{\mu t}

一番目の項を微分すると2b2(1−b2t2)2eμt\dfrac{2b^{2}}{(1-b^{2}t^{2})^{2}}e^{\mu t}とt=0t=0を代入したとき00になる項が出てくることがわかる。したがって、

m′′(0)=2b2(1−b2t2)2eμt+2μb2t(1−b2t2)2eμt+μ21−b2t2eμt∣t=0=2b2+μ2 \begin{align*} m^{\prime\prime}(0) &= \left. \dfrac{2b^{2}}{(1-b^{2}t^{2})^{2}}e^{\mu t} + \dfrac{2\mu b^{2}t}{(1 - b^{2}t^{2})^{2}} e^{\mu t} + \dfrac{\mu^{2}}{1 - b^{2}t^{2}} e^{\mu t} \right|_{t=0} \\ &= 2b^{2} + \mu^{2} \\ \end{align*}

よって分散は次のとおりである。

Var⁡(X)=E(X2)−(E(X))2=m′′(0)−(m′(0))2=2b2 \Var(X) = E(X^{2}) - (E(X))^{2} = m^{\prime\prime}(0) - (m^{\prime}(0))^{2} = 2b^{2}

■