plotseq.npEM | R Documentation |
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.
plot.npEM
, rnormmix
,
npEM
, spEMsymloc
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.