Question

假设我有一个矩阵列表：

matrix <- matrix(1:4, nrow = 2, ncol = 2)
list <- list(matrix, matrix, matrix)

由函数cbind()创建的矩阵：

long.matrix <- do.call(cbind, list)

      [,1] [,2] [,3] [,4] [,5] [,6]
[1,]    1    3    1    3    1    3
[2,]    2    4    2    4    2    4

我希望撤消从list获取long.matrix矩阵的过程。

我可以使用for循环手动执行此操作，但我正在搜索我认为应该存在的类似function(long.matrix, 3)的内容。有这样的事吗？

Answer 1

蛮力解决方案：

f <- function(long.matrix, num)
           lapply(split(long.matrix, 
                        rep(seq(num), each=(ncol(long.matrix)/num)*nrow(long.matrix))), 
                  function(x) matrix(x, nrow=nrow(long.matrix))
           )

f(long.matrix, 3)
## $`1`
##      [,1] [,2]
## [1,]    1    3
## [2,]    2    4
## 
## $`2`
##      [,1] [,2]
## [1,]    1    3
## [2,]    2    4
## 
## $`3`
##      [,1] [,2]
## [1,]    1    3
## [2,]    2    4

rep构建split的类别，以拆分数据。由于R是列专业，这里我们采用前四个，第二个四个，第三个四个条目。

填写示例long.matrix和3的当前维度的值，该函数会缩减为：

lapply(split(long.matrix, rep(seq(3), each=4)), function(x) matrix(x, nrow=2))

注意：

(r <- rep(seq(3), each=4) )
## [1] 1 1 1 1 2 2 2 2 3 3 3 3
split(long.matrix, r)
## $`1`
## [1] 1 2 3 4
## 
## $`2`
## [1] 1 2 3 4
## 
## $`3`
## [1] 1 2 3 4

然后将其中的每一个传递给matrix以获得所需的格式。

Answer 2

这样做：

listm=list()  #i=1
for(i in 1:3)listm[[i]]=long.matrix[,(2*i-1):(i*2)]

lapply版本

lapply(1:3,function(ii)long.matrix[,(2*ii-1):(ii*2)])

[[1]]
     [,1] [,2]
[1,]    1    3
[2,]    2    4

[[2]]
     [,1] [,2]
[1,]    1    3
[2,]    2    4

[[3]]
     [,1] [,2]
[1,]    1    3
[2,]    2    4

Answer 3

我更喜欢使用数组维度。然后，您可以为矩阵定义split方法：

split.matrix <- function(x, rslice = 1, cslice = 1) {
  if (ncol(x) %% cslice) stop("cslice not divisor of number of columns")
  if (nrow(x) %% rslice) stop("rslice not divisor of number of rows")

  x <- t(x)
  dim(x) <- c(dim(x)[1],
              dim(x)[2] / rslice,
              rslice)
  x <- lapply(seq_len(rslice), function(k, a) t(a[,,k]), a = x)

  if (cslice > 1) {
    x <- lapply(x, function(y, k) {

      dim(y) <- c(dim(y)[1],
                  dim(y)[2] / k,
                  k)
      y <- lapply(seq_len(k), function(k, a) a[,,k], a = y)
      y
    }, k = cslice)
  }
  if(length(x) == 1L) x <- x[[1]]

  x
}

split(long.matrix, 1, 3)
#[[1]]
#     [,1] [,2]
#[1,]    1    3
#[2,]    2    4
#
#[[2]]
#     [,1] [,2]
#[1,]    1    3
#[2,]    2    4
#
#[[3]]
#     [,1] [,2]
#[1,]    1    3
#[2,]    2    4

split(long.matrix, 1, 1)
#     [,1] [,2] [,3] [,4] [,5] [,6]
#[1,]    1    3    1    3    1    3
#[2,]    2    4    2    4    2    4

split(long.matrix, 2, 1)

#[[1]]
#     [,1] [,2] [,3] [,4] [,5] [,6]
#[1,]    1    3    1    3    1    3
#
#[[2]]
#     [,1] [,2] [,3] [,4] [,5] [,6]
#[1,]    2    4    2    4    2    4

split(long.matrix, 2, 3)
#[[1]]
#[[1]][[1]]
#[1] 1 3
#
#[[1]][[2]]
#[1] 1 3
#
#[[1]][[3]]
#[1] 1 3
#
#
#[[2]]
#[[2]][[1]]
#[1] 2 4
#
#[[2]][[2]]
#[1] 2 4
#
#[[2]][[3]]
#[1] 2 4

反向cbind（）函数R.

3 个答案: