从稀疏矩阵中获得上三角矩阵

时间:2013-02-01 19:53:38

标签: r cluster-analysis graph-theory linear-algebra sparse-matrix

我有一个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命令的建议。如果有必要,我可以尝试编写自己的功能,但我不想这样做。

0 个答案:

没有答案