有沒有匹配勛伯格算法的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)

會做的工作。