1 / 2
文档名称:

秩和比算法matlab程序,Matlab学习系32秩和比综合评价法.pdf

格式:pdf   大小:129KB   页数:2页
下载后只包含 1 个 PDF 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

秩和比算法matlab程序,Matlab学习系32秩和比综合评价法.pdf

上传人:鼠标 2023/6/8 文件大小:129 KB

下载得到文件列表

秩和比算法matlab程序,Matlab学习系32秩和比综合评价法.pdf

文档介绍

文档介绍:该【秩和比算法matlab程序,Matlab学习系32秩和比综合评价法】是由【鼠标】上传分享,文档一共【2】页,该文档可以免费在线阅读,需要了解更多关于【秩和比算法matlab程序,Matlab学习系32秩和比综合评价法】的内容,可以使用淘豆网的站内搜索功能,选择自己适合的文档,以下文字是截取该文章内的部分文字,如需要获得完整电子版,请下载此文档到您的设备,方便您编辑和打印。秩和?算法matlab程序,?.?秩和?综合评价法?、简单介绍秩和?法,是我国统计学家?凤调教授于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.?计算秩和?RSR(每个指标权重相同=1/m)若各评价指标的权重不同,则计算或加权秩和?WRSR其中,wj为第j个评价指标的权重。3.?计算概率单位编制RSR/WRSR频率分布表:将RSRi或WRSRi从?到?排列(值相同作为?组,或?定义分组),列出各组频数fi和各组累计频数∑fi;计算累积频率;再将pi转换为概率单位Probiti,其中,Probiti取标准正态分布的pi分位数+.?计算直线回归?程以累积频率所对应的概率单位Probiti为?变量,RSRi或WRSRi值为因变量,做线性回归:RSR/WRSR=a+b*Probit5.?分档排序按回归?程计算的RSR/WRSR估计值,对评价对象进?分档排序。三、Matlab实现例1?某市?民医院1983-1992年?作质量统计指标数据,如下表所?:其中,病死率x2,?平均住院?x6为负向指标。六项指标的权重依次为:=xlsread('');X=datas(:,2:end);w=[?????];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;?%计算加权秩和?:?加权,按?求和,再除以nfreq=tabulate(WRSR);?%统计WRSR的频数,频率,?freq的第3列为频率p=cumsum(freq(:,3))/100;?%计算累积频率p(end)=p(end)-1/(4*n);?%修正最后?个累积频率Probit=norminv(p,0,1)+5;?%计算标准正态分布的p分位数+5Probit=[ones(n,1),?Probit,?Probit.^2,?Probit.^3];[b,bint,r,rint,stats]=regress(WRSR,Probit);%三次多项式回归,??次或?次多项式回归R?和p值较差bstatsWRSRfit=Probit*b?%计算WRSR的估计值[s,ind]=sort(WRSRfit,'descend')?%对WRSR估计值从?到?排序运?结果:b?=??-??-?=?????=??=