有没有人知道替代scipy.stats.norm.pdf()函数的方法?我将我的python网站托管在Google App Engine上,而Google不支持SciPy。
我尝试了这个函数,但它并没有像scipy那样返回相同的结果:
def normpdf(x, mu, sigma):
u = (x-mu)/abs(sigma)
y = (1/(sqrt(2*pi)*abs(sigma)))*exp(-u*u/2)
return y
例如:
print scipy.stats.norm.pdf(20, 20, 10)
print normpdf(20, 20, 10)
print scipy.stats.norm.pdf(15, 20, 10)
print normpdf(15, 20, 10)
print scipy.stats.norm.pdf(10, 20, 10)
print normpdf(10, 20, 10)
返回以下数值:
0.0398942280401
0.0398942280401
0.0352065326764
0.0146762663174
0.0241970724519
0.0146762663174