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 2006

URL: http://cxc.harvard.edu/ciao3.4/guess.html
Hardcopy (PDF): A4 | Letter
AHELP for CIAO 3.4 guess Context: sherpa

Synopsis

Estimates initial parameter values and ranges, based on input data.

Syntax

sherpa> GUESS <arg> [# [ID]]

# specifies the dataset whose data are to be used to help estimated
parameter values. The ID modifier is used only for guesses involving
background data, and then if and only if the Sherpa state object
variable multiback is set to 1, i.e., if more than one background
dataset is to be associated with a single source dataset. The ID
modifier may be any unreserved string (e.g., A, foo, etc.), i.e., a
string that is not a parsable command.

Description

<arg> may be:

GUESS Command Arguments

Argument To perform an estimate for:
MODELS All current model components.
<sherpa_model_stack> All current model components in the named stack (SOURCE, BACKGROUND, PILEUP, [B]NOISE, or KERNEL).
<sherpa_modelname> The specified model component.
<modelname> The specified model component.

In addition, <arg> may be either ON or OFF, to turn Sherpa's parameter-value-guessing mechanism on and off respectively.

Note that:

  • When a name is assigned to a model component, the model component is established, and if a dataset has been input then estimates are automatically made for the initial parameter values and their ranges (i.e., the GUESS command is run without the user having to explicitly issue the command, unless the user has previous issued the command GUESS OFF, in which case the parameters take on default values).
  • If multiple datasets have been input, and a model component has been established without using the SOURCE command, such as in

sherpa> GAUSS[g]

then parameter estimates will be made using the first dataset. To refine the estimates for, e.g., dataset 3, do the following:

sherpa> SOURCE 3 = g
sherpa> GUESS SOURCE 3
  • If multiple datasets have been input, and a model component has been established using the SOURCE command, such as in

sherpa> SOURCE 4 = GAUSS[g]

then parameter estimates will be made using the numbered dataset, in this example, dataset number 4.

  • If a dataset has been input, and an RSP model defined, but INSTRUMENT has not been set, such as in

sherpa> RSP[a](example2.arf,example2.rmf)
sherpa> SOURCE = GAUSS[g]

then parameter estimates will be made using the counts data (divided by the observation time, if provided in the PHA header); i.e., they will not take into account the area information provided in the ARF. For good parameter estimates, always set the instrument!

sherpa> RSP[a](example2.arf,example2.rmf)
sherpa> INSTRUMENT = a
sherpa> SOURCE = GAUSS[g]
  • If a dataset has been input, and an RSP model defined with only an RMF (and no ARF), and the INSTRUMENT has been set, then the parameter estimates will be suboptimal, as they will not utilize effective area information that may have been included.
  • The guessing algorithm may not be optimal for complex model expressions.

See the CREATE command for more information and other examples.

Example 1

Estimate initial parameter values and ranges, based on input dataset number 1:

sherpa> PARAMPROMPT OFF
sherpa> POLY[modela]
sherpa> DATA example1.dat
sherpa> GUESS modela

This command estimates parameter values and ranges for component modela, based on dataset 1.

Example 2

Estimate initial parameter values and ranges, based on a specified input dataset:

sherpa> DATA 2 example2.dat
sherpa> GUESS MODELS 2

This command estimates parameter values and ranges for all current model components, based on dataset 2.

Example 3

Estimate initial parameter values and ranges for source models and their corresponding datasets, when multiple source models have been defined and multiple datasets have been input:

sherpa> ERASE ALL
sherpa> READ DATA 1 example1.dat
sherpa> READ DATA 2 example2.dat
sherpa> READ DATA 3 example3.dat
sherpa> POW[modelc]
sherpa> GAUSS[modelf]
sherpa> BBODY[modelh]
sherpa> SOURCE 1 = modelc
sherpa> SOURCE 2 = modelf
sherpa> SOURCE 3 = modelh
sherpa> GUESS SOURCE 1

This command estimates initial parameter values and ranges, based on dataset number 1, for the components of the source model for dataset number 1 (i.e., modelc). Note, however, that these estimates were already made automatically when component modelc was first established with the command POW[modelc].

sherpa> GUESS SOURCE 2

Estimates initial parameter values and ranges, based on dataset number 2, for the components of the source model for dataset number 2 (i.e., modelf). Note that until this command is issued, the initial parameter values for modelf are estimates based on dataset number 1.

sherpa> GUESS SOURCE 3

Estimates initial parameter values and ranges, based on dataset number 3, for the components of the source model for dataset number 3 (i.e., modelh). Note that until this command is issued, the initial parameter values for modelh are estimates based on dataset number 1.

Bugs

See the Sherpa bug pages online for an up-to-date listing of known bugs.

Hardcopy (PDF): A4 | Letter
Last modified: December 2006



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.