按日期范围创建虚拟变量

时间:2017-09-07 07:41:05

标签: r merge dummy-variable

我想创建虚拟变量,根据日期范围(即发布日期和下一版本发布日期之间)指示哪个产品版本生效,而不手动执行。我有几百个版本及其相应的发布日期,当下一个版本发布时,我将合并/加入这些数据与Purchases数据框。

目前我有两个数据框(版本和购买),如下所示:

View(Version)   
Type Version    Release_Date    Next_Release
A       1.2.3   2013-11-14     2014-01-11
B       1.3.1   2014-01-11     2014-02-20
A       1.5.1   2014-02-20     2014-03-08
A       1.5.2   2014-03-08     2014-04-06
B       1.5.3   2014-04-06     2014-04-12
A       1.5.4   2014-04-12     2014-04-15
B       1.5.5   2014-05-15     2014-05-20
B       1.6.1   2014-05-20     2014-06-26
A       1.6.2   2014-06-26     2014-07-14

View(Purchases)
TIMESTAMP   Amount 
2013-11-14   15.44
2013-11-14   13.39
2013-11-14   15.35
2013-11-15   86.43
2014-01-15   12.30
2014-01-17   23.55

我想为数据框中的每个版本创建虚拟变量'购买'根据版本生效的日期范围。即如果购买的时间戳在该版本生效的日期范围内,则该版本哑= 1,否则为0.

View(Purchases)
TIMESTAMP   Amount Version_1.2.3  Version_1.3.1  ....
2013-11-14   15.44      1              0
2013-11-14   13.39      1              0
2013-11-14   15.35      1              0
2013-11-15   86.43      1              0
2014-01-15   12.30      0              1
2014-01-17   23.55      0              1 
....

提前感谢您的任何建议或帮助。

1 个答案:

答案 0 :(得分:1)

任务可分为两个子任务:

  1. 在购买Release_Date之前找到最新的TIMESTAMP
  2. 从长格式转变为宽格式。
  3. 为了找到最新版本,data.table包提供了两种可能性,

    滚动加入

    library(data.table)
    setDT(Version)
    setDT(Purchases)
    
    Version[Purchases, on = .(Release_Date = TIMESTAMP), roll = TRUE,
            .(TIMESTAMP = i.TIMESTAMP, Amount, Version)]
    
        TIMESTAMP Amount Version
    1: 2013-11-14  15.44   1.2.3
    2: 2013-11-14  13.39   1.2.3
    3: 2013-11-14  15.35   1.2.3
    4: 2013-11-15  86.43   1.2.3
    5: 2014-01-15  12.30   1.3.1
    6: 2014-01-17  23.55   1.3.1
    

    非等于加入

    Version[Purchases, on = .(Release_Date <= TIMESTAMP), mult = "last",
            .(TIMESTAMP = i.TIMESTAMP, Amount, Version)]
    

    产生相同的结果。

    对于重塑,dcast()函数与length()一起用于聚合:

    # rolling join
    Version[Purchases, on = .(Release_Date = TIMESTAMP), roll = TRUE,
            .(TIMESTAMP = i.TIMESTAMP, Amount, Version)][
      , dcast(.SD, TIMESTAMP + Amount ~ Version, length)]
    
        TIMESTAMP Amount 1.2.3 1.3.1
    1: 2013-11-14  13.39     1     0
    2: 2013-11-14  15.35     1     0
    3: 2013-11-14  15.44     1     0
    4: 2013-11-15  86.43     1     0
    5: 2014-01-15  12.30     0     1
    6: 2014-01-17  23.55     0     1
    

    或者,如果在重塑期间重命名列

    # non-equi join
    Version[Purchases, on = .(Release_Date <= TIMESTAMP), mult = "last",
            .(TIMESTAMP = i.TIMESTAMP, Amount, Version)][
      , dcast(.SD, TIMESTAMP + Amount ~ paste0("Version_", Version), length)]
    
        TIMESTAMP Amount Version_1.2.3 Version_1.3.1
    1: 2013-11-14  13.39             1             0
    2: 2013-11-14  15.35             1             0
    3: 2013-11-14  15.44             1             0
    4: 2013-11-15  86.43             1             0
    5: 2014-01-15  12.30             0             1
    6: 2014-01-17  23.55             0             1