随机数的产生与模拟

随机数的产生与模拟1第一页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

用随机模拟方法解决实际问题时，首先要解决的是随机数的产生方法，或称随机变量的抽样方法。

本章目录2第二页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

伪随机数：在计算机上用数学方法产生均匀随机数是指按照一定的计算方法而产生的数列，它们具有类似于均匀随机变量的独立抽样序列的性质，这些数既然是依照确定算法产生的，便不可能是真正的随机数，因此常把用数学方法产生的随机数称为伪随机数。本章目录3第三页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

均匀分布随机数：本章目录4第四页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

均匀分布随机数：该定理说明了任意分布的随机数均可由均匀分布的随机数变换得到。常简称的随机数为均匀分布随机数。本章目录5第五页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

1均匀随机数的产生均匀随机数的产生：主要有线性同余法（LCG），组合同余法，反馈位移寄存器方法等

本章目录6第六页，共七十七页，编辑于2023年，星期四均匀随机数的产生：随机数的产生与模拟

1均匀随机数的产生本章目录线性同余法（LCG）的递推公式为：7第七页，共七十七页，编辑于2023年，星期四均匀随机数的产生：随机数的产生与模拟

1均匀随机数的产生本章目录当，上式称为混合同余发生器，当时，称为乘同余发生器，此时当模为素数时，称它为素数模乘同余发生器。

8第八页，共七十七页，编辑于2023年，星期四两个常用的混合式发生器：随机数的产生与模拟

1均匀随机数的产生本章目录9第九页，共七十七页，编辑于2023年，星期四常用的素数模乘同余发生器：随机数的产生与模拟

1均匀随机数的产生本章目录10第十页，共七十七页，编辑于2023年，星期四常用的素数模乘同余发生器：随机数的产生与模拟

1均匀随机数的产生本章目录11第十一页，共七十七页，编辑于2023年，星期四反馈位移寄存器法（FSR）：对寄存器中的二进制数码作递推运算，其中是给定的正整数，为给定的常数。取数列中连续的位构成一个位二进制整数，一直下去，一般地有

令

则即为FSR方法产生的均匀随机数列。随机数的产生与模拟

1均匀随机数的产生本章目录12第十二页，共七十七页，编辑于2023年，星期四组合发生器：先用一个随机数发生器产生的随机数列为基础，再用另一个发生器对随机数列进行重新排列得到的新数列作为实际使用的随机数。这种把多个独立的发生器以某种方式组合在一起作为实际使用的随机数，希望能够比任何一个单独的随机数发生器得到周期长、统计性质更优的随机数，即组合发生器。

随机数的产生与模拟

1均匀随机数的产生本章目录13第十三页，共七十七页，编辑于2023年，星期四组合发生器：

随机数的产生与模拟

1均匀随机数的产生本章目录Maclaren和Marsaglia在1965年提出的著名的组合发生器是组合同余发生器,该算法的具体步骤如下：14第十四页，共七十七页，编辑于2023年，星期四组合发生器：

1用第一个LCG产生个随机数，一般取。这个随机数被顺序地存放在矢量中。置;2用第二个LCG产生一个随机整数,要求;3令，然后再用第一个LCG产生一个随机数,令；置;4重复2~3，得随机数列，即为组合同余发生器产生的数列。若第一个LCG的模为，令，则为均匀随机数

随机数的产生与模拟

1均匀随机数的产生本章目录15第十五页，共七十七页，编辑于2023年，星期四由均匀分布随机数产生非均匀分布随机数的主要方法有：逆变换法，合成法和筛选法。

随机数的产生与模拟

2非均匀随机数的产生本章目录16第十六页，共七十七页，编辑于2023年，星期四1逆变换法：随机数的产生与模拟

2非均匀随机数的产生对任意分布函数,要产生服从该分布的随机数，由定理知其抽样步骤为：（1）由抽取;（2）

计算本章目录17第十七页，共七十七页，编辑于2023年，星期四1逆变换法：

随机数的产生与模拟

2非均匀随机数的产生本章目录例1已知(柯西分布)，试给出其抽样方法。18第十八页，共七十七页，编辑于2023年，星期四1逆变换法：随机数的产生与模拟

2非均匀随机数的产生解：设,则,因此其抽样步骤如下：（1）由抽取;（2）计算本章目录19第十九页，共七十七页，编辑于2023年，星期四1逆变换法：其SAS程序为（产生100个服从柯西分布的随机数）：dataex1;seed=678;doI=1to100;r=ranuni(seed);x=tan(3.14159*(r-0.5));output;end;run;随机数的产生与模拟

