我想从点(1,1)中仅过滤小于10个单位的行。我的数据框有两列,x和y。
这就是我的尝试:
filter(df, dist( rbind(c(1,2), c(x,y)) ) < 10 )
但是,这不起作用。它总是返回0行结果,虽然我知道它应该返回几行。我该怎么调试呢?我想在每次迭代中打印传递给x和y的每个值。
每个请求,这是dput(head(df))的输出:
structure(list(x = c(1, 2, 3, 4, 5), y = c(1, 1, 1, 1, 1)), .Names = c("x",
"y"), row.names = c(NA, 5L), class = "data.frame")
答案 0 :(得分:5)
我会使用您的数据但不受过滤器的影响。所以我会随机创建一些东西:
library(dplyr)
set.seed(42)
df <- data_frame(x = sample(20, size = 20, replace = TRUE),
y = sample(20, size = 20, replace = TRUE))
head(df)
# Source: local data frame [6 x 2]
# x y
# <int> <int>
# 1 19 19
# 2 19 3
# 3 6 20
# 4 17 19
# 5 13 2
# 6 11 11
问题在于dplyr::filter
需要逻辑向量。如果您手动检查dist(...)
的返回值,则会返回&#34; n-by-n&#34;阵列。目前尚不清楚filter
应该如何假设使用它。
如果您的数据实际上只是一个点(c(1, 2)
),那么您需要手动计算已知点与data.frame变量之间的距离,例如:
filter(df, sqrt( (x - 1)^2 + (y - 2)^2 ) < 10)
# Source: local data frame [2 x 2]
# x y
# <int> <int>
# 1 10 1
# 2 3 5
(我在这里假设欧几里德距离。)如果你有更多的尺寸和/或稍微不同的距离方程,应用程序应该是直截了当的。
如果您对df
中所有点之间的距离感兴趣(正如您对dist
所说的那样),那么您可能需要使用which(..., arr.ind = TRUE)
和一些技巧。或者也许在这些(df
)点和其他点之间进行外连接。