双向方差分析的 Levene 检验

Question

我想检查双向方差分析的方差假设的同质性。为此，我想使用 car 包中的函数 leveneTest()。 但是，它给了我错误：

Error in leveneTest.formula(formula, data, center = center) :  Levene's test is not appropriate with quantitative explanatory variables.

通过搜索错误，我发现解释变量的顺序很重要，该函数期望右侧是因子。我试过了，但没有用。

我也用相同的任务测试了一个教程：

library(car)

my_data <- ToothGrowth

leveneTest(len ~ supp*dose, data = my_data)

它给了我同样的错误，它应该给出结果：

Levene's Test for Homogeneity of Variance (center = median)
'  Df F value Pr(>F)'
'group  5  1.7086 0.1484'
      ' 54'  

是否有人遇到同样的错误，除了将数值变量转换为阶乘变量之外还有其他解决方案吗？

非常感谢您的帮助！

Answer 1

你必须这样做

leveneTest(len ~ factor(supp)*factor(dose), data = my_data)

而且我认为没有转换为因子就没有办法做到这一点，因为 Levene 的测试比较了组。组名是因子，而不是连续变量。