我试图采用中提出的这个解来确定一个简单正态分布的参数。尽管修改是很小的(基于维基百科),但结果却很差。有什么不对的地方吗?import numpy as npimport matplotlib.pyplot as plt
return 1./(math.sqrt(2.*math.pi)*sig)*np.exp(-np.power((x - mu)/sig, 2.)/2)
mu
我正在尝试用optim()实现混合高斯混合的MLE,使用R的本地数据集(来自质量的Geyser)。我的密码在下面。问题是optim工作得很好,但是,返回我传递给它的原始参数,并且还说它已经收敛。library(ggplot2)data("geyser")x.vector=externaldata
{ combined.vector=bigvec K