[R-br] Obtenção de intervalos de tempos
Éder Comunello
comunello.eder em gmail.com
Quinta Abril 23 10:29:39 BRT 2015
Bom dia,
Existe a função rle() que é útil nesses casos.
### <code r>
dat <- read.table(text='Horario Nivel Status
"01/01/2011 00:00" 29 0
"01/01/2011 00:01" 29 0
"01/01/2011 00:02" 30 1
"01/01/2011 00:03" 31 1
"01/01/2011 00:04" 29 0
"01/01/2011 00:05" 29 0
"01/01/2011 00:06" 29 0
"01/01/2011 00:07" 29 0
"01/01/2011 00:08" 30 1
"01/01/2011 00:09" 32 1
"01/01/2011 00:10" 34 1
"01/01/2011 00:11" 35 1
"01/01/2011 00:12" 35 1
"01/01/2011 00:13" 25 0', head=T, as.is=T)
dat$Horario <- as.POSIXlt(dat$Horario, format = '%d/%m/%Y %H:%M')
head(dat)
res <- do.call(data.frame, rle(dat$Status))
n <- nrow(res)
res$ini <- c(1, cumsum(res$lengths[-n])+1)
res$fim <- cumsum(res$lengths)
res$timeDiff <- with(res, dat[fim,1]-dat[ini,1])
res
# lengths values ini fim timeDiff
# 1 2 0 1 2 60 secs
# 2 2 1 3 4 60 secs
# 3 4 0 5 8 180 secs
# 4 5 1 9 13 240 secs
# 5 1 0 14 14 0 secs
### encontrar os intervalos de tempo em que a variavel "Status" permanence
no valor 1.
res[res$values==1,]
# lengths values ini fim timeDiff
# 2 2 1 3 4 60 secs
# 4 5 1 9 13 240 secs
### </code>
Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150423/98bf0e1c/attachment.html>
Mais detalhes sobre a lista de discussão R-br