Processing math: 100%

8 votos

Cálculo del valor esperado de la normal truncada

Utilizando el resultado de la relación de molinos, dejemos que XN(μ,σ2) entonces

E(X|X<α)=μσϕ(aμσ)Φ(aμσ)

Sin embargo, al calcularlo en R. no obtengo los resultados correctos como

> mu <- 1
> sigma <- 2 
> a <- 3 
> x <- rnorm(1000000, mu, sigma) 
> x <- x[x < a] 
> mean(x)
[1] 0.4254786
> 
> mu -  sigma * dnorm(a, mu, sigma) / pnorm(a, mu, sigma)
[1] 0.7124

¿Qué estoy haciendo mal?

15voto

Vitaly Zdanevich Puntos 95

La implementación de su fórmula es errónea porque, ϕ(xμσ)=12πe12(xμσ)2fX,μ,σ(x)=12πσe12(xμσ)2 Como puede ver, tenemos un extra σ en el denominador de fX,μ,σ(x) , lo que da como resultado: ϕ(xμσ)=σfX,μ,σ(x) dnorm le ofrece el método fX,μ,σ(x) donde hay que multiplicarlo por σ para obtener ϕ . Desde su σ=2 Esto puede hacerse prácticamente restando de nuevo el segundo término, que es 10.7124=0.2876 : 10.28760.2876=0.4247 que se acerca a su estimación.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X