在R中,我如何测试两个因素是否相等?

时间:2012-09-28 23:00:57

标签: r r-factor equivalence

我正在生成一个包含不同级别的因素列表,我希望能够检测到其中两个因素定义相同的分区。例如,我想检测以下所有内容是否相同:

x1 <- factor(c("a", "a", "b", "b", "c", "c", "a", "a"))
x2 <- factor(c("c", "c", "b", "b", "a", "a", "c", "c"))
x3 <- factor(c("x", "x", "y", "y", "z", "z", "x", "x"))
x4 <- factor(c("a", "a", "b", "b", "c", "c", "a", "a"), levels=c("b", "c", "a"))

这样做的最佳方式是什么?

1 个答案:

答案 0 :(得分:4)

我想您要确定双向制表具有与单向分类相同的填充级别数。 interaction中的默认设置是表示所有级别,即使未填充但设置drop = TRUE会更改它以适合您的目的:

> levels (interaction(x1,x2, drop=TRUE) )
[1] "c.a" "b.b" "a.c"
> length(levels(x1) ) == length(levels(interaction(x1,x2,drop=TRUE) ) )
[1] TRUE

概括性将考虑all( <the 3 necessary logical comparisons> )

 all( length(levels(x1) ) == length(levels(interaction(x1,x2,drop=TRUE) ) ),
      length(levels(x1) ) == length(levels(interaction(x1,x3,drop=TRUE) ) ),
      length(levels(x1) ) == length(levels(interaction(x1,x4,drop=TRUE) ) ) )
#[1] TRUE