目录

1、Mann-whitney 算法简介

2、定义

3、Mann-whitney 算法步骤

4、matlab函数

5、实例及matlab代码


独立双样本的非参数检验,不满足正态分布的小样本,秩和检验 X Y样本数量可以不相等

参考链接:https://blog.csdn.net/qq_34734303/article/details/80296316

https://blog.csdn.net/lhyer/article/details/50524174

1、Mann-whitney 算法简介

曼-惠特尼U检验又称“曼-惠特尼秩和检验”,是由H.B.Mann和D.R.Whitney于1947年提出的 [1] 。它假设两个样本分别来自除了总体均值以外完全相同的两个总体,目的是检验这两个总体的均值是否有显著的差别。

2、定义

曼-惠特尼U检验(Mann-Whitney U test)曼-惠特尼秩和检验可以看作是对两均值之差的参数检验方式的独立双样本T检验或相应的大样本正态检验的代用品。由于曼-惠特尼秩和检验明确地考虑了每一个样本中各测定值所排的秩,它比符号检验法使用了更多的信息。

3、Mann-whitney 算法步骤

具体步骤如下:

第一步: 将两组数据混合,并按照大小顺序编排等级。最小的数据等级为 1,第二小的数据 等级为 2,以此类推(若有数据相等的情形,则取这几个数据排序的平均值作为其等级)。

{1, 3, 3, 7}对应等级{1, 2.5, 2.5, 4}

第二步: 分别求出两个样本的等级和 W1、 W 2。

第三步: 计算曼 -惠特尼 U 检验统计量, n1 为第一个样本的量, n2 为第二个样本的量:

N=n1+n2

W1+W2=N(N+1)/2

So U1+U2=n1n2

选择 U1 和 U2 中最小者与临界值 Ua比较,当 U < Ua 时,拒绝 H0,接受 H1。在原假设为真的情况下, 随机变量 U 的均值和方差分别为,当 n1 和 n2都不小于 10 时,随机变量近似服从正态分布:

第四步: 作出判断。设第一个总体的均值为 μ 1,第二个总体的均值为 μ 2,则有:

4、matlab函数

MATLAB统计工具箱中提供了ranksum函数,用来做秩和检验,其调用格式如下:

 [p,h,stats]=ranksum(x,y,param1,val1,........)

根据样本观测值向量x和y做双侧秩和检验,原假设两独立样本x和y来自于具有相同中位数的连续分布,备择假设是x和y具有不同的中位数。此时可选成对出现的参数和参数值来控制计算结果。可选的参数名与参数值如下表:

参数名                  参数值及说明

'alpha'                检验的显著性水平,其取值介于0--1

默认值为0.05

‘method’         指定计算p值的方法,可能的取值情况如下

‘exact’:利用精确方法计算p值,适用于小样本(样本容量<10)情形

‘approximate’:利用正态近似计算p值,适用于大样本情形

输出参数分别为检验的p值、变量h和包含检验统计量信息的结构体变量stats,当p>a(显著性水平)或h=0时,接受原假设;当p<=a或h=1时,拒绝原假设。

5、实例及matlab代码

下面是两种不同加工方式的菜粕在黄牛瘤胃内培养16h的干物质降解率,用曼-惠特尼U检验比较其有无差异:

两种加工方式的菜粕瘤胃培养 16h的干物质降解率(%)

先按照大小顺序排列等级(见上表),而后计算 W1 = 38, W2 = 67, n1 = 6, n2 = 8。

假设两种菜粕的16h瘤胃干物质降解率除了平均水平以外在其它方面无差异,即检验:

H0:两种菜粕的16h瘤胃干物质降解率无差异;

H1:两种菜粕的16h瘤胃干物质降解率有差异。

计算U值:

%第一组
x=[39.33,44.10,35.89,43.35,47.61,43.71];
%第二组
y=[42.91,44.69,44.54,45.31,37.73,48.75,46.71,41.85];%调用ranksum进行检验
[p,h,stats]=ranksum(x,y,'method','exact')

p =

0.4136

h =

0

函数返回p=0.4136>0.05,U2值较小,选取 U2与 Uα(α=0.05)比较,通过查表(附表)可知 Uα = 8, U2 > Uα,即接受 H0,认为两种加工方式的菜粕瘤胃培养16h的干物质降解率无显著差异