2非均匀随机数的产生本章目录20第二十页，共七十七页，编辑于2023年，星期四2合成法：随机数的产生与模拟

2非均匀随机数的产生其想法是：如果X的密度难于抽样，而X关于Y的条件密度以及Y的密度函数均易于抽样，则X的随机数可如下产生：由Y的密度抽取y

由条件密度抽取x

则X服从本章目录21第二十一页，共七十七页，编辑于2023年，星期四2合成法：随机数的产生与模拟

2非均匀随机数的产生当为离散形式时，即,其中是密度函数，其抽样过程如下：1产生一个正的随机整数，使得，2产生分布为的随机数。本章目录22第二十二页，共七十七页，编辑于2023年，星期四2合成法：随机数的产生与模拟

2非均匀随机数的产生本章目录

设时梯形分布的密度函数为，试用合成法产生其随机数。

例223第二十三页，共七十七页，编辑于2023年，星期四2合成法：随机数的产生与模拟

2非均匀随机数的产生解：首先将进行分解,即，其中其抽样框图为本章目录24第二十四页，共七十七页，编辑于2023年，星期四2合成法

：随机数的产生与模拟

2非均匀随机数的产生产生产生令产生令输出YN本章目录25第二十五页，共七十七页，编辑于2023年，星期四2合成法：其SAS抽样程序如下（假若产生100个随机数,）：dataex2;seed=789;a=0.3;doI=1to100;r=ranuni(seed);r3=ranuni(seed);ifr1<=athendo;u=ranuni(seed);x=u;end;elsedo;u=ranuni(seed);v=ranuni(seed);x=max(u,v);end;output;end;run;随机数的产生与模拟

2非均匀随机数的产生本章目录26第二十六页，共七十七页，编辑于2023年，星期四3筛选抽样法：随机数的产生与模拟

2非均匀随机数的产生假设我们要从抽样，如果可将表示成，其中是一个密度函数且易于抽样，而，是常数，

本章目录27第二十七页，共七十七页，编辑于2023年，星期四3筛选抽样法：随机数的产生与模拟

2非均匀随机数的产生X的抽样可如下进行：1由抽取，由抽取

2如果，则;否则，转1则X的密度函数为本章目录28第二十八页，共七十七页，编辑于2023年，星期四3筛选抽样法：随机数的产生与模拟

2非均匀随机数的产生本章目录设，试用筛选法抽取其随机数。

例329第二十九页，共七十七页，编辑于2023年，星期四3筛选抽样法：随机数的产生与模拟

2非均匀随机数的产生解：因为：，即：则抽样框图如下：本章目录30第三十页，共七十七页，编辑于2023年，星期四3筛选抽样法：随机数的产生与模拟

2非均匀随机数的产生独立产生令NY本章目录31第三十一页，共七十七页，编辑于2023年，星期四3筛选抽样法：其SAS程序如下：dataex3;seed=789;doI=1to100;r1=ranuni(seed);r2=ranuni(seed);ifr1<=r2**3thendo;x=r2;output;end;end;run;随机数的产生与模拟

2非均匀随机数的产生本章目录32第三十二页，共七十七页，编辑于2023年，星期四蒙特卡罗（MonteCarlo）方法（即随机模拟方法）求解实际问题的基本步骤包括：1建模：对所求的问题构造一个简单而又便于实现的概率统计模型，使所求的解恰好是所建模型的参数或有关的特征量。2改进模型：根据概率统计模型的特点和计算实践的需要，尽量改进模型，以便减少误差和降低成本，提高计算效率。3模拟试验4求解：对模拟结果进行统计处理，给出所求问题的近似解。随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用本章目录33第三十三页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用计算定积分（1）随机投点法

①赋初值:试验次数n=0,成功次数m=0;规定投点试验的总次数N;②产生两个相互独立的均匀随机数置n=n+1;③判断n≤N是否成立,若成立转④,否则停止试验,转⑤;④判断条件是否成立,若成立置m=m+1,然后转②,否则转②;⑤计算m/N,则本章目录34第三十四页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用计算定积分（1）随机投点法

本章目录35第三十五页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用计算定积分(2)平均值估计法

平均值估计法的计算步骤：

①产生[0,1]区间的均匀随机数②计算③令=,则为积分值的近似解.本章目录36第三十六页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用计算定积分(3)重要抽样法

重要抽样法的计算步骤为：

①产生均匀随机数

