我有以下数据集(代码需要tslm
调用的forecast
包。
x <- rnorm(11, mean = 534, sd = 79)
y <- rnorm(9, mean = 800, sd = 56)
p <- list(x, y)
tsl <- list(); ts_trend <- list()
for(i in seq_along(p)) {
tsl[[i]] <- ts(p[[i]], start = c(2018, 1), frequency = 52)
}
for(i in seq_along(tsl)) {
ts_trend[[i]] <- tslm(tsl[[i]] ~ trend)
}
当我运行它时,我收到错误
Error in tsl[[i]] : subscript out of bounds
据我所知,下标显然不是越界的。 我在前面的循环中使用相同的引用,没有错误。
我不知道如何解决这个问题。 我错过了什么?
我们可以使用lapply
它会起作用
ts_trendN <- lapply(tsl, function(x) tslm(x ~ trend))
循环不工作的原因for
基于计算trend
的环境。 我们可以创建一个新环境,它会正常工作
for(i in seq_along(tsl)) {
ev1 <- new.env()
ev1$tsl1 <- tsl[[i]]
ts_trend[[i]] <- tslm(ev1$tsl1 ~ trend)
}
属性可能存在一些差异。 模型输出相同
library(broom)
identical(tidy(ts_trendN[[1]]), tidy(ts_trend[[1]]))
#[1] TRUE