Polynomial spline smoothing for nonlinear time series