概率统计Python计算:假设检验应用——基于成对数据的检验

2024-08-22 22:48

本文主要是介绍概率统计Python计算:假设检验应用——基于成对数据的检验,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

在这里插入图片描述
设有 n n n个相互独立的观测结果 ( X 1 , Y 1 ) (X_1,Y_1) (X1,Y1) ( X 2 , Y 2 ) (X_2,Y_2) (X2,Y2) ⋯ \cdots ( X n , Y n ) (X_n,Y_n) (Xn,Yn),诸对 X i X_i Xi Y i Y_i Yi受同一因素影响, D i = X i − Y i D_i=X_i-Y_i Di=XiYi~ N ( μ , σ 2 ) , i = 1 , 2 , ⋯ , n N(\mu,\sigma^2),i=1,2,\cdots,n N(μ,σ2),i=1,2,,n。其中 μ \mu μ σ 2 \sigma^2 σ2均未知。在指定显著水平 α \alpha α下,检验假设
H 0 : μ = μ 0 ( H 1 : μ ≠ μ 0 ) 或 H 0 : μ ≤ μ 0 ( H 1 : μ > μ 0 ) 或 H 0 : μ ≥ μ 0 ( H 1 : μ < μ 0 ) . H_0:\mu=\mu_0(H_1:\mu\not=\mu_0)\text{或}\\H_0:\mu\leq\mu_0(H_1:\mu>\mu_0)\text{或} \\H_0:\mu\geq\mu_0(H_1:\mu<\mu_0). H0:μ=μ0(H1:μ=μ0)H0:μμ0(H1:μ>μ0)H0:μμ0(H1:μ<μ0).
的问题,称为基于成对数据的检验问题。由于 D i D_i Di~ N ( μ , σ 2 ) , i = 1 , 2 , ⋯ , n N(\mu,\sigma^2),i=1,2,\cdots,n N(μ,σ2),i=1,2,,n,且 σ 2 \sigma^2 σ2未知,故可用 t t t检验法解决基于成对数据的检验问题。
由于 D i = X i − Y i D_i=X_i-Y_i Di=XiYi~ N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2) i = 1 , 2 , ⋯ , n i=1,2,\cdots,n i=1,2,,n,其中 σ 2 \sigma^2 σ2未知。在显著水平 α = 0.05 \alpha=0.05 α=0.05下,为检验假设 H 0 : μ = μ 0 H_0:\mu=\mu_0 H0:μ=μ0(或 H 0 : μ ≤ μ 0 H_0:\mu\leq\mu_0 H0:μμ0 H 0 : μ ≥ μ 0 H_0:\mu\geq\mu_0 H0:μμ0),scipy.stats包提供了函数
ttest_1samp(a, popmean, alternative=’two-sided’) \text{ttest\_1samp(a, popmean, alternative='two-sided')} ttest_1samp(a, popmean, alternative=’two-sided’)
其参数a表示序列 { d 1 = x 1 − y 1 , d 2 = x 2 − y 2 , ⋯ , d n = x n − y n } \{d_1=x_1-y_1,d_2=x_2-y_2,\cdots,d_n=x_n-y_n\} {d1=x1y1,d2=x2y2,,dn=xnyn},popmean表示 μ \mu μ的假设值 μ 0 \mu_0 μ0,alternative为三个选项之一’two-sided’,‘greater’或’less’,分别表示双侧假设、右侧假设及左侧假设,缺省值为表示双侧假设的’two-sided’。该函数的返回值包括两个数据:表示检验统计量值 d ‾ − μ 0 s / n \frac{\overline{d}-\mu_0}{s/\sqrt{n}} s/n dμ0的statistic和表示检验p值的pvalue。
例1将双胞胎分开来抚养,一个由父母亲自带大,另一个不是由父母亲自带大。现取14对双胞胎测试他们的智商,智商测试得分如下:

序号1234567891011121314
父母带大 X i X_i Xi2331251819252818252822143436
非父母带大 Y i Y_i Yi2231292428312715232726193028

希望比较两种不同的成长环境是否对孩子的智商有不同的影响。
解: 首先,注意到表中数据是成对出现的,14对双胞胎的智商 ( X i , Y i ) , i = 1 , 2 , ⋯ , 14 (X_i,Y_i),i=1,2,\cdots,14 (Xi,Yi),i=1,2,,14,任何两对双胞胎的成长过程互不影响,故可以认为是相互独立的。同一家庭的双胞胎孩子的智商 X i X_i Xi Y i Y_i Yi应当由某种联系。为考察一对双胞胎的不同成长环境对智商的影响,考虑两者的差 D i = X i − Y i , i = 1 , 2 , ⋯ , 14 D_i=X_i-Y_i,i=1,2,\cdots,14 Di=XiYi,i=1,2,,14。由于诸 D i D_i Di均受同样因素(孩子的成长环境)影响,故可认为具有相同的分布。假定 D i D_i Di~ N ( μ , σ 2 ) , i = 1 , 2 , ⋯ , 14 N(\mu,\sigma^2),i=1,2,\cdots,14 N(μ,σ2),i=1,2,,14(由诸 ( X i , Y i ) (X_i,Y_i) (Xi,Yi)的相互独立性,知诸 D i D_i Di也是相互独立的),本例即是需要在一定的显著水平 α \alpha α下检验假设
H 0 : μ = 0 ( H 1 : μ ≠ 0 ) . H_0:\mu=0(H_1:\mu\not=0). H0:μ=0(H1:μ=0).
下列代码完成例7-26中对双侧假设 H 0 H_0 H0在显著水平 α = 0.05 \alpha=0.05 α=0.05下的检验计算。

