Initial version. Work-in-progress
Work-in-progress (EM algorithm)