[英]Replace missing data with NA in R
我正在從存儲在 xts 對象列表中的幾個觀察構建一個時間序列 (xts)。 提取的數據有時會丟失,R報錯:
"Error in NextMethod(.Generic) : replacement has length zero"
我希望 R 改為報告 NA。 我想答案在於 tryCatch(),但我無法確定。
# Here is a MCVE:
Contract <- list(xts(1:12,order.by=Sys.Date()-1:12),
xts(1:10,order.by=Sys.Date()-1:10),
xts(1:8,order.by=Sys.Date()-2:9))
Vol <- xts(matrix(0, 12,3, byrow = FALSE),order.by=Sys.Date()-1:12)
for (A in 1:12){for (B in 1:3){
Vol[A,B] <- Contract[[B]][index(Vol)[A]]
}}
Vol
任何幫助將不勝感激。 (另外,如果有人像對雙循環進行矢量化一樣聰明......)
我們可以檢查將生成的觀察強制轉換為數字是否會產生numeric(0)
。
這是一個使用sapply()
的解決方案, sapply()
從所有日期的向量中創建矩陣。
all.dates <- Sys.Date() - 1:12
MX <- t(sapply(seq_along(all.dates), function(x)
sapply(seq_along(Contract), function(y) {
obs <- Contract[[y]][all.dates[x]]
if (identical(as.numeric(obs), numeric(0))) # reporting NA
xts(NA, all.dates[x])
else
obs
}
)))
Vol <- xts(MX, order.by=all.dates)
rm(MX) # clean up
> Vol
[,1] [,2] [,3]
2019-01-09 12 NA NA
2019-01-10 11 NA NA
2019-01-11 10 10 NA
2019-01-12 9 9 8
2019-01-13 8 8 7
2019-01-14 7 7 6
2019-01-15 6 6 5
2019-01-16 5 5 4
2019-01-17 4 4 3
2019-01-18 3 3 2
2019-01-19 2 2 1
2019-01-20 1 1 NA
您能否檢查一下這是否是您期望的結果?
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.