Sparse plus low rank matrix

Eventually this class will subclass Matrix objects, but for now this is a basic implementation that essentially only supports singular value decomposition.

Details

To learn more about S4 classes, please see https://adv-r.hadley.nz/s4.html.

Slots

sparse

sparseMatrix.

U

Matrix.

V

Matrix.

Examples

set.seed(528491)

n <- 50
m <- 40
k <- 3

A <- rsparsematrix(n, m, 0.1)

U <- Matrix(rnorm(n * k), nrow = n, ncol = k)
V <- Matrix(rnorm(m * k), nrow = m, ncol = k)

# construct the matrix, which represents A + U %*% t(V)
X <- sparseLRMatrix(sparse = A, U = U, V = V)

dim(X)
#> [1] 50 40

s <- svds(X, 5)  # efficient