convergencePlot: Show EM Algorithm Convergence

Description Usage Arguments Details Author(s) See Also Examples

View source: R/model_reporting.R

Description

Plot convergence of the log-likelihood to a maximum during execution of the EM algorithm for fitting mixture models. Optionally also show the Q-value, that is, the expected complete data log-likelihood.

Usage

1
convergencePlot(x, show.qval=FALSE, main="")

Arguments

x

an object of class mixmod or mdmixmod.

show.qval

logical; if TRUE, the Q-value will be plotted as well as the log-likelihood.

main

plot title.

Details

Log-likelihood is plotted vs. iteration in black. If show.qval is TRUE, Q-value is plotted in red.

Author(s)

Daniel Dvorkin

See Also

mixmod, mdmixmod, rocinfo; plot in package graphics.

Examples

1
2
3
4
5
6
## Not run: 
data(CiData)
fit <- mixmod(CiData$expression, 3)
convergencePlot(fit, TRUE, main="log-likelihood and Q-value")

## End(Not run)


lcmix documentation built on May 21, 2017, 4:05 a.m.
Search within the lcmix package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.