import numpy as np                      #导入numpy
from scipy.stats import ttest_1samp     #导入ttest_1samp
x=np.array([23, 31, 25, 18, 19, 25, 28, #设置样本数据18, 25, 28, 22, 14, 34, 36])
y=np.array([22, 31, 29, 24, 28, 31, 27,15, 23, 27, 26, 19, 30, 28])
alpha=0.05                              #显著水平
d=x-y                                   #计算di=xi-yi
_,pvalue=ttest_1samp(d, 0)              #计算检验p值
print('mu=0 is %s.'%(pvalue>=alpha))

程序的第3~7行按题面设置各项数据。第8行计算序列 d = { d 1 = x 1 − y 1 , d 2 = x 2 − y 2 , ⋯ , d n = x n − y n } d=\{d_1=x_1-y_1,d_2=x_2-y_2,\cdots,d_n=x_n-y_n\} d={d1=x1y1,d2=x2y2,,dn=xnyn},记为d。第9行调用函数ttest_1samp计算检验假设 H 0 H_0 H0的p值(由于此处不需要检验统计量值,故用下划线’_'将返回值中的statistic屏蔽掉),第10行计算检验并输出。

mu=0 is True.

表示接受假设 H 0 : μ = 0 H_0:\mu=0 H0:μ=0,即生长环境对孩子的智商没有显著影响。
写博不易,敬请支持:
如果阅读本文于您有所获,敬请点赞、评论、收藏,谢谢大家的支持!
代码诚可贵,原理价更高。若为AI学,读正版书好
返回《导引》

这篇关于概率统计Python计算:假设检验应用——基于成对数据的检验的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

Python中你不知道的gzip高级用法分享

《Python中你不知道的gzip高级用法分享》在当今大数据时代,数据存储和传输成本已成为每个开发者必须考虑的问题,Python内置的gzip模块提供了一种简单高效的解决方案,下面小编就来和大家详细讲... 目录前言:为什么数据压缩如此重要1. gzip 模块基础介绍2. 基本压缩与解压缩操作2.1 压缩文

MySQL 删除数据详解(最新整理)

《MySQL删除数据详解(最新整理)》:本文主要介绍MySQL删除数据的相关知识,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录一、前言二、mysql 中的三种删除方式1.DELETE语句✅ 基本语法: 示例:2.TRUNCATE语句✅ 基本语

Python设置Cookie永不超时的详细指南

《Python设置Cookie永不超时的详细指南》Cookie是一种存储在用户浏览器中的小型数据片段,用于记录用户的登录状态、偏好设置等信息,下面小编就来和大家详细讲讲Python如何设置Cookie... 目录一、Cookie的作用与重要性二、Cookie过期的原因三、实现Cookie永不超时的方法(一)

Python内置函数之classmethod函数使用详解

《Python内置函数之classmethod函数使用详解》:本文主要介绍Python内置函数之classmethod函数使用方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地... 目录1. 类方法定义与基本语法2. 类方法 vs 实例方法 vs 静态方法3. 核心特性与用法(1编程客

Python函数作用域示例详解

《Python函数作用域示例详解》本文介绍了Python中的LEGB作用域规则,详细解析了变量查找的四个层级,通过具体代码示例,展示了各层级的变量访问规则和特性,对python函数作用域相关知识感兴趣... 目录一、LEGB 规则二、作用域实例2.1 局部作用域(Local)2.2 闭包作用域(Enclos

Python实现对阿里云OSS对象存储的操作详解

《Python实现对阿里云OSS对象存储的操作详解》这篇文章主要为大家详细介绍了Python实现对阿里云OSS对象存储的操作相关知识,包括连接,上传,下载,列举等功能,感兴趣的小伙伴可以了解下... 目录一、直接使用代码二、详细使用1. 环境准备2. 初始化配置3. bucket配置创建4. 文件上传到os

使用Python实现可恢复式多线程下载器

《使用Python实现可恢复式多线程下载器》在数字时代,大文件下载已成为日常操作,本文将手把手教你用Python打造专业级下载器,实现断点续传,多线程加速,速度限制等功能,感兴趣的小伙伴可以了解下... 目录一、智能续传:从崩溃边缘抢救进度二、多线程加速:榨干网络带宽三、速度控制:做网络的好邻居四、终端交互

Python中注释使用方法举例详解

《Python中注释使用方法举例详解》在Python编程语言中注释是必不可少的一部分,它有助于提高代码的可读性和维护性,:本文主要介绍Python中注释使用方法的相关资料,需要的朋友可以参考下... 目录一、前言二、什么是注释?示例:三、单行注释语法:以 China编程# 开头,后面的内容为注释内容示例:示例:四

Python中win32包的安装及常见用途介绍

《Python中win32包的安装及常见用途介绍》在Windows环境下,PythonWin32模块通常随Python安装包一起安装,:本文主要介绍Python中win32包的安装及常见用途的相关... 目录前言主要组件安装方法常见用途1. 操作Windows注册表2. 操作Windows服务3. 窗口操作

Python中re模块结合正则表达式的实际应用案例

《Python中re模块结合正则表达式的实际应用案例》Python中的re模块是用于处理正则表达式的强大工具,正则表达式是一种用来匹配字符串的模式,它可以在文本中搜索和匹配特定的字符串模式,这篇文章主... 目录前言re模块常用函数一、查看文本中是否包含 A 或 B 字符串二、替换多个关键词为统一格式三、提