MissedEventsG¶
Implementation of recursion for exact missed-event G function.
Implements the exact-missed event probability calculations, as detailed in Hawkes, Jalali, and Colquhoun (1990). Specifically, this is equation 3.2.
Public FunctionsMissedEventsG(DeterminantEq const & _af, std::vector< Root > const & _roots_af, DeterminantEq const & _fa, std::vector< Root > const & _roots_fa, t_uint _nmax = 3)Initializes missed events G functor.
- Parameters
- _af -
Determinant equation for af.
- _fa -
Determinant equation for af.
- _roots_af -
Roots of determinant equation for af
- _roots_fa -
Roots of determinant equation for fa
- _nmax -
Switches to asymptotic values after t≥nmaxτ
MissedEventsG(QMatrix const & _qmatrix, t_real _tau, t_RootFinder const & _findroots, t_uint _nmax = 3)Initializes missed events functor.
Uses input root finding function to determine roots.
- Parameters
- _qmatrix -
Transition matrix
- _tau -
resolution/max length missed events
- _findroots -
A functor with which to find all roots. This function should take a DeterminantEq as its sole argument and return a std::vector<DCProgs::RootInterval>
- _nmax -
Switches to asymptotic values after t≥nmaxτ
Initializes missed-events functor.
- Parameters
- _qmatrix -
Transition matrix
- _tau -
resolution/max length missed events
- _nmax -
Switches to asymptotic values after t≥nmaxτ
- _xtol -
Tolerance for interval size
- _rtol -
Tolerance for interval size. The convergence criteria is an affine function of the root: xtol+rtolxcurrent=12|xa−xb|.
- _itermax -
maximum number of iterations for any of the three steps.
- _lowerbound -
Lower bound of the interval bracketing all roots. If None, the lower bound is obtained from find_lower_bound_for_roots().
- _upperbound -
Upper bound of the interval bracketing all roots. If None, the upper bound is obtained from find_upper_bound_for_roots().
MissedEventsG(MissedEventsG && _c)Move constructor.
Open to close transitions.
Close to open transitions.
Probability of no shut times detected between 0 and t.
Probability of no open times detected between 0 and t.
void set_nmax(t_uint _n)Sets t≥nmaxτ.
t_uint get_nmax()When to switch to asymptotic values.
t_real get_tau()Gets the value of missed event resolution;.
t_real get_tmax()tmax is the time after which approximate calculations are performed.
t_rmatrix const & get_af_factor()QAFe−QFFτ
t_rmatrix const & get_fa_factor()QFAe−QAAτ
Exact laplace of AF.
Exact laplace of FA.
QMatrix const & get_qmatrix()Returns current QMatrix.