我有一个简单的问题:假设我想估计一个缺少数据的(1,1,1)ARIMA模型。我为此使用forecast
。示例:
require(forecast)
set.seed(1)
Arima(c(NA,NA,rnorm(50),NA,NA),order = c(1,1,1)) # seems to simply ignore NAs
现在,我不确定这是一个面向编程还是面向统计的问题,因此我将其发布在这里。我的问题是Arima
如何处理NA。在典型的回归模型中,您有不同的选择,例如在lm
中:
行为
指示当数据出现时应该发生什么的功能 包含NA。默认值由选项的na.action设置来设置, 如果未设置,则为na.fail。 “新出厂”默认值为 na.omit。另一个可能的值为NULL,不执行任何操作。值不包括 可能有用。
由于ARIMA模型也是线性模型,因此为什么没有na.action选项?在文档中,我找不到任何东西。它确实具有na.interp
,但是我不确定何时在函数中使用它。链接到文档:https://cran.r-project.org/web/packages/forecast/forecast.pdf
任何提示表示赞赏。