Call: arima(x = StarSize.AllFrames, order = c(1, 0, 1), seasonal = list(order = c(1, 1, 1), period = 24)) Coefficients: ar1 ma1 sar1 sma1 0.9290 -0.7970 0.7530 -0.9713 s.e. 0.0093 0.0133 0.0141 0.0058 sigma^2 estimated as 753626609: log likelihood = -42693.82, aic = 85397.64