我试图将多个时间序列的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 (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发布于 2016-07-20 03:15:52
关于第三个循环有两个问题:
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)
}https://stackoverflow.com/questions/38462294
复制相似问题