基于指定间隔的地图值

时间:2016-08-31 04:03:58

标签: r data.table intervals

我们说我们有以下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

我想要的是将值dtnewdt中的区间匹配,以便我可以在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]

你能给我一些建议吗?

1 个答案:

答案 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)作为匹配行继续使用。