概率统计Python计算:离散型随机变量分布(bernoulli geom)

本文主要是介绍概率统计Python计算:离散型随机变量分布(bernoulli geom),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

在这里插入图片描述
Python的scipy.stats包中提供了各种随机变量的分布。每一种分布,其累积分布函数(分布函数)记为cdf。离散型变量分布的概率质量函数(分布律),记为pmf。除此之外,每个分布都有一个服从该分部变量发生器函数rvs,用来产生服从该分布的随机数。

1. bernoulli分布(0-1分布)

Python的scipy.stats包中,bernoulli类就是用来表示伯努利分布的。常用的三个函数说明见下表。

函数名参数功能
rvs(p, size)p:分布参数,size:产生的随机数个数,缺省值为1产生size个随机数
pmf(k, p)k:随机变量取值,p:与上同概率质量函数(分布律) P ( X = k ) P(X=k) P(X=k)
cdf(k, p)k:分布函数自变量,p:与上同累积概率函数(分布函数) F ( k ) F(k) F(k)
例1 下列代码利用bernoulli类对象的rvs函数模拟重复抛掷均匀分币试验。
from scipy.stats import bernoulli   #导入bernoulli
import numpy as np                  #导入numpy
x=bernoulli.rvs(p=1/2,size=500)     #产生500个服从p=1/2的0-1分布的随机数
hist, _=np.histogram(x, bins=2)     #统计取0、1的频数
hist/500                            #输出频率

其中的第3行调用bernoulli类对象的随机数发生函数rvs产生500个服从参数 p = 1 / 2 p=1/2 p=1/2的0-1分布(抛掷均匀分币0,1分别表示正面朝下和正面朝上)的随机数。第4行调用numpy的histogram函数统计500个数据中取0、1的频数。第5行输出频率。运行程序,输出

array([0.498, 0.502])

可见取0和1的频率分别为0.499和0.502,很好地模拟了抛掷均匀分币这一伯努利试验。

2. geom分布(几何分布)

scipy.stats包提供的geom类表示几何分布。常用的三个函数rvs、pmf和cdf的名称、参数和意义是完全一致的。这是因为0-1分布和几何分布均仅有一个表示一次试验成功概率的参数 p p p。由此可见,引入随机变量处理不同随机试验下的随机事件概率问题的形式是统一的。
例2 设灯泡在任意一天损坏的概率 p = 0.001 p=0.001 p=0.001,计算该灯泡的寿命至少为30天的概率。
:设灯泡的使用寿命(单位:天)为 X X X,则 X X X服从参数为 p = 0.001 p=0.001 p=0.001的几何分布。令 q = 1 − p = 0.999 q=1-p=0.999 q=1p=0.999,灯泡寿命至少为30天的概率
P ( X ≥ 30 ) = 1 − P ( X ≤ 29 ) = 1 − F ( 29 ) = 1 − ∑ k = 1 29 q 29 − k p = 1 − p 1 − q 29 1 − q = q 29 = 0.99 9 29 = 0.9714. P(X\geq30)=1-P(X\leq29)=1-F(29)\\ =1-\sum\limits_{k=1}^{29}q^{29-k}p=1-p\frac{1-q^{29}}{1-q}\\ =q^{29}=0.999^{29}=0.9714. P(X30)=1P(X29)=1F(29)=1k=129q29kp=1p1q1q29=q29=0.99929=0.9714.
下列代码验算本例中灯泡寿命至少为30天的概率。

from scipy.stats import geom        #导入geom
prob=1-geom.cdf(k=29,p=0.001)       #计算1-F(29)
print('P(X>=30)=1-F(29)=%.4f'%prob) #输出P(X>=30)

程序的第2行调用geom(第1行导入)的cdf函数,计算 1 − F ( 29 ) = 1 − ∑ k = 1 29 ( 1 − p ) k − 1 p 1-F(29)=1-\sum\limits_{k=1}^{29}(1-p)^{k-1}p 1F(29)=1k=129(1p)k1p。运行程序,输出

P(X>=30)=1-F(29)=0.9714

scipy.stats为每一种分布提供残存函数sf,该函数计算 P ( X > x ) = 1 − P ( X ≤ x ) = 1 − cdf ( x ) P(X>x)=1-P(X\leq x)=1-\text{cdf}(x) P(X>x)=1P(Xx)=1cdf(x)。例如,在上列程序中,将第2行代码换成
prob=geom.sf(k=29, p=0.01) \text{prob=geom.sf(k=29, p=0.01)} prob=geom.sf(k=29, p=0.01)
来计算服从参数为 p = 0.01 p=0.01 p=0.01的几何分布的随机变量 X X X的概率 P ( X ≥ 30 ) = P ( X > 29 ) P(X\geq30)=P(X>29) P(X30)=P(X>29)
写博不易,敬请支持:
如果阅读本文于您有所获,敬请点赞、评论、收藏,谢谢大家的支持!
代码诚可贵,原理价更高。若为AI学,读正版书好
返回《导引》

这篇关于概率统计Python计算:离散型随机变量分布(bernoulli geom)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1097641

相关文章

Python FastAPI+Celery+RabbitMQ实现分布式图片水印处理系统

《PythonFastAPI+Celery+RabbitMQ实现分布式图片水印处理系统》这篇文章主要为大家详细介绍了PythonFastAPI如何结合Celery以及RabbitMQ实现简单的分布式... 实现思路FastAPI 服务器Celery 任务队列RabbitMQ 作为消息代理定时任务处理完整

Python Websockets库的使用指南

《PythonWebsockets库的使用指南》pythonwebsockets库是一个用于创建WebSocket服务器和客户端的Python库,它提供了一种简单的方式来实现实时通信,支持异步和同步... 目录一、WebSocket 简介二、python 的 websockets 库安装三、完整代码示例1.

揭秘Python Socket网络编程的7种硬核用法

《揭秘PythonSocket网络编程的7种硬核用法》Socket不仅能做聊天室,还能干一大堆硬核操作,这篇文章就带大家看看Python网络编程的7种超实用玩法,感兴趣的小伙伴可以跟随小编一起... 目录1.端口扫描器:探测开放端口2.简易 HTTP 服务器:10 秒搭个网页3.局域网游戏:多人联机对战4.

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

Python使用自带的base64库进行base64编码和解码

《Python使用自带的base64库进行base64编码和解码》在Python中,处理数据的编码和解码是数据传输和存储中非常普遍的需求,其中,Base64是一种常用的编码方案,本文我将详细介绍如何使... 目录引言使用python的base64库进行编码和解码编码函数解码函数Base64编码的应用场景注意

Python基于wxPython和FFmpeg开发一个视频标签工具

《Python基于wxPython和FFmpeg开发一个视频标签工具》在当今数字媒体时代,视频内容的管理和标记变得越来越重要,无论是研究人员需要对实验视频进行时间点标记,还是个人用户希望对家庭视频进行... 目录引言1. 应用概述2. 技术栈分析2.1 核心库和模块2.2 wxpython作为GUI选择的优

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.

Python+PyQt5实现多屏幕协同播放功能

《Python+PyQt5实现多屏幕协同播放功能》在现代会议展示、数字广告、展览展示等场景中,多屏幕协同播放已成为刚需,下面我们就来看看如何利用Python和PyQt5开发一套功能强大的跨屏播控系统吧... 目录一、项目概述:突破传统播放限制二、核心技术解析2.1 多屏管理机制2.2 播放引擎设计2.3 专

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很