Utilizando el resultado de la relación de molinos, dejemos que X∼N(μ,σ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?