r - 如何为时间序列中的缺失数据创建 "NA"

我有几个看起来像这样的数据文件:

X code year month day pp
1 4515 1953     6   1  0
2 4515 1953     6   2  0
3 4515 1953     6   3  0
4 4515 1953     6   4  0
5 4515 1953     6   5  3.5

有时缺少数据，但是我没有NA，这些行根本不存在。数据丢失时，我需要创建NA。我虽然可以通过将其转换为Zoo对象并检查严格的规则性(我以前从未使用过Zoo)来确定何时发生，但我使用了以下代码:

z.date<-paste(CET$year, CET$month, CET$day, sep="/")
z <- read.zoo(CET,  order.by= z.date )
reg<-is.regular(z, strict = TRUE)

但是答案永远是正确的!

谁能告诉我为什么不起作用？甚至更好，请告诉我一种在数据丢失时创建NA的方法(带或不带Zoo包)？

谢谢

最佳答案

seq函数具有一些有趣的功能，可用于轻松生成完整的日期序列。例如，以下代码可用于生成从4月25日开始的日期序列:

编辑:此功能记录在?seq.Date中

start = as.Date("2011/04/25")
full <- seq(start, by='1 day', length=15)
full

 [1] "2011-04-25" "2011-04-26" "2011-04-27" "2011-04-28" "2011-04-29"
 [6] "2011-04-30" "2011-05-01" "2011-05-02" "2011-05-03" "2011-05-04"
[11] "2011-05-05" "2011-05-06" "2011-05-07" "2011-05-08" "2011-05-09"

现在，使用相同的原理通过每2天生成一次序列来生成一些“缺失”行的数据:

partial <- data.frame(
    date=seq(start, by='2 day', length=6),
    value=1:6
)
partial

        date value
1 2011-04-25     1
2 2011-04-27     2
3 2011-04-29     3
4 2011-05-01     4
5 2011-05-03     5
6 2011-05-05     6

要回答您的问题，可以使用向量下标或match函数创建具有NA的数据集:

with(partial, value[match(full, date)])
 [1]  1 NA  2 NA  3 NA  4 NA  5 NA  6 NA NA NA NA

要将结果与原始完整数据结合起来，请执行以下操作:

data.frame(Date=full, value=with(partial, value[match(full, date)]))
         Date value
1  2011-04-25     1
2  2011-04-26    NA
3  2011-04-27     2
4  2011-04-28    NA
5  2011-04-29     3
6  2011-04-30    NA
7  2011-05-01     4
8  2011-05-02    NA
9  2011-05-03     5
10 2011-05-04    NA
11 2011-05-05     6
12 2011-05-06    NA
13 2011-05-07    NA
14 2011-05-08    NA
15 2011-05-09    NA

关于r - 如何为时间序列中的缺失数据创建 "NA"，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/6058677/