leveneTest 在 2-levels 的 3-level 因子之间？

Question

到目前为止我做了什么：

我有一个 data.frame results 响应 Fail，以及三个因素 PREP、CLEAN 和 ADHES。 ADHES有3个等级：CrestCrystPoly

我计算了方差：

sigma..k=tapply(Fail,ADHES,var)
print(sqrt(sigma..k)):

Crest    Cryst     Poly 
17.56668 41.64679 39.42669

然后使用leveneTest 测试方差常数：

print(leveneTest(Fail~ADHES))

Levene's Test for Homogeneity of Variance (center = median)
      Df F value  Pr(>F)  
group  2   3.929 0.02588 *
      51                 

问题：
现在我想使用 Levene's 仅在因子 ADHES 的 Cryst 和 Poly 级别之间进行测试，但我无法在 R 中计算出执行此操作的语法。

Answer 1

感谢@PauloCardoso 给我的提示，我解决了：

    leveneTest(subset(results,ADHES == 'Cryst' | ADHES == 'Poly')[,5],
subset(results,ADHES == 'Cryst' | ADHES == 'Poly')[,3])

（'Fail' & 'ADHES' 分别是我的 data.frame 'results' 中的第 5 列和第 3 列）

奥布里加迪尼奥！