print("Datos")
print("Muestra de algodon")
Yn1 <- 25 
Yn1
print("Muestra de acetato")
Yn2 <- 25 
Yn2
print("Media")
print("Algodon")
Ymedia1 <- 20   
Ymedia1
("Acetato")
Ymedia2 <- 12    
Ymedia2
print("Desviacion")
print("Algodon")
Ysd1 <- 1.5     
Ysd1
print("Acetato")
Ysd2 <- 1.25     
Ysd2
print("Calculo de la varianza")
sp2 <- (((Yn1 - 1) * Ysd1^2 + (Yn2 - 1) * Ysd2^2) / (Yn1 + Yn2 - 2))
sp2
sp <- sqrt(sp2)
# Estadístico t
t <- (Ymedia1 - Ymedia2) / (sp * sqrt(1/Yn1 + 1/Yn2))
# Grados de libertad
gl <- Yn1 + Yn2 - 2
# Valor p unilateral (cola derecha)
p_value <- pt(t, df = gl, lower.tail = FALSE)
# Mostrar resultados
cat("Estadístico t =", round(t, 3), "\n")
cat("Grados de libertad =", gl, "\n")
cat("Valor p (unilateral) =", round(p_value, 4), "\n")