Vinicius,
Vinicius Magalhães wrote > > > + n <- length(*x.ts*) > + for (p in 0:maxord[1]) for (d in 0:maxord[2]) for (q in 0:maxord[3]) > + for (P in 0:maxord[4]) for (D in 0:maxord[5]) for (Q in 0:maxord[6]) > + { > + 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. Best regards. ----- Victor Delgado cedeplar.ufmg.br P.H.D. student www.fjp.mg.gov.br reseacher -- View this message in context: http://r.789695.n4.nabble.com/Error-in-solve-default-res-hessian-n-used-Lapack-routine-dgesv-system-is-exactly-singular-tp4428273p4435054.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.