R:创建行数未知的矩阵

Question

我已经编写了下面的代码来生成一个矩阵,其中包含对我来说相当复杂的模式.在这种情况下,我通过反复试验确定完成矩阵中有136行.

我可以编写一个函数来预先计算矩阵行的数量,但函数会有点复杂.在该示例中,矩阵中的行数=((4*3 + 1)+(3*3 + 1)+(2*3 + 1)+(1*3 + 1))*4.

有没有一种简单有效的方法在R中创建矩阵而不用硬连接矩阵语句中的行数？换句话说,是否有一种简单的方法让R在使用for循环时根据需要简单地向矩阵添加一行？

我提出了一个解决方案,每次通过循环使用rbind,但这似乎有点复杂,我想知道是否有一个更容易的解决方案.

对不起,如果这个问题与之前的问题有关.我无法使用本网站上的搜索功能或今天使用互联网搜索引擎找到类似的问题,尽管我认为我在过去的某处发现过类似的问题.

下面是2组示例代码,一组使用rbind,另一组使用试验和错误预先设置nrow = 136.

谢谢你的任何建议.

v1     <- 5
v2     <- 2
v3     <- 2
v4     <- (v1-1)

my.matrix <- matrix(0, nrow=136, ncol=(v1+4) )

i = 1

for(a in 1:v2) {
  for(b in 1:v3) {
    for(c in 1:v4) {
      for(d in (c+1):v1) {

        if(d == (c+1)) l.s = 4 
        else           l.s = 3

        for(e in 1:l.s) {

          my.matrix[i,c] = 1

            if(d == (c+1)) my.matrix[i,d]  = (e-1)
            else           my.matrix[i,d]  =  e

          my.matrix[i,(v1+1)] = a
          my.matrix[i,(v1+2)] = b
          my.matrix[i,(v1+3)] = c
          my.matrix[i,(v1+4)] = d

          i <- i + 1

        }
      }
    }
  }
}

my.matrix2 <- matrix(0, nrow=1, ncol=(v1+4) )
my.matrix3 <- matrix(0, nrow=1, ncol=(v1+4) )

i = 1

for(a in 1:v2) {
  for(b in 1:v3) {
    for(c in 1:v4) {
      for(d in (c+1):v1) {

        if(d == (c+1)) l.s = 4 
        else           l.s = 3

        for(e  in 1:l.s) {

          my.matrix2[1,c] = 1

          if(d == (c+1)) my.matrix2[1,d]  = (e-1)
          else           my.matrix2[1,d]  =  e

          my.matrix2[1,(v1+1)] = a
          my.matrix2[1,(v1+2)] = b
          my.matrix2[1,(v1+3)] = c
          my.matrix2[1,(v1+4)] = d

          i <- i+1

          if(i == 2) my.matrix3 <- my.matrix2
          else       my.matrix3 <- rbind(my.matrix3, my.matrix2)

          my.matrix2 <- matrix(0, nrow=1, ncol=(v1+4) )

        }
      }
    }
  }
}

all.equal(my.matrix, my.matrix3)

Answer 1

如果矩阵的大小有一些上限,则可以创建一个足以容纳所有数据的矩阵

my.matrix <- matrix(0, nrow=v1*v2*v3*v4*4, ncol=(v1+4) )

并在最后截断它.

my.matrix <- my.matrix[1:(i-1),]