data(stlouis)
s <- prelim.mix(stlouis,3) # preliminary manipulations
thetahat <- em.mix(s) # MLE under unrestricted general location model
loglik.mix(s, thetahat) # loglikelihood at thetahat
Run the code above in your browser using DataLab