概率统计Python计算:卡方分布分位点计算
nnn个相互独立,均服从N(0,1)N(0,1)N(0,1)的随机变量X1,X2⋯,XnX_1, X_2\cdots,X_nX1,X2⋯,Xn的平方和X12+X22+⋯+Xn2X_1^2+X_2^2+\cdots+X_n^2X12+X22+⋯+Xn2服从自由度为nnn的χ2\chi^2χ2分布,其密度函数为
f(x)={12n/2Γ(n/2)xn2−1e−x2x≥00x<0.{f(x)=} \begin{cases} \frac{1}{2^{n/2}\Gamma(n/2)}x^{\frac{n}{2}-1}e^{-\frac{x}{2}}&x\geq0\\ 0&x<0 \end{cases}.f(x)={2n/2Γ(n/2)1x2n−1e−2x0x≥0x<0.
下图展示了自由度nnn分别为2,4,6的f(x)f(x)f(x)的图像。
历史上均以残存函数S(x)=1−F(x)S(x)=1-F(x)S(x)=1−F(x)的反函数S−1(α)S^{-1}(\alpha)S−1(α)来表示单侧左分位点χα2(n)\chi^2_{\alpha}(n)χα2(n)的,即对显著水平α\alphaα
P(X≥χα2(n))=α.P(X\geq\chi^2_{\alpha}(n))=\alpha.P(X≥χα2(n))=α.
如下图所示:
将单侧右分位点表为χ1−α2(n)\chi^2_{1-\alpha}(n)χ1−α2(n),即1−P(X<χ1−α2(n))=α1-P(X<\chi^2_{1-\alpha}(n))=\alpha1−P(X<χ1−α2(n))=α,如下图所示。而1−P(X>χ1−α2(n))=P(X≤χ1−α2(n))=F(χ1−α2(n))1-P(X>\chi^2_{1-\alpha}(n))=P(X\leq\chi^2_{1-\alpha}(n))=F(\chi^2_{1-\alpha}(n))1−P(X>χ1−α2(n))=P(X≤χ1−α2(n))=F(χ1−α2(n))。
而双侧左右分位点表为满足P(χ1−α/22(n)≤X≤χα/22(n))=1−αP(\chi^2_{1-\alpha/2}(n)\leq X\leq\chi^2_{\alpha/2}(n))=1-\alphaP(χ1−α/22(n)≤X≤χα/22(n))=1−α的χ1−α/22(n)\chi^2_{1-\alpha/2}(n)χ1−α/22(n)和χα/22(n)\chi^2_{\alpha/2}(n)χα/22(n),如下图所示。
Python的scipy.stats包中,连续型分布类 rv_continuous的chi2对象表示χ2\chi^2χ2分布。χ2\chi^2χ2分布常用的计算分位点函数的调用接口见下表。
函数名 | 参数 | 意义 |
---|---|---|
ppf | q:表示显著水平α\alphaα,df:表示分布的自由度nnn | 左分位点χ1−α2(n)\chi_{1-\alpha}^2(n)χ1−α2(n) |
isf | q,df:与上同 | 右分位点χα2(n)\chi_{\alpha}^2(n)χα2(n) |
interval | alpha:表示置信水平1−α1-\alpha1−α,df:与上同 | 双侧分位点χ1−α/22(n)\chi_{1-\alpha/2}^2(n)χ1−α/22(n)和χα/22(n)\chi_{\alpha/2}^2(n)χα/22(n) |
例1 对显著水平α=0.05\alpha=0.05α=0.05,计算自由度n=24n=24n=24的χ2(n)\chi^2(n)χ2(n)分布的分位点。
解:下列代码完成本例计算。
from scipy.stats import chi2 #导入chi2
n=24 #设置自由度n
alpha=0.05 #设置alpha
a=chi2.ppf(q=alpha, df=n) #计算单侧左分位点
b=chi2.isf(q=alpha, df=n) #计算单侧右分位点
print('单侧左、右分位点:a=%.4f, b=%.4f'%(a, b))
a, b=chi2.interval(1-alpha, df=n) #计算双侧左分位点
print('双侧左、右分位点:a=%.4f, b=%.4f'%(a, b))
运行程序,输出
单侧左、右分位点:a=13.8484, b=36.4150
双侧左、右分位点:a=12.4012, b=39.3641
写博不易,敬请支持:
如果阅读本文于您有所获,敬请点赞、评论、收藏,谢谢大家的支持!
返回《导引》
概率统计Python计算:卡方分布分位点计算相关推荐
- 概率统计Python计算:标准正态分布分位点计算
标准正态分布对给定显著水平的分位点.设XXX~N(0,1)N(0,1)N(0,1),显著水平为α\alphaα.为计算右侧分位点zαz_{\alpha}zα(见下图),使得 P(X≤zα)=1−αP ...
- 概率统计Python计算:双正态总体方差比的双侧区间估计
计算两个正态总体方差比的区间估计涉及样本方差 s 1 2 s_1^2 s12, s 2 2 s_2^2 s22,样本容量 n 1 n_1 n1, n 2 n_2 n2和置信度 1 − α 1- ...
- 概率统计Python计算:条件概率和概率乘法公式
1. 古典概型中条件概率的计算 条件概率 P ( B ∣ A ) P(B|A) P(B∣A)是将样本空间限制在 A A A上, A ∩ B A\cap B A∩B的概率.因此,我们可以利用博文< ...
- 概率统计Python计算:F分布分位点计算
设XXX,YYY相互独立,且分别服从χ2(m)\chi^2(m)χ2(m)和χ2(n)\chi^2(n)χ2(n),则XY\frac{X}{Y}YX~F(m−1,n−1)F(m-1, n-1)F(m ...
- 概率统计Python计算:学生分布分位点计算
设XXX~N(0,1)N(0,1)N(0,1),YYY~χ2(n)\chi^2(n)χ2(n),且XXX与YYY相互独立,则XY/n\frac{X}{\sqrt{Y/n}}Y/nX~t(n)t(n ...
- 概率统计Python计算:单个正态总体均值双侧假设的Z检验
用p值法计算假设H0H_0H0的双侧检验,设aaa,bbb分别是检验统计量分布对应显著水平α\alphaα的左.右分位点,F(x)F(x)F(x)和S(x)S(x)S(x)分别为检验统计量的分布函数 ...
- 概率统计Python计算:单个正态总体均值单侧假设的Z检验
对正态总体参数的单侧假设检验,可以用如下的p值法进行.设显著水平为α\alphaα,考虑假设H0H_0H0的右侧检验.首先,注意到检验统计量的分布对应显著水平α\alphaα的右分位点bbb,实际上 ...
- python求上α分位点_什么是标准正态分布的上α分位点以及怎样求?
展开全部 标准正态分布的上α分位点:设X~N(0,1),对于任给的62616964757a686964616fe4b893e5b19e31333431363638α,(0Zα)= α的点Zα为标准正态 ...
- 概率统计Python计算:离散型随机变量分布(bernoulli geom)
Python的scipy.stats包中提供了各种随机变量的分布.每一种分布,其累积分布函数(分布函数)记为cdf.离散型变量分布的概率质量函数(分布律),记为pmf.除此之外,每个分布都有一个服从该 ...
最新文章
- access 查找工龄大于30_同济大学大学计算机access作业答案
- 信息学奥赛一本通 1324:【例6.6】整数区间
- 【环境搭建002】ubuntu 13 在vm 下的 NFS 搭建
- python零基础知识_python 零基础 基本知识
- 通过telnet自动下载cfg配置文件
- android 音乐播放器-------歌词同步 lrc
- 51Nod-2149 字符串出现位置【水题】
- VC6删除项目中的文件
- 使用MSChart控件绘制图表。
- java 导出word模板
- 股票实时数据如何分析?
- 苹果手机备份有必要吗_苹果手机备份最简单方式!原来大神都这样做
- 视频演示-Snapper快捷优秀的音频预览播放器演示
- 嘘,这份0项目经验应对面试的独门秘籍,千万别被你的面试官发现~
- 2020年河南省计算机对口升学分数线,2020年河南省单招分数线出来了吗?
- 文件(file)操作
- Hyperledger Fabric2中文文档-修改通道配置
- 关于一个学习计算机专业,迷茫的大一新生的看法和理解
- Qt 如何实现的 Meta Object
- dojo引用html模板,深入浅出dojo/request-本地HTML文件 带参数方案-遮罩层《一》_169IT.COM...
热门文章
- Scrum立会报告+燃尽图(十月十一日总第二次):需求分析
- 常见的WebShell管理工具(top10)
- 下载源码报错Cannot connect to the Maven process. Try again later. If the problem persists, check the Maven
- 快手光合计划完整版攻略
- google订阅消息推送
- 微信公众账号高级接口使用小结
- FPGA实现按键检测消抖程序
- PS 一键生成Android各个尺寸的图标
- 铜九铁路客运将于2008年9月1日正式开通
- 社群运营岗位部门团队KPI绩效sop工作流程计划方案