Question

假设我有一个索引列表，包括空的元素，如：

l <- list(c(1,2,3), c(1), c(1,5), NULL, c(2, 3, 5))

指定矩阵中的非零元素，如：

(m <- matrix(c(1,1,1,0,0, 1,0,0,0,0, 1,0,0,0,5, 0,0,0,0,0, 0,1,1,0,1), nrow=5, byrow=TRUE))
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    1    1    0    0
[2,]    1    0    0    0    0
[3,]    1    0    0    0    5
[4,]    0    0    0    0    0
[5,]    0    1    1    0    1

使用R从m创建l的最快方法是什么，假设矩阵非常大，比如50.000行和2000列？

Answer 1

您可以Filter来自＆＃39; l＆＃39;的非NULL列表元素并使用melt将其重塑为“data.frame＆＃39;用＆＃39;键/值＆＃39;列或`row / column＆＃39;索引列。

library(reshape2)
d2 <- melt(Filter(Negate(is.null), setNames(l, seq_along(l))))
Un1 <- unlist(l)
m1 <- matrix(0, nrow=length(l), ncol=max(Un1))
m1[cbind(as.numeric(d2[,2]), d2[,1])] <- 1
m1
#     [,1] [,2] [,3] [,4] [,5]
#[1,]    1    1    1    0    0
#[2,]    1    0    0    0    0
#[3,]    1    0    0    0    1
#[4,]    0    0    0    0    0
#[5,]    0    1    1    0    1

或者

 library(Matrix)
 as.matrix(sparseMatrix(as.numeric(d2[,2]), d2[,1], x=1))
 #     [,1] [,2] [,3] [,4] [,5]
#[1,]    1    1    1    0    0
#[2,]    1    0    0    0    0
#[3,]    1    0    0    0    1
#[4,]    0    0    0    0    0
#[5,]    0    1    1    0    1

Answer 2

你可以这样做：

do.call(rbind, lapply(l, function(x) (1:max(unlist(l)) %in% x)+0L))

#     [,1] [,2] [,3] [,4] [,5]
#[1,]    1    1    1    0    0
#[2,]    1    0    0    0    0
#[3,]    1    0    0    0    1
#[4,]    0    0    0    0    0
#[5,]    0    1    1    0    1

即使akrun解决方案应该更快！

R：索引列表，包括空索引到二进制矩阵

2 个答案: