grupo_A <- c(70, 75, 80, 85, 90, 95)
grupo_B <- c(85, 87, 89, 91, 93, 95)
fisher_test <- var.test(grupo_A, grupo_B)
print(fisher_test)
##
## F test to compare two variances
##
## data: grupo_A and grupo_B
## F = 6.25, num df = 5, denom df = 5, p-value = 0.06573
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
## 0.8745684 44.6648864
## sample estimates:
## ratio of variances
## 6.25
if (fisher_test$p.value > 0.05) {
cat("Las varianzas son iguales, se usa una prueba t con varianzas iguales.\n")
t_test <- t.test(grupo_A, grupo_B, var.equal = TRUE)
} else {
cat("Las varianzas no son iguales, se usa una prueba t con varianzas desiguales.\n")
t_test <- t.test(grupo_A, grupo_B, var.equal = FALSE)
}
## Las varianzas son iguales, se usa una prueba t con varianzas iguales.
print(t_test)
##
## Two Sample t-test
##
## data: grupo_A and grupo_B
## t = -1.8235, df = 10, p-value = 0.09822
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -16.664307 1.664307
## sample estimates:
## mean of x mean of y
## 82.5 90.0
cat("Intervalo de confianza del 95% para la diferencia de medias:\n")
## Intervalo de confianza del 95% para la diferencia de medias:
print(t_test$conf.int)
## [1] -16.664307 1.664307
## attr(,"conf.level")
## [1] 0.95