IdealG¶

class dcprogs.likelihood.IdealG(*args)[source]¶

Ideal Likelihood.

This object can be instantiated one of several way:

With a QMatrix

>>> matrix = QMatrix(array([...]), 2)
>>> idealg = IdealG(matrix)

af(self, t) → DCProgs::t_rmatrix[source]¶

AF transitions with respect to time

Implements the ideal likelihood:

\[e^{t\mathcal{Q}_{FF}}\mathcal{Q}_{FA}.\]

Parameters:	t – A scalar or something to a numpy array. In the latter case, the return is an array of matrices.

fa(self, t) → DCProgs::t_rmatrix[source]¶

FA transitions with respect to time

Implements the ideal likelihood:

\[e^{t\mathcal{Q}_{AA}}\mathcal{Q}_{AF}.\]

Parameters:	t – A scalar or something to a numpy array. In the latter case, the return is an array of matrices.

final_occupancies¶

Equilibrium occupancies for final states.

Computes the right eigenvector of \(\mathcal{G}_{FA}\mathcal{G}_{AF}\), where \(\mathcal{G}_{FA}\) is the laplacian for \(s=0\) of the likelihood.

initial_occupancies¶

Equilibrium occupancies for initial states.

Computes the left eigenvector of \(\mathcal{G}_{AF}\mathcal{G}_{FA}\), where \(\mathcal{G}_{AF}\) is the laplacian for \(s=0\) of the likelihood.

laplace_af(self, s) → DCProgs::t_rmatrix[source]¶

AF transitions with respect to scale

Implements the laplace transform of the likelihood:

\[(sI - \mathcal{Q}_{AA})^{-1}\mathcal{Q}_{AF}.\]

Parameters:	s – A scalar or something to a numpy array. In the latter case, the return is an array of matrices.

laplace_fa(self, s) → DCProgs::t_rmatrix[source]¶

FA transitions with respect to scale

Implements the laplace transform of the likelihood:

\[(sI - \mathcal{Q}_{FF})^{-1}\mathcal{Q}_{FA}.\]

Parameters:	s – A scalar or something to a numpy array. In the latter case, the return is an array of matrices.