I would like to model some data that exhibits strong hourly seasonality.
Additionally the data depends on the ambient temperature and if it is a holiday or workday.
@formula(y ~ 1 + T_amb + T_amb^2 + hour_of_day(t))
One option would be to encode the hour of day in a set of dummy variables or to use a Fourier approach.
Additionally I want to switch the hour of day model if it is a holiday or workday.
@formula(y ~ 1 + T_amb + T_amb^2 + hour_of_day_holiday(t) + hour_of_day_workday(t))
How would I do that?
Is there an easier way to include such a switch?