如何用Matlab计算lsqcurvefit的95％置信区间？

Question

由于Matlab中有一些固定参数的问题，我不得不从std切换。对fit的{​​{1}}命令。

对于普通lsqcurvefit命令，其中一个输出参数为fit，我可以从中计算每个参数的+/-和r ^ 2值。

gof也应该可以。但我不认为它是输出参数之一。

或者用我的问题换句话说：如何从lsqcurvefit计算一个拟合参数的+/-？ 有人可以帮我吗？

谢谢，Niko

Answer 1

是的。获取lsqcurvefit的所有输出参数，并在nlparci中使用它们，如下所示：

[x,resnorm,residual,exitflag,output,lambda,jacobian] =...
                               lsqcurvefit(@myfun,x0,xdata,ydata);
conf = nlparci(x,residual,'jacobian',jacobian)

现在conf包含N x 2拟合参数的N矩阵。 conf的每一行给出相应参数的上限和下限95％置信区间。