Попробуйте эту функцию:
type_anova <- function(df, feature, type) {
b=aov(reformulate(type, feature), df)
F=summary(b)[[1]][["Pr(>F)"]][1]
p=summary(b)[[1]][["F value"]][1]
if (p <= 0.05) message="Reject"
else message="Accept"
cat("\nF Statistics: ", F," p-value: ", p," Null Hypothesis: ",message)
}
Использование ее на mtcars
наборе данных
type_anova(mtcars, "am", "cyl")
#F Statistics: 0.0022 p-value: 11 Null Hypothesis: Accept