For a PfPR time series, \(x\), fit
a spline function. To fit only a subset of the interpolation
points, define options$shock_ix
Usage
fit_shock(xds_obj, options = list())
For a PfPR time series, \(x\), fit
a spline function. To fit only a subset of the interpolation
points, define options$shock_ix
fit_shock(xds_obj, options = list())