我有一个0和1的对称但稀疏的矩阵,表示两个对象之间是否存在边缘或链接---邻接矩阵。我试图在R中的集群包中使用命令PAM,这需要不相似矩阵(特别是该矩阵的上三角形)或距离。
下面找到一个有代表性的矩阵(虽然我的矩阵有更多的边和节点,不像这个矩阵实际上是稀疏的)。
0 1 0 0
1 0 0 0
0 0 0 1
0 0 1 0
我需要我的矩阵的上三角形版本,这意味着我需要获得这个矩阵。我的问题是,我的矩阵至少是单数的,因此矩阵命令,如常规LU分解或Choleski分解不适用于它。
这是我在尝试64 * 64稀疏矩阵时得到的错误消息:
Warning message:
In .local(x, ...) :
Exact singularity detected during LU decomposition: U[i,i]=0, i=20.
我在稀疏矩阵上阅读了几篇帖子,包括“Large Sparse Matrix to Triangular Matrix R”
但它们似乎并不完全符合我的要求。
我想知道是否有人有关于我应该尝试获取所需矩阵的R命令的建议。如果有必要,我可以尝试编写自己的功能,但我不想这样做。