Re: [问题] 一笔资料转多笔

楼主: cywhale (cywhale)   2016-12-01 23:49:59
手痒也试了一下,就用data.table::foverlaps做
## data.table::foverlaps, reference: #1N_NLOIo
library(data.table)
library(magrittr)
dt <- fread('id start_y start_s end_y end_s
1 100 1 102 2
2 101 2 103 1
3 101 2 101 2') %>%
.[,{.(id=id,
start=as.numeric(paste(start_y,start_s,sep=".")),
end = as.numeric(paste(end_y,end_s,sep=".")))}] %>% setkey(start,end)
gt <- CJ(start_y=100:103, s=1:2) %>%
.[,{.(start=as.numeric(paste(start_y,s,sep=".")))}] %>% unique() %>%
.[,end:=start]
gx <- foverlaps(gt, dt, type="within", which=TRUE) %>%
.[which(!is.na(yid)),]
cbind(dt[gx$yid,.(id)], gt[gx$xid,.(start)]) %>% .[order(id),] %>%
.[,{.(id=id,
year=substr(start,1,nchar(start)-2),
s=substr(start,nchar(start),nchar(start))
)}]
# id year s
# 1: 1 100 1
# 2: 1 100 2
# 3: 1 101 1
# 4: 1 101 2
# 5: 1 102 1
# 6: 1 102 2
# 7: 2 101 2
# 8: 2 102 1
# 9: 2 102 2
#10: 2 103 1
#11: 3 101 2
※ 引述《criky (立业成家)》之铭言:
: [问题类型]:
:
: 程式咨询(我想用R 做某件事情,但是我不知道要怎么用R 写出来)
: [软件熟悉度]:
: 新手(没写过程式,R 是我的第一次)
: [问题叙述]:
: 若我有资料字段如下:
:
: id start_y start_s end_y end_s
: 1 100 1 102 2
: 2 101 2 103 1
: 3 101 2 101 2
: year:
: 如何转成下面的样子:
: id year s
: 1 100 1
: 1 100 2
: 1 101 1
: 1 101 2
: 1 102 1
: 1 102 2
: 2 101 2
: 2 102 1
: 2 102 2
: 2 103 1
: 3 101 2 (只有一笔)
: 谢谢回答~
: [程式范例]:
:
:
: [环境叙述]:
:
: 请提供 sessionInfo() 的输出结果,
: 里面含有所有你使用的作业系统、R 的版本和套件版本资讯,
: 让版友更容易找出错误
:
: [关键字]:
:
: 选择性,也许未来有用
:
作者: celestialgod (天)   2016-12-01 23:56:00
最后一行应该用tstrsplit就好CJ 好酷!!CJ在某种意义上是expand.grid,真是有趣的函数
作者: criky (2501-2)   2016-12-02 00:05:00
谢谢cy大,这篇解答好多括号= =
楼主: cywhale (cywhale)   2016-12-02 00:07:00
忘了还有tstrsplit 感谢提醒~~ ^^
作者: celestialgod (天)   2016-12-02 00:07:00
我帮忙排成简单易懂的好了@@算了,跟我习惯不同XDD
楼主: cywhale (cywhale)   2016-12-02 00:17:00
XDD 没关系 括号Rstudio会take care 我也常看不懂别人的排版<-其实是自己怪,另推data.table J, SJ, CJ好用
作者: Edster (Edster)   2016-12-02 11:37:00
你的排版真坚持

Links booklink

Contact Us: admin [ a t ] ucptt.com