R创建一个字符向量列表

Question

我正在尝试将邻接矩阵处理为社区列表，输出应为向量列表（字符类型）

例如

l1[1] = c("a","b","c")
l1[2] = c("d")
l1[3] = c("e","f")

我正在尝试在for循环中实现这一点：

l1=vector("list")

for(kn in 1:nrow(adjFinal)){
  temp = c()

  for(tm in 1:ncol(adjFinal))
    {
    if(adjFinal[kn,tm]==1)
      temp=c(temp,colnames(adjFinal)[tm])
  }

  l1[kn] = temp
}

但是我不断收到“要替换的项目数不是替换长度的倍数”的警告，并且列出了长度为1的3个向量，例如：

l1[1] = "a"
l1[2] = "d"
l1[3] = "e"

更新：这是adjFinal：

3 x 6 sparse Matrix of class "dgCMatrix"
  a b c d e f
b 1 1 1 . . .
d . . . 1 . .
e . . . . 1 1

和dput（adjFinal）：

new("dgCMatrix"
    , i = c(0L, 0L, 0L, 1L, 2L, 2L)
    , p = 0:6
    , Dim = c(3L, 6L)
    , Dimnames = list(c("b", "d", "e"), c("a", "b", "c", "d", "e", "f"))
    , x = c(1, 1, 1, 1, 1, 1)
    )

先感谢您。

Answer 1

你可以试试：

apply(adjFinal,1,function(x) colnames(adjFinal)[which(x==1)])
#$b
#[1] "a" "b" "c"
#$d
#[1] "d"
#$e
#[1] "e" "f"

考虑到您有一个稀疏矩阵，因此如果您有一个较大的矩阵，则上述方法可能效率不高，因为apply会将稀疏矩阵强制为常规矩阵。 为避免这种情况，您还可以尝试：

adj2<-as(adjFinal,"dgTMatrix")
split(colnames(adjFinal)[adj2@j+1],adj2@i)

这将获得相同的结果，但效率更高。