このR関数を組み込み関数で再表現しますか？

Question

[i,j]（y1と呼ぶ）と値i:jのベクトル（y2と呼ぶ）の一連のセンチネル値があるとします。たとえば：

y1 <- c(649, 652, 656, 658, 660, 661, 669, 671) 
y2 <- 649:671 

それとも

y1 <- c(719, 725, 726, 728, 729, 731, 734, 740) 
y2 <- 718:740 

ありませんが、この機能を考えてみます。

Cody5 <- function(y1, y2){ 
    y3 <- rep(NA, length(y2)) 
    j <- 1 
    i <- 0 
    while(j < length(y1)){ 
     i <- i + 1 
     if(y2[i]<y1[j]){ 
      y3[i]<-y1[j] 
     } else { 
      y3[i]<-y1[j<-j+1] 

     } 
    } 
    if(i < length(y3)){ 
     for(l in i:length(y3)){ 
      y3[l] <- y1[j] 
     } 
    } 
    y3 
} 

私はcut()を使用してCody5を交換しようとしているが、私はそれが正しいように見えることはできません。

編集：

私はそれを書いておく必要があります。番兵（y1）と元のベクトルのエントリがにソートされています。

Answer 1

あなたの機能のより効率的な実装は次のようになります。例えば、データと

y1[findInterval(y2, y1, rightmost.closed=TRUE)+1] 

：

y1 <- c(649, 652, 656, 658, 660, 661, 669, 671) 
y2 <- 649:671 
Cody5(y1, y2) 
# [1] 652 652 652 656 656 656 656 658 658 660 660 661 669 669 669 669 669 669 669 669 671 671 671 
y1[findInterval(y2, y1, rightmost.closed=TRUE)+1] 
# [1] 652 652 652 656 656 656 656 658 658 660 660 661 669 669 669 669 669 669 669 669 671 671 671 

Answer 2

mycut <- function(x, y) { 
    x[cut(y, c(-Inf, x), 
      labels = FALSE, right = FALSE, 
      include.lowest = TRUE)] 
} 

y1 <- c(649, 652, 656, 658, 660, 661, 669, 671) 
y2 <- 649:671 

all.equal(Cody5(y1,y2), 
      mycut(y1,y2)) 
#[1] TRUE 

y1 <- c(719, 725, 726, 728, 729, 731, 734, 740) 
y2 <- 718:740 

all.equal(Cody5(y1,y2), 
      mycut(y1,y2)) 
#[1] TRUE