我有一些包含入院和出院日期的数据,我需要(对于每个人和承认 ID #)创建一个从宽到长的数据框,其中包含人员和承认 ID,以及承认和出院之间的顺序指示介于入院和出院之间(包括入院和出院)的每一天的值。
我偶然发现了一种非常骇人听闻的方法,它利用了 dlply 如何命名它创建的列表中的项目;但是,我担心它可能容易出错。无论哪种方式,它都感觉很笨拙,我想知道是否有一种使用更少代码的更简洁的方法。
person <- c(1, 2, 3, 3)
admit <- c(1, 1, 1, 2)
admit.date <- as.Date(c("1/1/2010", "1/1/2010", "1/1/2010", "2/1/2010"), "%m/%d/%Y")
discharge.date <- as.Date(c("1/1/2010", "1/1/2010", "1/1/2010", "2/1/2010"), "%m/%d/%Y") + c(1,2,2,2)
df1 <- data.frame(person, admit, admit.date, discharge.date)
df1 ## where I start
library(plyr)
los_seq <- function(df) { seq(df$admit.date, df$discharge.date, 1)}
lst1 <- dlply(df1, .(person, admit), los_seq)
vec1 <- unlist(lst1)
## now it gets really hackish
df2 <- data.frame(v1 = paste(names(vec1), vec1, sep="__"))
df2$person <- substr(df2$v1, 1, regexpr("\\.", df2$v1)-1)
df2$admit <- substr(df2$v1, regexpr("\\.", df2$v1)+1, regexpr("\\.", df2$v1)+1)
df2$date <- as.Date(as.numeric(substr(df2$v1, regexpr("__", df2$v1)+2, nchar(df2$v1))), origin="1970-01-01")
df2[,-1] ## this is how I need the result to look