Wilcoxon符号秩检验
2024-06-01 14:05:33
Wilcoxon符号秩检验
它适用于T检验中的成对比较，但并不要求成对数据之差服从正态分布，只要求对称分布即可。检验成对观测数据之差是否来自均值为0的总体(产生数据的总体是否具有相同的均值)。
在Matlab中，秩和检验由函数ranksum实现。
命令为： [p,h]=ranksum(x,y,alpha)
其中x，y可为不等长向量，alpha为给定的显著水平，它必须为0和1之间的数量。p返回 产生两独立样本的总体是否相同的显著性概率，h返回假设检验的结果。如果x和y的总 体差别不显著，则h为零；如果x和y的总体差别显著，则h为1。如果p接近于零，则可对 原假设质疑。
clear;
clc;
path1 = 'C:\Users\Administrator\Desktop\重新整理血管网络\真实数据.xls';
path2 = 'C:\Users\Administrator\Desktop\重新整理血管网络\真实数据1.xls';
path3 = 'C:\Users\Administrator\Desktop\重新整理血管网络\生成数据9.xls';
path4 = 'C:\Users\Administrator\Desktop\重新整理血管网络\生成数据1.xls';
path5 = 'C:\Users\Administrator\Desktop\重新整理血管网络\wgan-gp1.xls';
path6 = 'C:\Users\Administrator\Desktop\重新整理血管网络\wgan-gp200.xls';mypath1 = path1;%对比两个结果
mypath2 = path2;
A = 'A2:A500';
B = 'B2:B500';
C = 'C2:C500';
D = 'D2:D500';
E = 'E2:E500';
F = 'F2:F500';
G = 'G2:G500';
H = 'H2:H500';
I = 'I2:I500';
J = 'J2:J500';
K = 'K2:K500';
L = 'L2:L500';
M = 'M2:M500';
N = 'N2:N500';
O = 'O2:O500';
disp('主分支长度')
[num1] = xlsread(mypath1,A);
[num2] = xlsread(mypath2,A);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('左分支长度')
[num1] = xlsread(mypath1,B);
[num2] = xlsread(mypath2,B);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('右分支长度')
[num1] = xlsread(mypath1,C);
[num2] = xlsread(mypath2,C);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('主左长度比')
[num1] = xlsread(mypath1,D);
[num2] = xlsread(mypath2,D);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('主右长度比')
[num1] = xlsread(mypath1,E);
[num2] = xlsread(mypath2,E);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('主分叉卷曲度')
[num1] = xlsread(mypath1,F);
[num2] = xlsread(mypath2,F);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('左分叉卷曲度')
[num1] = xlsread(mypath1,G);
[num2] = xlsread(mypath2,G);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('右分叉卷曲度')
[num1] = xlsread(mypath1,H);
[num2] = xlsread(mypath2,H);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('主左角度')
[num1] = xlsread(mypath1,I);
[num2] = xlsread(mypath2,I);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('主右角度')
[num1] = xlsread(mypath1,J);
[num2] = xlsread(mypath2,J);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('左右角度')
[num1] = xlsread(mypath1,K);
[num2] = xlsread(mypath2,K);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('主端点长度')
[num1] = xlsread(mypath1,L);
[num2] = xlsread(mypath2,L);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p);disp('左端点长度')
[num1] = xlsread(mypath1,M);
[num2] = xlsread(mypath2,M);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p)disp('右端点长度')
[num1] = xlsread(mypath1,N);
[num2] = xlsread(mypath2,N);
mean1 = mean(num1); %平均值
mean2 = mean(num2);
std1 = std(num1); %标准差
std2 = std(num2);
var1 = var(num1); %方差
var2 = var(num2);
[p,h] = ranksum(num1,num2);
fprintf('原始数据1--平均值:%f,标准差:%f,方差:%f\n',mean1,std1,var1);
fprintf('原始数据2--平均值:%f,标准差:%f,方差:%f\n',mean2,std2,var2);
fprintf('h值:%d,p值:%f\n\n',h,p)Wilcoxon符号秩检验相关推荐
- R语言编写自定义函数使用Wilcoxon符号秩检验（Wilcoxon signed rank）实现多分组非参数成对检验（pairwise）、并使用p.adjust函数调整概率值
R语言编写自定义函数使用Wilcoxon符号秩检验(Wilcoxon signed rank)实现多分组非参数成对检验(Nonparametric pairwise multiple comparis ...
- Wilcoxon 符号秩检验（python-matlab-spss）
Wilcoxon 符号秩检验 威尔科克森符号秩检验亦称威尔科克伦代符号的等级检验,是由威尔科克森(F·Wilcoxon)于1945年提出的.该方法是在成对观测数据的符号检验基础上发展起来的,比传统的单 ...
- R语言中的Wilcoxon符号秩检验与配对学生t检验
在这篇文章中,我们将探索比较两组依赖(即配对)定量数据的检验:Wilcoxon符号秩检验和配对学生t检验.这些检验之间的关键区别在于Wilcoxon的检验是非参数检验,而t检验是参数检验.在下文中,我 ...
- wilcoxon符号秩检验matlab,Wilcoxon符号秩统计量的性质.ppt
Wilcoxon符号秩统计量的性质.ppt 第三章 单一样本的推断问题 主要内容 第一节 符号检验和分位数推断 例3.1. 假设某地16座预出售的楼盘均价,单位(百元/平方米)如下表所示: 36 32 ...
- Wilcoxon符号秩检验详解
Wilcoxon符号秩检验详解 Wilcoxon符号秩检验是一种常用的非参检验方法,用于比较两个相关样本之间的差异.本文将对Wilcoxon符号秩检验进行详细介绍. 基本原理 Wilcoxon符号秩检 ...
- 非参数检验之符号检验、Wilcoxon符号秩检验、游程检验
目录 一.符号检验 例2.1下面是世界上71个大城市的花费指数(包括租金)按递增次序排列如下(这里上海是44位,其指数为63.5): R代码: 二.Wilcoxon符号秩检验 例2.3下面是10个欧洲 ...
- t检验、Wilcoxon 符号秩检验、Mann-Whitney 秩和检验 Python实现
@创建于:20211216 @修改于:20211216 文章目录 1 理论知识 2 代码实现 3 统计库函数 1 理论知识 知乎:一文详解t检验 知乎:t检验还是Mann-Whitney U检验 电子 ...
- 【非参数统计02】单一样本的位置推断问题：符号检验、分位数推断、Cox-Staut趋势存在性检验、随机游程检验、Wilcoxon符号秩检验、位置参数区间估计、正态记分、分布一致性检验、稳健性评价
目录导引 2 单一样本的位置推断 2.1 符号检验 2.1.2 大样本场合 2.1.3 配对样本 2.1.4 分位数检验 2.2 Cox-Staut 趋势存在性检验 2.2.1 最优权重 2.2.2 ...
- Wilcoxon符号秩检验-MATLAB实现
Wilcoxon符号秩检验 非参数检验 非参数意味着总体数据没有正态分布 如果数据对之间的差异是非正态分布的,则应使用Wilcoxon Signed-Rank检验 两个稍微不同的测试版本 Wilcox ...
最新文章
- centos 安装 NTFS支持
- SpringBoot第十二篇：springboot集成apidoc
- css折叠样式（1）——使用css样式的三种方式
- Linux防火墙配置—访问外网WEB
- 【EventBus】EventBus 源码解析 ( EventBus 构建 | EventBus 单例获取 | EventBus 构造函数 | EventBus 构建者 )
- NOIP2007 树网的核 [BZOJ2282][Sdoi2011]消防
- 用几条shell命令快速去重10G数据
- CSS text-indent 属性
- macOS的Fn键实用技巧
- WIN10更新后超慢问题的解决
- 《大数据技术原理与应用》—— 提纲
- 计算机毕业设计springboot+vue餐厅点餐系统源码
- GTK使用cairo绘图教程
- 公司组织管理与权限管理的设计原则
- ios系统越狱才可以享受的100项操作
- Revit插件中的“喷淋对齐”“链接CAD”功能操作
- 谈判技巧——执行合同
- NCBI|宏基因组原始数据上传
- 【Luogu5348】密码解锁（莫比乌斯反演，数论）
- Lect5_Model_free_Control