4

我正在生成一个包含不同级别的因素的大列表,并且我希望能够检测其中两个何时定义相同的分区。例如,我想将以下所有内容检测为彼此等效:

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"))

做这个的最好方式是什么?

4

1 回答 1

5

我想您想确定双向制表与单向分类具有相同数量的填充级别。中的默认设置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
于 2012-09-28T23:23:25.680 回答