如何在R中创建自定义矩阵,其中前几列为零,其余为Identity Matrix?

时间:2015-06-26 12:59:18

标签: arrays r

假设我们有一个订单A的矩阵k2*k,其中k=k1+k2k1=k-k2。矩阵A的特征是第一k1列为零,剩余k2*k2为单位矩阵。如何在R中创建这样的自定义矩阵?

注意:

对于k1k2的较小尺寸,这很容易。但我正在寻找一个自动命令来处理更大尺寸的k1k2

1 个答案:

答案 0 :(得分:2)

k <- 20
k2 <- 15
k1 <- k - k2

diagonal <- diag(k2)
zeros <- matrix(0, nrow = k2, ncol = k1)

result <- cbind(zeros, diagonal)