Как создать разреженную диагональную матрицу?

Question

Это может быть очень глупый вопрос, но я не могу найти правильный способ создания разреженной диагональной матрицы в R. Я нашел функции:

• diag.spam()
• spdiags()

и использовал их с загруженной библиотекой Matrix и пакетом spam, но R, похоже, не распознал эти функции. Кто-нибудь знает функцию или библиотеку, которую мне нужно скачать?

Мне это нужно, потому что я хочу создать диагональные матрицы больше 256 на 256.

Answer 1

Функция Diagonal() в пакете Matrix. (Matrix - это «рекомендуемый» пакет, что означает, что он автоматически доступен при установке R).

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

Diagonal(n) создает идентификационную матрицу nxn. Если вы хотите создать диагональную матрицу с указанной диагональю x, используйте Diagonal(x=<your vector>)

Answer 2

Используйте bandSparse библиотеки Matrix.

, чтобы получить матрицу n-на-n с m по диагонали, напишите:

bandSparse(n,n,0,list(rep(m, n+1)))