lazymatrix provides a framework for working with
transformed sparse data matrices using lazy
evaluation. This approach is particularly useful for handling large
sparse matrices where transformations (e.g., centering and scaling)
would otherwise result in dense matrices, consuming significant memory
and computational resources.
You can install the development version of lazymatrix from GitHub with:
# install.packages("pak")
pak::pak("vsegersall/lazymatrix")Note: Once the package is accepted to CRAN, you can also install it via:
install.packages("lazymatrix")This is a basic example which shows you how to solve a common problem:
library(lazymatrix)
#>
#> Attaching package: 'lazymatrix'
#> The following object is masked from 'package:base':
#>
#> norm
set.seed(123)
sparse_matrix <- Matrix::Matrix(0, 5, 3)
sparse_matrix[sample(length(sparse_matrix), 5)] <- rnorm(5)
b <- rnorm(3)
lazy_matrix <- LazyMatrix(sparse_matrix, scale = "sd", location = "mean")
print(lazy_matrix)
#> An object of class "LazyMatrix"
#> Slot "data":
#> 5 x 3 sparse Matrix of class "dgCMatrix"
#>
#> [1,] . 1.55870831 .
#> [2,] . . .
#> [3,] . . .
#> [4,] -0.5604756 0.07050839 -0.2301775
#> [5,] . 0.12928774 .
#>
#> Slot "col_scales":
#> [1] 0.2506523 0.6769030 0.1029385
#>
#> Slot "row_scales":
#> [1] 0.89992066 0.00000000 0.00000000 0.31560781 0.07464431
#>
#> Slot "col_locations":
#> [1] -0.1120951 0.3517009 -0.0460355
#>
#> Slot "row_locations":
#> [1] 0.51956944 0.00000000 0.00000000 -0.24004825 0.04309591
lazy_matrix %*% b
#> 5 x 1 Matrix of class "dgeMatrix"
#> [,1]
#> [1,] -1.751397
#> [2,] 2.225999
#> [3,] 2.225999
#> [4,] -4.596694
#> [5,] 1.896092This package is currently in active development. While the core functionality is stable, the API may evolve as the project matures. Contributions are currently closed while the package is prepared for its first CRAN release. However, bug reports and feature requests are welcome via the GitHub Issue Tracker.