Hola Miguel Angel, pues no sé si estás viendo algo que yo no veo, pero en teoria os resultados del modelo con interacción dan unos residuos que no son normales, por lo que no veo donde dices que puede cambiar el resultado. Te copio lo que me salió en el proceso, está en el bloque nuevo y lo he vuelto a repetir. Si me puedes decir qué estás viendo, te lo agradezco
res.aov2 <- aov(DIff_ODI2 ~ NHD * Grupo, data = df)
> summary(res.aov2)
Df Sum Sq Mean Sq F value Pr(>F)
NHD 2 793 397 8.972 0.000188 ***
Grupo 1 31273 31273 707.590 < 2e-16 ***
NHD:Grupo 2 86 43 0.973 0.379603
Residuals 194 8574 44
—
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> aov_residuos <- residuals(object = res.aov2 )
> shapiro.test(x = aov_residuos )
Shapiro-Wilk normality test
data: aov_residuos
W = 0.97572, p-value = 0.001529
