Really good code, thank you John.
just one question, I want to have a constrained spline between a max and min values but with 'maxvalue' and 'minvalue' you say (and I try it and it happens): "This constraint is only a necessary constraint. It is not sufficient. In some circumstances the spline may pass slightly above this maximum value". Is there a way to avoid that?
And why is not exactly? is it a fomulation problem or a numerical? because if is numerical i think i can work with that.
Thanks so much for the code. For my data works fine. Some questions though :
1. How can I retrieve the standard error (a posteriori sigma zero ) for the overall result?
2. How can I retrieve the standard errors (or Covariance matrix) for each coefficient
I need these std errors for evaluating the results.