机器学习算法与Python实战 | 常见统计概率分布实现（内含python代码）

最新推荐文章于 2024-10-10 17:34:13 发布

双木的木

最新推荐文章于 2024-10-10 17:34:13 发布

阅读量1k

点赞数 24

分类专栏： python拓展学习统计学文章标签： python 人工智能开发语言概率论数据分析

本文链接：https://blog.csdn.net/csdn_xmj/article/details/136130105

版权

python拓展学习同时被 2 个专栏收录

155 篇文章 2 订阅

订阅专栏

统计学

17 篇文章 0 订阅

订阅专栏

本文来源公众号“机器学习算法与Python实战”，仅用于学术分享，侵权删，干货满满。

原文链接：https://mp.weixin.qq.com/s/0Lgmdvey70wXQcP1XQvylQ

在平时的科研中，我们经常使用统计概率的相关知识来帮助我们进行城市研究。因此，掌握一定的统计概率相关知识非常有必要。

本文涉及的概念包括：

随机变量(Random Variable)
密度函数(Density Functions)
伯努利分布(Bernoulli Distribution)
二项式分布(Binomial Distribution)
均匀分布(Uniform Distribution)
泊松分布(Poisson Distribution)
正态分布(Normal Distribution)
长尾分布(Long-Tailed Distribution)
学生 t 检验分布(Student’s t-test Distribution)
对数正态分布(Lognormal Distribution)
指数分布(Exponential Distribution)
威布尔分布(Weibull Distribution)
伽马分布(Gamma Distribution)
卡方分布(Chi-square Distribution)
中心极限定理(Central Limit Theorem)

1 随机变量

1.1 离散随机变量

随机实验的所有可能结果都是随机变量。一个随机变量集合用 X 表示。

如果实验可能的结果是可数的，那么它被称为离散随机变量。例如，如果你抛硬币 10 次，你能得到的正面数可以用一个数字表示。或者篮子里有多少苹果仍然是可数的。

1.2 连续随机变量

还有一些是不能以离散方式表示的值。例如，一个人可能有 1.7 米高，1米 80 厘米，1.6666666...米高等等。

2 密度函数

我们使用密度函数来描述随机变量 X 的概率分布。

2.1 PMF：概率质量函数

返回离散随机变量 X 等于 x 的值的概率。所有值的总和等于 1。PMF 只能用于离散变量。

PMF。来源：https://en.wikipedia.org/wiki/Probability_mass_function

2.2 PDF：概率密度函数

它类似于连续变量的 PMF 版本。返回连续随机变量 X 在某个范围内的概率。

PDF。来源：https://byjus.com/maths/probability-density-function/

2.3 CDF：累积分布函数

返回随机变量 X 取小于或等于 x 的值的概率。

CDF（指数分布的累积分布函数）。来源：https://en.wikipedia.org/wiki/Cumulative_distribution_function

3 离散分布

3.1 伯努利分布

我们只有一个试验（只有一个观察结果）和两个可能的结果。例如，抛硬币试验。

我们有一个真的(1)的结果和一个假的(0)的结果。假设我们接受正面为真（我们可以选择正面为真或成功）。那么，如果正面朝上的概率是 p ，相反情况的概率就是 1-p 。

import seaborn as sns
from scipy.stats import bernoulli
# 单一观察值 
# 生成数据 (1000 points, possible outs: 1 or 0, probability: 50% for each)
data = bernoulli.rvs(size=1000,p=0.5)
# 绘制图形
ax = sns.distplot(data_bern,kde=False,hist_kws={"linewidth": 10,'alpha':1})
ax.set(xlabel='Bernouli', ylabel='freq')

3.2 二项式分布

伯努利分布是针对单个观测结果的。多个伯努利观测结果会产生二项式分布。例如，连续抛掷硬币。

试验是相互独立的。一个尝试的结果不会影响下一个。

二项式分布可以表示为B(n, p)。 n 是试验次数， p 是成功的概率。

让我们进行一个实验，我们连续抛掷一枚公平的硬币 20 次。

