Class N_state_model

Class containing functions for the N-state model.

Instance Methods

__init__(self)
Initialise the class by placing API_common methods into the API.

_CoM(self, pivot_point=None, centre=None)
Centre of mass analysis.

_cone_pdb(self, cone_type=None, scale=1.0, file=None, dir=None, force=False)
Create a PDB file containing a geometric object representing the various cone models.

_minimise_bc_data(self, model)
Extract and unpack the back calculated data.

numpy rank-3 array, numpy rank-1 array. _minimise_setup_atomic_pos(self, sim_index=None)
Set up the atomic position data structures for optimisation using PCSs and PREs as base data sets.

tuple of (list, numpy rank-1 array, numpy rank-1 array, numpy rank-1 array) _minimise_setup_tensors(self, sim_index=None)
Set up the data structures for optimisation using alignment tensors as base data sets.

numpy rank-1 array. _minimise_setup_fixed_tensors(self)
Set up the data structures for the fixed alignment tensors.

_target_fn_setup(self, sim_index=None, scaling=True)
Initialise the target function for optimisation or direct calculation.

list of [SpinContainer instance, str, str] or [InteratomContainer instance, str, str] base_data_loop(self)
Loop over the base data of the spins - RDCs, PCSs, and NOESY data.

calculate(self, spin_id=None, verbosity=1, sim_index=None)
Calculation function.

list of floats create_mc_data(self, data_id=None)
Create the Monte Carlo data by back-calculation.

float default_value(self, param)
The default N-state model parameter values.

grid_search(self, lower=None, upper=None, inc=None, constraints=False, verbosity=0, sim_index=None)
The grid search function.

bool is_spin_param(self, name)
Determine whether the given parameter is spin specific.

list of float map_bounds(self, param, spin_id=None)
Create bounds for the OpenDX mapping function.

minimise(self, min_algor=None, min_options=None, func_tol=None, grad_tol=None, max_iterations=None, constraints=False, scaling=True, verbosity=0, sim_index=None, lower=None, upper=None, inc=None)
Minimisation function.

tuple of (int, int, float) model_statistics(self, model_info=None, spin_id=None, global_stats=None)
Return the k, n, and chi2 model statistics.

list of (float or None) return_data(self, data_id)
Return the base data for the given data ID.

str return_data_name(self, param)
Return a unique identifying string for the N-state model parameter.

list of float return_error(self, data_id=None)
Create and return the spin specific Monte Carlo Ri error structure.

str return_grace_string(self, param)
Return the Grace string representation of the parameter.

str return_units(self, param)
Return a string representing the parameters units.

set_error(self, model_info, index, error)
Set the parameter errors.

set_param_values(self, param=None, value=None, spin_id=None, error=False, force=True)
Set the N-state model parameter values.

sim_init_values(self)
Initialise the Monte Carlo parameter values.

sim_pack_data(self, data_id, sim_data)
Pack the Monte Carlo simulation data.

list of float sim_return_param(self, model_info, index)
Return the array of simulation parameter values.

Inherited from api_base.API_base: back_calc_ri, bmrb_read, bmrb_write, constraint_algorithm, data_init, data_names, data_type, deselect, duplicate_data, eliminate, get_param_names, get_param_values, has_errors, model_desc, model_loop, model_type, molmol_macro, num_instances, overfit_deselect, pymol_macro, read_columnar_results, return_conversion_factor, return_data_desc, return_value, set_selected_sim, set_update, sim_return_chi2, sim_return_selected, skip_function, test_grid_ops

Inherited from object: __delattr__, __format__, __getattribute__, __hash__, __new__, __reduce__, __reduce_ex__, __repr__, __setattr__, __sizeof__, __str__, __subclasshook__

