Plotting sequences of estimates from non- or semiparametric EM-like Algorithm
Returns plots of the sequences of scalar parameter
estimates along iterations from an object of class npEM.
## S3 method for class 'npEM' plotseq(x, ...)
x |
an object of class |
... |
further parameters that are passed to |
plotseq.npEM returns a figure with one plot for each component
proportion, and, in the case of spEMsymloc, one plot for each
component mean.
Benaglia, T., Chauveau, D., and Hunter, D. R. (2009), An EM-like algorithm for semi- and non-parametric estimation in multivariate mixtures, Journal of Computational and Graphical Statistics (to appear).
Bordes, L., Chauveau, D., and Vandekerkhove, P. (2007), An EM algorithm for a semiparametric mixture model, Computational Statistics and Data Analysis, 51: 5429-5443.
## Example from a normal location mixture n <- 200 set.seed(100) lambda <- c(1/3,2/3) mu <- c(0, 4); sigma<-rep(1, 2) x <- rnormmix(n, lambda, mu, sigma) b <- spEMsymloc(x, mu0=c(-1, 2), stochastic=FALSE) plotseq(b) bst <- spEMsymloc(x, mu0=c(-1, 2), stochastic=TRUE) plotseq(bst)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.