isItStillActive: Function to determine if a feature is still active.

In HCBravoLab/metagenomeSeq: Statistical analysis for sparse high-throughput sequencing

Function to determine if a feature is still active.

Description

In the Expectation Maximization routine features posterior probabilities routinely converge based on a tolerance threshold. This function checks whether or not the feature's negative log-likelihood (measure of the fit) has changed or not.

Usage

isItStillActive(eps, tol, stillActive, stillActiveNLL, nll)

Arguments

eps	Vector of size M (features) representing the relative difference between the new nll and old nll.
tol	The threshold tolerance for the difference
stillActive	A vector of size M booleans saying if a feature is still active or not.
stillActiveNLL	A vector of size M recording the negative log-likelihoods of the various features, updated for those still active.
nll	Vector of size M with the current negative log-likelihoods.

Value

None.

See Also

fitZig

HCBravoLab/metagenomeSeq documentation built on March 17, 2024, 3:21 p.m.