[英]What is the fastest way to construct this sparse matrix in R
那这个呢?
M <- matrix(0,nrow = 5,ncol = 9)
M[1,1:3] <- 1
M[2,1:3] <- -1
diag(M[3:5,7:9]) <- 1
如果你想要一个单线你可以这样做:
t(`[<-`(`[<-`(`[<-`(matrix(0, 9, 5), 1:3, 1), 10:12, -1), 7:9, 3:5, diag(3)))
#> [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
#> [1,] 1 1 1 0 0 0 0 0 0
#> [2,] -1 -1 -1 0 0 0 0 0 0
#> [3,] 0 0 0 0 0 0 1 0 0
#> [4,] 0 0 0 0 0 0 0 1 0
#> [5,] 0 0 0 0 0 0 0 0 1
或者如果你想要一些真正的代码高尔夫,
`[<-`(`[<-`(matrix(0,5,9),c(1+0:2*5,0:2*6+33),1),2+0:2*5,-1)
我们可以使用bdiag
在一行中做到这一点
library(Matrix)
as.matrix(bdiag(cbind(rbind(rep(1, 3), rep(-1, 3)), 0, 0), diag(3)))
-输出
# [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
#[1,] 1 1 1 0 0 0 0 0
#[2,] -1 -1 -1 0 0 0 0 0
#[3,] 0 0 0 0 0 1 0 0
#[4,] 0 0 0 0 0 0 1 0
#[5,] 0 0 0 0 0 0 0 1
矩阵只是一个带有dim
属性的向量。 matrix(c(rep(c(1,-1,0,0,0),3), rep(0,17), 1, rep(c(rep(0,5), 1), 2)), ncol = 9)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.