help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: rcond = 0 when using leasqr. How to check result?


From: oxyopes
Subject: Re: rcond = 0 when using leasqr. How to check result?
Date: Thu, 23 Aug 2007 20:19:36 +0200

Oh, I think leasqr is well written, and i dont want to edit it.
I prefere to optionally define a third function in my script by adding
the exponential and the offset. Then I'll use the tips from the thread
''defining a function from 2 other ones'' in this list.
Thanks again ...

> > As you see, the term p(3) in the second is not used.
> That is probably the reason for the singularity. When 'leasqr'
> estimates the gradient (or Jacobian), nothing 'happens' for that
> parameter and the respective component(s) will be zero.
>
> As Sörens reply shows, the error (or warning) happens after converging
> (or not converging --- check [..., ..., _kvg_, ...] = leasqr (...) ),
> when additional information is computed, which you possibly do not
> need, the optimization is already completed.
>
> I would guess that the unused parameter does not change the result of
> optimization, once 'leasqr' succeeds. The reason for succeeding anyway
> is probably that in optimization the pseudoinverse is used, not the
> inverse. Still I would advise you not to pass p(3) as an optimization
> parameter; put it into a global variable or construct a suitable
> anonymous function.
>
> The problem may still occur sometimes, probably if 'leasqr' converges
> to regions with very flat gradients. Then, if the warning annoys you,
> you can simply edit 'leasqr' and check for singularity.



reply via email to

[Prev in Thread] Current Thread [Next in Thread]