在 R 中构造这个稀疏矩阵的最快方法是什么

Question

在 R 中构建上述矩阵的最快方法是什么？不知何故，我觉得必须有比下面更好的方法。

M <- t(matrix(c(1,1,1,-1,-1,-1),nrow=3))
M <- rbind(M, matrix(rep(0,9), nrow=3))
M <- cbind(M, matrix(rep(0,5*3), ncol=3))
M <- cbind(M,rbind(matrix(rep(0,2*3),ncol=3),diag(3)))

Answer 1

那这个呢？

M <- matrix(0,nrow = 5,ncol = 9)
M[1,1:3] <- 1
M[2,1:3] <- -1
diag(M[3:5,7:9]) <- 1