n <- 25 # Tamaño de la muestra
x_barra <- 2800 # Consumo promedio de la muestra
sigma_muestral <- 280 # Desviación estándar muestral
mu <- 2700 # Valor para el cual se quiere calcular la probabilidad
sigma <- sigma_muestral / sqrt(n)
Z <- (mu - x_barra) / sigma
probabilidad <- pnorm(Z)
print(probabilidad)
## [1] 0.03707277
consumo <- c(710, 732, 654, 760, 774, 675, 757, 718, 814, 755, 736, 774, 768, 770, 677, 747, 823, 660, 797, 796)
media <- mean(consumo)
desviacion <- sd(consumo)
n <- length(consumo)
valor <- 745
Z <- (valor - media) / (desviacion / sqrt(n))
probabilidad <- 1 - pnorm(Z)
print(probabilidad)
## [1] 0.4946017