print("Datos")
print("cantidad 1")
Yn1 <- 25
Yn1
print("cantidad 2")
Yn2 <- 36
Yn2
print("media")
Ymedia1 <- 81
Ymedia1
Ymedia2 <- 76
Ymedia2
print("desviacion")
Ysigma1 <- 5.2
Ysigma1
Ysigma2 <- 3.4
Ysigma2
print("Estadistico Z")
z <- (Ymedia1 - Ymedia2) / sqrt((Ysigma1^2 / Yn1) + (Ysigma2^2 / Yn2))
z
print("Valor p bilateral")
p_value <- 2 * pnorm(-abs(z))
p_value
# Mostrar resultados
cat("Estadístico z =", round(z, 3), "\n")
cat("Valor p (bilateral) =", round(p_value, 4), "\n")