②用直接抽样法产生随机数,即由计算则③计算=,则是的估计量.本章目录37第三十七页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用计算定积分(4)分层抽样法分层抽样法的计算步骤如下：

本章目录38第三十八页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用(1)随机投点法多重积分随机投点法计算步骤为:计算多重定积分本章目录39第三十九页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用(1)随机投点法计算多重定积分本章目录40第四十页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用(2)平均值估计法计算多重定积分多重积分的平均值法计算步骤为

本章目录41第四十一页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用(2)平均值估计法计算多重定积分本章目录42第四十二页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用计算多重定积分用蒙特卡罗方法计算积分值时，误差的阶数为,它与多重积分的重数k无关，而用其他数值方法计算多重积分时，其误差与重数k是有关的，可见当k≥3时，使用蒙特卡罗方法计算多重积分将显现出很大的优越性．

本章目录43第四十三页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用例4：用上述四种方法计算

(1)随机投点法

dataE1;Dok=1to1000;m=0;Doh=1to1000;a=ranuni(32789);b=ranuni(32789);ifb<=(exp(a)-1)/(exp(1)-1)thenm=m+1;end;I1=m/1000*(exp(1)-1)+1;output;E1=abs(I1-(exp(1)-1));end;run;procmeansdata=e1MeanVar;varI1;run;本章目录44第四十四页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用例4：用上述四种方法计算

(2)平均值估计法

dataE2;Dok=1to1000;s=0;Doi=1to1000;x=ranuni(32789);fx=exp(x);s=s+fx;end;I2=s/1000;output;E2=abs(I2-(exp(1)-1));end;run;procmeansdata=e2MeanVar;varI2;run;

本章目录45第四十五页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用例4：用上述四种方法计算

(3)重要抽样法dataE3;dok=1to1000;s=0;Doi=1to1000;r=ranuni(32789);x=(3*r+1)**(1/2)-1;s=s+exp(x)/(1+x);end;I3=3/(2*1000)*s;output;E3=abs(I3-(exp(1)-1));End;run;procmeansdata=e3MeanVar;varI3;run;本章目录46第四十六页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用例4：用上述四种方法计算

(4)分层抽样法dataE4;Dok=1to1000;s1=0;s2=0;Doi=1to400;ri=ranuni(32789);r1=0.5*ri;f1=exp(r1);s1=s1+f1;end;Doj=1to600;rj=ranuni(32789);r2=0.5+0.5*rj;f2=exp(r2);s2=s2+f2;end;I4=s1*(1/800)+s2*(1/1200);output;E4=abs(I4-(exp(1)-1));end;run;procmeansdata=e4MeanVar;varI4;run;

本章目录47第四十七页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用例4：用上述四种方法计算

结果模拟方法均值方差随机投点法（I1）1.71768340.000728670平均值估计法(I2)1.71801360.00027492重要抽样法(I3)1.71818630.000024254分层抽样法(I4)1.71812820.000062622

=e-1=1.71828,这些方法的I值与真实值很接近，而方差也都比较小，同时看出，这次模拟其方差有以下关系：Var(I4)<Var(I3)<Var(I2)<Var(I1)。本章目录48第四十八页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用(1)随机投点法例5：计算二重积分

Datae5;Dok=1to1000;m=0;h=0;Doh=1to1000;a1=ranuni(32789);a2=ranuni(32789);b=ranuni(32789);ifb<=(exp(a1+a2)-1)/(exp(2)-1)thenm=m+1;end;I5=(exp(2)-1)*(m/1000)+1;output;E5=abs(I5-(exp(1)-1)**2);end;run;procmeansdata=e5MeanVar;varI5;run;

本章目录49第四十九页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用(2)平均值估计法例5：计算二重积分

dataE6;Dok=1to1000;m=0;s=0;Doh=1to1000;a1=ranuni(32789);a2=ranuni(32789);if0<=a1<=1and0<=a2<=1thendo;fx=exp(a1)*exp(a2);s=s+fx;end;end;I6=s/1000;output;E6=abs(I6-(exp(1)-1)**2);end;run;procmeansdata=e6MeanVar;varI6;run;本章目录50第五十页，共七十七页，编辑于2023年，星期四应用实例随机数的产生与模拟

3MonteCarlo方法在解确定性问题中的应用例5：计算二重积分

