가장 빠른 방법은 각 간격에서 함수의 최대 값을 찾는 것입니다.

Question

저는 벡터가 v이고 간격 벡터가 w입니다. 각각의 간격에서 함수의 최대 값, f(x)을 찾고 싶습니다. 결과를 찾는 코드를 따르는 것보다 빠른 방법이 있습니까? 예를 들면 :

v = c(3.5, 2.5, 4, 6.5, 10, 2.3, 1.8, 4.7, 12, 11.5) 
w = c(0, 5, 15, 20) 
f = function(x){x^2} 
> max = unlist(list(sapply(split(v, cut(v, w),drop = TRUE), 
       function(v) v[which.max(f(v))])), use.names = FALSE) 
> max 
[1] 4.7 12.0 

Answer 1

무엇 findInterval 및 tapply에 대해. findInterval는 cut 비슷하지만 최대 값

tapply(f(v),findInterval(v,w),max) 
# 1  2 
# 22.09 144.00 

을 원하는 아니면 사실을 사용할 수 있다면 요인

tapply(v,findInterval(v,w),function(x)x[which.max(f(x))]) 
# 1 2 
# 4.7 12.0 

또는 다음으로 변환하는 오버 헤드없이 기능은 단조 모든 양의 값을 증가하고 있음 하고. 당신이 경계에서 발생하는 지정해야합니다

f(tapply(v,findInterval(v,w),max)) 

참고 w`는 3 개 세트 (하나가 비어)를 정의`것을

library(microbenchmark) 
    microbenchmark(
     mnel = tapply(v,findInterval(v,w),max), 
     flodel = unname(vapply(split(f(v), cut(v, w), drop = TRUE), max, numeric(1L))), 
     flodel2 = unname(vapply(split(seq_along(v), findInterval(v, w)), function(i, v, fv)v[i][which.max(fv[i])], numeric(1L), v, f(v)))) 
# Unit: microseconds 
# expr  min  lq median  uq  max neval 
# mnel 260.945 262.9155 264.2265 276.0645 458.670 100 
# flodel 331.218 334.3585 336.0580 351.1985 694.715 100 
#flodel2 124.998 127.3230 128.5170 137.0505 354.545 100