在R中展开矩阵

Question

我想将线性程序的约束放在矩阵的行中。 我从列表的元素中得到约束，但是列表中的几个元素由零向量组成。 在我的情况下，相应的约束并不重要，也不必在约束矩阵中。

有没有一种方法可以预先定义矩阵，这样您就可以继续向矩阵中添加行，而不必担心矩阵的大小。 您只想继续向矩阵添加行，最后您将看到矩阵有多少行。

这是到目前为止的代码：

B<-matrix(HERE TO DO THE TRICK SUCH THAT THE FOR LOOP BELOW GIVES NO ERROR)

for (i in 1:(length(cars))){
  if (sum(cars[[i]])) > 0{
      B[i,]<-cars[[i]]
  }
}

cars是一个在每个元素中都有一个向量的列表。 我不想取消列出，因为我的列表太大了，因此我想摆脱具有由零组成的向量的元素。

Answer 1

我将创建NA的矩阵

B <- matrix(NA, ncol = length(cars[[1]]), nrow = (length(cars)))

然后运行循环以填充矩阵

for (i in 1:(length(cars))){
  if (sum(cars[[i]])) > 0{
      B[i,]<-cars[[i]]
  }
}

并在此之后排除所有带有NA行

B <- B[-which(is.na(B[ ,1])), ]

Answer 2

我推荐rbind() 。 它合并了多个数据帧（或矩阵）的行。 请参阅cbind()合并列。

cars <- 1:5
B <- NULL
for(i in 1:length(cars)){
  B <- rbind(B, mtcars[i, ])
}
B
##                    mpg cyl disp  hp drat    wt  qsec vs am gear carb
## Mazda RX4         21.0   6  160 110 3.90 2.620 16.46  0  1    4    4
## Mazda RX4 Wag     21.0   6  160 110 3.90 2.875 17.02  0  1    4    4
## Datsun 710        22.8   4  108  93 3.85 2.320 18.61  1  1    4    1
## Hornet 4 Drive    21.4   6  258 110 3.08 3.215 19.44  1  0    3    1
## Hornet Sportabout 18.7   8  360 175 3.15 3.440 17.02  0  0    3    2

要组合多个多维数组，请参见abind包。