秩和比算法matlab程序,Matlab学习系32. 秩和比综合评价法.docx
32. 秩和比综合评价法
一、简单介绍
秩和比法,是我国统计学家田凤调教授于1988年提出的一种综合评价方法,是利用秩和比(RSR, Rank-sum ratio)进行统计分析的一种方法,该方法在医疗卫生等领域的多指标综合评价、统计预测预报、统计质量控制等方面已得到广泛的应用。
基本原理是一种将样本多项指标值通过秩变换,得到无量纲统计量统计量RSR(介于0~1之间),再运用参数统计分析方法研究RSR的分布。不论所分析的问题是什么,计算的RSR越大越好,以RSR值对评价对象的优劣直接排序或分档排序,从而对评价对象做出综合评价。
优点:是非参数统计分析,对指标的选择无特殊要求,适于各种评价对象;由于计算用的数值是秩次,可以消除异常值的干扰,它融合了参数分析的方法,结果比单纯采用非参数法更为精确,既可以直接排序,又可以分档排序,使用范围广泛。
缺点:是排序的主要依据是利用原始数据的秩次,最终算得的RSR值反映的是综合秩次的差距,而与原始数据的顺位间的差距程度大小无关,这样在指标转化为秩次是会失去一些原始数据的信息,如原始数据的大小差别等。
二、算法步骤
1. 编秩
设有n个评价对象,m个评价指标的样本数据(n行m列),分别对每个指标列的数据编秩:正向指标(值越大越好)从小到大编秩,负向指标(值越小越好)从大到小编秩,值相同的数据编平均秩。得到秩矩阵R=(rij)n×m
注:编秩即对数据排序,其顺序号作为秩。
例如,有5人的数学成绩:甲80分, 乙73分, 丙65分, 丁92分, 戊73分。成绩是正向指标,从小到大排序:丙
编秩:丙=1;甲=4;丁=5;排第2和第3位的乙和戊成绩相同取平均秩:乙=戊=(2+3)/2=2.5。
2. 计算秩和比RSR(每个指标权重相同=1/m)
若各评价指标的权重不同,则计算或加权秩和比WRSR
其中,wj为第j个评价指标的权重。
3. 计算概率单位
编制RSR/WRSR频率分布表:将RSRi或WRSRi从小到大排列(值相同作为一组,或自定义分组),列出各组频数fi和各组累计频数∑fi;计算累积频率;再将pi转换为概率单位Probiti,其中,Probiti取标准正态分布的pi分位数+5.
4. 计算直线回归方程
以累积频率所对应的概率单位Probiti为自变量,RSRi或WRSRi值为因变量,做线性回归:RSR/WRSR=a+b*Probit
5. 分档排序
按回归方程计算的RSR/WRSR估计值,对评价对象进行分档排序。
三、Matlab实现
例1 某市人民医院1983-1992年工作质量统计指标数据,如下表所示:
其中,病死率x2, 平均住院日x6为负向指标。
六项指标的权重依次为:0.0930.4180.1320.1000.0980.159
datas=xlsread('RSR.xlsx');
X=datas(:,2:end);
w=[0.093 0.418 0.132 0.100 0.098 0.159];
X(:,[2 6])=-X(:,[2 6]); %负向指标转换为正向指标
R=tiedrank(X); %对X的各列分别编秩
[n,m]=size(R);
W=repmat(w,n,1);
WRSR=sum(W.*R,2)/n; %计算加权秩和比: 加权,按行求和,再除以n
freq=tabulate(WRSR); %统计WRSR的频数,频率, freq的第3列为频率
p=cumsum(freq(:,3))/100; %计算累积频率
p(end)=p(end)-1/(4*n); %修正最后一个累积频率
Probit=norminv(p,0,1)+5; %计算标准正态分布的p分位数+5
Probit=[ones(n,1), Probit, Probit.^2, Probit.^3];
[b,bint,r,rint,stats]=regress(WRSR,Probit);
%三次多项式回归, 一次或二次多项式回归R方和p值较差
b
stats
WRSRfit=Probit*b %计算WRSR的估计值
[s,ind]=sort(WRSRfit,'descend') %对WRSR估计值从大到小排序
运行结果:
b = 10.7950 -6.4953 1.3137 -0.0851
stats = 0.7932 7.6735 0.0178 0.0066
WRSRfit = 0.4335
0.3851
0.4133
0.4655
0.5280
0.5953
0.6634
0.7243
0.7463
0.5454
s =
秩和比算法matlab程序,Matlab学习系32. 秩和比综合评价法.docx相关推荐
- 影像重采样matlab程序,matlab重采样代码
本设计通过编 写 matlab 程序,借助对声音音色和音调的双重复合改变,实现输出声音的改变.共振 峰频 率的改变是基 于重采样 实现的,从重 采样原理 知道,这也同...... 本设计通过编 写 m ...
- 标号法matlab程序,matlab程序大全
上例 Matlab 程序如下: %使用表8.4 给出的小波基 ld=[0 0.0... matlab经典代码大全_计算机软件及应用_IT/计算机_专业资料.哈哈哈 MATLAB 显示正炫余炫图:plo ...
- 投影仪标定matlab程序,matlab投影仪标定
matlab投影仪标定 matlab中的投影仪标定方法属于投射标定图案法,投射的图案只能是棋盘格,该方法将投影仪当作逆向的相机,它具有映射图像二维点到三维坐标的功能,因此可以使用张正友标定方法来对它进 ...
- rbpf粒子滤波slam matlab程序_学习笔记(优达学城)- 车辆定位之粒子滤波器(整合版)...
1.代码传送门 首先,一如既往的,打开传送门! Fred159/CarND-Kidnapped-Vehicle-Projectgithub.com 代码, 很重要,但更重要的是从代码的行与行之间探索 ...
- vof模型matlab程序,Fluent学习笔记(11)----VOF模型
喷水现象,属于流体力学中的射流问题.fluent来对该问题进行研究. 实例介绍: 尺寸:edge4=0.8m,edge1=0.09m,edge2=0.02m,edge5=0.2m, 1. gambit ...
- 小世界网络matlab程序,matlab的dijkstra算法问题(小世界网络)
本帖最后由 枫叶舞流星 于 2020-12-19 21:18 编辑 运行WS_net.m文件 得到小世界网络图的平均路径长度为0,不知道哪里错了,求大神指点 好像是dijkstra.m文件的错误,di ...
- 哈密顿圈 MATLAB程序,matlab求最佳哈密顿圈遇到错误
问题描述: matlab求最佳哈密顿圈遇到错误 求最佳哈密顿圈遇到错误 程序如下 %用矩阵翻转方法来实现二边逐次修**过程,求最佳哈密尔顿圈(H圈) clc clear load('zd.mat'); ...
- 龙格库塔法解微分方程组的matlab程序,MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc...
MATLAB实例源码教程:龙格库塔法求解微分方程组源代码实例.doc MATLAB实例源码教程龙格库塔法求解微分方程组源代码实例题目用经典 Runge-Kutta方法求下列一阶微分方程组的近似解y1 ...
- dmc matlab程序,matlab编的DMC程序.doc
matlab编的DMC程序 clear all; % close all; %系统模型建立 num=[0.8]; den=[225 1]; [a,b,c,d]=tf2ss(num,den); % st ...
- 约束优化的matlab程序,Matlab多约束最优化问题请高手求助,我编的程序
问题补充,这是我自己编的, 我刚学matlab没有多长时间,请大家帮忙看下程序有什么不对的地方! x1==csvread('x1.txt'); %x1 是300*10的矩阵 [x1,y,z,v ...
最新文章
- mongoose 实用 API 总结
- linux 查看java进程_Linux进程查看及管理工具(ps, vmstat, dstat, glances等)
- React中的路由react-router
- 【我拼搏的2016】-苦逼运维如何变身为SRE成长经历
- IIS服务器绑定域名问题
- golang 所有进程休眠_Golang-进程无故消失?
- php怎么求最小公倍数,PHP编程求最大公约数与最小公倍数的方法示例
- Java编码规范注意
- 作为测试人员,不要说有错误,要说有什么错误
- ubuntu删除OpenCV
- 下载webpack依赖报错“Refusing to install package with name “webpack“ under a package”
- Android柠檬水面试题,柠檬水的饮料摊(面试案例)
- Android Studio 之万恶 gradle
- 看了几个技术入股的帖子,忍不住写个自己亲身经历吧
- 用C#开发的一个通用的地铁换乘查询工具
- 实现数组的两两分组,任意倍数分组
- stm32 串口2空闲中断死机_STM32串口之空闲中断
- Java实现交叉表_Sql交叉表简单实现 | 学步园
- 推荐一款小巧而强大的截屏软件 FastStone Capture
- 推荐一招适合新手做网站引流量的方法,方法简单但实际有效