data(iris)
irisMatrix <- as.matrix(iris[,1:4])
irisClass <- iris[,5]
msEst <- mstep(modelName = "EEE", data = irisMatrix,
z = unmap(irisClass))
names(msEst)
em(modelName = msEst$modelName, data = irisMatrix,
mu = msEst$mu, Sigma = msEst$Sigma, pro = msEst$pro)
do.call("em", c(list(data = irisMatrix), msEst)) ## alternative call
Run the code above in your browser using DataLab