目录

  • 1 单样本和两样本均值的t检验
  • 2 单样本和两样本关于中位数的非参数检验
  • 3 拟合优度的X2X^2X2检验
  • 4 列联表的X2X^2X2检验及Fisher检验

1 单样本和两样本均值的t检验

(1)单样本

举例:H0:μ=3.1⇌H1:μ≠3.1H_0:\mu= 3.1 \rightleftharpoons H_1:\mu \neq 3.1H0​:μ=3.1⇌H1​:μ​=3.1

import scipy.stats as stats
import numpy as npnp.random.seed(1010)
x = np.random.normal(3,1,500)
stat,p_value = stats.ttest_1samp(x,3,1)
print(stat,p_value)

(2)两样本

举例:H0:μ1=μ2⇌H1:μ≠μ2H_0:\mu_ 1 = \mu_2 \rightleftharpoons H_1:\mu \neq \mu_2H0​:μ1​=μ2​⇌H1​:μ​=μ2​

import scipy.stats as stats
import numpy as npnp.random.seed(1010)
x = np.random.normal(3,1,500)
y = np.random.normal(3.2,1,500)
stat,p_value = stats.ttest_ind(x,y,equal_var=False)#不假定方差相等
print(stat,p_value)

2 单样本和两样本关于中位数的非参数检验

(1)单样本之Wilcoxon符号秩检验

H0:M=3.1⇌H1:M≠3.1H_0:M = 3.1 \rightleftharpoons H_1:M \neq 3.1H0​:M=3.1⇌H1​:M​=3.1

import scipy.stats as stats
import numpy as npnp.random.seed(1010)
x = np.random.normal(3,1,500)
stat,p_value = stats.wilcoxon(x-3.1,correction=True)
print(stat,p_value)

(2)两样本之Wilcoxon秩和检验

H0:M1=M2⇌H1:M1<M2H_0:M_1 = M_2 \rightleftharpoons H_1:M_1 < M_2H0​:M1​=M2​⇌H1​:M1​<M2​

import scipy.stats as stats
import numpy as npnp.random.seed(1010)
x = np.random.normal(3,1,500)
y = np.random.normal(3.2,1,500)stat,p_value = stats.mannwhitneyu(x,y,alternative='less')
print(stat,p_value)

3 拟合优度的X2X^2X2检验

(1)X2X^2X2检验

拟合优度检验是检验一组样本观测值是不是来自一个已知分布的总体。拟合优度的x2x^2x2检验的原理为:根据假设的分布可以在一些任意划分的k个范围计算出期望值E1,E2,...,EkE_1,E2,...,E_kE1​,E2,...,Ek​,而根据在这些区域所得到的观测值O1,O2,...,OkO_1,O_2,...,O_kO1​,O2​,...,Ok​可以得到检验统计量X2=∼i=1k(Oi−Ei)2EiX^2 = \sim_{i=1}^k\frac{(O_i-E_i)^2}{E_i}X2=∼i=1k​Ei​(Oi​−Ei​)2​,渐渐服从x2x^2x2分布。

