About Chandra Archive Proposer Instruments & Calibration Newsletters Data Analysis HelpDesk Calibration Database NASA Archives & Centers Chandra Science Links

Skip the navigation links
Last modified: December 2007

URL: http://cxc.harvard.edu/sherpabeta/projection.sl.html
AHELP for Sherpa 4.0 Beta projection Context: sl.sherpa

Synopsis

Estimate confidence intervals for selected thawed parameters

Syntax

proj(; [id=1], [[otherids], [model parameters]] );

Description

Confidence interval bounds are determined for each selected parameter in turn. A given parameter's value is varied along a grid of values while the values of all the other nominally thawed parameters are allowed to float to new best-fit values (compare to covar, where the values of all the other nominally thawed parameters remain fixed to their best-fit values). This method of estimating confidence interval bounds gives truly accurate results only in special cases:

An estimated confidence interval is accurate if and only if:

  • the chi^2 or logL surface in parameter space is approximately shaped like a multi-dimensional paraboloid, and
  • the best-fit point is sufficiently far from parameter space boundaries.

One may determine if these conditions hold, for example, by plotting the fit statistic as a function of each parameter's values (the curve should approximate a parabola) and by examining contour plots of the fit statistics made by varying the values of two parameters at a time (the contours should be elliptical, and parameter space boundaries should be no closer than approximately 3 sigma from the best-fit point).

If no arguments are given this function, it is assumed that the data id is the default data id, and that limits should be computed for all thawed parameters.

If arguments are given, each argument is examined to see if it is a Sherpa model parameter. If model parameters are given, it is assumed that limits for only those parameters should be computed. Any argument that is not a model parameter is assumed to be a data id.

Example

sherpa> proj();
sherpa> proj(1);
sherpa> proj(1,2,g.fwhm,g.ampl);

The proj function can be called with no arguments, or with arguments that are either data set IDs or model parameters. If no arguments are model parameters, then limits for all thawed parameters are computed. If any arguments are model parameters, then limits are computed for only those parameters.

Bugs

See the bugs pages on the Sherpa Beta website for an up-to-date listing of known bugs.

Last modified: December 2007



The Chandra X-Ray Center (CXC) is operated for NASA by the Smithsonian Astrophysical Observatory.
60 Garden Street, Cambridge, MA 02138 USA.    Email: cxcweb@head.cfa.harvard.edu
Smithsonian Institution, Copyright © 1998-2004. All rights reserved.