R语言如何按指定步长从向量生成固定增量元素列表

针对你需要按100个基因的步长、逐次生成从100个元素到全长16000个元素的向量列表的需求，不需要调整accumulate()的逐元素运算逻辑，直接按索引切片的方式实现效率最高，完全适配后续逐次聚类的计算场景。

基础R实现（无额外依赖，性能最优） #

直接遍历步长序列，每次截取对应长度的向量片段即可，1.6w个基因的场景下运行耗时几乎可以忽略：

# 替换为你自己的基因向量即可，此处为示例数据
gene_vec <- paste0("gene_", 1:16000)
step_size <- 100
# 生成步长序列：100、200、300……直到覆盖向量总长度
length_seq <- seq(step_size, length(gene_vec), by = step_size)
# 按索引切片生成目标列表
gene_list <- lapply(length_seq, function(n) gene_vec[1:n])

# 如果需要从16000个基因逐次递减100个的序列，直接反转步长序列即可
# length_seq_rev <- seq(length(gene_vec), step_size, by = -step_size)
# gene_list_rev <- lapply(length_seq_rev, function(n) gene_vec[1:n])

生成的列表一共包含160个元素，第i个元素的向量长度恰好为i*100，完全匹配你的需求。

purrr 包实现（适配tidyverse使用习惯） #

如果你平时习惯用purrr系列函数，不需要硬套accumulate()的逐元素累加逻辑，直接用map()做切片即可，逻辑和基础R实现一致：

library(purrr)
gene_list_purrr <- map(
  .x = seq(100, length(gene_vec), by = 100),
  .f = ~gene_vec[1:.x]
)

如果一定要用accumulate()实现效果也可以，只是多了一步提前拆分的操作，运行效率略低于直接切片：提前把原始基因向量按每100个元素拆分成组，再让accumulate()每次拼接一整组元素即可。

# accumulate 实现方式（非必要不推荐）
gene_groups <- split(gene_vec, ceiling(seq_along(gene_vec)/100))
gene_list_acc <- accumulate(gene_groups, c)

聚类分析注意事项 #

如果你的实验设计是随机剔除100个基因后再聚类，不要直接按原始向量的固定顺序截取片段，不然会因为基因在向量中的固有排序引入偏差。提前打乱基因顺序再做切片即可，记得固定随机种子保证结果可复现：

set.seed(123) # 可替换为任意整数，固定后随机结果可重复
gene_shuffled <- sample(gene_vec)
gene_list_random <- lapply(seq(100, 16000, 100), function(n) gene_shuffled[1:n])

内容的提问来源于stack exchange，提问作者SHADJI