21

在 Mathematica 中，有Clip[x, {min, max}] 给出xfor min<=x<=max、minforx<min和maxfor的命令x>max，请参见

http://reference.wolfram.com/mathematica/ref/Clip.html （镜像）

在 R 中实现这一目标的最快方法是什么？理想情况下，它应该是一个可列出的函数，并且理想情况下应该在单个值、向量、矩阵或数据帧上工作......

4

4 回答 4

25

Rcpp为此clamp：

cppFunction('NumericVector rcpp_clip( NumericVector x, double a, double b){
    return clamp( a, x, b ) ;
}')

这是一个快速基准测试，显示了它与讨论的其他方法相比如何执行：

pmin_pmax_clip <- function(x, a, b) pmax(a, pmin(x, b) )
ifelse_clip <- function(x, a, b) {
  ifelse(x <= a,  a, ifelse(x >= b, b, x))
}
operations_clip <- function(x, a, b) {
  a + (x-a > 0)*(x-a) - (x-b > 0)*(x-b)
}
x <- rnorm( 10000 )
require(microbenchmark)

microbenchmark( 
  pmin_pmax_clip( x, -2, 2 ), 
  rcpp_clip( x, -2, 2 ), 
  ifelse_clip( x, -2, 2 ), 
  operations_clip( x, -2, 2 )
)
# Unit: microseconds
#                        expr      min        lq   median        uq       max
# 1     ifelse_clip(x, -2, 2) 2809.211 3812.7350 3911.461 4481.0790 43244.543
# 2 operations_clip(x, -2, 2)  228.282  248.2500  266.605 1120.8855 40703.937
# 3  pmin_pmax_clip(x, -2, 2)  260.630  284.0985  308.426  336.9280  1353.721
# 4       rcpp_clip(x, -2, 2)   65.413   70.7120   84.568   92.2875  1097.039
于 2012-12-13T23:26:54.757 回答
20

这是一个嵌套pmin和pmax设置边界的方法：

 fenced.var <- pmax( LB, pmin( var, UB))

很难找到更快的方法。包装在默认范围为 3 和 7 的函数中：

fence <- function(vec, UB=7, LB=3) pmax( LB, pmin( vec, UB))

> fence(1:10)
 [1] 3 3 3 4 5 6 7 7 7 7
于 2012-12-13T22:45:35.180 回答
9

这是一个适用于向量和矩阵的函数。

myClip <- function(x, a, b) {
    ifelse(x <= a,  a, ifelse(x >= b, b, x))
}

myClip(x = 0:10, a = 3,b = 7)
#  [1] 3 3 3 3 4 5 6 7 7 7 7

myClip(x = matrix(1:12/10, ncol=4), a=.2, b=0.7)
# myClip(x = matrix(1:12/10, ncol=4), a=.2, b=0.7)
#      [,1] [,2] [,3] [,4]
# [1,]  0.2  0.4  0.7  0.7
# [2,]  0.2  0.5  0.7  0.7
# [3,]  0.3  0.6  0.7  0.7

这是另一个：

myClip2 <- function(x, a, b) {
    a + (x-a > 0)*(x-a) - (x-b > 0)*(x-b)
}

myClip2(-10:10, 0, 4)
# [1] 0 0 0 0 0 0 0 0 0 0 0 1 2 3 4 4 4 4 4 4 4
于 2012-12-13T21:51:19.153 回答
3

我相信那将clamp()来自raster包装。

library(raster)
clamp(x, lower=-Inf, upper=Inf, ...)
于 2018-01-10T03:19:44.693 回答