2013-01-15 74 views
0

通过时间间隔

在这里变换数据设定为具有一排是用于个别的一些数据与id = 1

id time status 
-------------- 
1 t status 

t是时间的一些事件,并status是要么如果发生事件则为1或如果没有发生则为0(在这种情况下,t是研究的持续时间)。

假设t位于a2a3之间。

我的目标是我的数据转换成如下:

id period start stop status 
--------------------------- 
1 1  0  a1 0  
1 2  a1 a2 0  
1 3  a2 t status 

个人1的总时间分为三个区间,其中存在(0, a1)和任何情况下(a1, a2)

问题

你能想到一种有效的方法来编写一个R函数,它输入一个数据集和一个向量,并输出转换后的数据集?


编辑

第1部分 我一直在问一个具体的例子。这里是一个:

id time status 
    -------------- 
    1 5 1 

a1=1a2=3a3=7

第2部分我也被要求展示我的尝试。这是

> data <- data.frame(id=1, time=5, status=1) 
> a <- c(1, 3, 7) 
> N <- nrow(data) 
> data$period <- ifelse(data$time < a[1], 1, 
+      ifelse(data$time < a[2], 2, 
+        ifelse(data$time < a[3], 3, 4))) 
> 
> 
> dataTemp1 <- data.frame(matrix(nrow=N, ncol=ncol(data))) 
> names(dataTemp1) <- names(data) 
> dataTemp2 <- data.frame(matrix(nrow=N, ncol=ncol(data))) 
> names(dataTemp2) <- names(data) 
> dataTemp3 <- data.frame(matrix(nrow=N, ncol=ncol(data))) 
> names(dataTemp3) <- names(data) 
> dataTemp4 <- data.frame(matrix(nrow=N, ncol=ncol(data))) 
> names(dataTemp4) <- names(data) 
> 
> for(j in 1:N) 
+ { 
+ if(data[j, "period"] == 1){ 
+  data[j, "start"] <- 0 
+  data[j, "stop"] <- data[j, "time"] 
+ } else if(data[j, "period"] == 2){ 
+  dataTemp1[j, c("id", "time", "period")] <- 
+  data[j, c("id", "time", "period")] 
+  dataTemp1[j, "start"] <- 0 
+  dataTemp1[j, "stop"] <- a[1] 
+  dataTemp1[j, "status"] <- 0 
+  
+  data[j, "start"] <- a[1] 
+  data[j, "stop"] <- data[j, "time"] 
+ } else if(data[j, "period"] == 3){ 
+  dataTemp1[j, c("id", "time", "period")] <- 
+  data[j, c("id", "time", "period")] 
+  dataTemp1[j, "start"] <- 0 
+  dataTemp1[j, "stop"] <- a[1] 
+  dataTemp1[j, "status"] <- 0 
+  
+  dataTemp2[j, c("id", "time", "period")] <- 
+  data[j, c("id", "time", "period")] 
+  dataTemp2[j, "start"] <- a[1] 
+  dataTemp2[j, "stop"] <- a[2] 
+  dataTemp2[j, "status"] <- 0 
+  
+  data[j, "start"] <- a[2] 
+  data[j, "stop"] <- data[j, "time"]  
+ } else if(data[j, "period"] == 4){ 
+  dataTemp1[j, c("id", "time", "period")] <- 
+  data[j, c("id", "time", "period")] 
+  dataTemp1[j, "start"] <- 0 
+  dataTemp1[j, "stop"] <- a[1] 
+  dataTemp1[j, "status"] <- 0 
+  
+  dataTemp2[j, c("id", "time", "period")] <- 
+  data[j, c("id", "time", "period")] 
+  dataTemp2[j, "start"] <- a[1] 
+  dataTemp2[j, "stop"] <- a[2] 
+  dataTemp2[j, "status"] <- 0 
+  
+  dataTemp3[j, c("id", "time", "period")] <- 
+  data[j, c("id", "time", "period")] 
+  dataTemp3[j, "start"] <- a[2] 
+  dataTemp3[j, "stop"] <- a[3] 
+  dataTemp3[j, "status"] <- 0 
+  
+  data[j, "start"] <- a[3] 
+  data[j, "stop"] <- data[j, "time"] 
+ } 
+ } 
> 
> dataTemp1 <- dataTemp1[complete.cases(dataTemp1), ] 
> dataTemp2 <- dataTemp2[complete.cases(dataTemp2), ] 
> dataTemp3 <- dataTemp3[complete.cases(dataTemp3), ] 
> dataTemp4 <- dataTemp4[complete.cases(dataTemp4), ] 
> 
> data <- rbind(data, dataTemp1, dataTemp2, dataTemp3, dataTemp4) 
> data[, "period"] <- ifelse(data[, "start"] == 0, 1, 
+       ifelse(data[, "start"] == a[1], 2, 
+         ifelse(data[, "start"] == a[2], 3, 
+           ifelse(data[, "start"] == a[3], 4, 
+             5)))) 
> data <- data[order(data$id, data$start), 
+    c("id", "period", "start", "stop", "status")] 
> data 
    id period start stop status 
2 1  1  0 1  0 
3 1  2  1 3  0 
1 1  3  3 5  1 
+0

你应该提供一个可重复的例子。 ai是什么日期?为什么不提供一些数值,而不仅仅是符号,还可以显示你尝试过什么? – agstudy

+0

@agstudy:我做了编辑。但是,我想要一个功能而不是一个只适用于一个例子的程序。 – user7064

+0

@Arun:Wahou,thx!如果你让它成为答案,我会接受它! – user7064

回答

0

我将它写成一个适当的重复性解决方案:

df <- data.frame(id=1, time=5, status=2) 
a <- c(1, 3, 7) 

res.fn <- function(df, a) { 
    id <- rep(1, length(a)) 
    period <- 1:length(a) 
    start <- c(0, a[1:(length(a)-1)]) 
    stop <- c(a[1:(length(a)-1)], df$time) 
    status <- c(rep(0, length(a)-1), df$status) 
    data.frame(id, period, start, stop, status) 
} 
> res.fn(df, a) 

    id period start stop status 
1 1  1  0 1  0 
2 1  2  1 3  0 
3 1  3  3 5  2 
+0

这很完美! – user7064

+0

是的,当然,我会的。谢谢 – user7064