Prueba de Levene con anova bidireccional

Question

Estoy intentando realizar un levenetest en mi anova de dos vías.

leveneTest(Score ~ Party + Group, data = Data) 

Pero R sólo me grita diciendo

Error in leveneTest.formula(formula(y), data = model.frame(y), ...) : 
  Model must be completely crossed formula only.

Answer 1

Creo que está obteniendo un error ya que su modelo no incluye un término de interacción. Yo intentaría primero crear un ajuste del modelo utilizando la función lm función:

fit <- lm(Score ~ Party + Group, data = Data)

y, a continuación, utilizando el levene.test en la función s20x ya que es un poco menos quisquilloso.

levene.test(fit)

Dicho todo esto, las pruebas de Levene son en realidad una forma muy pobre de comprobar la igualdad de varianza (EoV). Son demasiado sensibles a las desviaciones, mientras que los ANOVA son bastante robustos. Es mejor limitarse a observar un gráfico (por ejemplo plot(fit, which =1) ) y evaluarlo visualmente para EoV. Si es para una tarea que lo requiere, entonces ignore este último consejo. Lo mismo vale para la normalidad y la prueba de Shapiro-Wilks ....

Answer 2

Tienes que utilizar un diseño completamente cruzado para la fórmula, lo que significa que tienes que utilizar * en lugar de + en tu fórmula. Porque + no implica el término de interacción.

Haz esto:

leveneTest(Score ~ Party * Group, data = Data)

Y debería funcionar.