Call:
arima(x = Educational_IntroductionForJohnCambell, order = c(1, 0, 1), seasonal = list(order = c(1, 
    1, 1), period = 24))

Coefficients:
         ar1      ma1     sar1    sma1
      0.8794  -0.6982  -0.9583  0.6986
s.e.  0.0146   0.0202   0.0070  0.0221

sigma^2 estimated as 430863550:  log likelihood = -35021.9,  aic = 70053.8