import matplotlib.pyplot as plt
from scipy.stats import binom
n = 20  
# 实验次数
p = 0.5 
# 成功的概率
r = list(range(n + 1)) 
# the number of success
# pmf值
pmf_list = [binom.pmf(r_i, n, p) for r_i in r ]
# 绘图
plt.bar(r, pmf_list)plt.show()

它看起来像正态分布，但请记住这些值是离散的。

现在这次，你有一枚欺诈硬币。你知道这个硬币正面向上的概率是 0.7。因此，p = 0.7。

带有偏差硬币的二项式分布

该分布显示出成功结果数量增加的概率增加。

p : 成功的概率 n : 实验次数 q=1-p : 失败的概率

3.3 均匀分布

所有结果成功的概率相同。掷骰子，1 到 6。

掷 6 次。

data = np.random.uniform(1, 6, 6000)

掷 6000 次。

3.4 Poisson 分布

它是与事件在给定时间间隔内发生频率相关的分布。

Po(λ) ， λ是在指定时间间隔内预期发生的事件次数。它是在该时间间隔内发生的事件的已知平均值。 X 是事件在指定时间间隔内发生的次数。如果事件遵循泊松分布，则：X ~ Po(λ)。

在泊松分布中，事件彼此独立。事件可以发生任意次数。两个事件不能同时发生。

如每 60 分钟接到 4 个电话。这意味着 60 分钟内通话的平均次数为 4。让我们绘制在 60 分钟内接到 0 到 10 个电话的概率。

import matplotlib.pyplot as plt
from scipy.stats 
import poisson
r = range(0,11) 
# 呼叫次数
lambda_val = 4 
# 均值
# 概率值
data = poisson.pmf(r, lambda_val)
# 绘图
fig, ax = plt.subplots(1, 1, figsize=(8, 6))
ax.plot(r, data, 'bo', ms=8, label='poisson')
plt.ylabel("Probability", fontsize="12")
plt.xlabel("# Calls", fontsize="12")
plt.title("Poisson Distribution", fontsize="16")
ax.vlines(r, 0, data, colors='r', lw=5, alpha=0.5)

结果如下：

4 连续分布

4.1 正态分布

最著名和最常见的分布（也称为高斯分布），是一种钟形曲线。它可以通过均值和标准差定义。正态分布的期望值是均值。

曲线对称。均值、中位数和众数相等。曲线下总面积为 1。

大约 68%的值落在一个标准差范围内。~95% 落在两个标准差范围内，~98.7% 落在三个标准差范围内。

import scipy
mean = 0
standard_deviation = 5
x_values = np. arange(-30, 30, 0.1)
y_values = scipy.stats.norm(mean, standard_deviation)
plt.plot(x_values, y_values. pdf(x_values))

正态分布的概率密度函数为：

QQ 图

我们可以使用 QQ 图来直观地检查样本与正态分布的接近程度。

计算每个数据点的 z 分数并对其进行排序，然后在 y 轴上表示它们。X 轴表示值的排名的分位数。

这个图上的点越接近对角线，分布就越接近正态分布。

import numpy as np
import statsmodels.api as sm
points = np.random.normal(0, 1, 1000)
fig = sm.qqplot(points, line ='45')
plt.show()

4.2 长尾分布

尾巴是分布的长而窄的部分，离群值就位于其中。当一侧尾巴不同于另一侧时，就称为偏斜。下图是长尾分布的 QQ 图。

左偏分布的代码：

import matplotlib.pyplot as plt
from scipy.stats import skewnorm
def generate_skew_data(n: int, max_val: int, skewness: int):    
  # Skewnorm function    
  random = skewnorm.rvs(a = skewness,loc=max_val, size=n)    
  plt.hist(random,30,density=True, color = 'red', alpha=0.1)    
  plt.show()

generate_skew_data(1000, 100, -5) # negative (-5)-> 左偏分布

右偏分布的代码：

generate_skew_data(1000, 100, 5) # positive (5)-> 右偏分布

学生 t 检验分布

