Re: Error in solve.default(res$hessian * n.used) :Lapack routine dgesv: system is exactly singular
Vinicius Magalhães wrote
+ n <- length(x.ts)
+ for (p in 0:maxord) for (d in 0:maxord) for (q in 0:maxord)
+ for (P in 0:maxord) for (D in 0:maxord) for (Q in 0:maxord)
+ fit <- arima(x.ts, order=c(p,d,q),
+ seas = list(order=c(P,D,Q),
+ frequency(x.ts)), method = "CSS")
+ fit.aic <- -2 * fit$loglik + (log(n) + 1) * length(fit$coef)
I think you are using 'x.ts' instead 's.ts'. Is it correct?!
Maybe you can also put here objects likely 'x', 'x.ts' to help us give feedback to you routine.
Professor in department of Economics,
UFOP - Univ. Federal de Ouro Preto, Brazil