Как сделать ANOVA тест нескольких параметров в R

Question

У меня есть структура данных, организованная как на следующем рисунке (две популяции и некоторые параметры). Я запутался с кодированием. Какой способ сделать тест ANOVA, который понижает пригодность (значение F, p) каждого параметра, пожалуйста?

Answer 1

Вы можете попробовать что-то вроде этого, но, вероятно, необходимо учитывать несколько сравнений:

library(data.table)

DT <- data.table("GROUP" = c("GR1", "GR1", "GR1", "GR2", "GR2")
               , "Weight" = c(78, 85, 84.3, 70, 67)
               , "BloodPres" = c(11, 14, 13, 12, 12)
               , "Heart" = c(140, 130, 142, 135, 120)
               , "Glucose" = c(80, 110, 95, 97, 105))

Outcomes <- c("Weight", "BloodPres", "Heart", "Glucose")

sapply(Outcomes, function(my) {
       f <- as.formula(paste(my, "~GROUP", sep=""))
       summary(aov(f, data=DT))
})

Это работает в одностороннем порядке, но вам нужно настроить его для двустороннего (см. Правильное использование sapply с Anova для нескольких подмножеств в R ).