내가 어떻게 계산

가정하자 및 밀도 함수 표준 정규 분포의 분포 함수이다. $\phi(\cdot)$ $\Phi(\cdot)$

적분을 어떻게 계산할 수 있습니까?

\int_{- \infty}^{\infty} Φ (\frac{w - a}{b}) ϕ (w) d w

$\int^{\infty}_{-\infty}\Phi\left(\frac{w-a}{b}\right)\phi(w)\,\mathrm dw$

mathematical-statistics normal-distribution integral

— 하디 산지
소스

이건 다 괜찮아 이것을 포함하는보다 일반적인 결과에 대한 초기 언급은 Ellison (1964, J.Am.Stat.Assoc, 59, 89-95); 정리 2의 목록 1을 참조하십시오.

답변:

보다 일반적인 표기법은

y (μ, σ) = \int Φ (\frac{x - μ}{σ}) ϕ (x) d x = Φ (\frac{- μ}{\sqrt{1 + σ^{2}}}) .

$y(\mu, \sigma) = \int\Phi\left(\frac{x-\mu}{\sigma}\right)\phi(x) dx = \Phi\left(\frac{-\mu}{\sqrt{1+\sigma^2}}\right).$

이것은 와 대해 적분을 차별화하여 닫힌 형태로 표현 될 수있는 기본 적분을 생성 함으로써 찾을 수 있습니다. $\mu$ $\sigma$

\frac{\partial y}{\partial μ} (μ, σ) = - \frac{1}{\sqrt{2 π} \sqrt{σ^{2} + 1}} e^{- \frac{1}{2} \frac{μ^{2}}{σ^{2} + 1}},

$\frac{\partial y}{\partial \mu}(\mu, \sigma) = -\frac{1}{\sqrt{2 \pi } \sqrt{\sigma ^2+1}}e^{-\frac{1}{2}\frac{\mu ^2}{\sigma ^2+1}},$

\frac{\partial y}{\partial σ} (μ, σ) = \frac{μ σ}{\sqrt{2 π} {(σ^{2} + 1)}^{3 / 2}} e^{- \frac{1}{2} \frac{μ^{2}}{σ^{2} + 1}} .

$\frac{\partial y}{\partial \sigma}(\mu, \sigma) = \frac{\mu\sigma }{\sqrt{2 \pi } \left(\sigma ^2+1\right)^{3/2}}e^{-\frac{1}{2}\frac{\mu ^2}{\sigma ^2+1}}.$

이 시스템은 초기 상태로 시작 통합 될 수 = = , (쉽게 분화에 의해 확인된다) 소정의 용액을 얻었다. $y(0,1)$ $\int\Phi(x)\phi(x)dx$ $1/2$

— 우버
소스

나는 숫자 통합을 통해 답을 재확인와의 비율을 윤곽

계약이 범위에서 열한 유효 숫자였다.

- 2 \leq μ \leq 2

$-2 \le \mu \le 2$

0 < σ \leq 2

$0 \lt \sigma \le 2$

— whuber

와우, 영리한 해결책.

— Cam.Davidson.Pilon 2016

나는이 검사가 거의 수행 될 수 있다고 생각합니다. 적분 아래의 첫 번째 항은 균일 한 [0,1] 랜덤 변수입니다. 일반 pdf는 대칭이므로 적분은

이어야합니다.

\frac{1}{2}

$1 \over 2$

— soakley

@soakley 귀하의 접근 방식은

에 적용되지만 다른

인수에 어떻게 적용되는지는 확실하지 않습니다 .

y (0, 1)

$y(0,1)$

y

$y$

— whuber

@whuber 이해하지 못해서 미안하지만, 미분과 초기 조건에 대한 두 개의 닫힌 양식이 있으면 거기에서 최종 솔루션으로 어떻게 이동합니까? 다시 말해, 미분과 초기 조건에 대한 닫힌 양식 표현식으로 무엇을 했습니까?

— user106860

하자 및 독립적 정상 랜덤 변수 일 및 표준 정규 확률 변수. 그런 다음 $X$ $Y$ $X \sim N(a,b^2)$ $Y$ 따라서 총 확률 법칙을 사용하여

P {X \leq Y ∣ Y = w} = P {X \leq w} = Φ (\frac{w - a}{b}) .

