代码之家 › 专栏 › 技术社区 › Omry Atia

提取向量中与索引最近的相邻索引

vectorization r

Omry Atia · 技术社区 · 7 年前

N ,和一个随机索引 ind 在这个向量中。我想以最对称的方式提取11个最接近 ind公司 在最有效的方法中,比如不超过原始向量的长度,也不会得到负指数。

如果 N = 200 ,和 ind = 199

ind = 3 ,结果应为1:11。

如果 ind = 195 ,结果应为190:200。

我可以和你一起做 if

ff = function(ind, N){
  vv = sapply(ind, function(x) c((x-5):(x+5)))
  if (max(vv) > N) {
    vv = vv - max(vv) + N
  } else if (min(vv) < 6) {
    vv = vv - min(vv) + 1
  }
  return(vv)
}

有没有更有效的方法(例如子集)来实现这一点?

1 回复 | 直到 7 年前

chinsoon12 7 年前

以下是几种方法:

1) 使用 order

sort(head(order(abs(1:N - ind)), 11L))

2) 排序索引和 ind 然后进入前11个位置

sort((1:N)[as.integer(names(head(sort(abs(setNames(1:N, 1:N) - ind)), 11L)))])

3) 使用基数排序:

sort(head(sort(abs(1:N - ind), index.return=TRUE)$ix, 11L))

推荐文章

Costinean Sebastian · 如何优化使用for循环处理大型数组的MATLAB脚本的执行时间?"

1 年前

Gyro · 为什么我的%xmm3寄存器在vbroadcastsd中使用第一个参数,而不是第四个参数?

1 年前

aeiou · 对日志表进行矢量化处理,以确定最新的可用性

1 年前

IchKenneDeinenNamen · 当将向量存储为类属性时,为什么NumPy的矢量化评估速度较慢?

2 年前

roshoka · 如何在不使用for循环的情况下将函数应用于林空间的每个元素

2 年前

Xavier · 变长矢量的矢量化和

2 年前

bbbb · 矢量化或在大型数据帧上使用多重处理[重复]

2 年前

Sterling Butters · NumPy 1D近邻中循环的避免

2 年前

AnthonyML · 修改2d numpy数组中分隔符后的行中的所有元素

2 年前

Tim · np.vectorize和relativelta返回“relativeltaonlydiffs datetime/date”

2 年前