y <- msts(dat.day.ts$DailyEnergy, seasonal.periods=c(7,30.4,182.4),start=c(2007,1,1)) fit <- tbats(y) plot(fit, main="Multiple Season Decomposition"