我有以下索引向量
V_ind = cumsum(c(10,9,8,7,6,5,4,3,2,1))
[1] 10 19 27 34 40 45 49 52 54 55
我创建了以下FOR LOOP
k=1
for(ind in V_ind){
if(ind<=10){
print("ok")
}else{
print(c(V_ind[1:k]))
k = k + 1
}
}
结果
[1] "ok"
[1] 10
[1] 10 19
[1] 10 19 27
[1] 10 19 27 34
[1] 10 19 27 34 40
[1] 10 19 27 34 40 45
[1] 10 19 27 34 40 45 49
[1] 10 19 27 34 40 45 49 52
[1] 10 19 27 34 40 45 49 52 54
然而,我试图实现的是以下结果
[1] "ok"
[1] 10
[1] 9 10 19
[1] 8 9 10 18 19 27
[1] 7 8 9 10 17 18 19 26 27 34
[1] 6 7 8 9 10 16 17 18 19 25 26 27 33 34 40
[1] 5 6 7 8 9 10 15 16 17 18 19 24 25 26 27 32 33 34 39 40 45
[1] 4 5 6 7 8 9 10 14 15 16 17 18 19 23 24 25 26 27 31 32 33 34 38 39 40 44 45 49
[1] 3 4 5 6 7 8 9 10 13 14 15 16 17 18 19 22 23 24 25 26 27 30 31 32 33 34 37 38 39 40 43 44 45 48 49 52
[1] 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17 18 19 21 22 23 24 25 26 27 29 30 31 32 33 34 36 37 38 39 40 42 43 44 45 47 48 49 51 52 54
结果如下:
在第一次迭代中,我们只打印OK
在第二次迭代中,我们提取向量
V_ind
的第一个元素在第三次迭代中,我们提取向量
V_ind
的第一个和第二个元素,以及向量V_ind
- 1的第一个元素,即数字9
。在第四次迭代中,我们提取向量
V_ind
的第一,第二和第三个元素,以及第一个元素减1,即9
,第一个元素减2,即8
,第二个元素减1,即。18
.在第五次迭代中,我们提取向量
V_ind
的第一、二、三、四元素,第一个元素分别减去1、2、3,即7,8,9
,第二个元素减去1、2,即17,18
,第三个元素减去1,即26
。
这个过程一直持续到FOR LOOP结束。这有可能在R中以一种通用的方式完成吗?
使用purrr
的一个选项可以是:
map(.x = accumulate(V_ind, c),
~ map2(.x,
rev(seq_along(.x) - 1),
function(y, z) seq(y - z, y, 1)) %>%
reduce(c))
[[1]]
[1] 10
[[2]]
[1] 9 10 19
[[3]]
[1] 8 9 10 18 19 27
[[4]]
[1] 7 8 9 10 17 18 19 26 27 34
[[5]]
[1] 6 7 8 9 10 16 17 18 19 25 26 27 33 34 40
[[6]]
[1] 5 6 7 8 9 10 15 16 17 18 19 24 25 26 27 32 33 34 39 40 45
[[7]]
[1] 4 5 6 7 8 9 10 14 15 16 17 18 19 23 24 25 26 27 31 32 33 34 38 39 40 44 45 49
[[8]]
[1] 3 4 5 6 7 8 9 10 13 14 15 16 17 18 19 22 23 24 25 26 27 30 31 32 33 34 37 38 39 40 43 44 45 48 49 52
[[9]]
[1] 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17 18 19 21 22 23 24 25 26 27 29 30 31 32 33 34 36 37 38 39 40 42 43 44 45 47 48
[42] 49 51 52 54
[[10]]
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41
[42] 42 43 44 45 46 47 48 49 50 51 52 53 54 55
如果它很重要,你可以简单地加上"OK"迭代retrospesctively:
append("OK",
map(.x = accumulate(V_ind, c),
~ map2(.x,
rev(seq_along(.x) - 1),
function(y, z) seq(y - z, y, 1)) %>%
reduce(c)))
同样,如果您需要从原始向量中省略最后一个数字:
append("OK",
map(.x = accumulate(head(V_ind, -1), c),
~ map2(.x,
rev(seq_along(.x) - 1),
function(y, z) seq(y - z, y, 1)) %>%
reduce(c)))
for (i in seq_along(V_ind)) {
if (i == 1) {
print("ok")
} else if (i == 2) {
print(V_ind[1])
} else {
out_vector <- V_ind[seq(i - 1)]
max_minus <- i - 2
minus_indices <- rep(seq(max_minus), rev(seq(max_minus)) + 1)
minus_vector <- c()
for (j in rev(seq(max_minus))) {
minus_vector <- c(minus_vector, rev(seq(0, j)))
}
out_vector <- numeric(length(minus_vector))
for (k in seq_along(out_vector)) {
out_vector[k] <- V_ind[minus_indices[k]] - minus_vector[k]
}
print(c(out_vector, V_ind[i - 1]))
}
}
[1] "ok"
[1] 10
[1] 9 10 19
[1] 8 9 10 18 19 27
[1] 7 8 9 10 17 18 19 26 27 34
[1] 6 7 8 9 10 16 17 18 19 25 26 27 33 34 40
[1] 5 6 7 8 9 10 15 16 17 18 19 24 25 26 27 32 33 34 39 40 45
[1] 4 5 6 7 8 9 10 14 15 16 17 18 19 23 24 25 26 27 31 32 33 34 38 39 40 44 45 49
[1] 3 4 5 6 7 8 9 10 13 14 15 16 17 18 19 22 23 24 25 26 27 30 31 32 33 34 37 38 39 40 43 44 45 48 49 52
[1] 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17 18 19 21 22 23 24 25 26 27 29 30 31 32 33 34 36 37 38 39 40 42 43 44 45 47 48 49
[43] 51 52 54
可以显式地定义要减去的下标,以及要减去多少。(将+1加到seq中以减去0),然后只需添加最后一项(V_ind[i -1]),其中不对向量
进行减法操作。sequence
扮演关键角色的另一个选项
lapply(seq_along(x), function(n){
x[rep(1:n, n:1)] - rev(sequence(1:n) - 1)
})
# [[1]]
# [1] 10
#
# [[2]]
# [1] 9 10 19
#
# [[3]]
# [1] 8 9 10 18 19 27
#
# [[4]]
# [1] 7 8 9 10 17 18 19 26 27 34
其中x
是向量的子集:
x = cumsum(10:7)
如果需要,只需c
"ok">