到目前为止我做了什么:
我有一个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
,还有三个因素Fail
,PREP
&amp; CLEAN
。
ADHES
有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中找到执行此操作的语法。
答案 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!