《EXCEL在医学统计学中的应用.doc》由会员分享,可在线阅读,更多相关《EXCEL在医学统计学中的应用.doc(28页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、【精品文档】如有侵权,请联系网站删除,仅供学习与交流EXCEL在医学统计学中的应用.精品文档.资料和方法 1. 资料 主要根据高等医药院校教材卫生统计学3、上海第一医学院卫生统计教研组编写的医学统计方法4、中国医学百科全书医学统计学5、医用数理统计方法6和临床统计学7寿命表及其应用8等书中有关统计公式和例题进行编程。 2. 编程软件 应用Microsoft Excel 7.0 电子表格软件。编程的函数主要有: 2.1统计函数示例详见统计函数解释工作簿(请下载)。1)求集中趋势的统计函数:AVERAGE(number1,number2,)等,详细请见统计函数解释工作簿的集中趋势的统计函数工作表。
2、2)求离散趋势的统计函数:AVEDEV(number1,number2, )等,详细请见统计函数解释工作簿中的离散趋势的统计函数工作表。3)求频数分布的统计函数:FREQUENCY(data_array,bins_array)等,详细请见统计函数解释工作簿中的频数分布的统计函数工作表。4)正态分布的统计函数:NORMDIST(x,mean,standard_dev,cumulative)等,详细请见统计函数解释工作簿中的正态分布的统计函数工作表。5)总体均数估计和假设检验的统计函数:TTEST(array1,array2,tails, type)等,详细请见统计函数解释工作簿中的总体均数估计和
3、假设检验的统计函数工作表。6)方差分析的统计函数:FTEST(array1,array2)等,详细请见统计函数解释工作簿中的方差分析的统计函数工作表。7)X2检验的统计函数: CHIDIST(x,degrees_freedom)等,详细请见统计函数解释工作簿中的卡方检验的统计函数工作表。8)秩和检验的统计函数: RANK(number,ref,order)等,详细请见统计函数解释工作簿中的秩和检验的统计函数工作表。9)回归与相关:CORREL(array1,array2)等,详细请见统计函数解释工作簿中的回归与相关工作表。2.2常用数学、逻辑和查找函数省略,见EXCEL使用指南,示例详见一般函
4、数解释工作簿(请下载)。3应用EXCEL进行统计分析的方法用EXCEL进行统计分析有四种方法来实现:1) 直接应用EXCEL统计函数公式;2) 应用EXCEL分析工具库;3) 应用EXCEL电子表格和函数编制应用程序;4) 应用EXCEL统计宏。其中统计宏的编制需要一定软件基础,难以展开讨论,以下对其他三种方法结合实例分析。 3.1应用EXCEL统计函数公式进行统计计算 EXCEL函数公式不但有强大的科学运算函数和财务分析函数,而且还有丰富的统计函数,这些统计函数可直接对电子表格上指定区域内数据进行统计计算。以两个样本的均数T检验为例:3.1.1TTEST 求两个样本的均数t 检验TTEST判
5、断两个样本是否可能来自两个具有相同均值的总体。 语法TTEST (array1,array2,tails,type)Array1为第一个数据集;Array2为第二个数据集。 Tails指明单尾或双尾检验。如果 tails = 1,函数 TTEST 使用单尾检验。如果 tails = 2,函数 TTEST 使用双尾检验。 Type为 t 检验的类型。 1成对检验方法;2等方差双样本检验;3异方差双样本检验。 例1. 某单位研究饮食中缺乏维生素E与维生素A含量的关系,将同种属的大白鼠按性别相同,年龄、体重相近者配成对字,共8对,并将每对中的两头动物随机分到正常饲料组和维生素E缺乏组,过一定时期将大
6、白鼠杀死,测得其肝中维生素A的含量,见表,问不同饲料的大白鼠肝中维生素A含量有无差别? (点击图片可放大)操作步骤:(1)在J44:J51区域内输入正常饲料组维生素A含量,(2)在K44:K51区域内输入维生素E缺乏组维生素A含量,(3)在I44格中输入“TTEST(J44:J51,K44:K51,2,1)”。(4)输入回车键后,I44格就显示P值结果:“0.004001”。 推断分析:因为P值0.05,故可认为两种饲料组所得肝中维生素A含量的差别有显著性意义,即维生素E缺乏对大白鼠肝中维生素A含量有影响。详细请见统计函数解释工作簿中总体均数估计和假设检验的统计函数工作表。 例2. 某克山病区
7、测得11例克山病患者与13名健康人的血磷值(mmol/L)如下,问该地急性克山病患者与健康人的血磷值是否不同? 患者X1:0.84 1.05 1.20 1.20 1.39 1.53 1.67 1.80 1.87 2.07 2.11 健康人X2:0.54 0.64 0.64 0.75 0.76 0.81 1.16 1.20 1.34 1.35 1.45 1.87 操作步骤:(1)在J28:J38区域内输入患者血磷值;(2)在K28:K40区域内输入健康人血磷值;(3)在I28格中输入TTEST(J28:J38,K28:K40,2,2);(4)输入回车键后,I28格就显示P值结果:0.019337
8、。 推断分析:因为P值0.05,故两样本的均数差别不显著。详细请见统计函数解释工作簿中总体均数估计和假设检验的统计函数工作表。3.1.2卡方检验 也是常用的统计方法,可用CHITEST函数直接计算,举例说明如下: 语法CHITEST (actual_range, expected_range)Actual_range用以检验期望值的观测值区域。 Expected_range指包含行列汇总的乘积除以总和后的数据区域。 该函数公式为 CHITEST=p(Xc2)。 其中:Aij = i行j列中的实际值;Eij = i行j列中的期望值。r = 行数;c = 列数。 CHITEST 返回 chi-sq
9、uared 统计的概率和自由度df,其中 df=(r-1)(c-1)。 例4. 在二乙基亚硝胺诱发大白鼠鼻咽癌的实验中,一组单纯用亚硝胺向鼻腔滴注(鼻注组),另一组在鼻注基础上加肌注维生素B12,实验结果见表,问两组发癌率有无差别? 操作步骤:1)在J3:K4区域内输入实际频数数据;2)在J6:K7区域内输入理论频数数据;3)在I3格中输入CHITEST(J3:K4,J6:K7), 4)输入回车键后,I3格就显示P值: 0.01092。 推断分析:卡方检验结果显示P值=0.01092,0.05,故可认为两组发癌率有差别,说明增加肌注维生素B12有可能提高大白鼠的鼻咽癌发生率。详细请见统计函数解
10、释工作簿中卡方检验的统计函数工作表。 3.2应用EXCEL分析工具库进行统计分析 详细请见统计分析工具解释工作簿(请下载)。在EXCEL中还提供统计分析工具库有如下:t-检验:成对双样本均值分析工具、双样本等方差假设分析工具,双样本异方差假设分析工具及其公式;F - 检验:双样本方差分析分析工具,Anova:单因素方差分析工具,Anova:可重复双因素分析工具,Anova:无重复双因素分析工具,相关系数分析工具及其公式,协方差分析工具,直方图分析工具,双样本均值分析,抽样分析工具,相关和回归分析工具,描述统计分析工具,排位和百分比排位分析工具,随机数发生器分析工具,指数平滑分析工具,移动平均分
11、析工具。这些分析工具库使用时都需要对选项进行选择,下面就单因素方差分析、无重复双因素分析、可重复双因素分析工具、相关和回归分析工具的应用,举例说明如下: 3.2.1Anova:单因素方差分析工具 此分析工具通过简单的方差分析(anova),对两个以上样本均值进行相等性假设检验。 例5. 某职业病防治院对31名石棉矿工中的石棉肺患者、可疑患者及非患者进行了用力肺活量(L)测定,结果见表5.1,问三组石棉矿工的用力肺活量有无差别? (点击图片可放大)分析提示:此为成组设计的定量资料,是单因素三个样本均数的比较,可用成组方差分析。 操作步骤: 1)设置数据区域,并输入数据。如本例J6:J17为石棉肺
12、患者数据区,K6:K15为可疑患者数据区,L6:L17为非患者数据区。各列第一行为标志位,即石棉肺患者、可疑患者和非患者。 2)选取工具栏中数据分析。 3)选定Anova:单因素方差分析。 4)选择确定按钮。 5)在输入区域 框输入J6:L17。 6)在分组方式框选择逐列,即样本数据按列分组。 7)打开标志位于第一行复选框。 8)在显著性检验水平框输入0.05 9)在输出区域框输入A14。 10)再选择确定按钮。 11)有关单因素方差分析检验结果显示如下:详细请见统计分析工具解释工作簿中单因素方差分析工作表。 (点击图片可放大)3.2.2Anova:“无重复双因素分析分析工具此分析工具通过双因
13、素 anova 分析(但每组数据只包含一个样本),对两个以上样本均值进行相等性假设检验。 例6. 为研究雌激素对子宫发育的作用,以四个种类的未成年雌性大白鼠各3只,每只按一种剂量注射雌激素,经一定时期取出子宫,称重,结果见表。试比较雌激素的作用在三种剂量间的差别,同时比较四个种系大白鼠之间的差别。 (点击图片可放大)分析提示:此是计量资料,按雌激素剂量分为三组,这是一个因素,按大白鼠的不同种系分组,又是一个因素,该试验属随机分组设计,故应用两因素方差分析。 操作步骤: 1)设置数据区域,在垂直和水平两个方向上,进行分组分层设计数据区域并输入数据。如本例、B1、C1、D1为垂直分组,A2、A3、
14、A4、A5为水平分层,将样本数据输入到相应的各组各层区域内,即甲种系大白鼠按雌激素剂量分别将106、116、145(g)输入到B2、C2、D2单元格,乙、丙和丁种系大白鼠子宫重量依此类推输入到B3、C3和D3, B4、C4和D4,B5、C5和D5。 2)选取工具栏中数据分析。 3)选定Anova:无重复双因素分析。 4)选择确定按钮。 5)在输入区域 框输入A1:D5。 6)打开标志复选框。 7)在显著性检验水平框输入0.05。 8)在输出区域框输入A20。 9)再选择确定按钮。 10)有关Anova:无重复双因素分析检验结果显示如下:详细请见统计分析工具解释工作簿中无重复双因素方差分析工作表
15、。(点击图片可放大)检验结果:P0.01,在水准拒绝H0,接受H1, 推断分析:故认为雌激素不同剂量组间的差别有显著性意义,大白鼠不同种系间差别亦有显著性意义。详细请见统计分析工具解释工作簿中无重复双因素方差分析工作表。3.2.3Anova:可重复双因素分析工具 此分析工具是对单因素 anova 分析的扩展,即每一组数据包含不止一个样本。 例6. 为探讨某化学反应中温度和催化剂对收率的影响,有人选了4种温度和三种不同的催化剂甲,乙,丙,对所有可能的组合在相同条件都重复2次试验,所得数据见表7.15,试判断温度、催化剂的作用以及它们之间的交互作用是否显著。 (点击图片可放大)解题思路:此题是某化
16、学反应中温度和催化剂二因素作用下重复2次试验的计量资料,考虑温度和催化剂两个因素是否确系影响收率的重要因素,进一步考虑这两个因素有无交互作用,根据题意和数据特点可用重复试验的两因素方差分析。 操作步骤:1)设置数据区域,在垂直和水平两个方向上,进行分组分层设计数据区域并输入数据。 如本例B1、C1、D1、E1为不同温度分组,A2和A3、A4和A5、A6和A7为甲,乙,丙催化剂分层,其中A2和A3为甲催化剂同一层次重复数据区域,A4和A5,A6和A7为乙,丙催化剂二个层次的重复数据区域。将样本数据输入到相应的各组各层区域内。如重复数据为3次,则每一层应有3个数据区域,A2,A3和A4;A5,A6
17、和A7;A8,A9和A10各为一层。重复次数增加,相应以此类推。 2)选取工具栏中数据分析。 3)选定Anova:可重复双因素分析。 4)选择确定按钮。显示方差分析:可重复双因素分析对话框。 5)在输入区域 框输入A1:E7。 6)在每个样本的行数框内输入2,意味重复试验2次的数据。 7)在显著性检验水平框输入0.05。 8)在输出区域框输入A20。 9)再选择确定按钮。 10)有关方差分析:可重复双因素分析检验结果显示如下: 检验结果:1)温度因素的F值(98.30363454),催化剂因素的F值(63.14869238),温度和催化剂交互作用的F值(1.767047357);2)温度因素的
18、P值1.03007E-08(P0.01); 催化剂因素的P值4.2678E-07(P0.05);推断分析:温度和催化剂交互作用不显著;温度因素对收率有极显著影响;催化剂因素对收率也有极显著影响。详细请见统计分析工具解释工作簿中可重复双因素方差分析工作表。 3.2.4相关系数分析工具 此分析工具及其公式可用于判断两组数据集(可以使用不同的度量单位)之间的关系。总体相关性计算的返回值为两组数据集的协方差除以它们标准偏差的乘积: 注释如果需要返回两个单元格区域的相关系数,可直接使用CORREL工作表函数。 例7. 某地一年级12名女大学生的体重与肺活量数据如下,试求肺活量(L)Y对体重(kg)X的直
19、线回归方程。 体 重(kg)42424646465050505252X5858肺活量(L)2.552.202.752.402.802.813.413.103.462.85Y3.503.00分析提示:此为求肺活量(L)Y对体重(kg)X的直线回归方程,可用直线回归与相关分析。 操作步骤: 1)设置数据区域,并输入数据。如本例J4:J16为体重(kg)X自变量数据区,K4:K16为肺活量(L)Y因变量数据区。 2)选取“工具”栏中“数据分析”。 3)选定“相关系数”。 4)选择“确定”按钮。显示“相关系数”对话框。 5)在“输入区域”框输入J4:K16。 6)在“分组方式”框选择“逐列”,即自变量
20、和因变量数据按列分组。 7)打开“标志位于第一行”复选框。 8)在“输出区域”框输入A14。 9)再选择“确定”按钮。 10)有关相关系数分析结果显示如下: 体重(kg)肺活量(L)体重(kg)1肺活量(L)0.7494821 检验结果:X相关系数r值(0.749482342)。详细请见统计分析工具解释工作簿中相关与回归工作表。 3.2.5回归分析工具 例7. 同上 操作步骤: 1)设置数据区域,并输入数据。如本例J4:J16为体重(kg)X自变量数据区,K4:K16为肺活量(L)Y因变量数据区 。 2)选取工具栏中数据分析。 3)选定回归分析。 4)选择确定按钮。显示回归分析对话框。 5)在
21、Y 值输入区域 框输入K4:K16。 6)在X值输入区域 框输入J4:J16。 7)在分组方式框选择逐列,即自变量和因变量数据按列分组。 8)打开标志位于第一行复选框。 9)在输出区域框输入A34。 10)再选择确定按钮。 11)有关回归分析结果显示如下: (点击图片可放大)检验结果:回归系数b值(0.058826087),截距 a值(0.000413043);直线回归方程(Y=0.000413+0.058826X)。 推断分析:今1=1, 2=10,查附表F界值表,得P0.01。按=0.05水准拒绝H0,接受H1,故可认为一年级女大学生肺活量与体重之间有直线关系。详细请见统计分析工具解释工作
22、簿中相关与回归工作表。 3.3应用EXCEL电子表格和函数编制应用程序进行统计分析 根据统计的原理和计算的步骤,应用上述EXCEL函数,就可以编写成常用的统计程序,实际编程,举例如下: 3.3.1统计程序设计方法 示例详细请见统计程序设计解释工作簿(请下载)。在EXSEL电子表上编写医学统计应用程序,尽管与真正意义上程序设计有很大差别,但是在程序设计某些基本方面还是很相似,如根据需要解决的具体问题画出流程图,按流程图应用计算机语句(在此应用数学公式、数据管理命令和统计函数)一步一步地解决流程中每一个方面小问题,并把数据、公式和函数写在电子表格的一定区域和单元格,最终得到我们需要解决问题的最后结
23、果。这一过程与EXCEL的宏的制作过程相似,只不过没有记录成一个宏的名字,其解决问题全过程都罗列在电子表上,比较直观,容易发现问题。最大的好处在于,在解决课题具体统计问题时,EXCEL等没有现成程序和软件,自己动手编程,能恰到好处解决问题,为科研带来很大方便。特别在反复计算或数据量大时更是如此。为了能清楚地说明编程方法,采用举例说明,便于读者理解。 3.3.2计量资料的统计程序设计以几何均数t检验和单因素方差分析及二二比较为例。几何均数t检验程序流程框图: 程序设计过程: 根据程序流程框图的要求,在电子表格上输入EXCEL函数,为了更简洁表达意思,把在C3单元格中输入=IF(B3=0,LOG1
24、0(B3),简化成C3=IF(B3=0,LOG10(B3)。详细情况可打开成组设计的两样本几何均数的t检验程序观看,过程说明如下: 1)设立数据区并且转换数据先B3:B22作为样本甲的数据输入区,C3:C22作为对样本甲的数据转换成对数形式的区域,C3=IF(B3D33,0.01,IF(AND(B30B33),0.05)双侧:B37=IF(B30E33,0.01,IF(AND(B30C33),0.05)5)分析、推断和结论:省略4.2.2程序详细请见统计程序设计解释工作簿中几何均数t检验工作表。例8选甲型流感病毒血凝抑制抗体滴度(倒数)5者24人,随机分为两组,每组12人,用甲型流感病毒活疫苗
25、进行免疫,一组用气雾法,另一组用鼻腔喷雾法。免疫后一月采血,分别测定血凝抑制抗体滴度,结果如下,问两法免疫的效果有无差别? 气 雾 组,X1:40 20 30 25 10 15 25 30 40 10 15 30 鼻腔喷雾组,X2:50 40 30 35 60 70 30 20 25 70 35 25计算结果: t值(2.93443),B35为单侧检验结果(0.01),E35为单侧检验P值(0.00384);B37为双侧检验结果(9.21,方差明显不齐,IF(AND(B225.99),方差齐性不好,方差齐性好)10.设置方差分析表。 求组内、组间和总变异:B27=E19-E22;B28=(B1
26、62/B18+C162/C18+D162/D18)-E22;B29=B27-B28 11.求组内、组间和总自由度:C27=E18-1;C28=3-1;C29=E18-3 12.求均方:D28=B28/C28;D29=B29/C29 13.求F值:E28=D28/D29 14.查F界值:B31=VLOOKUP(D30,A40:Y90,B30);E31=VLOOKUP(D30,A100:Y149,B30)。 15.求P值:B32=IF(E28E31,0.01,IF(AND(E28B31),0.05)。 16.作两两比较。设置两两比较表。将各组均数排序并自动填入表中:A46=LARGE(B17:D1
27、7,1)A47=LARGE(B17:D17,2)A48=LARGE(B17:D17,3) 17.根据排序情况自动将各组样本数填入表中:B46=HLOOKUP(A46,B17:D18,2)B47=HLOOKUP(A47,B17:D18,2);B48=HLOOKUP(A48,B17:D18,2)。18.求对比组均数差值:(Xa-Xb):D46=A46-A48;D47=A46-A47;D48=A47-A48 19.求各对比组q值:F46= D46/SQRT(D29/2)*(1/B46+1/B48);F47=D47/SQRT(D29/2)*(1/B46+1/B47);F48= D48/SQRT(D29
28、/2)*(1/B47+1/B48)。 20.设置对比组组距和自由度:E46=B49,B49=3;E47=D49,D49=2;E48=F49,F49=2;H48=C29 21.查找q界值(0.05),并填入表中:C51=VLOOKUP(H49,AA62:AJ77,B49);C52=VLOOKUP(H49,AA62:AJ77,D49);C53=VLOOKUP(H49,AA62:AJ77,F49);G46=C51;G47=C52;G48=C53。 22.查找q界值(0.01), 并填入表中:E51=VLOOKUP(H49,AA82:AJ97,B49);E52=VLOOKUP(H49,AA82:AJ9
29、7,D49);E53=VLOOKUP(H49,AA82:AJ97,F49);H46=E51;H47=E52;H48=E53。 23.求P值:C54=IF(F46E51,0.01,IF(AND(F46C51),0.05);C55=IF(F47E52,0.01,IF(AND(F47C52),0.05);C56=IF(F48E53,0.01,IF(AND(F48C53),0.05);I46=C54;I47=C55;I48=C56。 24.设置F0.05界值表区域:A61:Y111区域为0.05F界值表;B61:Y61区域为分子自由度区域。A62:A111区域为分母自由度区域。输入F0.05界值:B6
30、2:Y111区域填入F0.05界值。 25. 设置F0.01界值表区域:A120:Y170区域为0.01界值表;B120:Y120区域为分子自由度区域。A121:A170区域为分母自由度区域。输入F0.01界值:B121:Y170区域填入F0.01界值。 26.设置Q界值表:确定Q0.05界值表区域:AA61:AJ77为Q0.05界值表区域;AB61:AJ61为组距区域;AA62:AA77为自由度区域;输入Q0.05界值:AA62:AJ77区域填入Q0.05界值。 确定Q0.01界值表区域;AA81:AJ97为Q0.01界值表区域;AB81:AJ81为组距区域;AA82:AA97为自由度区域;
31、输入Q0.01界值:AA82:AJ97区域填入Q0.01界值。27程序:详细请见统计程序设计解释工作簿中3样本成组方差分析工作表。示例请见EXCEL实用临床统计程序.xls。例9例题和结果同例5。 操作调试:(1)按照上述方法在工作表中组成一个三样本成组方差分析电子表格。 (2)分别将三样本的原始数据输入到B3:B15,C3:C15和D3:D15区域中。(3)A25:G29为方差分析表,B27为总方差(10.8),B28为组间方差(9.265656566),B29为组内方差(1.534343434);C27为总自由度(30),C28为组间自由度(2),C29为组内自由度(28);D28为组间均
32、方差(4.632828283),D29为组内均方差(0.05479798);E28为方差值(84.5437788),F28为P值大于或小于0.05或0.01的情况(0.01),G28为P的具体值(1.36453E-12);B31为f值0.05界值(3.340389299),E31为f值0.01界值(5.452932328)。 (4)A44:I48为二二比较的分析表。检验结果:F46、F47和F48分别为1组/3组、1组/2组和2组/3组对比组的q值,I46、I47和I48分别为1组/3组、1组/2组和2组/3组对比组的P值。 3.3.4 计数资料的统计程序设计以有序数据的卡方检验为例。程序流程
33、图: 程序设计过程: 根据程序流程框图的要求,另选一张工作表在电子表格上开始输入EXCEL函数。 1)确定行列输入数据区域:B5:E6 有序数据的卡方数据输入区域。 2)求各列的和:B7=SUM(B5:B6);C7=SUM(C5:C6);D7=SUM(D5:D6);E7=(E5:E6)。 算中间值: 4)求ac:B9 =B5*B4;C9=C5*C4;D9=D5*D4;E9=E5*E4。F9 =SUM(B9:E9)。 求nc:B10 =B7*B4;C10=C7*C4;D10=D7*D4;E10=E7*E4。F10 =SUM(B10:E10)。 求nc2:B11 =B7*B42;C11 =C7*C
34、42;D11 =D7*D42;E11 =E7*E42。F10 =SUM(B11:E11)。 计算第一行所占比率:P=A/N= B13 =F5/F7 计算卡方值:B14 =(F9-(F5*F10)/F7)2)/(F11-(F102)/F7)*B13*(1-B13)5)求卡方的概率P值:B15 =CHIDIST(B14,1)。 6)根据已知的P值,作推断和结论:B16=IF(B15=0.01,B150.05),两药有显著性差异,两药没有显著性差异 )。7)程序:详细请见统计程序设计解释工作簿中有序数据的卡方检验工作表。例10用某药治疗胃癌,对溃疡型和梗阻型胃癌的治疗结果见表7,试分析此药对两型胃癌
35、的疗效。 3.3.5非参数资料的统计程序设计以配对设计的秩和检验为例。配对设计的秩和检验的程序框图程序设计:1)数据输入及排序区域: (1)设B3:B12为甲样本数据区域,C3:C13为乙样本数据区域。 (2)计算差值:F3=B3-C3, F4=B4-C4, F5F12依次类推,以下用表示。(3)符号标志:D3=IF(F30,1, ),。 (4)取绝对值:E3=ABS(F3),。 (5)说明:在数据输入后,秩和检验前,须对数据资料以差值的绝对值为关键词,进行递增排序。该区域排序可用Excel工作区顶部的菜单栏中数据排序命令,通过排序对话框来实现。2)编秩、秩次校正和确定检验统计量 (1)符号标志:B15=IF(D3= ,0,D3),。(2)自小而大排列绝对值: C15=SMALL(E3:E12,1),C16=SMALL(E3:E12,2),。 (3)复制输入区域绝对值:D15=E3,。 (4)对排序后的数据统一编秩次:E15=RANK(D15,D15:D24,1),。 (5)去零校正:F15=E15-E13。 (6)去零秩次:G15=IF(F150,F15, ),。去除差值为零的秩次。 (7)同号处理: H15=IF(G15= , ,1); H16=IF(AND(G16 ,AND(G16=G15,B16*B150),G16+1,G16)