Call: arima(x = Music_MariahCareyObsessed, order = c(1, 0, 1), seasonal = list(order = c(1, 1, 1), period = 24)) Coefficients: ar1 ma1 sar1 sma1 0.9488 -0.7279 0.1368 -0.9118 s.e. 0.0083 0.0180 0.0205 0.0089 sigma^2 estimated as 560077195: log likelihood = -35432.92, aic = 70875.84