SciPy

scipy.interpolate.LSQUnivariateSpline.get_residual

LSQUnivariateSpline.get_residual(self)[source]

Return weighted sum of squared residuals of the spline approximation.

This is equivalent to:

sum((w[i] * (y[i]-spl(x[i])))**2, axis=0)