# \donttest{
n <- nrow(iris)
d <- ncol(iris)-1
G <- 3
emEst <- me(modelName="VVI", data=iris[,-5], unmap(iris[,5]))
names(emEst)
args(bic)
bic(modelName="VVI", loglik=emEst$loglik, n=n, d=d, G=G)
# do.call("bic", emEst) ## alternative call
# }
Run the code above in your browser using DataLab