以下是创建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
我旁边的编码伙伴让他的两个人都去工作。