How to find uncertainties of estimated parameters in Levenberg- Marquardt algorithm

Hello I have used Levenberg-Marquardt algorithm in 3 parameter problem but I unable to estimate uncertainties associated with it, Please help me in this regard

[Chris Martin]

I need the formula

Answer 1

[Torsten]

If you have the statistics toolbox, use nlinfit together with nlparci and nlpredci.

Best wishes

Torsten.

[Chris Martin]

No I have used my own matlab code

[Torsten]

Then you will have to consult the literature cited for nlparci and nlpredci on how confidence intervals for the parameters are to be computed.

See e.g.

http://onlinelibrary.wiley.com/book/10.1002/0471725315

Best wishes

Torsten.

[Chris Martin]

There is nothing in that link

[Torsten]

Just wanted to give a link to this book:

http://www.amazon.com/Nonlinear-Regression-George-F-Seber/dp/0471471356/ref=sr_1_2?ie=UTF8&qid=1446721830&sr=8-2&keywords=Seber+wild

Best wishes

Torsten.

Answer 2

[Star Strider]

Finding the covariance matrix of a nonlinear regression requires that you calculate the Jacobian matrix at the convergence (the last parameter estimates calculated). In the context of the Wikipedia notation, the partial derivatives are with respect to each parameter at each data point, so the partial derivative of f(1) with respect to x(1) refers to the value of the partial derivative of the function at the first value of the independent variable with respect to the first parameter. It continues row-wise down the values of the independent variable, and column-wise across the parameters.

In a nonlinear regression, some or all of the partial derivatives of the function with respect to each parameter are by definition functions of themselves or of other parameters, so expect that in your result. (This is known as the model being ‘nonlinear in the parameters’.) This is easier if you have the Symbolic Math Toolbox, since it can calculate the symbolic Jacobian for you.

After that, you can plug your evaluated Jacobian matrix into the design matrix for a linear problem to calculate the covariance matrix, as explained in Least squares, regression analysis, and statistics. You can calculate the parameter confidence intervals from the diagonals of the covariance matrix and the t-distribution.

[Chris Martin]

Give me the formula.......

[Star Strider]

No. This is your project. You do the maths and write the code. If you have problems getting it to run correctly, I will help you.

[Chris Martin]

I got it thank you very much.......

[Star Strider]

My pleasure.