Непараметрический тест, эквивалент n-образного дисперсионного анализа.

Я хотел протестировать модель со следующим:

summary(aov(dep~ind.1*ind.2))

Но значения p не взаимозаменяемы, поскольку не соблюдаются предположения о нормальности и гомоскедастичности. Я ищу непараметрический тест, который мог бы заменить эту двустороннюю Anova (и, в более общем плане, n-way Anova).

Тест Дурбина-Ватсона - хорошее решение?

Я пытаюсь провести тест Дарбина-Ватсона, но мне это не удается!

require(lmtest)
dwtest(dep~ind.1*ind.2) # Fail
dwtest(lm(dep~ind.1*ind.2)) # I get only one p.value instead of the three I expected

Чтобы мой вопрос можно было воспроизвести, вот некоторые данные:

set.seed(34)
dep = runif(24,0,1)
ind.1 = rep(c(1,2),12)
ind.2 = rep(c(1,2),each=12)