Question

我试图通过使用我自己的函数“foo”来估计ARIMA模型。错误代码找不到R中的函数图无处可去。真的困扰我。知道怎么回事吗？我试过stats:::plot仍然无法正常工作。

```{r warning=FALSE, message=FALSE}
require(forecast)
foo <- function(N, y1=1, p=0.8, q=1-p, seed=NULL, stepwise=TRUE, verbose=0) {
if (!is.null(seed)) set.seed(seed)
y <- c(y1, rep(NA,N-1))
for (i in 2:N) y[i] <- p*y[i-1] + q*rnorm(1)
if (verbose) plot(y, type='l', xlab='time')
return (list(y=y, aa=auto.arima(y, stepwise=stepwise)))
}
v <- foo(N=100, seed=10, verbose=1); v$aa$coef
```

Answer 1

require(forecast)
Loading required package: forecast
> foo <- function(N, y1=1, p=0.8, q=1-p, seed=NULL, stepwise=TRUE, verbose=0) {
+     if (!is.null(seed)) set.seed(seed)
+     y <- c(y1, rep(NA,N-1))
+     for (i in 2:N) y[i] <- p*y[i-1] + q*rnorm(1)
+     if (verbose) plot(y, type='l', xlab='time')
+     return (list(y=y, aa=auto.arima(y, stepwise=stepwise)))
+ }
> v <- foo(N=100, seed=10, verbose=1); v$aa$coef
  ar1 
0.8946974

你的代码运作良好..它甚至绘制了一个情节。所以也许问题是加载预测包。或者尝试使用stats :: plot之前的情节不是三个冒号而是两个冒号。

错误：在R中找不到功能图

1 个答案: