leveneTest在2級3級因子之間？

Question

到目前為止，我所做的是：

我有一個響應Fail的data.frame results ，以及三個因素PREP ， CLEAN和ADHES 。 ADHES具有3個級別： Crest Cryst Poly

我計算了方差：

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

Crest    Cryst     Poly 
17.56668 41.64679 39.42669

然后使用leveneTest測試方差的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來僅在因子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列）

Obrigadinho！