5

在 SO ( LINK ) 上的一个问题中,一位发帖者提出了一个问题,我给出了一个有效的答案,但有一部分让我感到困扰,list从向量创建一个作为索引列表传递。所以说我有这个向量:

n <- 1:10
#> n
# [1]  1  2  3  4  5  6  7  8  9 10

假设我想将它分解成一个向量列表,每个向量的长度为 3。完成此任务的最佳(最短代码量或最快)方法是什么?我们想扔掉第 10 项,因为10 %% 310/3 ( length(n) - 10 %% 3) 中还有 1 ( ) 的余数。

这是期望的结果

list(1:3, 4:6, 7:9)

这将为我们提供那些不能组成三个组的索引:

(length(n) + 1 - 10 %% 3):length(n)

编辑

这是 Wojciech Sobala 在与之链接的另一个线程上发布的一个有趣的方法(我让他们在这里回答,如果他们回答,我将删除此编辑)

n <- 100
l <- 3
n2 <- n - (n %% l)
split(1:n2, rep(1:n2, each=l, length=n2))

作为一个函数:

indices <- function(n, l){
    if(n > l) stop("n needs to be smaller than or equal to l")
    n2 <- n - (n %% l)
    cat("numbers", (n + 1 - n %% l):n, "did not make an index of length", l)
    split(1:n2, rep(1:n2, each=l, length=n2))
}
4

3 回答 3

5

不确定这是否有效?

x = function(x, n){ 
    if(n > x) stop("n needs to be smaller than or equal to x")
    output = matrix(1:(x-x%%n), ncol=(x-x%%n)/n, byrow=FALSE)
    output
}

编辑:将输出更改为列表

x = function(x, n){ 
    if(n > x) stop("n needs to be smaller than or equal to x")
    output = matrix(1:(x-x%%n), ncol=(x-x%%n)/n, byrow=TRUE)
    split(output, 1:nrow(output))
}

Example:
x(10, 3)
$`1`
[1] 1 2 3

$`2`
[1] 4 5 6

$`3`
[1] 7 8 9
于 2012-05-19T03:02:55.743 回答
4
xx <- 1:10
xxr <- rle(0:(length(1:10)-1) %/% 3)  # creates an rle object
fac3 <- rep( xxr$values[xxr$lengths == 3], each=3)  #selects the one of length 3
                                     # and recreates the shortened grouping vector
tapply(xx[ 1:length(fac3)],          # a shortened original vector
                       fac3, list)   # split into little lists
$`0`                                # Hope you don't mind having names on your list
[1] 1 2 3

$`1`
[1] 4 5 6

$`2`
[1] 7 8 9
于 2012-05-19T03:12:59.823 回答
3

这不是最短的,但这里有一个小递归版本:

wrap <- function(n,x,lx,y) {
    if (lx < n) return (y)
    z <- x[-(1:n)]
    wrap(n, z, length(z), c(y, list(x[1:n])))
}

wrapit <- function(x,n) {
    wrap(n,x,length(x),list())
}

> wrapit(1:10,3)
[[1]]
[1] 1 2 3

[[2]]
[1] 4 5 6

[[3]]
[1] 7 8 9
于 2012-05-19T03:46:57.697 回答