Mann-Whitney检验(曼-惠特尼秩和检验)及matlab代码相关推荐

  1. 曼惠特尼U检验与特征筛选

    曼惠特尼U检验与特征筛选 曼-惠特尼U检验(Mann-Whitney U test)曼-惠特尼秩和检验可以看作是对两均值之差的参数检验方式的T检验或相应的大样本正态检验的代用品.由于曼-惠特尼秩和检验 ...

  2. 曼-惠特尼U检验Mann–Whitney U Test

    sklearn实战-乳腺癌细胞数据挖掘(博主亲自录制视频教程) https://study.163.com/course/introduction.htm?courseId=1005269003&am ...

  3. 曼-惠特尼U检验与威尔科克森符号秩检验

    文章目录 1. Mann–Whitney *U*-test (曼-惠特尼U检验) 1.1 曼-惠特尼U检验定义 1.2 曼-惠特尼U检验实现 2. Wilcoxon Signed-Rank Test( ...

  4. 【SPSS】两独立样本的曼-惠特尼U检验详细操作教程(附案例实战)

  5. 【杂纪】从ROC曲线到AUC值,再到Mann–Whitney U统计量

    统计检验中的两类错误 在进行假设检验时,分别提出原假设(Null Hypothesis)和备择假设(Alternative Hypothesis),检验结果可能出现的两类错误: 原假设实际上是正确的, ...

  6. AUC(一):AUC与Mann–Whitney U test

    在算法面试中,常常会被一个问题AUC的物理含义是什么.本文参考论文1,介绍一下AUC的物理含义. 概念 在二分类模型中,预测值ppp表示事件发生的概率.对于分类任务,需要给出一个阈值(threshol ...

  7. stata代码命令全集:跟检验do文件、stata常用命令、代码大合集、面板熵值法

     一.常用模型代码整理 1.数据来源: 2.时间跨度:无 3.区域范围:无 4.指标说明: 包含如下模型代码: l OLS模型 l Heckman两阶段模型 l PSM+DID模型 l 固定效应模型( ...

  8. 基于霍夫曼(Huffman)图像编码的图像压缩和重建-含Matlab代码

    目录 一.引言 二.霍夫曼Huffman编码 2.1 霍夫曼编码流程 2.2 输入数据的编码 三.霍夫曼解码 四.实验结果 五.参考文献 六.Matlab代码(GUI界面)获取 一.引言 随着通信与信 ...

  9. 格子玻尔兹曼流体代码_格子玻尔兹曼方法(LBM)学习:对流-扩散问题(附MATLAB代码)...

    (๑❛ᴗ❛๑) 麻烦各位读者收藏之余点个喜欢或赞呢,咱也更有干劲了~ OrzSunspot:格子玻尔兹曼方法(LBM)学习:等温不可压缩流体流动问题(附MATLAB代码)​zhuanlan.zhihu ...

最新文章

  1. Git相关二三事(git reflog 和彩色branch)【转】
  2. C++多线程的简单程序
  3. 27/100. Intersection of Two Linked Lists
  4. 谷歌在华遭遇首例关键词官司
  5. 2018数据技术嘉年华-金融峰会·重庆站即将起航!
  6. 【Kafka】kafka Failed to acquire lock on file .lock in /data/kafak-logs a kafka instance
  7. 业务线开发流程图(四)
  8. 集成灶带给我的是无尽烦恼,大家的集成灶用得如何?
  9. Ubuntu 安装 中文输入法(Google 拼音)
  10. 一个文科小白的数据分析师之路
  11. scrapy[skp]快速入门
  12. 远程连接SQL server 数据库
  13. CSS: 前端的视差滚动,写出看起来高逼格的网页效果
  14. linux权限777什么意思,chmod 权限777是什么意思
  15. EditText焦点
  16. 阶段1 - 03. 常用API第一部分 - 04. ArrayList集合
  17. from PyQt6 import QtCore, QtGui, QtWidgets ImportError: DLL load failed while importing QtGui:
  18. ASEMI代理ADI(亚德诺)AD5934YRSZ-REEL7车规级芯片
  19. Git教程 | (9) 自定义Git和使用SourceTree
  20. 微信小程序开发者工具-使用技巧(更新中~)

热门文章

  1. 最好的医生是自己:把吃出来的病吃回去
  2. 无线网手动添加服务器,无线网络手动设置的问题解决v
  3. android检测usb存储,USBdevice识别为存储设备并找到路径
  4. Windows Server 2012 R2部署OWA服务(三)安装部署OWA服务
  5. hjr-MUD游戏(三):android原生游戏编写
  6. 记一个简单的Android计时器制作过程。
  7. 大谷新作,AI复原泰戈尔:世界以痛吻我,要我报之以歌
  8. MD5中使用16进制消息摘要
  9. float在python中的意思_Python中的float()
  10. linux下进行Oracle自动备份配置