Левен пост-специальный тест в R

Question

У меня есть проблема, связанная с анализом данных в R. Одна из моих гипотез в основном состоит в том, что мои группы будут отличаться в плане разброса оценок, что указывает на разницу в конечностях между группами.

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

Это код, который я пробовал до сих пор:

leveneTest(SS_mean~RA01, DF)
DF$residuals <- abs(DF$SS_mean - DF$SS_mean_big) #SS_mean = Participants score, 
#SS_mean_big = mean for each group.

---

Мой тестовый и специальный тест выглядит следующим образом:

levene.anova<-aov(residuals~RA01, DF) #RA01 is the groups. Four in total
summary(levene.anova)
TukeyHSD(levene.anova)

---

ANOVA по остаткамтакже оказалось значимым, но значение p изменилось с 0,04 (тест Левена) до 0,01 (ANOVA по остаткам).Когда я читал об этом, казалось, что тест Левена - это просто ANOVA на resiudals, и поэтому он должен дать мне те же результаты.И я также не уверен, какой тест я должен использовать.Я подумал о Даннетте, так как он включает в себя базовый уровень, который соответствует одной из моих групп.

Наконец, я сделал leveneTest для остатков, а также "leveneTest (residuals ~ RA01)", что оказалось значительным.Для меня лучше использовать непараметрический тест, например, h-тест Крускала-Уоллиса, и вместо этого провести специальный тест на моем тесте Крускала-Уоллиса?И если это так, какой будет соответствующий тест?Должен ли я использовать попарный u-тест Манна-Уитни или тест Данна?

Поскольку я впервые делаю что-то подобное, я не уверен в том, что это законный анализ, я был бы очень признателен за вашу помощьили введите!

Answer 1

Тест Левена действительно должен дать то же значение p, что и ANOVA для остатков.

См., Например, этот код:

data("mtcars")
mtcars$cyl <- as.factor(mtcars$cyl)

# Calculate means and add them to data
cyl_means <- aggregate(disp ~ cyl, data = mtcars, FUN = mean)
colnames(cyl_means)[2] <- "disp_mean"
mtcars2 <- merge(mtcars, cyl_means, by = "cyl")

# Residuals and anova
mtcars2$residuals <- abs(mtcars2$disp - mtcars2$disp_mean)
res.aov <- aov(residuals ~ cyl, data = mtcars2)
summary(res.aov)

# Levene's test
lawstat::levene.test(mtcars$disp, mtcars$cyl, location = "mean")

Возможно, вы случайно запустили тест Брауна-Форсайта, который используется по умолчанию в lawstat::levene.test и использует медиану вместо среднего для вычисления остатков.

Используйте Dunnett, если вас интересует только сравнение групп с одной базовой группой. Используйте TukeyHSD, если вы хотите все парные сравнения между группами.