En R esto debería funcionar.
X <- rnorm(100000,-3,4)
X1<- abs(x)
mean(X1<1)
.14854
Esto crea 100.000 variables aleatorias normales con media -3 y sd= 4. (R utiliza sd como segundo parámetro. Asegúrese de utilizar 2 si su varianza media es 4 en su ejemplo). La segunda línea de código sólo toma el valor absoluto de las variables, a continuación, la línea tres sólo comprueba si el porcentaje es menor que 1.
R tiene una función ( pnorm
) que es la función de densidad acumulativa para la distribución normal (también existe para muchas otras distribuciones distintas de la normal). pnorm
básicamente dice qué parte del pdf es menor o igual que un determinado número. Para tu ejemplo, estás buscando números entre -1 y 1 (es decir, el valor absoluto de 1), así que pnorm(1,-3,4)
le dirá $P(x<1)$ . Ahora que le dará más de lo que quieres, porque usted no quiere que los números más pequeños que $(x<-1)$ . Así que tienes que cortar esa parte de la distribución.
Escribir el código pnorm(1,mean=-3,sd=4) - pnorm(-1,mean=-3, sd=4)
debería darle la parte de la distribución a la izquierda del valor $1$ pero a la derecha del valor $-1$ . Recibí una respuesta de .14998
cuando ejecuté este código. Espero que esto ayude.