----------------------------------------------------------------
iteration number    1, loglikelihood = -1223.3493 
iteration number    2, loglikelihood = -860.4353