开发Python函数normdist(x,mu,sigma)

时间:2016-08-15 08:35:16

标签: python math statistics

我不明白我收到的这个家庭作业问题,我们的任务是开发Python函数normdist(x,mu,sigma),它评估多变量 k维向量x的高斯概率密度函数,均值向量μ和协方差 矩阵Σ。在k = 1的特殊情况下,该函数评估单变量高斯概率 标量x的密度函数,平均值μ和标准偏差σ。

我的尝试如下:

def normcdf(x, mu, sigma):
    t = x-mu;
    y = 0.5*erfcc(-t/(sigma*sqrt(2.0)));
    if y>1.0:
        y = 1.0;
    return y

def normpdf(x, mu, sigma):
    u = (x-mu)/abs(sigma)
    y = (1/(sqrt(2*pi)**k*abs(sigma)))*exp(-u*u/2)
    return y

def normdist(x, mu, sigma, k):
    if k:
        y = normcdf(x,mu,sigma)
    else:
        y = normpdf(x,mu,sigma)
    return y

以上代码记入Cerin

我如何处理k = 1的情况?

0 个答案:

没有答案