如何将行转换为R中的因子?

RNA*_*RNA 1 r matrix r-factor

例如,我有一个矩阵:

> a = cbind(sample(c(0,1),6,replace=T), sample(c(0,1),6,replace=T))
> a
     [,1] [,2]
[1,]    0    0
[2,]    0    0
[3,]    0    1
[4,]    1    0
[5,]    1    0
[6,]    1    1
Run Code Online (Sandbox Code Playgroud)

我想制作一个对象bab是一个因素,每个级别代表 中的不同行a。在这种情况下,b将是:

> b
[1] 1 1 2 3 3 4
Levels: 1 2 3 4
Run Code Online (Sandbox Code Playgroud)

我可以用一种肮脏的方式来做,但我想知道是否有一个优雅的解决方案?

A5C*_*2T1 5

不知道您当前的“肮脏”方式是什么,这是一个可能的解决方案:

> aFac <- interaction(data.frame(a), lex.order=TRUE)
> factor(aFac, levels = levels(aFac), labels = seq_along(levels(aFac)))
[1] 1 1 2 3 3 4
Levels: 1 2 3 4
Run Code Online (Sandbox Code Playgroud)

在哪里:

a <- structure(c(0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 1L), 
               .Dim = c(6L, 2L), .Dimnames = list(NULL, NULL))
Run Code Online (Sandbox Code Playgroud)

我使用的唯一原因lex.order = TRUE是为了匹配您的特定输出。


另一种可能性是:

> aFac <- interaction(data.frame(a), lex.order=TRUE, drop = TRUE)
> factor(as.numeric(aFac))
[1] 1 1 2 3 3 4
Levels: 1 2 3 4
Run Code Online (Sandbox Code Playgroud)

就是drop = TRUE从 中删除任何未使用的级别interaction,正如我们在下面的注释中的示例中得到的那样。

为了证明 的影响drop = TRUE,请考虑以下事项,并注意所产生的因子水平:

> b <- structure(c(1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1), .Dim = c(6L,2L))
> bFac1 <- interaction(data.frame(b), lex.order=TRUE)
> bFac2 <- interaction(data.frame(b), lex.order=TRUE, drop=TRUE)
> factor(as.numeric(bFac1))
[1] 3 4 3 2 2 4
Levels: 2 3 4
> factor(as.numeric(bFac2))
[1] 2 3 2 1 1 3
Levels: 1 2 3
Run Code Online (Sandbox Code Playgroud)


dro*_*opy 5

一个可能的解决方案:

 b <- apply(a, 1, paste, collapse="_")
 b <- factor(b, levels=unique(b), labels=1:length(unique(b)))
Run Code Online (Sandbox Code Playgroud)