r - 按列从矩阵中提取 block

Question

假设我有一个包含 1000 列的矩阵。我想每隔一个 n 创建一个新矩阵来自原始矩阵的列，从列 i 开始.

所以可以说 n=3和 i=5 ，那么旧矩阵中我需要的列是 5,6,7,11,12,13,17,18,19等等。

Answer 1

使用两个 seq() s 创建开始和停止边界，然后使用 mapply()建立真正的列索引间隔。然后只是普通的括号符号从你的矩阵中提取。

set.seed(1)
# using 67342343's test case
M <- matrix(runif(100^2), ncol = 100)
n <- 3
i <- 5

starts <- seq(i, ncol(M), n*2)
stops <- seq(i+(n-1), ncol(M), n*2)
col_index <- c(mapply(seq, starts, stops)) # thanks Jaap and Sotos

col_index
[1]  5  6  7 11 12 13 17 18 19 23 24 25 29 30 31 35 36 37 41 42 43 47 48 49 53 54 55 59 60 61 65 66 67 71 72 73 77 78
[39] 79 83 84 85 89 90 91 95 96 97

M[, col_index]