getZ: Calculate the current Z estimate responsibilities (posterior...

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

Calculate the current Z estimate responsibilities (posterior probabilities)

Description

Calculate the current Z estimate responsibilities (posterior probabilities)

Usage

getZ(z, zUsed, stillActive, nll, nllUSED)

Arguments

z	Matrix (m x n) of estimate responsibilities (probabilities that a count comes from a spike distribution at 0).
zUsed	Matrix (m x n) of estimate responsibilities (probabilities that a count comes from a spike distribution at 0) that are actually used (following convergence).
stillActive	A vector of size M booleans saying if a feature is still active or not.
nll	Vector of size M with the current negative log-likelihoods.
nllUSED	Vector of size M with the converged negative log-likelihoods.

Value

A list of updated zUsed and nllUSED.

See Also

fitZig

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