我将在R中使用相当大(7 e6 x 4.5 e3)但非常稀疏的矩阵工作.所以我试图理解如何有效地处理稀疏矩阵.我有两个相关的问题.
第一:我已经了解该Matrix包自动链接到LAPACK和SuiteSparse编译的dll.(我在Windows中工作.)我的印象是,与使用LAPACK套件使用密集矩阵相比,使用SuiteSparse例程可缩短执行时间.但低于运行测试表明,在运行时为稀疏的版本矩阵的是多较密的版本慢.
> library(Matrix)
> sparse <- sparseMatrix(1:4, 1:4, x=rnorm(4))
> dense <- as.matrix(sparse)
> x <- 1:4
> system.time(for (i in 1:10000) sparse %*% x)
user system elapsed
0.23 0.00 0.23
> system.time(for (i in 1:10000) dense %*% x)
user system elapsed
0 0 0
> system.time(for (i in 1:1000) solve(sparse))
user system elapsed
3.94 0.00 3.94
> system.time(for (i in 1:1000) solve(dense))
user system elapsed
0.05 0.00 0.05
Run Code Online (Sandbox Code Playgroud)
a)我是否正确Matrix …
当我使用corrplot :: corrplot()绘制相关矩阵时,y轴上的底部标签(1)是半截断的,因为绘图的底部位于绘图区域的最底部,并且1以底部轴为中心.我想用这个情节出版.如何在底部留出更多空间,以便不切断此底部y轴标签?在此先感谢您的情节和帮助.除上述问题外,这是一个非常好的情节.拉里亨斯克