我们说我们有以下data.table
:
library(data.table)
dt <- data.table(x=c(-0.01, -0.001, 0, 0.01,0.02,0.03,1,3,4,10,20,25), value=c(rep(1,3), rep(2, 3), rep(3, 2), rep(4, 2), rep(5, 2)))
然后,我使用以下内容指定间隔:
library(Hmisc)
dt[, int:=cut2(x, cuts=c(-Inf, 0.01, 1, 4, 20, Inf))]
返回,
x value int
1: -0.010 1 [ -Inf, 0.01)
2: -0.001 1 [ -Inf, 0.01)
3: 0.000 1 [ -Inf, 0.01)
4: 0.010 2 [ 0.01, 1.00)
5: 0.020 2 [ 0.01, 1.00)
6: 0.030 2 [ 0.01, 1.00)
7: 1.000 3 [ 1.00, 4.00)
8: 3.000 3 [ 1.00, 4.00)
9: 4.000 4 [ 4.00,20.00)
10: 10.000 4 [ 4.00,20.00)
11: 20.000 5 [20.00, Inf]
12: 25.000 5 [20.00, Inf]
假设有新的data.table
dtnew <- data.table(x=c(-0.001, 0.4, 0.3, 5, 25))
x
1: -0.001
2: 0.400
3: 0.300
4: 5.000
5: 25.000
我想要的是将值dtnew
与dt
中的区间匹配,以便我可以在value
中获得dt
。预期的回报是
x value int
1: -0.001 1 [ -Inf, 0.01)
2: 0.400 2 [ 0.01, 1.00)
3: 0.300 2 [ 0.01, 1.00)
4: 5.000 4 [ 4.00,20.00)
5: 25.000 5 [20.00, Inf]
我尝试将时间间隔分配给dtnew
,以便merge(dtnew, dt, by = 'int')
使用以下内容:
dtnew[, int:=cut2(y, cuts=c(-Inf, 0.01, 1, 4, 20, Inf))]
的时间间隔与dt
中使用的时间间隔相同。但是,第1行的间隔与dt
不同,如下所示:
1: -0.001 -Inf
2: 0.400 [ 0.01, 1.00)
3: 0.300 [ 0.01, 1.00)
4: 5.000 [ 4.00,20.00)
5: 25.000 [20.00, Inf]
你能给我一些建议吗?
答案 0 :(得分:4)
这是一个简单的滚动连接:
dt[dtnew, on="x", roll=Inf]
# x value int
# 1: -0.001 1 [ -Inf, 0.01)
# 2: 0.400 2 [ 0.01, 1.00)
# 3: 0.300 2 [ 0.01, 1.00)
# 4: 5.000 4 [ 4.00,20.00)
# 5: 25.000 5 [20.00, Inf]
对于dtnew
中的每一行,通过查找dt
找到x
中的匹配行。如果来自dtnew
的特定值不匹配,而是属于间隙,即dt
中的两个值之间,则转发上一个/上一个值
,对于dtnew$x = 0.4
,它位于第6行0.03
和第7行1.00
之间。因此前一行(= 6)作为匹配行继续使用。