Parameter Errors from lsqcurvefit

7 次查看(过去 30 天)
Unkle77
Unkle77 2015-6-18
评论: Unkle77 2015-6-19
Hi Everybody,
I've been using the lsqcurvefit function to fit a Weibull function to a small set of data (4 to 10 data points) I've been using lsqcurvefit as it allows me to constrain the parameters as otherwise the fits become unrealistic with near right-angle turns at a single point. However I need to calculate the errors of the Weibull parameters, especially the error in ymax as this is needed for another graph. I've tried the code below which has been suggested in other threads:
[Q,R] = qr(J,0);
mse = sum(abs(r).^2)/(size(J,1)-size(J,2));
Rinv = inv®;
Sigma = Rinv*Rinv'*mse;
se = sqrt(diag(Sigma));
but this has not produced any numbers for my data.
any suggestions?

回答(1 个)

Torsten
Torsten 2015-6-19
I guess that the formulae above work for unconstrained fitting, but are no longer valid if you constrain your parameters.
By the way: "inv" is a "registered trademark" ?
Best wishes
Torsten.
  1 个评论
Unkle77
Unkle77 2015-6-19
I don't understand why a constrained fitting should result in it not working, surely it would just result in a larger error?
Yeah, no idea why (R) became a registered trade mark character, just for clarity it should be ( R ).

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Linear Least Squares 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by