Synopsis
Load a set of templates and use it as a model component.
Syntax
load_template_model(modelname, templatefile, dstype=Data1D, sep=' ', comment='#', method=linear_interp, template_interpolator_name='default') modelname - str templatefile - str dstype - data class to use, optional sep - str, optional comment - str, optional method - func template_interpolator_name - str
Description
A template model can be considered to be an extension of the table model supported by `load_table_model` . In the template case, a set of models (the "templates") are read in and then compared to the data, with the best-fit being used to return a set of parameters.
Examples
Example 1
Load in the templates from the file "index.tmpl" as the model component "kerr", and set it as the source model for the default data set. The optimisation method is switched to use a grid search for the parameters of this model.
>>> load_template_model("kerr", "index.tmpl")
>>> set_source(kerr)
>>> set_method('gridsearch')
>>> set_method_opt('sequence', kerr.parvals)
>>> fit()Example 2
Fit a constant plus the templates, using the neville scheme for integrating the template onto the data grid. The Monte-Carlo based optimiser is used.
>>> load_template_model('tbl', 'table.lis',
... sherpa.utils.neville)
>>> set_source(tbl + const1d.bgnd)
>>> set_method('moncar')PARAMETERS
The parameters for this function are:
| Parameter | Definition |
|---|---|
| modelname | The identifier for this table model. |
| templatefile | The name of the file to read in. This file lists the template data files. |
| dstype | What type of data is to be used. Supported values include `Data1D` (the default) and `Data1DInt` . |
| sep | The separator character. The default is ' ' . |
| comment | The comment character. The default is '#' . |
| method | The interpolation method to use to map the input data onto the coordinate grid of the data set. Linear, nearest-neighbor, and polynomial schemes are provided in the sherpa.utils module. |
| template_interpolator_name | The method used to interpolate within the set of templates. The default is default . A value of none turns off the interpolation; in this case the grid-search optimiser must be used to fit the data. |
Notes
Examples of interpolation schemes provided by `sherpa.utils` are: `linear_interp` , `nearest_interp` , and `neville` .
The template index file is the argument to `load_template_model` , and is used to list the data files. It is an ASCII file with one line per template, and each line containing the model parameters (numeric values), followed by the MODELFLAG column and then the file name for the data file (its name must begin with FILE). The MODELFLAG column is used to indicate whether a file should be used or not; a value of 1 means that the file should be used, and a value of 0 that the line should be ignored. The parameter names are set by the column names.
The data file - the last column of the template index file - is read in and the first two columns used to set up the x and y values ( `Data1D` ) or xlo, xhi, and y values ( `Data1DInt` ). These files must be in ASCII format.
The method parameter determines how the template data values are interpolated onto the source data grid.
The template_interpolator_name parameter determines how the dependent axis (Y) values are interpolated when the parameter values are varied. This interpolation can be turned off by using a value of none , in which case the grid-search optimiser must be used. See `load_template_interpolator` for how to create a valid interpolator. The "default" interpolator uses `sherpa.models.KNNInterpolator` with k=2 and order=2.
Bugs
See the bugs pages on the Sherpa website for an up-to-date listing of known bugs.
See Also
- data
- dataspace1d, dataspace2d, datastack, fake, get_bkg_chisqr_plot, get_bkg_delchi_plot, get_bkg_fit_plot, get_bkg_model_plot, get_bkg_plot, get_bkg_ratio_plot, get_bkg_resid_plot, get_bkg_source_plot, load_arf, load_arrays, load_ascii, load_bkg, load_bkg_arf, load_bkg_rmf, load_data, load_grouping, load_image, load_multi_arfs, load_multi_rmfs, load_pha, load_quality, load_rmf, load_staterror, load_syserror, load_table, pack_image, pack_pha, pack_table, unpack_arf, unpack_arrays, unpack_ascii, unpack_bkg, unpack_data, unpack_image, unpack_pha, unpack_rmf, unpack_table
- filtering
- load_filter
- info
- get_default_id, list_bkg_ids, list_data_ids, list_model_ids, show_bkg_model, show_bkg_source
- modeling
- add_model, add_user_pars, clean, create_model_component, delete_bkg_model, delete_model, delete_model_component, get_model, get_model_autoassign_func, get_model_component, get_model_component_image, get_model_component_plot, get_model_plot, get_num_par, get_num_par_frozen, get_num_par_thawed, get_order_plot, get_par, get_pileup_model, get_source, get_source_component_image, get_source_component_plot, get_source_contour, get_source_image, get_source_plot, get_xsabund, get_xscosmo, get_xsxsect, get_xsxset, image_model, image_model_component, image_source, image_source_component, integrate, link, load_table_model, load_template_interpolator, load_user_model, normal_sample, reset, save_model, save_source, set_bkg_model, set_bkg_source, set_full_model, set_model, set_model_autoassign_func, set_pileup_model, set_source, set_xsabund, set_xscosmo, set_xsxsect, set_xsxset, t_sample, uniform_sample
- plotting
- get_cdf_plot, get_pdf_plot, get_pvalue_plot, get_pvalue_results, plot_cdf, plot_model, plot_model_component, plot_pdf, plot_pvalue, plot_scatter, plot_source, plot_source_component, plot_trace
- psfs
- delete_psf, load_conv
- saving
- save_arrays, save_data, save_delchi, save_error, save_filter, save_grouping, save_image, save_pha, save_quality, save_resid, save_staterror, save_syserror, save_table
- statistics
- load_user_stat
- utilities
- calc_chisqr, calc_energy_flux, calc_model_sum, calc_photon_flux, calc_source_sum, calc_stat, eqwidth
- visualization
- contour_model, contour_ratio, contour_resid