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

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)