Description Usage Arguments Value See Also

View source: R/getNegativeLogLikelihoods.R

Maximum-likelihood estimates are approximated using the EM algorithm where we treat mixture membership $delta_ij$ = 1 if $y_ij$ is generated from the zero point mass as latent indicator variables. The log-likelihood in this extended model is $(1-delta_ij) log f_count(y;mu_i,sigma_i^2 )+delta_ij log pi_j(s_j)+(1-delta_ij)log (1-pi_j (sj))$. The responsibilities are defined as $z_ij = pr(delta_ij=1 | data and current values)$.

1 | ```
getNegativeLogLikelihoods(z, countResiduals, zeroResiduals)
``` |

`z` |
Matrix (m x n) of estimate responsibilities (probabilities that a count comes from a spike distribution at 0). |

`countResiduals` |
Residuals from the count model. |

`zeroResiduals` |
Residuals from the zero model. |

Vector of size M of the negative log-likelihoods for the various features.

metagenomeSeq documentation built on May 31, 2017, 11:53 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.