2

我正在使用一个函数来识别一个序列,然后以分钟为单位计算序列的持续时间。当我在最后阶段将结果与数据 cbind 时,将返回持续时间,但返回的相邻列带有“NA”,而不是这些列中最初的值

d<-read.table(text='Date.Time Aerial
794  "2012-10-01 08:18:00"      1
795  "2012-10-01 08:34:00"      1
796  "2012-10-01 08:39:00"      1
797  "2012-10-01 08:42:00"      1
798  "2012-10-01 08:48:00"      1
799  "2012-10-01 08:54:00"      1
800  "2012-10-01 08:58:00"      1
801  "2012-10-01 09:04:00"      1
802  "2012-10-01 09:05:00"      1
803  "2012-10-01 09:11:00"      1
1576 "2012-10-01 09:17:00"      2
1577 "2012-10-01 09:18:00"      2
804  "2012-10-01 09:19:00"      1
805  "2012-10-01 09:20:00"      1
1580 "2012-10-01 09:21:00"      2
1581 "2012-10-01 09:23:00"      2
806  "2012-10-01 09:25:00"      1
807  "2012-10-01 09:32:00"      1
808  "2012-10-01 09:37:00"      1
809  "2012-10-01 09:43:00"      1', header=TRUE, stringsAsFactors=FALSE, row.names=1)
#Give correct data type
d$Aerial<- as.numeric(d$Aerial)
d$Date.Time<- as.POSIXct(d$Date.Time)

功能(识别天线 2 重复的序列和序列的持续时间):

fun1 <- function(data,aerial){
  data_above <- 1L*(data$Aerial == aerial)
  id_start <- paste(data$Date.Time[which(diff(c(0L,data_above))==1)])
  id_end <- paste(data$Date.Time[which(diff(c(data_above,0L))== -1)])
  res <- cbind(data[id_start,1:1],Duration=difftime(id_end,id_start, units='mins'))
  return(res)
}
fun1(d,2)

回报:

        Duration
[1,] NA        1
[2,] NA        2

持续时间是正确的,但是我希望它返回应该在相关列中的数据:

     Date.Time                     Duration
[1,] 2012-10-01 09:11:00            1
[2,] 2012-10-01 09:21:00            2

我的实际 data.frame 有很多列而不仅仅是 Date.Time 并且它仍然为所有这些返回 NA

4

1 回答 1

2

我会这样做:

fun1 <- function(data,aerial) {
    data_above <- 1L * (data$Aerial == aerial)
    id_start <- data$Date.Time[which(diff(c(0L,data_above)) == 1)]
    id_end <- data$Date.Time[which(diff(c(data_above, 0L)) == -1)]
    res <- cbind(data[data$Date.Time %in% id_start, 1, drop=FALSE], 
                  Duration = difftime(id_end,id_start, units='mins'))
   return(res)
}
fun1(d,2)

#                Date.Time Duration
# 1576 2012-10-01 09:17:00   1 mins
# 1580 2012-10-01 09:21:00   2 mins

这里要注意的点:

  • 当您对 a 进行子集化data.frame并且它只返回 1 个元素时,这样做df[, 1]会产生一个向量。使用是安全的df[, 1, drop = FALSE

  • 传递非data.frame参数(意思是,所有参数cbind都不是 data.frames)将导致输出为matrix. 它期望至少有一个参数是data.frame。所以,如果你不使用drop = FALSE并且输出结果在 1 行,那么它将是一个向量,结果将是一个矩阵(见第一点)

  • 我认为您不需要在paste这里使用id_startand id_end

  • cbind您访问的第一个参数data.frame data是不正确的。您必须id_start在 的所有值中查询Date.Time。可以使用%in%如图所示完成。

希望这可以帮助。

于 2013-03-29T17:14:01.327 回答