Inherited from api_common.API_common (private): _base_data_loop_spin, _create_mc_relax_data, _data_init_dummy, _data_init_spin, _eliminate_false, _has_errors_spin, _is_spin_param_true, _model_loop_single_global, _model_loop_spin, _model_type_global, _model_type_local, _num_instances_spin, _overfit_deselect_dummy, _return_data_relax_data, _return_error_relax_data, _return_no_conversion_factor, _return_value_general, _set_error_spin, _set_param_values_global, _set_param_values_spin, _set_selected_sim_global, _set_selected_sim_spin, _set_update, _sim_init_values_spin, _sim_pack_relax_data, _sim_return_chi2_spin, _sim_return_param_spin, _sim_return_selected_global, _sim_return_selected_spin, _test_grid_ops_general

Class Variables

default_value_doc = Desc_container("N-state model default valu...

return_data_name_doc = Desc_container("N-state model data type...

_table = uf_tables.add_table(label= "table: N-state data type ...

set_doc = Desc_container("N-state model set details")

Inherited from api_base.API_base: eliminate_doc, write_doc

Properties

Inherited from object: __class__

Method Details

__init__(self)
(Constructor)

Initialise the class by placing API_common methods into the API.

Overrides: object.__init__

_CoM(self, pivot_point=None, centre=None)

Centre of mass analysis.

This function does an analysis of the centre of mass (CoM) of the N states. This includes calculating the order parameter associated with the pivot-CoM vector, and the associated cone of motions. The pivot_point argument must be supplied. If centre is None, then the CoM will be calculated from the selected parts of the loaded structure. Otherwise it will be set to the centre arg.

Parameters:

• pivot_point (list of float of length 3) - The pivot point in the structural file(s).
• centre (list of float of length 3) - The optional centre of mass vector.

_cone_pdb(self, cone_type=None, scale=1.0, file=None, dir=None, force=False)

Create a PDB file containing a geometric object representing the various cone models.

Currently the only cone types supported are 'diff in cone' and 'diff on cone'.

Parameters:

• cone_type (str) - The type of cone model to represent.
• scale (float) - The size of the geometric object is eqaul to the average pivot-CoM vector length multiplied by this scaling factor.
• file (str) - The name of the PDB file to create.
• dir (str) - The name of the directory to place the PDB file into.
• force (int) - Flag which if set to True will cause any pre-existing file to be overwritten.

_minimise_bc_data(self, model)

Extract and unpack the back calculated data.

Parameters:

• model (class instance) - The instantiated class containing the target function.

_minimise_setup_atomic_pos(self, sim_index=None)

Set up the atomic position data structures for optimisation using PCSs and PREs as base data sets.

Parameters:

• sim_index (None or int) - The index of the simulation to optimise. This should be None if normal optimisation is desired.

Returns: numpy rank-3 array, numpy rank-1 array.

The atomic positions (the first index is the spins, the second is the structures, and the third is the atomic coordinates) and the paramagnetic centre.

_minimise_setup_tensors(self, sim_index=None)

Set up the data structures for optimisation using alignment tensors as base data sets.

Parameters:

• sim_index (None or int) - The index of the simulation to optimise. This should be None if normal optimisation is desired.

Returns: tuple of (list, numpy rank-1 array, numpy rank-1 array, numpy rank-1 array)

The assembled data structures for using alignment tensors as the base data for optimisation. These include:

• full_tensors, the data of the full alignment tensors.
• red_tensor_elem, the tensors as concatenated rank-1 5D arrays.
• red_tensor_err, the tensor errors as concatenated rank-1 5D arrays.
• full_in_ref_frame, flags specifying if the tensor in the reference frame is the full or reduced tensor.

_minimise_setup_fixed_tensors(self)

Set up the data structures for the fixed alignment tensors.

Returns: numpy rank-1 array.

The assembled data structures for the fixed alignment tensors.

_target_fn_setup(self, sim_index=None, scaling=True)

Initialise the target function for optimisation or direct calculation.

Parameters:

• sim_index (None or int) - The index of the simulation to optimise. This should be None if normal optimisation is desired.
• scaling (bool) - If True, diagonal scaling is enabled during optimisation to allow the problem to be better conditioned.

base_data_loop(self)

Loop over the base data of the spins - RDCs, PCSs, and NOESY data.

This loop iterates for each data point (RDC, PCS, NOESY) for each spin or interatomic data container, returning the identification information.

Returns: list of [SpinContainer instance, str, str] or [InteratomContainer instance, str, str]

A list of the spin or interatomic data container, the data type ('rdc', 'pcs', 'noesy'), and the alignment ID if required.

Overrides: api_base.API_base.base_data_loop

calculate(self, spin_id=None, verbosity=1, sim_index=None)

Calculation function.

Currently this function simply calculates the NOESY flat-bottom quadratic energy potential, if NOE restraints are available.

Parameters:

• spin_id (None or str) - The spin identification string (unused).
• verbosity (int) - The amount of information to print. The higher the value, the greater the verbosity.
• sim_index (None) - The MC simulation index (unused).

Overrides: api_base.API_base.calculate

create_mc_data(self, data_id=None)

Create the Monte Carlo data by back-calculation.

Parameters:

• data_id (str) - The list of spin ID, data type, and alignment ID, as yielded by the base_data_loop() generator method.

Returns: list of floats

The Monte Carlo Ri data.

Overrides: api_base.API_base.create_mc_data

default_value(self, param)

The default N-state model parameter values.

Parameters:

• param (str) - The N-state model parameter.

Returns: float

The default value.

Overrides: api_base.API_base.default_value

grid_search(self, lower=None, upper=None, inc=None, constraints=False, verbosity=0, sim_index=None)

The grid search function.

Parameters:

• lower (array of numbers) - The lower bounds of the grid search which must be equal to the number of parameters in the model.
• upper (array of numbers) - The upper bounds of the grid search which must be equal to the number of parameters in the model.
• inc (array of int) - The increments for each dimension of the space for the grid search. The number of elements in the array must equal to the number of parameters in the model.
• constraints (bool) - If True, constraints are applied during the grid search (elinating parts of the grid). If False, no constraints are used.
• verbosity (int) - A flag specifying the amount of information to print. The higher the value, the greater the verbosity.

Overrides: api_base.API_base.grid_search

is_spin_param(self, name)

Determine whether the given parameter is spin specific.

Parameters:

• name (str) - The name of the parameter.

Returns: bool

False

Overrides: api_base.API_base.is_spin_param

map_bounds(self, param, spin_id=None)

Create bounds for the OpenDX mapping function.

Parameters:

• param (str) - The name of the parameter to return the lower and upper bounds of.
• spin_id (None) - The spin identification string (unused).

Returns: list of float

The upper and lower bounds of the parameter.

Overrides: api_base.API_base.map_bounds

minimise(self, min_algor=None, min_options=None, func_tol=None, grad_tol=None, max_iterations=None, constraints=False, scaling=True, verbosity=0, sim_index=None, lower=None, upper=None, inc=None)

Minimisation function.

Parameters:

• min_algor (str) - The minimisation algorithm to use.
• min_options (array of str) - An array of options to be used by the minimisation algorithm.
• func_tol (None or float) - The function tolerance which, when reached, terminates optimisation. Setting this to None turns of the check.
• grad_tol (None or float) - The gradient tolerance which, when reached, terminates optimisation. Setting this to None turns of the check.
• max_iterations (int) - The maximum number of iterations for the algorithm.
• constraints (bool) - If True, constraints are used during optimisation.
• scaling (bool) - If True, diagonal scaling is enabled during optimisation to allow the problem to be better conditioned.
• verbosity (int) - A flag specifying the amount of information to print. The higher the value, the greater the verbosity.
• sim_index (None or int) - The index of the simulation to optimise. This should be None if normal optimisation is desired.
• lower (array of numbers) - The lower bounds of the grid search which must be equal to the number of parameters in the model. This optional argument is only used when doing a grid search.
• upper (array of numbers) - The upper bounds of the grid search which must be equal to the number of parameters in the model. This optional argument is only used when doing a grid search.
• inc (array of int) - The increments for each dimension of the space for the grid search. The number of elements in the array must equal to the number of parameters in the model. This argument is only used when doing a grid search.

Overrides: api_base.API_base.minimise

model_statistics(self, model_info=None, spin_id=None, global_stats=None)

Return the k, n, and chi2 model statistics.

k - number of parameters. n - number of data points. chi2 - the chi-squared value.

Parameters:

• model_info (None) - The data returned from model_loop() (unused).
• spin_id (None or str) - The spin identification string. This is ignored in the N-state model.
• global_stats (None or bool) - A parameter which determines if global or local statistics are returned. For the N-state model, this argument is ignored.

Returns: tuple of (int, int, float)

The optimisation statistics, in tuple format, of the number of parameters (k), the number of data points (n), and the chi-squared value (chi2).

Overrides: api_base.API_base.model_statistics

return_data(self, data_id)

Return the base data for the given data ID.

Parameters:

• data_id (list of str) - The list of spin ID, data type, and alignment ID, as yielded by the base_data_loop() generator method.

Returns: list of (float or None)

The base data.

Overrides: api_base.API_base.return_data

return_data_name(self, param)

Return a unique identifying string for the N-state model parameter.

Parameters:

• param (str) - The N-state model parameter.

Returns: str

The unique parameter identifying string.

Overrides: api_base.API_base.return_data_name

return_error(self, data_id=None)

Create and return the spin specific Monte Carlo Ri error structure.

Parameters:

• data_id (str) - The list of spin ID, data type, and alignment ID, as yielded by the base_data_loop() generator method.

Returns: list of float

The Monte Carlo simulation data errors.

Overrides: api_base.API_base.return_error

return_grace_string(self, param)

Return the Grace string representation of the parameter.

This is used for axis labelling.

Parameters:

• param (str) - The specific analysis parameter.

Returns: str

The Grace string representation of the parameter.

Overrides: api_base.API_base.return_grace_string

return_units(self, param)

Return a string representing the parameters units.

Parameters:

• param (str) - The name of the parameter to return the units string for.

Returns: str

The parameter units string.

Overrides: api_base.API_base.return_units

set_error(self, model_info, index, error)

Set the parameter errors.

Parameters:

• model_info (int) - The global model index originating from model_loop().
• index (int) - The index of the parameter to set the errors for.
• error (float) - The error value.

Overrides: api_base.API_base.set_error

set_param_values(self, param=None, value=None, spin_id=None, error=False, force=True)

Set the N-state model parameter values.

Parameters:

• param (list of str) - The parameter name list.
• value (list) - The parameter value list.
• spin_id (None) - The spin identification string (unused).
• error (bool) - A flag which if True will allow the parameter errors to be set instead of the values.
• force (bool) - A flag which if True will cause current values to be overwritten. If False, a RelaxError will raised if the parameter value is already set.

Overrides: api_base.API_base.set_param_values

sim_init_values(self)

Initialise the Monte Carlo parameter values.

Overrides: api_base.API_base.sim_init_values

sim_pack_data(self, data_id, sim_data)

Pack the Monte Carlo simulation data.

Parameters:

• data_id (list of str) - The list of spin ID, data type, and alignment ID, as yielded by the base_data_loop() generator method.
• sim_data (list of float) - The Monte Carlo simulation data.

Overrides: api_base.API_base.sim_pack_data

sim_return_param(self, model_info, index)

Return the array of simulation parameter values.

Parameters:

• model_info (int) - The global model index originating from model_loop().
• index (int) - The index of the parameter to return the array of values for.

Returns: list of float

The array of simulation parameter values.

Overrides: api_base.API_base.sim_return_param

Class Variable Details

default_value_doc

Value:

Desc_container("N-state model default values")

return_data_name_doc

Value:

Desc_container("N-state model data type string matching patterns")

_table

Value:

uf_tables.add_table(label= "table: N-state data type patterns", captio
n= "N-state model data type string matching patterns.")