Iterative dimension estimation works with function mlcc.reps
For a time being mlcc.kmeans is exported.
Iteratively adjusting penalty is implemented in adjusted.cluster.BIC
adding TODO list
extending iterative estimation for mlcc.bic
changing plot and print for mlcc.fit object