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

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

设置循环范围时，需要将冒号设置为开头。因此length(ozone)应为1:length(ozone)。否则，您只迭代最后一个臭氧项目。
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)
}

for循环不迭代向量中的所有值

1 个答案: