我如何在lsqcurvefit中计算95％的置信区间？

Question

由于Matlab中固定参数的一些问题，我不得不从标准转换。 fit lsqcurvefit命令。

对于正常的fit命令，输出参数之一是gof ，从中我可以计算出每个参数的+/-和r ^ 2值。

对于lsqcurvefit也应该是可能的。 但是我没有把它作为输出参数之一。

或者换句话说：我如何从lsqcurvefit计算出fitparamter的+/-？ 有人可以帮我吗？

谢谢，尼可

Answer 1

是的。 获取lsqcurvefit所有输出参数，并像下面这样在nlparci使用它们：

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

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