当前位置: 首页 > news >正文

python统计分析——单样本均值检验

参考资料:python统计分析【托马斯】

1、单样本均值的t检验

        检验一个正态分布数据的均值和一个参考值的差异,我们一般使用单样本t检验,该检验基于t分布。

        如果我们知道一个正态分布总体的均值和标准差,那么我们可以计算对应的标准误,并用来自正态分布的值来判断得到一个特定值的可能性大小。然而,在实际中,我们必须从样本中估计均值和标准差;并且用来描述来自正态分布的样本均值的t分布。

        我们从均值为7、标准差为3的正态分布中抽取100个数据。

# 导入库
import numpy as np
from scipy import stats
import random# 设置正态分布
normal=stats.norm(7,3)
# 设置随机种子,用于复现结果
random.seed(1234)
# 从正态分布中抽取100个数据
sample=normal.rvs(100)# 输出样本均值和标准差
print("样本均值:",sample.mean())
print("样本标准差:",sample.std(ddof=1))

检验样本均值6.83是否与总体均值7是否有显著差异,使用如下函数:

t,pVal=stats.ttest_1sample(data,checkVale)

t,pVal=stats.ttest_1samp(sample,7)
print("t统计量:",t)
print("p值",pVal)

结论:p值大于0.05,样本均值与总体均值差异不显著。

可参考:python统计分析——单样本t检验_单样本t分布检验python代码-CSDN博客

2、单样本均值的Wilcoxon符号秩和检验

        在非正态分布的数据中,我们不应该使用单样本t检验。相反我们必须使用均值的非参数检验方法。我们可以进行Wilcoxon符号秩和检验。

wil,pVal=stats.wilcoxon(data-checkValue)

该方法有3个步骤:

        ①计算每个观测值和感兴趣的值的差异;

        ②忽略差异的符号,将它们按照大小排序;

        ③将所有负(或正)秩次的秩次加起来,也就是哪些低于(或高于)选定的假设值的秩次。

案例参考:

# 输入案例数据
sample=np.array([5260,5470,5640,6180,6380,6515,6805,7515,7515,8230,8770])
# 进行Wilcoxon秩和检验
wil,pVal=stats.wilcoxon(sample-7515)
print("统计量:",wil)
print("p值:",pVal)

结论:样本的中位数与7515差异不显著。

相关文章:

  • Linux根据时间删除文件或目录
  • nginx 设置跨域
  • C# 多线程编程:线程锁与无锁并发
  • Qt常用容器之:QVector
  • 【vue核心技术实战精讲】1.9 Vue指令之v-model双向数据绑定
  • 华为云使用指南02
  • vue 列表渲染
  • k8s 如何获取加入节点命名
  • 浅谈iOS开发中的自动引用计数ARC
  • 使用llamafile 构建本地大模型运用
  • spring boot的返回值里面含有net.sf.json.JSONObject 报错net.sf.json.JSONNull[“empty“])]
  • <深度学习入门学习笔记P1>——《深度学习》
  • 后端返回文件流pdf 下载
  • 【AIGC调研系列】Starling-LM-7B模型与其他模型相比的优势和劣势
  • 更高效稳定 | 基于ACM32 MCU的编程直流电源应用方案
  • 【跃迁之路】【735天】程序员高效学习方法论探索系列(实验阶段492-2019.2.25)...
  • 4. 路由到控制器 - Laravel从零开始教程
  • fetch 从初识到应用
  • Java IO学习笔记一
  • JavaScript新鲜事·第5期
  • Markdown 语法简单说明
  • oschina
  • Python打包系统简单入门
  • rc-form之最单纯情况
  • Spring Cloud(3) - 服务治理: Spring Cloud Eureka
  • vagrant 添加本地 box 安装 laravel homestead
  • 大整数乘法-表格法
  • 技术胖1-4季视频复习— (看视频笔记)
  • 你不可错过的前端面试题(一)
  • 如何邀请好友注册您的网站(模拟百度网盘)
  • 实习面试笔记
  • 试着探索高并发下的系统架构面貌
  • ​configparser --- 配置文件解析器​
  • #mysql 8.0 踩坑日记
  • (2.2w字)前端单元测试之Jest详解篇
  • (3)nginx 配置(nginx.conf)
  • (poj1.3.2)1791(构造法模拟)
  • (附源码)springboot码头作业管理系统 毕业设计 341654
  • (六)软件测试分工
  • (七)理解angular中的module和injector,即依赖注入
  • (未解决)macOS matplotlib 中文是方框
  • (一)kafka实战——kafka源码编译启动
  • (一)硬件制作--从零开始自制linux掌上电脑(F1C200S) <嵌入式项目>
  • ****** 二 ******、软设笔记【数据结构】-KMP算法、树、二叉树
  • .bat批处理出现中文乱码的情况
  • .cn根服务器被攻击之后
  • .NET Core 版本不支持的问题
  • .NET Core日志内容详解,详解不同日志级别的区别和有关日志记录的实用工具和第三方库详解与示例
  • .NET gRPC 和RESTful简单对比
  • .NET NPOI导出Excel详解
  • .Net程序猿乐Android发展---(10)框架布局FrameLayout
  • .net打印*三角形
  • .net开发时的诡异问题,button的onclick事件无效
  • .net利用SQLBulkCopy进行数据库之间的大批量数据传递
  • .sh文件怎么运行_创建优化的Go镜像文件以及踩过的坑