Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.graphics.apps.gnuplot > #328

Re: Wrong fit direction / zrange / unreadable graph file

Show all headers | View raw

On 19.05.2011 11:07, Ingo Thies wrote:

> Hmm, I do not see why
>
> d/da (a+b*x)
>
> should cause troubles at a=0, unless there is some kind of scaling the
> differences (rather than differentials in numerics) da by a, so that da
> may become zero if a is.

... which is exactly what we do.  The derivative in question is usually 
computed as (f(a*1.001,b)-f(a,b))/(a*0.001).  As that would yield 
nonsense for a==zero, we replace all values in the region of 
-1e30<a<+1e-30 by +1e-30.

>> The second part of the problem is that we don't normalize parameters to
>> a common scale, so parameters of wildly different magnitude cause the
>> fit to grind to a halt.

> In the given example the number of iterations is pretty small, so that I
> won't assume any kind of numerical stall here.

Actually, in the problematic cases reported here, the reason the number 
of iterations was so small was that the fit progress was so slow that 
the algorithm mistook it for convergence (relative change of chisquare 
below the limit).  Note how 'n' got stuck at basically the first 
non-zero value it acquired (1e-30), and never changed again.  That's a 
typical symptom of parameter scale mismatch.

Back to comp.graphics.apps.gnuplot | Previous | Next — Previous in thread | Next in thread | Find similar

Thread

Wrong fit direction / zrange / unreadable graph file Katharina Miwege <k.miwege@googlemail.com> - 2011-05-15 23:59 -0700
  Re: Wrong fit direction / zrange / unreadable graph file Christoph Bersch <usenet@bersch.net> - 2011-05-16 11:04 +0200
    Re: Wrong fit direction / zrange / unreadable graph file Ingo Thies <ingo.thies@gmx.de> - 2011-05-16 11:34 +0200
      Re: Wrong fit direction / zrange / unreadable graph file Christoph Bersch <usenet@bersch.net> - 2011-05-16 11:52 +0200
        Re: Wrong fit direction / zrange / unreadable graph file Richard Owlett <rowlett@pcnetinc.com> - 2011-05-16 06:25 -0500
          Re: Wrong fit direction / zrange / unreadable graph file Christoph Bersch <usenet@bersch.net> - 2011-05-16 14:44 +0200
        Re: Wrong fit direction / zrange / unreadable graph file Ingo Thies <ingo.thies@gmx.de> - 2011-05-16 16:31 +0200
          Re: Wrong fit direction / zrange / unreadable graph file Katharina Miwege <k.miwege@googlemail.com> - 2011-05-16 09:01 -0700
            Re: Wrong fit direction / zrange / unreadable graph file Christoph Bersch <usenet@bersch.net> - 2011-05-17 08:46 +0200
              Re: Wrong fit direction / zrange / unreadable graph file Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-05-19 01:34 +0200
          Re: Wrong fit direction / zrange / unreadable graph file Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-05-19 01:49 +0200
            Re: Wrong fit direction / zrange / unreadable graph file Ingo Thies <ingo.thies@gmx.de> - 2011-05-19 11:07 +0200
              Re: Wrong fit direction / zrange / unreadable graph file Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-05-19 14:02 +0200
      Re: Wrong fit direction / zrange / unreadable graph file Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-05-19 01:39 +0200
        Re: Wrong fit direction / zrange / unreadable graph file Ingo Thies <ingo.thies@gmx.de> - 2011-05-19 10:58 +0200

csiph-web