MA(q)模型抛出错误

时间:2015-09-07 22:41:06

标签: r time-series

以下是创建MA(q)模型的代码,该模型包含值' thetas' sigsq'以及' T'。'。

masim <- function(thetas, sigsq, T)
{
  q <- length(thetas)

  noise <- rnorm(T+q, sd=sqrt(sigsq))
  x <- c(noise[1:q], rep(0,T))
  for( i in (q+1):(T+q)){
    x[i] <- thetas %*% x[i-(1:q)] + noise[i]
  }
x <- x[(q+1):(T+q)]
x
}

现在,如果我在masim上运行这样的自相关,例如:

acf(masim(0.5,1,10000))

它的工作正常,有很好的ACF情节。但是如果我运行这个例子:

acf(masim(2,1,10000))

然后我收到了这个错误。

Error in plot.window(...) : need finite 'ylim' values
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf

我旁边的编码伙伴让他的两个人都去工作。

0 个答案:

没有答案