r - 如何创建稀疏对角矩阵？

Question

这可能是一个非常愚蠢的问题，但我似乎找不到在 R 中创建稀疏对角矩阵的合适方法。我找到了函数:

• diag.spam()
• spdiags()

并将它们与下载的库 Matrix 和包 spam 一起使用，但 R 似乎无法识别这些函数。有人知道我需要下载的函数或库吗？

我需要它，因为我想创建大于 256 x 256 的对角矩阵。

Answer 1

Diagonal() Matrix 中的函数包裹。 (Matrix 是“推荐”包，这意味着它会在您安装 R 时自动可用。)

library(Matrix)
m <- Diagonal(500)
image(m)

Diagonal(n)创建一个 n x n 单位矩阵。如果要创建具有指定对角线的对角矩阵 x , 使用 Diagonal(x=<your vector>)