r - 在 R 中重塑数组时智能地填充给定维度

Question

假设我对网格n上的某个字段进行了观察100x100，并且这些数据存储为单个 vector obs。我想在100x100xn不知道什么n是先进的情况下将其重塑为一个数组。在matlab中我可以使用

reshape(obs,100,100,[]);

或在蟒蛇

np.reshape(obs,(100,100,-1))

R中是否有类似的功能，还是我必须手动计算最后一个索引的大小？

Answer 1

试试这个（我相信你可以适应你的进一步需求）：

shaper <- function(obs, a, b) {
 array(obs, dim=c(a, b, length(obs)/a/b))
}
shaper(obs, 100, 100)

如果您不确定所需的尺寸是否正确，您可以检查是否有剩菜，如下所示：

shaper <- function(obs, a, b) {
  dimension <- length(obs)/a/b 
  if (dimension %% 1 != 0) { 
   stop("not correctly divisible")
  }
  else {
   return(array(obs, dim=c(a, b, dimension)))
  }
}
shaper(obs, 100, 100)

还添加了将任意数量的维度作为输入的功能，它会尝试将其扩展 1。

shaper <- function(obs, ...) {
 len.remaining <- length(obs)
 for (i in c(...)) {
   len.remaining <- len.remaining / i
 }
 if (len.remaining %% 1 != 0) { 
  stop("not correctly divisible")
 }
 else {
  return(array(obs, dim=c(..., len.remaining)))
 } 
}

现在可以使用：

obs <- rep(1, 100 * 100 * 5)
> res <- shaper(obs, 100, 100)
> dim(res) 
[1] 100 100 5
> res <- shaper(obs, 10, 10, 100)
> dim(res)
[1] 10 10 100 5

Answer 2

这是你想要的吗？

 dim(obs) <- c(100,100,length(v)/100/100)

例如：

v <- seq(2*2*3)
dim(v) <- c(2,2,length(v)/2/2)

, , 1

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

, , 2

     [,1] [,2]
[1,]    5    7
[2,]    6    8

, , 3

     [,1] [,2]
[1,]    9   11
[2,]   10   12

请注意，对于

  v <- seq(2*2*3+1)

你会得到一个错误：

Error in dim(v) <- c(2, 2, length(v)/2/2) : 
  dims [product 12] do not match the length of object [13]

Answer 3

处理上面对偶性的代码，我实现了一个 R 函数，如 np.reshape 语法

npreshape <- function(x,shape) {
    x<- array(x)
    ind <- which(shape ==  -1)
    shape[ind] = prod(dim(x))/prod(shape[-ind])
    return(array(x,shape))
}

这是它在控制台中的工作方式

> npreshape(1:10,c(-1,5))
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    3    5    7    9
[2,]    2    4    6    8   10