模拟方法均值方差随机投点法(I5)2.94815100.0865186平均值估计法(I6)2.95222550.0015766对于多元积分也有Var(I6)<Var(I5)本章目录51第五十一页，共七十七页，编辑于2023年，星期四随机服务系统研究的对象是服务系统，如到理发店理发，理发师与顾客构成了一个服务系统；到商店买东西，售货员与顾客就构成了一个服务系统。随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用本章目录52第五十二页，共七十七页，编辑于2023年，星期四随机服务系统一般具有三要素，顾客、排队规则和窗口

随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用本章目录53第五十三页，共七十七页，编辑于2023年，星期四

1顾客顾客到达排队系统的过程也称为输入过程。顾客的来源和到达排队系统的情况是多种多样的。顾客来源可能是有限的，也可能是无限的。顾客到达方式可能是连续的，也可能离散；可能是一个一个的，也可能是成批的或大量的；顾客相继到达的间隔时间可以是确定型的，也可以是随机型的；顾客的到达可以相互独立，也可以是相互关联的。如果描述顾客相继到达的间隔时间分布和所含参数（如期望值、方差等）都与时间无关，则称为平稳（Stationary）输入过程，否则称之为非平稳输入过程。

随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用本章目录54第五十四页，共七十七页，编辑于2023年，星期四2排队规则

常用的规则有：损失制（LossingSystem）：顾客到达时，若所有服务台均被占，该顾客就自动消失。如通常使用的损失制电话系统。等待制（WatingSystem）：顾客到达时，若所有服务台均被占，他们就排成队伍，等待服务。服务次序可采用以下各种规则：

先到先服务：即按到达的次序接受服务。

后到先服务：即后到的顾客、先接受服务。如在有的流水装配线上，后到的零件先装配；在通讯系统中，最后到达的信息一般最有价值。

随机的服务：当服务机构得空时，在等待顾客中、随机地选取一名进行服务，也即每一等待的顾客被选到的概率相同。

优先权服务：如医院对重患或急诊患者予以优先治疗、重要电话先接通等。

多个服务台：当顾客到达时可以按如下规则在每个服务台前排成一个队：第1，n+1,2n+1,…个顾客排入第一队；第2，n+2,2n+2，…个顾客排入第二队等等。或者排成一个公共的队，当一个服务台得空时，队首顾客进入服务。队列数目排队队列有单列和多列之分。顾客排队后由于等待时间过长而中途离队，但也有不允许中途离队的情况，这种情况必须坚持到服务完为止。在多队列排队情况下，各队列之间的顾客有的可以相互转移，有的不允许转移。

随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用本章目录55第五十五页，共七十七页，编辑于2023年，星期四

3窗口：服务台的个数可以是一个或几个，可以是单个服务，也可以是成批服务。

随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用本章目录56第五十六页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用例6得意系统可靠性估计。设系统由顺次连接的两个元件组成。两个元件中，任何一个元件发生故障系统就停止工作。第一个元件有两个组成部分A,B（它们并联）。第二个元件有一个部件C组成。试用MonteCarlo法求:1估计系统工作的概率，已知组成部件的工作概率分别为：2绝对误差，其中为系统的可靠性。可用分析的方法获得。进行50次试验。本章目录57第五十七页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用解：易知。其SAS程序如下：dataex6;seed=12345;doI=1to50;a=ranuni(seed);b=ranuni(seed);c=ranuni(seed);if(a<0.8orb<0.85)andc<0.6thendo;d=1;output;end;elsedo;d=0;output;end;end;run;procmeansdata=ex6;vard;outputout=resultsum=dsumn=n;run;dataout;setresult;p=dsum/n;pdiff=0.582-p;run;procprintdata=out;varppdiff;run;本章目录58第五十八页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用其输出结果为：

OBSPPDIFF10.580.002本章目录59第五十九页，共七十七页，编辑于2023年，星期四例7某服务系统有三个服务员，输入流为泊松流，服务时间为指数分布。每个顾客服务时间等于0.5分钟。求在T=4分钟内被服务顾客的数学期望。随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用解：设表示第个顾客到达的时刻，表示第个顾客到达的时间间隔。则有。其SAS程序如下：

本章目录60第六十页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用dataex7;s1=0;s2=0;s3=0;T=0;Dountil(T>4);R=ranuni(-1);P=0.2*(-log(R));T=T+p;Ss1=s1;ss2=s2;ss3=s3;If(T<=ss1)and(T<=ss2)and(T<=ss3)thend=0;Elsed=1;If(T>ss1)thendo;s1=T+0.5;end;If(T<=ss1)and(T<=ss2)and(T>ss3)thendo;s3=T+0.5;end;Output;End;Run;Procmeansdata=ex7;Vard;Outputout=resultsum=dsum;Run;Procprintdata=result;Run;本章目录61第六十一页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

