Как сделать одностороннюю анову со всеми столбцами данных славы?

Question

Я хочу сделать тест ANOVA с моими данными об изобилии разных станций, я хочу знать, как правильно написать скрипт, потому что все примеры с функцией aov связаны с формулой, но я не уверен, что эта формула учитывает все столбцы, сравнивая их средства.

Следующие данные являются одним из примеров того, как я использую этот тест. Но в столбцах нет категориальной переменной, поэтому я уверен, что этот путь неверен.

Спасибо за вашу помощь

set.seed(200)
D <- data.frame(a=sample(15),b=sample(15), c=sample(15))

A<-aov(a~c, data = D)

Answer 1

Вы должны сначала изменить data.frame от широкого до длинного формата.

Следующее использует внешний пакет, reshape2, чтобы изменить данные из широкого в длинный формат.

molten <- reshape2::melt(D)
head(molten)

model <- lm(value ~ variable, data = molten)
anova(model)