QA/QC ルーチンを実行しようとしている 1 時間ごとの降水量データの時系列があります。私がやりたいことの 1 つは、データが欠落している時間の長さでグループ化された、欠落しているデータのある間隔のカウントのヒストグラムを作成することです。つまり、1 時間のデータが欠落している期間がいくつあるか、2 時間連続して欠落している期間がいくつあるか、 3 つの場合はいくつになるかなど...ネストされたループでこれを行うことができると思いますが、もっと良い方法があるかどうか疑問に思います。
時系列は連続しています (すべての時間が表されます)。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")
出力では、1 時間の期間が 1 つ、2 時間の期間が 1 つ、3 時間の期間が 1 つ、欠損データとして認識されます。ありがとう！