daily.fit <- arima(sd_l_daily_adj$Price,
                   order=c(8,0,0),
                   seasonal=list(order=c(5,0,1), period=7),
                   xreg = sd_l_daily_adj$Holiday,
                   include.mean=FALSE)