$P\{X \leq Y \mid Y = w\} = P\{X \leq w\} = \Phi\left(\frac{w-a}{b}\right).$

이제,

으로 표현 될 수

에 주목하여 그

P {X \leq Y} = \int_{- \infty}^{\infty} P {X \leq Y ∣ Y = w} ϕ (w) d w = \int_{- \infty}^{\infty} Φ (\frac{w - a}{b}) ϕ (w) d w .

$P\{X \leq Y\} = \int_{-\infty}^\infty P\{X \leq Y \mid Y = w\}\phi(w)\,\mathrm dw = \int_{-\infty}^\infty \Phi\left(\frac{w-a}{b}\right)\phi(w)\,\mathrm dw.$

P {X \leq Y} = P {X - Y \leq 0}

$P\{X \leq Y\} = P\{X-Y \leq 0\}$

Φ (\cdot)

$\Phi(\cdot)$

, 따라서 우리는 얻을

X - Y \sim N (a, b^{2} + 1)

$X-Y \sim N(a,b^2+1)$

whuber의 답변 결과와 동일합니다.

\int_{- \infty}^{\infty} Φ (\frac{w - a}{b}) ϕ (w) d w = Φ (\frac{- a}{\sqrt{b^{2} + 1}})

$\int_{-\infty}^\infty \Phi\left(\frac{w-a}{b}\right)\phi(w)\,\mathrm dw = \Phi\left(\frac{-a}{\sqrt{b^2+1}}\right)$

— 디립 사르 베이트
소스

또 다른 해결책은 다음과 같습니다. 우리는 정의합니다.

\begin{aligned} I (γ) & = \int_{- \infty}^{\infty} Φ (ξ x + γ) N (x | 0, σ^{2}) d x, \end{aligned}

$\begin{align*} I(\gamma) & =\int_{-\infty}^{\infty}\Phi(\xi x+\gamma)\mathcal{N}(x|0,\sigma^{2})dx, \end{align*}$

γ = - ξ μ

$\gamma=-\xi\mu$

I (γ)

$I(\gamma)$

I (0) = 0

$I(0)=0$

γ

$\gamma$

\begin{aligned} \frac{d I}{d γ} & = \int_{- \infty}^{\infty} N ((ξ x + γ) | 0, 1) N (x | 0, σ^{2}) d x \\ = \int_{- \infty}^{\infty} \frac{1}{\sqrt{2 π}} \exp (- \frac{1}{2} {(ξ x + γ)}^{2}) \frac{1}{\sqrt{2 π σ^{2}}} \exp (- \frac{x^{2}}{2 σ^{2}}) d x . \end{aligned}

$\begin{align*} \frac{dI}{d\gamma} & =\int_{-\infty}^{\infty}\mathcal{N}((\xi x+\gamma)|0,1)\mathcal{N}(x|0,\sigma^{2})dx\\ & =\int_{-\infty}^{\infty}\frac{1}{\sqrt{2\pi}}\exp\left(-\frac{1}{2}\left(\xi x+\gamma\right)^{2}\right)\frac{1}{\sqrt{2\pi\sigma^{2}}}\exp\left(-\frac{x^{2}}{2\sigma^{2}}\right)dx. \end{align*}$

\begin{aligned} {(ξ x + γ)}^{2} + \frac{x^{2}}{σ^{2}} & = \underset{= a}{\underset{⏟}{(ξ^{2} + σ^{- 2})}} x^{2} + \underset{= b}{\underset{⏟}{- 2 γ ξ}} x + \underset{= c}{\underset{⏟}{γ^{2}}} \\ = a {(x - \frac{b}{2 a})}^{2} + (c - \frac{b^{2}}{4 a}) (c - \frac{b^{2}}{4 a}) \\ = γ^{2} - \frac{4 γ^{2} ξ^{2}}{4 (ξ^{2} + σ^{- 2})} \\ = γ^{2} (1 - \frac{ξ^{2}}{ξ^{2} + σ^{- 2}}) \\ = γ^{2} (\frac{1}{1 + ξ^{2} σ^{2}}) \end{aligned}

