确定超标的持续时间

作者: Colin
发布时间: 2024-08-17 09:45:44 (10天前)
转自：

2 条回复

0#
回复此人
尼拉。 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 如果您的时间序列是1分钟的时间序列（即：1分钟时间段的时间序列），则很容易使用超过某个阈值的间隔长度 <code> rle </code> ： </p> <pre> <code> xx = rle(turb >1.5) sum(xx$values==TRUE & xx$lengths >=15) </code> </pre> <P> 所以这里为了得到这个时间系列，一个解决方案是近似它以创建一个更准确的新时间系列。 </p> <pre> <code> library(xts) xx = xts(turb,time) yy = na.approx(merge(xts(,seq.POSIXt(min(time),max(time),by=1)), xx)) ## optional plot the new and the old time series plot(x = yy, xlab = "Time", minor.ticks = FALSE, col = "red") points(x = xx, col = "darkgreen",pch=20) </code> </pre> <P> <img src =“https://i.stack.imgur.com/ZNXKB.png”alt =“在此处输入图片说明”/> </p> <P> 然后我计算间隔的数量，如上所述： </p> <pre> <code> xx = rle(as.vector(coredata(yy>1.5))) sum(xx$values==TRUE & xx$lengths >=15) [1] 6 </code> </pre> <P> 注意：这里我发现只有6个间隔.. </p> </DIV>

编辑

登录后才能参与评论