選択された行が一定の最小時間だけ異なるように時系列をサブセット化

Question

私は時系列を保存するためにRのdata.tableを使用しています。私は、選択された時間の連続する行が、選択された最後の行から少なくともN秒離れているようなサブセットを返したいとします。私は

library(data.table) 
x <- data.table(t=c(0,1,3,4,5,6,7,10,16,17,18,20,21), v=1:13) 
x 
    t v 
1: 0 1 
2: 1 2 
3: 3 3 
4: 4 4 
5: 5 5 
6: 6 6 
7: 7 7 
8: 10 8 
9: 16 9 
10: 17 10 
11: 18 11 
12: 20 12 
13: 21 13 

を持っていると私は、最初の行から始まる離れて少なくとも5秒は、ある行をサンプリングしたい場合は、私は時間/値ペアでdata.tableを取得する必要があります：

y <- x[...something...] 
y 
    t v 
1: 0 1 
2: 5 5 
3: 10 8 
4: 16 9 
5: 21 13 

時間のサンプルは規則的に間隔を置く必要はないので、私はすべてのM行をとることはできません。もちろん、私はdata.tableの行を手動でループすることでこれを行うことができますが、data.tablesのインデックスを使用してこれを表現するより便利な方法があるのだろうかと思います。ここで

Answer 1

は、ローリングは、あなたのサブセットで、行、wのセットを見つけるために参加し、使用するカップルの方法があります。

t_plus = 5 

# one join per row visited 
w <- c() 
nxt <- 1L 
while(!is.na(nxt)){ 
    w <- c(w, nxt) 
    nxt <- x[.(t[nxt]+t_plus), on=.(t), roll=-Inf, which=TRUE] 
} 

# join once on all rows 
w0 <- x[.(t+5), on=.(t), roll=-Inf, which=TRUE] 

w <- c() 
nxt <- 1L 
while (!is.na(nxt)){ 
    w <- c(w, nxt) 
    nxt <- w0[nxt] 
} 

次にあなたがx[w]のようにサブセットすることができます。

---

コメント原則として

、OPの条件を満たし、他のサブセット「離れて、少なくとも5秒」があるかもしれません。これは最初の行から順方向に進むことによって見つけられたものです。

第2の方法は、上記のHenrikとリンクしたQ &の@DavidArenburg's answerに基づいています。質問は同じように見えますが、私はそのアプローチをここで完全に働かせることはできませんでした。

一般的に、私はwでやっているように、Rのループで物事を成長させるのは悪い考えです。パフォーマンスの問題に遭遇している場合は、このコードで改善するのがよいかもしれません。