我如何在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％置信區間。