SplinePCFit function

Smoothing spline Anova method