R 汇总POSIX时间序列中缺失数据的间隔

我有一个小时降水量数据的时间序列，我正试图通过QA/QC例行程序运行该数据。我想做的一件事是创建一个缺失数据间隔计数的柱状图，该柱状图根据数据缺失的时间长度分组，即缺失一小时数据的时段有多少，连续两小时的时段有多少，三小时的时段有多少，等等……我可能可以用一些嵌套循环来实现这一点，但我想知道是否有更好的方法

时间序列是连续的（表示所有小时）。datetime为POSIXct，数据为数字，缺失数据为NA。可以使用以下内容创建一个简短示例：

precip <- structure(list(datetime = structure(c(1114905600, 1114909200, 
1114912800, 1114916400, 1114920000, 1114923600, 1114927200, 1114930800, 
1114934400, 1114938000, 1114941600, 1114945200, 1114948800), class = c("POSIXct", 
"POSIXt"), tzone = "UTC"), precip = c(1.1, NA, 2, 0, NA, NA, 
NA, 0, 0, NA, NA, 0.5, 0.3)), .Names = c("datetime", "precip"
), row.names = c(NA, -13L), class = "data.frame")

precip使用
rle
（行程编码）：

R>rle\u res rle\u na表格（rle\u na$length）
1 2 3 
1 1 1 

R> rle_res <- as.data.frame(unclass(rle(is.na(precip$precip))))
R> rle_na <- subset(rle_res, values==TRUE)  # filter NA values
R> table(rle_na$length)
1 2 3 
1 1 1