print("Datos")
print("muestra")
Yn <- 20     
Yn
print("promedio mestral")
Ymedia_muestra <- 33.1  
Ymedia_muestra
print("media bajo Ho")
Ymedia_hipotetica <- 35  
Ymedia_hipotetica
print("desviacion muestral")
Ysd_muestra <- 4.3       
Ysd_muestra
print("Estadistico t")
t <- (Ymedia_muestra - Ymedia_hipotetica) / (Ysd_muestra / sqrt(Yn))
t
print("Grados de libertad")
gl <- Yn - 1
gl
# Valor p para prueba unilateral (cola izquierda)
p_value <- pt(t, df = gl, lower.tail = TRUE)

print("Resultados")
cat("Estadístico t =", round(t, 3), "\n")
cat("Grados de libertad =", gl, "\n")
cat("Valor p (unilateral) =", round(p_value, 4), "\n")