2

次の問題があります。私がベクトルを持っているとしましょう：

v = [1,2,3,4,5,1,2,3,4,...]

以前にサンプリングされたポイントからのしきい値よりも絶対的なマグニチュードの差が大きい、ベクトルからポイントを順番にサンプリングしたいと思います。つまり、私のしきい値は2だとしましょう。

インデックス1から開始し、最初のポイント1をサンプリングします。次にv [3]で条件が満たされ、3をサンプリングします（3-1> = 2であるため）。次に3、新しいサンプリングされたポイントが参照になり、それをチェックします。次のサンプリングされたポイントは5で、これはv [5]（5-3> = 2）です。次に、次のポイントはv [6]（abs（1-5）> = 2）である1です。

残念ながら、Rの私のコードは時間がかかりすぎています。基本的に、私はアレイを繰り返しスキャンして、一致するものを探しています。しかし、このアプローチは素朴だと思います。このタスクは、アレイを1回通過するだけで実行できると感じています。でも方法はわかりません。助けていただければ幸いです。私が直面している問題は、次のサンプルポイントの場所が配列内のどこにでもある可能性があることであり、現在のポイントから最後まで配列をスキャンして見つける必要があると思います。

ありがとう。

4

2 に答える 2

12

ループなしでこれを行う方法がわかりません。そのため、次の方法があります。

my.sample <- function(x, thresh) {
   out <- x
   i <- 1
   for (j in seq_along(x)[-1]) {
      if (abs(x[i]-x[j]) >= thresh) {
         i <- j
      } else {
         out[j] <- NA
      }
   }
   out[!is.na(out)]
}

my.sample(x = c(1:5,1:4), thresh = 2)
# [1] 1 3 5 1 3
于 2012-11-21T02:00:40.143 に答える
1

少しの再帰を使用して、ループなしでこれを行うことができます。

vsearch = function(v, x, fun=NULL) {
  # v: input vector
  # x: threshold level

  if (!length(v) > 0) return(NULL)

  y = v-rep(v[1], times=length(v))
  if (!is.null(fun)) y = fun(y)

  i = which(y >= x)

  if (!length(i) > 0) return(NULL)
  i = i[1]

  return(c(v[i], vsearch(v[-(1:(i-1))], x, fun=fun)))
}

上記のベクトルで：

> vsearch(c(1,2,3,4,5,1,2,3,4), 2, abs)
[1] 3 5 1 3
于 2012-11-21T08:22:33.513 に答える