# Two seasonal periods in ARIMA using R

I'm currently using R to predict a time series with these instructions:

X <- ts(datas, frequency=24)
X.arima <- Arima(X, order=c(2,1,0), seasonal=c(1,1,1))
plot.ts(pred\$pred)


As you can see I've data each hour, and I chose the seasonal period of 24 (one day).

I would like to improve my forecasting using an additional seasonal period in order to include the seasonal component of the week (seasonal length of 7*24=168 data)

Is there any method for this? How do you do it?

UPDATE: I've read this (your) blog page, maybe can I use the external regressors to simulate a second seasonal period?

-
Yes, you could use some Fourier terms as regressors and deal with the seasonality that way. –  Rob Hyndman Jan 15 '13 at 22:53

There are no R packages that handle multiple seasonality for ARIMA models as far as I know. You could try the forecast package which implements multiple seasonality using models based on exponential smoothing. The dshw, bats and tbats functions will all handle data with two seasonal periods.

-

I am going through similar situation,Tyranitar. However, I found this paper - http://www.research.att.com/techdocs/TD_100381.pdf

The paper is about predicting the mobile network traffic prediction using the automatic double seasonal ARIMA. As it is a research paper, it has clearly described the algorithm that one can adopt to adopt multi-seasonal ARIMA prediction. So far, it has given me enough background to proceed further with my research. Please have a look, you can relate many formulation with your requirement.

-
Welcome to the site, @Neo182. Would you mind giving an overview of the paper so that users can decide it it's what they're looking for? It might also be nice to have an official citation in case of linkrot. Since you're new here, you should read our FAQ, which contains info about the site like this. –  gung May 8 '13 at 22:15