旅鼠 感谢,但是lower=c(0.0001,-Inf,0.0001,-Inf,-Inf,0.0001),应该是lower=c(0.0001,-Inf,0.0001,-Inf,0.0001),此处是笔误。
nan.xiao 这里的初始值就是通过先画图再一拍大腿估计得到的,俗称 eyeball estimation。 不过说句心里话,既然用了这个例子,文章的重点应该放在 mixture model 上,而不是推导求解需要用到的 MLE ……