正态但有尾（更厚、更长）。

t 分布和 z 分布。来源：https://www.geeksforgeeks.org/students-t-distribution-in-statistics/

t 分布是具有较厚尾部的正态分布。如果可用数据较少（约 30 个），则使用 t 分布代替正态分布。

在 t 分布中，自由度变量也被考虑在内。根据自由度和置信水平在 t 分布表中找到关键的 t 值。这些值用于假设检验。

t 分布表可以移步：https://www.sjsu.edu/faculty/gerstman/StatPrimer/t-table.pdf。

4.3 对数正态分布

随机变量 X 的对数服从正态分布的分布。

import numpy as np
import matplotlib.pyplot as plt
from scipy import stats
X = np.linspace(0, 6, 1500)
std = 1
mean = 0
lognorm_distribution = stats.lognorm([std], loc=mean)
lognorm_distribution_pdf = lognorm_distribution.pdf(X)
fig, ax = plt.subplots(figsize=(8, 5))
plt.plot(X, lognorm_distribution_pdf, label="μ=0, σ=1")
ax.set_xticks(np.arange(min(X), max(X)))
plt.title("Lognormal Distribution")
plt.legend()plt.show()

结果如下：

4.4 指数分布

我们在 Poisson 分布中研究了在一定时间间隔内发生的事件。在指数分布中，我们关注的是两个事件之间经过的时间。如果我们把上面的例子倒过来，那么两个电话之间需要多长时间？

因此，如果 X 是一个随机变量，遵循指数分布，则累积分布函数为：

代码如下：

from scipy.stats import expon
import matplotlib.pyplot as plt
x = expon.rvs(scale=2, size=10000) # 2 calls
# 绘图
plt.hist(x, density=True, edgecolor='black')

x 轴表示时间间隔的百分比

4.5 韦伯分布

它是指时间间隔是可变的而不是固定的情况下使用的指数分布的扩展。在 Weibull 分布中，时间间隔被允许动态变化。

β 是形状参数，如果是正值，则事件发生的概率随时间而增加，反之亦然。 η 是尺度参数。

代码如下：


import matplotlib.pyplot as plt
x = np.arange(1,100.)/50.
def weib(x,n,a):    
  return (a / n) * (x / n)**(a - 1) * np.exp(-(x / n)**a)

count, bins, ignored = plt.hist(np.random.weibull(5.,1000))
x = np.arange(1,100.)/50.
scale = count.max()/weib(x, 1., 5.).max()
plt.plot(x, weib(x, 1., 5.)*scale)
plt.show()

结果如下：

4.6 Gamma 分布

指与第 n 个事件发生所需的时间有关的分布，而指数分布则与首次事件发生的时间有关。

代码如下：


import numpy as np
import scipy.stats as stats
import matplotlib.pyplot as plt
#Gamma distributions
x = np.linspace(0, 60, 1000)
y1 = stats.gamma.pdf(x, a=5, scale=3)
y2 = stats.gamma.pdf(x, a=2, scale=5)
y3 = stats.gamma.pdf(x, a=4, scale=2)
# plots
plt.plot(x, y1, label='shape=5, scale=3')
plt.plot(x, y2, label='shape=2, scale=5')
plt.plot(x, y3, label='shape=4, scale=2')
#add legend
plt.legend()
#display 
plotplt.show()

Gamma 分布。X 轴表示随机变量 X 可能取到的潜在值，Y 轴表示分布的概率密度函数（PDF）值

4.7 卡方分布

它用于统计检验。这通常在实际分布中不会出现。

# x轴范围0-10，步长0.25
X = np.arange(0, 10, 0.25)
plt.subplots(figsize=(8, 5))
plt.plot(X, stats.chi2.pdf(X, df=1), label="1 dof")
plt.plot(X, stats.chi2.pdf(X, df=2), label="2 dof")
plt.plot(X, stats.chi2.pdf(X, df=3), label="3 dof")
plt.title("Chi-squared Distribution")
plt.legend()
plt.show()

结果如下：