IGRAPH,POSIX,和data.table(igraph, POSIX, and data.t

2019-09-27 06:28发布

在较早前的问题 ,我了解到,图表来折叠这些数据非常有用

require(data.table)    
set.seed(333)
t <- data.table(old=1002:2001, dif=sample(1:10,1000, replace=TRUE))
t$new <- t$old + t$dif; t$foo <- rnorm(1000); t$dif <- NULL

> head(t)
    old  new        foo
1: 1002 1007 -0.7889534
2: 1003 1004  0.3901869
3: 1004 1014  0.7907947
4: 1005 1011  2.0964612
5: 1006 1007  1.1834171
6: 1007 1015  1.1397910

以获得仅那些行,使得new[i] = old[i-1] 那么结果可能会被加入到一个表,谁都有自己的出发点用户

i <- data.table(id=1:3, start=sample(1000:1990,3))

    > i
   id start
1:  1  1002
2:  2  1744
3:  3  1656

具体地,仅在第一当n=3 ,计算步骤,将溶液

> library(igraph)
> i[, t[old %in% subcomponent(g, start, "out")[1:n]], by=.(id)]
   id  old  new        foo
1:  1 1002 1007 -0.7889534
2:  1 1007 1015  1.1397910
3:  1 1015 1022 -1.2193666
4:  2 1744 1750 -0.1368320
5:  2 1750 1758  0.3331686
6:  2 1758 1763  1.3040357
7:  3 1656 1659 -0.1556208
8:  3 1659 1663  0.1663042
9:  3 1663 1669  0.3781835

在实施这个时候的设置是相同的,但newold ,并start有POSIXct类,

set.seed(333)
u <- data.table(old=seq(from=as.POSIXct("2013-01-01"), 
                          to=as.POSIXct("2013-01-02"), by="15 mins"), 
                dif=as.difftime(sample(seq(15,120,15),97,replace=TRUE),units="mins"))
u$new <- u$old + u$dif; u$foo <- rnorm(97); u$dif <- NULL
j <- data.table(id=1:3, start=sample(seq(from=as.POSIXct("2013-01-01"), 
                                       to=as.POSIXct("2013-01-01 22:00:00"), by="15 mins"),3))

> head(u)
                   old                 new        foo
1: 2013-01-01 00:00:00 2013-01-01 01:00:00 -1.5434407
2: 2013-01-01 00:15:00 2013-01-01 00:30:00 -0.2753971
3: 2013-01-01 00:30:00 2013-01-01 02:30:00 -1.5986916
4: 2013-01-01 00:45:00 2013-01-01 02:00:00 -0.6288528
5: 2013-01-01 01:00:00 2013-01-01 01:15:00 -0.8967041
6: 2013-01-01 01:15:00 2013-01-01 02:45:00 -1.2145590

> j
   id               start
1:  1 2013-01-01 22:00:00
2:  2 2013-01-01 21:00:00
3:  3 2013-01-01 13:30:00

命令

> j[, u[old %in% subcomponent(h, V(h)$name %in% as.character(start), "out")[1:n]], by=.(id)]
Empty data.table (0 rows and 4 cols): id,old,new,foo

返回一个空载体,这似乎是由于内部部分u[...] 我不太看到那里的问题是在这种情况下,不知道是否有人发现一个错误。

文章来源: igraph, POSIX, and data.table