For循环未遍历向量中的所有值

Question

我试图将多个时间序列的每年NA值相加。 我创建了一个仅包含日期的向量，然后将其与日期和臭氧值的向量进行比较。 我首先创建了一个for循环，该循环遍历我的日期向量的所有日期以仅提取年份。 然后，我正在将年向量与臭氧和日期值向量中的年进行比较（这是我的问题所在），然后我希望将每年放入向量中并求和NA值。 当我尝试遍历臭氧矢量的年份时，我要么创建了一个无限循环，要么仅提取了第一个变量。

我尝试了几种不同的方法来迭代日期和时间的向量，但是非常困难。 下面列出了我尝试过的其他内容以及我的代码：

totalNA <- c()
for (i in npsfiles[1]) {
  timeozone <- import(i, date="DATE", date.format = "%m/%d/%Y", header=TRUE, na.strings="-999")
  O3 <- c(timeozone$O3)
  tyr <- substr(timeozone$date, 1, 10)
  ozone <- O3
  names(ozone) <- tyr
  for (d in as.character(tyr)) {
    theyears <- format.Date(d, "%Y")
    print(theyears)
    for (y in length(ozone)) {
      if (names(ozone[y]) == theyears) {
      totalNA <- c(totalNA, ozone[y])}
      print(y)
      }
  }   
} 

我的问题在第三个for循环中。 我也尝试过：

for (y in seq(1:length(ozone))
for (y in ozone)
for (y in names(ozone))  

这是一些数据：

ABBR            DATE         O3
BADL-VC 10/01/1987 02   -999
BADL-VC 10/01/1987 03   -999
BADL-VC 01/01/1989 10   12
BADL-VC 01/01/1989 11   13
BADL-VC 01/01/1989 12   14

Answer 1

关于第三循环有两个问题：

1. 设置循环范围时，需要将开始设置为以冒号结尾。 所以length(ozone)应该是1:length(ozone) 。 否则，您仅迭代最后一个臭氧项。
2. names(ozone)返回完整日期而不是年份。 因此，您需要格式化这些元素，就像处理theyears 。

考虑对第三循环进行以下调整：

for (y in 1:length(ozone)) {
    ozoneyear <- format.Date(names(ozone[y]), "%Y")
    if (ozoneyear == theyears) {
      totalNA <- c(totalNA, ozone[y])
    }
    print(y)
}