convergencePlot: Show EM Algorithm Convergence

Description Usage Arguments Details Author(s) See Also Examples

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 2, 2019, 6:49 p.m.