Question

我想获取时间序列列表的百分比差异，但是当列表包含空list（）元素时，我遇到了错误。

例如，给定一个列表myList，如下所示：

set.seed(230)
zoo.Date <- as.Date("2003-02-01") + c(1, 3, 7, 9, 14) - 1
myList <- list(zoo(c(NA,rnorm(4)), zoo.Date), 
           zoo(c(rnorm(4), NA), zoo.Date), 
           list(),
           zoo(c(rnorm(2), NA, rnorm(2)), zoo.Date))

我可以将每个元素与myList %>% lapply(function(x) diff(x))

区分开来

但是一旦我尝试除以前一个元素myList %>% lapply(function(x) diff(x) / lag(x, -1))，我就会收到错误Error in attr(x, "tsp") <- c(1, NROW(x), 1) : cannot assign 'tsp' to zero-length vector。

问题是list（）元素。我可以摆脱它，但我想确保myList中的第四个对象对应于结果（滞后）列表的第四个对象。换句话说，我希望输出看起来像：

[[1]]
2003-02-03 2003-02-07 2003-02-09 2003-02-14 
        NA  -1.431227   2.245691   0.100597 

[[2]]
2003-02-03 2003-02-07 2003-02-09 2003-02-14 
-0.9078077 12.8523189 -1.2523403         NA 

[[3]]
NA

[[4]]
2003-02-03 2003-02-07 2003-02-09 2003-02-14 
 -1.177138         NA         NA  -2.203719

或[[3]]为list()或NULL或其他内容为空。

我尝试使用myList %>% unlist(recursive=F) ...和myList %>% lapply(function(x) ifelse(length(x) == 0, NA, x)) %>% lapply(function(x) diff(x) / lag(x, -1))之类的内容，但都没有效果。

Answer 1

解决此问题的一种方法是使用tryCatch来捕获空列表元素引发的错误并返回NA（或其他任何内容）：

lapply(myList, function(x) tryCatch(
  { diff(x) / lag(x, -1) }, 
  error=function(e){ NA }
  ))

# [[1]]
# 2003-02-03 2003-02-07 2003-02-09 2003-02-14 
#         NA  -1.431227   2.245691   0.100597 
# 
# [[2]]
# 2003-02-03 2003-02-07 2003-02-09 2003-02-14 
# -0.9078077 12.8523189 -1.2523403         NA 
# 
# [[3]]
# [1] NA
# 
# [[4]]
# 2003-02-03 2003-02-07 2003-02-09 2003-02-14 
#  -1.177138         NA         NA  -2.203719

时间序列列表与空元素的百分比差异

1 个答案: