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

Fitting: How does gnuplot calculate the covariance matrix?

From	Ingo Thies <ingo.thies@gmx.de>
Newsgroups	comp.graphics.apps.gnuplot
Subject	Fitting: How does gnuplot calculate the covariance matrix?
Date	Fri, 08 Apr 2011 13:05:01 +0200
Lines	33
Message-ID	<9088euFi3iU1@mid.individual.net> (permalink)
Mime-Version	1.0
Content-Type	text/plain; charset=ISO-8859-1
Content-Transfer-Encoding	8bit
X-Trace	individual.net dSPCvd7zvK4o54HDp7nD1Qi9glz5He6pZRRmVKi6XkOjSxmgMR
Cancel-Lock	sha1:LdMeDASQb4cct1AH/FIh2mnnPQo=
User-Agent	Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.9.2.14) Gecko/20110223 Thunderbird/3.1.8
Path	csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!news.stben.net!border3.nntp.ams.giganews.com!border1.nntp.ams.giganews.com!nntp.giganews.com!de-l.enfer-du-nord.net!feeder1.enfer-du-nord.net!fu-berlin.de!uni-berlin.de!individual.net!not-for-mail
Xref	x330-a1.tempe.blueboxinc.net comp.graphics.apps.gnuplot:223

Show key headers only | View raw

Hi all,

I am trying to find out how exactly gnuplot calculates the covariance
matrix for a fitted function.

The background: I am still mistrusting the error ellipses one gets from
the eigenvalues. As already mentioned last year (aroung May or so), I
suspect that the resulting error contours are underestimated since the
best-fit chi^2 is assumed to be zero there, but it isn't.

I have compared the error ellipses with the true error contours from an
independent chi^2-2D-mapping routine, and found that the latter may be
considerably larger than the ones one would get from gnuplot's cov matrix.

However, basic information about covariance (e.g. from Wikipedia) are
misleading here, since they only deal with the covariance of measured
data themselves, not the errors. Also I didn't manage to get a clear
image from the sourcecode yet.

For the long term, it would be helpful to include the basic formulas and
the most important caveats and pitfalls (e.g. how may I trust the
results or not) into the gnuplot.pdf.

In our institute, there are people working with gnuplot and its fitting
algorihm, but lack of important information about how they are working
in detail.

Many thanks for any helpful reply!

-- 
Gruß,
      Ingo

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

Thread

Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-08 13:05 +0200
  Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-08 17:18 +0200
    Re: Fitting: How does gnuplot calculate the covariance matrix? Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-04-10 19:31 +0200
      Re: Fitting: How does gnuplot calculate the covariance matrix? Charles Allen <ca137tmp@earthlink.net> - 2011-04-10 20:14 -0500
        Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-14 15:41 +0200
      Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-11 21:42 +0200
        Re: Fitting: How does gnuplot calculate the covariance matrix? Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-04-11 23:27 +0200
          Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-12 11:28 +0200
            Re: Fitting: How does gnuplot calculate the covariance matrix? Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-04-16 00:03 +0200
              Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-16 18:06 +0200
      Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-14 15:07 +0200
      Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-15 13:47 +0200
  Re: Fitting: How does gnuplot calculate the covariance matrix? Hans-Bernhard Bröker <HBBroeker@t-online.de> - 2011-04-10 19:04 +0200
    Re: Fitting: How does gnuplot calculate the covariance matrix? Ingo Thies <ingo.thies@gmx.de> - 2011-04-11 22:13 +0200

csiph-web