keywords p-splines; mixed models; times series forecasting; varying-coefficient models; harmonic regression