leveneTest介于2级3级因子之间?

时间:2015-10-19 09:20:11

标签: r

到目前为止我做了什么:

我有一个data.frame <?php global $wpdb; $result = $wpdb->get_results("SELECT *, field1, field2 case when field1+field2 > some_value then 1 else 2 end as my_alias FROM my_table ORDER BY my_alias ASC"); print_r($result); ?> ,回复为results,还有三个因素FailPREP&amp; CLEANADHES有3个级别:ADHES Crest Cryst

我计算了差异:

Poly

然后使用sigma..k=tapply(Fail,ADHES,var) print(sqrt(sigma..k)): Crest Cryst Poly 17.56668 41.64679 39.42669 来测试方差的常数:

leveneTest

问题:
现在我想用Levene来测试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 &amp;因子Cryst的{​​{1}}级别,但我无法在R中找到执行此操作的语法。

1 个答案:

答案 0 :(得分:1)

感谢提示@PauloCardoso给了我解决的问题:

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

(&#39;失败&#39;&amp;&#39; ADHES&#39;分别是我的data.frame&#39;结果&#39;中的第5列和第3列)

Obrigadinho!