plotConvergence.mnem {mnem} | R Documentation |
This function plots the convergence of the different EM iterations (four figures, e.g. par(mfrow=(2,2))).
## S3 method for class 'mnem' plotConvergence(x, col = NULL, type = "b", convergence = 0.1, ...)
x |
mnem object |
col |
vector of colors for the iterations |
type |
see ?plot.default |
convergence |
difference of when two log likelihoods are considered equal; see also convergence for the function mnem() |
... |
additional parameters for the plots/lines functions |
plot of EM convergence
Martin Pirkl
sim <- simData(Sgenes = 3, Egenes = 2, Nems = 2, mw = c(0.4,0.6)) data <- (sim$data - 0.5)/0.5 data <- data + rnorm(length(data), 0, 1) result <- mnem(data, k = 2, starts = 1) par(mfrow=c(2,2)) plotConvergence(result)