我试图将邻接矩阵处理成一个社区列表, 输出应该是矢量列表(字符类型)
例如
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
}
但我一直警告"要更换的物品数量不是更换长度的倍数" ,以及3个长度为1的向量的列表 像:
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)
)
提前谢谢。
答案 0 :(得分:3)
您可以尝试:
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)
将获得相同的结果,但效率更高。