import scipy.stats as stats
import numpy as np# 初始化举例的数据
size = 5000
x= stats.poisson.rvs(10,size = size,random_state=1010)
fig = plt.figure(figsize=(12,5))
H = plt.hist(x,20,color='white','hatch='/')
E = np.diff(stats.poisson.cdf(H[1],10))*size
h_0 = H[0]
# 卡方检验
chisq,p = stats.chisquare(f_obs=h_0,f_exp=E)
print(chisq,p)

拟合优度检验还可以用Kolmogorov-Smirnov检验

stat,p_value = stats.kstest(x,lambda x:stats.poisson.cdf(x,10))
print(stat,p_value)

4 列联表的X2X^2X2检验及Fisher检验

(1)Pearson X2X^2X2检验检验

import scipy.stats as stats
import numpy as np#列联表
X =np.array([[10,20,4],[20,9,8]])#2×3的列联表chi2,p,df,exp = stats.chi2_contingency(X)
print(p,df)

(2)似然比X2X^2X2检验

import scipy.stats as stats
import numpy as np#列联表
X =np.array([[10,20,4],[20,9,8]])#2×3的列联表chi2,p,df,exp = stats.chi2_contingency(X,lambda_='log_likelihood')
print(p,df)

(3)Fisher X2X^2X2检验

import scipy.stats as stats
import numpy as np#列联表
X1 =np.array([[10,20],[20,9]])#2×2的列联表
oddsratio,p = stats.fisher_exact(X1)
print(oddsratio,p)

【python 数据分析】不同情况下的t检验、Wilcoxon符号秩检验、Wilcoxon秩和检验、卡方检验、Fisher检验相关推荐

  1. 如何打开python的交互窗口-Python多版本情况下四种快速进入交互式命令行的操作技巧...

    原标题:Python多版本情况下四种快速进入交互式命令行的操作技巧 因为工作需求或者学习需要等原因,部分小伙伴的电脑中同时安装了Python2和Python3,相信在Python多版本的切换中常常会遇 ...

  2. python实现双屏情况下,第二个屏播发视频

    python实现双屏情况下,第二个屏播发视频 业务情况 环境准备 业务情况 需要在系统启动的情况下,在第二个显示屏上播发MP4的视频, 环境准备 安装QQ影音 将video.mp4设置成QQ影音默认播 ...

  3. R语言中的Wilcoxon符号秩检验与配对学生t检验

    在这篇文章中,我们将探索比较两组依赖(即配对)定量数据的检验:Wilcoxon符号秩检验和配对学生t检验.这些检验之间的关键区别在于Wilcoxon的检验是非参数检验,而t检验是参数检验.在下文中,我 ...

  4. python logging默认情况下打印_python logging日志打印过程解析

    一. 基础使用 1.1 logging使用场景 日志是什么?这个不用多解释.百分之九十的程序都需要提供日志功能.Python内置的logging模块,为我们提供了现成的高效好用的日志解决方案.但是,不 ...

  5. Python数据分析笔记(下)

    数据的描述 1.数据的集中趋势 (1)平均值 import pandas as pd # 1.算数平均值 cars_score=pd.read_csv(r'D:\Projects\Python\Doi ...

  6. python交互式窗口怎么换行_Python多版本情况下四种快速进入交互式命令行的操作技巧

    原标题:Python多版本情况下四种快速进入交互式命令行的操作技巧 因为工作需求或者学习需要等原因,部分小伙伴的电脑中同时安装了Python2和Python3,相信在Python多版本的切换中常常会遇 ...

  7. [译] Web 爬虫下的 Python 数据分析:中情局全球概况图解

    原文地址:Data Analytics with Python by Web scraping: Illustration with CIA World Factbook 原文作者:Tirthajyo ...

  8. 我是什么情况下开始学python,现在能用python做什么

    目录 一.首先一对一回答题主的问题 二.详细说说我和python的缘分始终 2.1.爬虫篇 2.2.数据分析篇 2.3.数据可视化篇 2.4.web开发篇 2.5.总结 三.关于python副业兼职 ...

  9. 《Python数据分析与挖掘实战》第10章(下)——DNN2 筛选得“候选洗浴事件”3 构建模型

    本文是基于<Python数据分析与挖掘实战>的实战部分的第10章的数据--<家用电器用户行为分析与事件识别> 做的分析. 接着前一篇文章的内容,本篇博文重点是处理用水事件中的属 ...

  10. 干货 | 50题带你入门Python数据分析(下)

    25 显示列唯一值 df[''].unique() ⚠️这里显示具体的唯一值 26 显示列唯一值个数 df[''].nunique() ⚠️这里显示的是唯一值个数 27 返回列最大/小值的索引 df[ ...

最新文章

  1. NK3C框架(MyBatis、Durid)连接SQL Server
  2. android studio 首字母提示 设置 大小写敏感
  3. SpringCloud教程-注册中心(Consul)(SpringCloud版本Greenwich.SR4)
  4. windows下安装RabbitMQ消息服务器 + 读写队列
  5. Stream流中的常用方法_skip
  6. ajax 接收json数据的进一步了解
  7. pl sql 连接mysql_PL/SQL 连接mysql步骤
  8. 【AMAD】django-activity-stream
  9. 小程序组件库开发之抽奖游戏组件
  10. 记录——《C Primer Plus (第五版)》第十章编程练习第四题
  11. php sockets有什么用,PHP中Sockets与流有什么关系啊!!!!!
  12. 推荐一个Oracle数据库学习网站
  13. 【Matlab】求解函数导数
  14. QQ出现大规模盗号,qq被盗发布不良信息怎么办
  15. Win10 + Ubuntu20.04 双系统+双硬盘安装
  16. 多旋翼无人机ROSC++开发例程(一):环境配置
  17. rgb加色模式与cmyk减色模式
  18. Maya 2013 下载以及破解教程
  19. NAR | 深圳先进院马迎飞组发布底盘噬菌体的高通量制备方法
  20. 95后创业风口已来临?牛逼又低调的95后们都在干什么

热门文章

  1. JWT全面解读、使用步骤
  2. macd金叉股共振_股票MACD金叉什么意思
  3. 徐思 201771010132
  4. RF- BuiltIn_Run Keyword关键字系列
  5. 偏差-方差分解简要推导
  6. 网站设计| 网站设计模板素材
  7. Racket Cheat Sheet
  8. Exp4 恶意代码分析 20154328 常城
  9. 诸神战纪堕落天使java_诸神战纪-堕落天使路西法
  10. C#实现office文档转换为PDF或xps