R : Copyright 2003, The R Development Core Team Version 1.7.1 (2003-06-16) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative project with many contributors. Type `contributors()' for more information. Type `demo()' for some demos, `help()' for on-line help, or `help.start()' for a HTML browser interface to help. Type `q()' to quit R. [Previously saved workspace restored] > deaths <- ts(data= scan(file="deaths.tsm"), + start=c(1973, 1), frequency=12) Read 72 items > plot(deaths) > help(arima) > arima(deaths, order = c(0,1,1), seasonal = list(order=c(0,1,1))) Call: arima(x = deaths, order = c(0, 1, 1), seasonal = list(order = c(0, 1, 1))) Coefficients: ma1 sma1 -0.4264 -0.5584 s.e. 0.1226 0.1787 sigma^2 estimated as 99481: log likelihood = -425.53, aic = 857.06