MATLAB Answers

0

Confidence intervals estimates for cubic smoothing spline (csaps)

Is there a built in function that will return the confidence bands for csaps (or the output from ppval)? Extra information that might be useful:
  1. I have data that is smoothed both with and without weights.
  2. I estimated the smoothing parameter using cross-validation using both k-fold and leave-one-out. Aside: why doesn't MATLAB do this automatically?
I found a paper, Smoothing Noisy Data with Spline Functions by M. F. Hutchinson and F. R. De Hoog that describes a method for estimating the confidence intervals, but I'd rather not code it myself if a script already exists.
Edit: The title was poorly worded and it doesn't seem I can change it now. I meant confidence bands, not confidence interval.

  0 Comments

Sign in to comment.

2 Answers

Answer by Matthew Taliaferro on 29 May 2017
 Accepted Answer

I wrote my own code and posted it here. Hopefully it's useful for someone.

  0 Comments

Sign in to comment.


Answer by David Ding
on 26 May 2017

Hi Mark,
Unfortunately the output of "csaps" function requires further processing to obtain the confidence bands. The following detailed explanations on the Curving Fitting Toolbox might be useful:
Thanks,
David

  1 Comment

Thank you for taking the time to answer my question. The doc page for confint() and predint() both say you cannot use them for calculating confidence or prediction intervals for non-parametric regression models. I wish it had been that easy :(
Again, thanks for the help.

Sign in to comment.