pyemma.msm.ChapmanKolmogorovValidator

class pyemma.msm.ChapmanKolmogorovValidator(test_model, test_estimator, memberships, mlags=None, conf=0.95, err_est=False, n_jobs=None, show_progress=True)
__init__(test_model, test_estimator, memberships, mlags=None, conf=0.95, err_est=False, n_jobs=None, show_progress=True)

Initialize self. See help(type(self)) for accurate signature.

Methods

__init__(test_model, test_estimator, memberships) Initialize self.
estimate(X, **params) Estimates the model given the data X
fit(X[, y]) Estimates parameters - for compatibility with sklearn.
get_params([deep]) Get parameters for this estimator.
load(file_name[, model_name]) Loads a previously saved PyEMMA object from disk.
save(file_name[, model_name, overwrite, …]) saves the current state of this object to given file and name.
set_params(**params) Set the parameters of this estimator.

Attributes

estimates Returns estimates at different lagtimes
estimates_conf Returns the confidence intervals of the estimates at different lagtimes (if available).
lagtimes
logger The logger for this class instance
memberships
model The model estimated by this Estimator
n_jobs Returns number of jobs/threads to use during assignment of data.
name The name of this instance
predictions Returns tested model predictions at different lagtimes
predictions_conf Returns the confidence intervals of the estimates at different lagtimes (if available)
show_progress whether to show the progress of heavy calculations on this object.
test_estimator
test_model
estimate(X, **params)

Estimates the model given the data X

Parameters:
  • X (object) – A reference to the data from which the model will be estimated
  • params (dict) – New estimation parameter values. The parameters must that have been announced in the __init__ method of this estimator. The present settings will overwrite the settings of parameters given in the __init__ method, i.e. the parameter values after this call will be those that have been used for this estimation. Use this option if only one or a few parameters change with respect to the __init__ settings for this run, and if you don’t need to remember the original settings of these changed parameters.
Returns:

estimator – The estimated estimator with the model being available.

Return type:

object

estimates

Returns estimates at different lagtimes

Returns:Y – each row contains the n observables computed at one of the T lag t imes.
Return type:ndarray(T, n)
estimates_conf

Returns the confidence intervals of the estimates at different lagtimes (if available).

If not available, returns None.

Returns:
  • L (ndarray(T, n)) – each row contains the lower confidence bound of n observables computed at one of the T lag times.
  • R (ndarray(T, n)) – each row contains the upper confidence bound of n observables computed at one of the T lag times.
fit(X, y=None)

Estimates parameters - for compatibility with sklearn.

Parameters:X (object) – A reference to the data from which the model will be estimated
Returns:estimator – The estimator (self) with estimated model.
Return type:object
get_params(deep=True)

Get parameters for this estimator.

Parameters:deep (boolean, optional) – If True, will return the parameters for this estimator and contained subobjects that are estimators.
Returns:params – Parameter names mapped to their values.
Return type:mapping of string to any
classmethod load(file_name, model_name='default')

Loads a previously saved PyEMMA object from disk.

Parameters:
  • file_name (str or file like object (has to provide read method)) – The file like object tried to be read for a serialized object.
  • model_name (str, default='default') – if multiple models are contained in the file, these can be accessed by their name. Use pyemma.list_models() to get a representation of all stored models.
Returns:

obj

Return type:

the de-serialized object

logger

The logger for this class instance

model

The model estimated by this Estimator

n_jobs

Returns number of jobs/threads to use during assignment of data.

Returns:
  • If None it will return the setting of ‘PYEMMA_NJOBS’ or
  • ’SLURM_CPUS_ON_NODE’ environment variable. If none of these environment variables exist,
  • the number of processors /or cores is returned.

Notes

This setting will effectively be multiplied by the the number of threads used by NumPy for algorithms which use multiple processes. So take care if you choose this manually.

name

The name of this instance

predictions

Returns tested model predictions at different lagtimes

Returns:Y – each row contains the n observables predicted at one of the T lag times by the tested model.
Return type:ndarray(T, n)
predictions_conf

Returns the confidence intervals of the estimates at different lagtimes (if available)

If not available, returns None.

Returns:
  • L (ndarray(T, n)) – each row contains the lower confidence bound of n observables computed at one of the T lag times.
  • R (ndarray(T, n)) – each row contains the upper confidence bound of n observables computed at one of the T lag times.
save(file_name, model_name='default', overwrite=False, save_streaming_chain=False)

saves the current state of this object to given file and name.

Parameters:
  • file_name (str) – path to desired output file
  • model_name (str, default='default') – creates a group named ‘model_name’ in the given file, which will contain all of the data. If the name already exists, and overwrite is False (default) will raise a RuntimeError.
  • overwrite (bool, default=False) – Should overwrite existing model names?
  • save_streaming_chain (boolean, default=False) – if True, the data_producer(s) of this object will also be saved in the given file.

Examples

>>> import pyemma, numpy as np
>>> from pyemma.util.contexts import named_temporary_file
>>> m = pyemma.msm.MSM(P=np.array([[0.1, 0.9], [0.9, 0.1]]))
>>> with named_temporary_file() as file: # doctest: +SKIP
...    m.save(file, 'simple') # doctest: +SKIP
...    inst_restored = pyemma.load(file, 'simple') # doctest: +SKIP
>>> np.testing.assert_equal(m.P, inst_restored.P) # doctest: +SKIP
set_params(**params)

Set the parameters of this estimator. The method works on simple estimators as well as on nested objects (such as pipelines). The former have parameters of the form <component>__<parameter> so that it’s possible to update each component of a nested object. :returns: :rtype: self

show_progress

whether to show the progress of heavy calculations on this object.