CALL: FIT_POLYNOMIAL
The histogram is referenced via a global variable that has usually
been returned by a call to BOOK_HISTOGRAM or GET_HISTOGRAM. This
argument may also be an expression that results in a histogram.
The argument is not modified on return.
A 1-dimensional array that contains the ordinates of the
points to be fitted.
The vectors x, y and err_y must have the same length,
unless err_y is given as a scalar.
Not changed on exit.
A 1-dimensional array that contains the coordinates of the
points to be fitted.
The vectors x, y and err_y must have the same length,
unless err_y is given as a scalar.
Not changed on exit.
A 1-dimensional array that contains the errors on the
coordinates of the points to be fitted.
This argument can also be a scalar, which is then used as
error for all points. If this argument is a vector, then
it must have the same length as x and y.
Not changed on exit.
Terms in the polynomial expansion:
polynomial = a0 + a1*x + a2*x**2 + a3*x**2 + ... + an*x**n
The input values for the terms are not used in the fitting
procedure. Initial values are obtained from an unweighted
maximum likelihood fit.
This argument must be modifiable, i.e. it can not be a constant
nor an expression.
The errors returned are the square roots of the diagonal elements
of the covariance matrix, multiplied by a correction factor equal
to:
sqrt(chi**2/(#points - #parameters))
This correction factor is equal to 1 if the errors to the data
points have been assigned correctly.
These arguments must be modifiable, that is, they can neither be
constants nor expressions.
The options string may contain the keywords listed in the table
below. The keywords should be separated by blanks or commas.
Additional text may be added - no error message will be issued
if unrecognised keywords are used.
---------------------------------------------------------------------
Option Meaning Defaults
---------------------------------------------------------------------
EQUAL Assume equal weight of all points Not default
NOPLOT Don't plot the histogram and fit Default
NOPRINT Don't print details during the fit Default
PLOT Plot the histogram with fit Not default
POISSON Assume Poisson errors Default
PRINT Print details during the fit Not default
---------------------------------------------------------------------
Note: the EQUAL and POISSON options are ignored for matrix fits since
the errors on the points are explicitely provided in that case.
Keyword index.
Formatted on 10/11/98.