Processing math: 100%

MissedEventsG

class dcprogs.likelihood.MissedEventsG(*args, **kwargs)[source]

Computes missed-events likelihood.

Exact calculations take place for times smaller than nmaxτ. Asymptotic calculations take over for larger times.

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

Likelihood of an observed open time of length t

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

Factor accounting for minimum shut time

It is the likelihood QAFeQFFτ of a shut time of length τ.

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

Likelihood of a shut time of length t

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

Factor accounting for minimum open time

It is the likelihood QFAeQAAτ of an open time of length τ.

final_CHS_occupancies(*args)[source]

Computes final CHS occupancies.

final_occupancies

Equilibrium occupancies for final states.

Computes the right eigenvector of eGFAeGAF, where eGFA is the laplacian for s=0 of the likelihood.

initial_CHS_occupancies(*args)[source]

Computes initial CHS occupancies.

initial_occupancies

Equilibrium occupancies for initial states.

Computes the left eigenvector of eGAFeGFA, where eGAF is the laplacian for s=0 of the likelihood.

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

Exact missed-events G in Laplace space.

The exact expression is eGAF(s)=AR(s)esτQAFeQFFτ, with AR(s)=[sIQAAQAFτ0esteQFFttQFA]1.

Parameters:s – The laplace scale. A real scalar or something convertible to a numpy array.
Returns:A matrix if the input is scalar, an array of matrices otherwise, with the shape of the input.
laplace_fa(self, s) → DCProgs::t_rmatrix[source]

Exact missed-events G in Laplace space.

The exact expression is eGFA(s)=FR(s)esτQFAeQAAτ, with FR(s)=[sIQFFQFAτ0esteQAAttQAF]1.

Parameters:s – The laplace scale. A real scalar or something convertible to a numpy array.
Returns:A matrix if the input is scalar, an array of matrices otherwise, with the shape of the input.
nmax

Cut-off time of exact calculations in units of tau.

nopen

Number of open-states.

nshut

Number of shut-states.

tau

Resolution or maximum length of the missed events.

tmax

Cut-off time of exact calculations tmax=(nmax1)tau.

For practical reasons, the minimum observation time has alreadybeen removed here.