python简单实现最大似然估计&scipy库的使用

python简单实现最大似然估计

1、scipy库的安装
wim+R输入cmd，然后cd到python的pip路径，即安装：pip install scipy即可

2、导入scipy库
from scipy.sats import norm
导入scipy.sats中的norm
3、案例分析

from scipy.stats import norm
import matplotlib.pyplot as plt
import numpy as np

'''
norm.cdf 返回对应的累计分布函数值
norm.pdf 返回对应的概率密度函数值
norm.rvs 产生指定参数的随机变量
norm.fit 返回给定数据下，各参数的最大似然估计（MLE）值
'''
x_norm = norm.rvs(size=200)
#在这组数据下，正态分布参数的最大似然估计值
x_mean, x_std = norm.fit(x_norm)
print ('mean, ', x_mean)
print ('x_std, ', x_std)
plt.hist(x_norm, normed=True, bins=15)#归一化直方图（用出现频率代替次数），将划分区间变为 20（默认 10）
x = np.linspace(-3,3,50)#在在(-3,3)之间返回均匀间隔的50个数字。
plt.plot(x, norm.pdf(x), 'r-')
plt.show()

运行如下：

本文是通过官网地址：https://nbviewer.jupyter.org/github/lijin-thu/notes-python/blob/master/04-scipy/04.03-statistics-with-scipy.ipynb进行学习

**感悟： 有一次去听有关于大数据到人工智能（Big Data To AI）的讲座，谢老师说：“数据是呈现一种现象，状态的一种指标”老师用这一简短的一句解释了什么是数据。**