1

假设我想根据向量中近似x出现的值来过滤向量:y

x <- c(1.123456789, 2.123456789, 3.123456789)
y <- c(1.12345, 2.12345)

如果我不想进行近似比较,我会使用%in%

x %in% y
[1] FALSE FALSE FALSE

我需要的结果是:

# something like: x %near_in% y
[1] TRUE TRUE FALSE

dplyr::near(x, y, tol)函数的帮助文件中提到了“x, y: Numeric vectors to compare”,但这当然不完全正确,y必须是x' 长度或单个值,因为near()所做的就是使用该abs()函数:

near <- function (x, y, tol = .Machine$double.eps^0.5) 
{
    abs(x - y) < tol
}

如果我们这样做,我们会看到abs()takey的值并重复它们直到它需要(不是没有警告),我们得到:

abs(x - y)
[1] 0.000006789 0.000006789 2.000006789
Warning message:
In x - y : longer object length is not a multiple of shorter object length

我目前的解决方案是使用sapply()ony的元素来创建一个n x m矩阵(此处为 3 x 2),然后用于apply()查看其中是否any()有行(的值xTRUE

apply(sapply(y, function(y_val) near(x, y_val, 0.0001)), 1, any)
[1] TRUE TRUE FALSE

但这似乎很麻烦!如果我有数千个值怎么办,y我不会创建一个包含数千行的临时矩阵吗?有更好的办法吗?

4

2 回答 2

2

您可以floorround值:

tol <- 1e-5
floor(x/tol)
#> [1] 112345 212345 312345
floor(y/tol)
#> [1] 112345 212345
floor(x/tol) %in% floor(y/tol)
#> [1]  TRUE  TRUE FALSE
于 2021-09-05T06:20:28.703 回答
0

我不太确定。但请查看DescTools库,您可以在其中找到最接近的值:

library(DescTools)

y %in% sapply(x, function(i) Closest(y, i))

[1] TRUE TRUE
于 2021-09-05T06:34:21.923 回答