AEST Actividad

# Definir los valores iniciales
SST <- 10800  # Suma de cuadrados total
SSA <- 4560   # Suma de cuadrados entre tratamientos
N <- 30       # Número total de observaciones
k <- 3        # Número de tratamientos

# Calcular SSE
SSE <- SST - SSA

# Grados de libertad
df_total <- N - 1
df_trat <- k - 1
df_error <- N - k

# Calcular los cuadrados medios
MSA <- SSA / df_trat
MSE <- SSE / df_error

# Estadístico F
F_calculado <- MSA / MSE

# Valor crítico F para α = 0.05
alpha <- 0.05
F_critico <- qf(1 - alpha, df_trat, df_error)

# Resultados
cat("Resultados del ANOVA:\n")

## Resultados del ANOVA:

cat("SSE:", SSE, "\n")

## SSE: 6240

cat("MSA (Cuadrado medio entre tratamientos):", MSA, "\n")

## MSA (Cuadrado medio entre tratamientos): 2280

cat("MSE (Cuadrado medio del error):", MSE, "\n")

## MSE (Cuadrado medio del error): 231.1111

cat("F calculado:", F_calculado, "\n")

## F calculado: 9.865385

cat("F crítico:", F_critico, "\n")

## F crítico: 3.354131

# Conclusión
if (F_calculado > F_critico) {
  cat("Conclusión: Rechazamos la hipótesis nula. Hay diferencias significativas entre las medias.\n")
} else {
  cat("Conclusión: No se rechaza la hipótesis nula. No hay evidencia suficiente de diferencias entre las medias.\n")
}

## Conclusión: Rechazamos la hipótesis nula. Hay diferencias significativas entre las medias.