双向方差分析的 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 的检验比较了组。组名是因子，不是连续变量。