MoE_plotLogLik function

Plot the Log-Likelihood of a MoEClust Mixture Model

Plots the log-likelihood at every iteration of the EM/CEM algorithm used to fit a MoEClust mixture model.

MoE_plotLogLik(res,
               type = "l",
               xlab = "Iteration",
               ylab = "Log-Likelihood",
               xaxt = "n",
               ...)

Arguments

• res: An object of class "MoEClust" generated by MoE_clust, or an object of class "MoECompare" generated by MoE_compare. Models with a noise component are facilitated here too.
• type, xlab, ylab, xaxt: These graphical parameters retain their usual definitions from plot.
• ...: Catches unused arguments, or additional arguments to be passed to plot.

Returns

A plot of the log-likelihood versus the number EM iterations. A list with the vector of log-likelihood values and the final value at convergence can also be returned invisibly.

Note

plot.MoEClust is a wrapper to MoE_plotLogLik which accepts the default arguments, and also produces other types of plots.

res$LOGLIK can also be plotted, to compare maximal log-likelihood values for all fitted models.

Examples

data(ais)
res <- MoE_clust(ais[,3:7], gating= ~ BMI, expert= ~ sex, tau0=0.1,
                 G=2, modelNames="EVE", network.data=ais)
(ll <- MoE_plotLogLik(res))

See Also

MoE_clust, plot.MoEClust,

Author(s)

Keefe Murphy - <keefe.murphy@mu.ie >