有没有匹配勋伯格算法的R样条？

Question

是否有 R 算法可以在最小化 L 的同时拟合平滑样条，

L = ρ ∑ (i from 0 to n-1) wi(yi-Si(xi))² + (1 - ρ) ∫ (x from 0 to x_(n-1)) (S''(x)) ² dx

也许 smooth.spline 是可能的，但我没有成功找到好的参数。

（这里的方程可以看的更清楚： https://www.iro.umontreal.ca/~simardr/ssj/doc/html/umontreal/iro/lecuyer/functionfit/SmoothingCubicSpline.html ）

Answer 1

pspline::smooth.Pspline(  x = input_x,
                          y = input_y, 
                          norder = 2, 
                          method = 1,
                          spar = rho)

会做的工作。