4随机模拟方法在随机服务系统中的应用此过程共进行6次模拟，可得其平均值为

=16，即在4分钟内平均服务了16个顾客。本章目录62第六十二页，共七十七页，编辑于2023年，星期四随机模拟方法不仅在求解确定性和随机性复杂系统的问题，它在理论研究方面也大为可有。比如有些问题从理论上已经得出了圆满的结论，但因没有经过实践验证比较，暂时没有被应用。这时若使用随机模拟方法先反复加以比较验证，再用于实践中就更可靠了。还有些问题，从理论上证明很困难，而科学家从其他方面的知识及经验，对所研究问题有某些猜想，这时随机模拟方法就是一个有效可行的方法。下面仅举例说明用随机模拟的方法在比较系统聚类方法上的应用随机数的产生与模拟

5随机模拟方法在理论研究中的应用本章目录63第六十三页，共七十七页，编辑于2023年，星期四例8随机数的产生与模拟

5随机模拟方法在理论研究中的应用假设数据来自和的总体，用SAS来计算，比较系统聚类法的八种常用方法在分类时之间的分类效果的好坏。本章目录64第六十四页，共七十七页，编辑于2023年，星期四例8解：基本思想为：随机数的产生与模拟

5随机模拟方法在理论研究中的应用本章目录65第六十五页，共七十七页，编辑于2023年，星期四例8解：基本思想为：随机数的产生与模拟

5随机模拟方法在理论研究中的应用2.用八种常用的系统聚类方法对容量为2n个样品的数据进行聚类，计算各种聚类方法的错分率(即判错个数所占的比例)(j=1,2,…,8);本章目录66第六十六页，共七十七页，编辑于2023年，星期四例8解：基本思想为：随机数的产生与模拟

5随机模拟方法在理论研究中的应用3.重复以上两步N次，得(j=1,

…,8;i=1,…,N),计算平均错分率：

(j=1,2,…,8)。本章目录67第六十七页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

5随机模拟方法在理论研究中的应用其SAS程序为(以average为例)：%macrocreatedata(mdata=,leixing=,mv1=,mv2=,mvar1=,mvar2=,mvar3=);data&mdata;dropiu1u2;fenlei=&leixing.;doi=1to50;u1=rannor(0);u2=rannor(0);x1=&mv1.+sqrt(&mvar1.)*u1;x2=&mv2.+(&mvar2.*u1+sqrt(&mvar1.*&mvar3.-&mvar2.*&mvar2.)*u2)/sqrt(&mvar1.);output;end;/*产生来自两元正态总体的随机数据*/run;%mendcreatedata;

本章目录68第六十八页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

5随机模拟方法在理论研究中的应用其SAS程序为(以average为例)：%macrodatacluster(mdata=,method=);data&mdata.;setab;run;procclusterdata=&mdata.method=&method.outtree=cnoprint;varx1x2;copyfenleix1x2;run;/*对两个来自不同两元正态总体的随机数据进行聚类*/proctreedata=cout=abcncl=2noprint;copyfenleix1x2;run;dataresult1;setabc;result=0;iffenlei^=clusterthenresult=1;run;procsortdata=result1;byfenlei;run;procmeansdata=result1noprint;varresult;byfenlei;outputout=resultsum=errorsum;run;/*计算出错分的个数*/procappendbase=r_resultdata=result;run;%menddatacluster;

本章目录69第六十九页，共七十七页，编辑于2023年，星期四随机数的产生与模拟

5随机模拟方法在理论研究中的应用其SAS程序为(以average为例)：%macro

analyze;%doi=1

%to

50;%createdata(mdata=a,leixing=1,mv1=0,mv2=0,mvar1=1,mvar2=0,mvar3=1);%createdata(mdata=b,leixing=2,mv1=3,mv2=3,mvar1=1,mvar2=0,mvar3=1);%datacluster(mdata=ab,method=average);%end;%mendanalyze;%analyze;datarr;setr_result;iferrorsum>25

thenerrorsum=50-errorsum;errorratio=errorsum/50;run;/*计算错分率*/proc

sort

data=rr;byfenlei;run;proc

means

data=rrnoprint;output

out=rmean=err_ratio;varerrorratio;byfenlei;run;/*计算平均错分率*/proc

print

data=r;varfenleierr_ratio;title

'总错判率:';r