sseq <- c(1, seq(5, 120, by = 5)) for(i in 1:length(sseq)){ assign(paste("arima", i, sep=""), arima0(data.ts[sseq[i]:(sseq[i]+115)], order=c(1,1,1))) } pred1 = predict(arima1, n.ahead = 5, se.fit = TRUE)$pred how do I traverse the arima models so I repeat the above prediction procedure(bold) on all arima models. Also, how do i automatically create one huge vector of the predictions from all the arima models? -- View this message in context: http://r.789695.n4.nabble.com/outputting-arima-models-tp2535695p2535695.html Sent from the R help mailing list archive at Nabble.com.