$\begin{align*} \left(\xi x+\gamma\right)^{2}+\frac{x^{2}}{\sigma^{2}} & =\underbrace{\left(\xi^{2}+\sigma^{-2}\right)}_{=a}x^{2}+\underbrace{-2\gamma\xi}_{=b}x+\underbrace{\gamma^{2}}_{=c} \\ &=a\left(x-\frac{b}{2a}\right)^{2}+\left(c-\frac{b^{2}}{4a}\right) \left(c-\frac{b^{2}}{4a}\right)\\ & =\gamma^{2}-\frac{4\gamma^{2}\xi^{2}}{4\left(\xi^{2}+\sigma^{-2}\right)}\\ &=\gamma^{2}\left(1-\frac{\xi^{2}}{\xi^{2}+\sigma^{-2}}\right)\\ &=\gamma^{2}\left(\frac{1}{1+\xi^{2}\sigma^{2}}\right) \end{align*}$

\begin{aligned} \frac{d I}{d γ} & = \frac{1}{2 π σ} \exp (- \frac{1}{2} (c - \frac{b^{2}}{4 a})) \sqrt{\frac{2 π}{a}} \int_{- \infty}^{\infty} \sqrt{\frac{a}{2 π}} \exp (- \frac{1}{2} a {(x - \frac{b}{2 a})}^{2}) d x \\ = \frac{1}{2 π σ} \exp (- \frac{1}{2} (c - \frac{b^{2}}{4 a})) \sqrt{\frac{2 π}{a}} \\ = \frac{1}{\sqrt{2 π σ^{2} a}} \exp (- \frac{1}{2} (c - \frac{b^{2}}{4 a})) \\ = \frac{1}{\sqrt{2 π (1 + σ^{2} ξ^{2})}} \exp (- \frac{1}{2} \frac{γ^{2}}{1 + ξ^{2} σ^{2}}) \end{aligned}

$\begin{align*} \frac{dI}{d\gamma} & =\frac{1}{2\pi\sigma}\exp\left(-\frac{1}{2}\left(c-\frac{b^{2}}{4a}\right)\right)\sqrt{\frac{2\pi}{a}}\int_{-\infty}^{\infty}\sqrt{\frac{a}{2\pi}}\exp\left(-\frac{1}{2}a\left(x-\frac{b}{2a}\right)^{2}\right)dx\\ & =\frac{1}{2\pi\sigma}\exp\left(-\frac{1}{2}\left(c-\frac{b^{2}}{4a}\right)\right)\sqrt{\frac{2\pi}{a}}\\ &=\frac{1}{\sqrt{2\pi\sigma^{2}a}}\exp\left(-\frac{1}{2}\left(c-\frac{b^{2}}{4a}\right)\right)\\ & =\frac{1}{\sqrt{2\pi\left(1+\sigma^{2}\xi^{2}\right)}}\exp\left(-\frac{1}{2}\frac{\gamma^{2}}{1+\xi^{2}\sigma^{2}}\right) \end{align*}$

\begin{aligned} I (γ) & = \int_{- \infty}^{γ} \frac{1}{\sqrt{2 π (1 + σ^{2} ξ^{2})}} \exp (- \frac{1}{2} \frac{z^{2}}{1 + ξ^{2} σ^{2}}) d z \\ = Φ (\frac{γ}{\sqrt{1 + ξ^{2} σ^{2}}}) \end{aligned}

$\begin{align*} I(\gamma) &=\int_{-\infty}^{\gamma}\frac{1}{\sqrt{2\pi\left(1+\sigma^{2}\xi^{2}\right)}}\exp\left(-\frac{1}{2}\frac{z^{2}}{1+\xi^{2}\sigma^{2}}\right)dz\\ &=\Phi\left(\frac{\gamma}{\sqrt{1+\xi^{2}\sigma^{2}}}\right) \end{align*}$

이것은 암시

\begin{aligned} \int_{- \infty}^{\infty} Φ (ξ x) N (x | μ, σ^{2}) d x & = I (ξ μ) \\ = Φ (\frac{ξ μ}{\sqrt{1 + ξ^{2} σ^{2}}}) . \end{aligned}

$\begin{align*} \int_{-\infty}^{\infty}\Phi(\xi x)\mathcal{N}(x|\mu,\sigma^{2})dx &=I(\xi\mu)\\ &=\Phi\left(\frac{\xi\mu}{\sqrt{1+\xi^{2}\sigma^{2}}}\right). \end{align*}$

— 제니 라이닝 거
소스