小编Lar*_*ker的帖子

密集与稀疏矩阵代数的速度

我将在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 …

r matrix sparse-matrix

4
推荐指数
1
解决办法
653
查看次数

R corrplot作物底轴标签

当我使用corrplot :: corrplot()绘制相关矩阵时,y轴上的底部标签(1)是半截断的,因为绘图的底部位于绘图区域的最底部,并且1以底部轴为中心.我想用这个情节出版.如何在底部留出更多空间,以便不切断此底部y轴标签?在此先感谢您的情节和帮助.除上述问题外,这是一个非常好的情节.拉里亨斯克

r r-corrplot

3
推荐指数
1
解决办法
2065
查看次数

标签 统计

r ×2

matrix ×1

r-corrplot ×1

sparse-matrix ×1