2016-07-28 7 views
0

私はこのようになります秒ごとの秒の値をCSV、持っている:R:欠落データのヒストグラム

"x","timestamp","value" 
"1",2016-01-01 00:00:00,124 
"2",2016-01-01 00:00:01,121 
"3",2016-01-01 00:00:02,NA 
"4",2016-01-01 00:00:03,NA 
"5",2016-01-01 00:00:04,NA 
"6",2016-01-01 00:00:05,123 
"7",2016-01-01 00:00:06,122 
"8",2016-01-01 00:00:07,124 
"9",2016-01-01 00:00:08,NA 
"10",2016-01-01 00:00:09,124 

だから、いくつかのデータが欠落しているとNAとしてマークされています。今私は欠けているデータブロックの長さのヒストグラムを作りたいと思う。与えられた例では、長さが1 sec (1)であり、2 sec (0)であり、残りが3 sec (1)であるなど、欠落しているデータブロックの数をカウントすることを意味します。私の実際のデータで

は、ビン/間隔は少し違うだろうし、私はこれらの8つのカテゴリを考える:

= 1 sec 
2 to 5 sec 
6 to 10 sec 
11 to 30 sec 
31 to 300 sec 
301 to 3600 sec 
3600 to 86400 sec 
> 86400 sec 

だから私の考えは、CSVのすべてのラインを介して実行RコードをできるようにしましたNAの値を検出するたびに、実際の値を再び見つけるまで行数をカウントします。 8つのカテゴリは、NA値のフィッティングブロックが検出されるたびに+1とカウントアップされる整数変数です。

完全なR-noobとして、私はそれをどうやって行うのか分かりません。ヘルプは高く評価されるだろう:)

+0

ランレングス符号化を使用してください。 '?rle'のヘルプを参照してください。 – Andrie

答えて

0

を私は(繰り返し可能な乱数を生成するset.seedを使用して)が時系列ソリューションでなければなりませんが、あなたが始めるために確信しています:

set.seed(42) 

# Create some sample data 
df <- data.frame(x = 1:100, 
       timestamp = seq(from = Sys.time() - 99, to = Sys.time(), by = "secs"), 
       value = sample(c(NA, 1:3), 100, replace = TRUE)) 

# Runs of identical data 
runs <- rle(is.na(df$value)) 

# Those that are missing 
missing <- which(runs$values) 

# The end positions in the sequence that are missing 
positions <- cumsum(runs$lengths) 

# The start times 
start <- df$timestamp[positions[missing] - runs$lengths[missing] + 1] 
end <- df$timestamp[positions[missing]] 

# Time difference 
delta <- difftime(end, start, "seconds") 

# Combine in a usable data.frame 
output <- data.frame(StartRow = positions[missing] - runs$lengths[missing] + 1, 
        EndRow = positions[missing], 
        StartTime = start, 
        EndTime = end, 
        Duration = delta) 
+0

すばらしい、ありがとう!それは完璧なベクトルを生産し続けています:) – Bene

0

多分これは便利

です
temp <- rle(diff(c(0,cumsum(is.na(df1$value))))) 
runs <- temp$lengths[temp$values==1] 
table(cut(runs,breaks = c(0,1,5,10,30,300,3600,86400,Inf),right = T)) 
hist(runs,breaks = c(1,5,10,30,300,3600,86400))