SAS软件与数学建模.pptx

《SAS软件与数学建模.pptx》由会员分享，可在线阅读，更多相关《SAS软件与数学建模.pptx（129页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、统计统计生物物理化学经济保险心理教育管理第1页/共129页内容简介模块简介模块简介SAS与与Excel的通讯的通讯常用过程常用过程SAS的学习的学习曲线拟合非线性规划神经网络多元统计矩阵运算回归分析第2页/共129页模块简介模块简介SAS系统的组成SAS数据库部分：其模块为SAS/BASE。它也是其它模块的基础，即其它模块是建立在其基础之上的SAS分析核心：这一部分包括了许多模块：SAS/STAT：统计分析模块：回归分析、方差分析、属性数据分析、多变量分析、判别和聚类分析、残存分析、心理测验分析和非参数分析等8类40多个过程。SAS/ETS：经济预测或时间序列分析模块。如实用预测（逐步回归、指

2、数平滑等）序列相关校正回归、分布滞后回归、ARIMA模型、状态空间方法、谱分析和互谱分析等。SAS/OR：运筹学和工程管理模块：可进行线性和非线性规划，还包括项目管理，时间安排和资源分配等问题的一整套方法。SAS/QC：质量控制和试验设计模块。SAS/IML：矩阵运算模块。SAS/LAB：菜单驱动的面向任务的解释引导式数据分析模块。SAS/INSIGHT：可视化数据探索工具模块。SAS/SPECTRAVIEW：多维数据观测、分析、研究的交互式立体可视化工具模块第3页/共129页模块简介模块简介SAS系统的组成SAS开发及呈现工具：SAS/AF：应用开发工具。采用面向对象的技术，开发用户自己的图

3、形用户界面（GUI）的应用系统。SAS/EIS：行政管理系统或个人的信息系统SAS/GRAPH：图形软件包SAS/GIS：集地理信息系统功能与空间数据的显示分析于一体的软件SAS对分布处理模式的支持及其数据仓库设计：SAS/ACCESS：与外部数据库文件的接口模块。SAS/CONNECT：在网络环境下，使各平台上的SAS系统建立内在联系模块。实现分布处理，从而有效地利用各平台数据和机器资源SAS/SHARE：实行SAS系统中数据库的并发控制的模块第4页/共129页SAS与与Excel的通讯的通讯SAS中数据Excel第5页/共129页SAS系统基础1 1 概述概述1.1SAS1.1SAS显示管

4、理系统显示管理系统1.2SAS1.2SAS数据库与数据集数据库与数据集2 SAS2 SAS编程基础编程基础2.12.1一个简单一个简单SASSAS程序程序2.2SAS2.2SAS变量变量2.3SAS2.3SAS表达式表达式2.4SAS2.4SAS函数函数2.4.1SAS2.4.1SAS函数概述函数概述2.4.22.4.2常用常用SASSAS函数函数3 DATA3 DATA步步3.1DATA3.1DATA步流程步流程3.2INPUT3.2INPUT语句语句3.33.3数据集的整理数据集的整理3.43.4数据集的连接与合并数据集的连接与合并3.5SAS3.5SAS与外部数据的交换与外部数据的交换第

5、6页/共129页StatisticalAnalysisSystem(缩写为SAS)：是用于数据分析与决策支持的大型集成式模块化软件系统在数据处理和统计分析领域，SAS系统被誉为国际上的标准软件系统第7页/共129页四大任务：数据访问数据访问访问任何形式和来源的数据访问任何形式和来源的数据数据管理数据管理将数据转变为可用形式将数据转变为可用形式数据呈现数据呈现以文字、图形和表格的形式展现数据信息以文字、图形和表格的形式展现数据信息数据分析数据分析获取有意义的信息获取有意义的信息第8页/共129页1.SAS1.SAS显示管理系统显示管理系统 SASSAS系统的启动系统的启动点击图标点击图标，启动，

6、启动SASSAS系统系统命令窗口命令窗口菜单栏菜单栏工具栏工具栏LOGLOG窗窗程序编辑窗程序编辑窗状态栏状态栏结果窗结果窗资源管理器资源管理器结果输出窗结果输出窗演示演示第9页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据库数据库:主要用于存放数据集和目录主要用于存放数据集和目录.SASSAS系统系统库库操作系统操作系统文件夹文件夹第10页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据库数据库:主要用于存放数据集和目录主要用于存放数据集和目录.临时库临时库(WORK):SAS(WORK):SAS系统的临时工作目录系统的临时工作目录,

7、退出系统时退出系统时,其中的内容全部删除其中的内容全部删除永久库永久库:与临时库相反与临时库相反.SAS.SAS提供提供SASUSER,SASHELPSASUSER,SASHELP第11页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据库数据库:主要用于存放数据集和目录主要用于存放数据集和目录.创建用户永久库创建用户永久库vv 用操作系统中的资源管理器用操作系统中的资源管理器,建立一文件夹建立一文件夹程序方式程序方式:libname:libname库名库名文件夹所在路径文件夹所在路径 菜单方式菜单方式:点击图标点击图标,然后按要求操作然后按要求操作.vv在在SASS

8、AS系统中建库系统中建库第12页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据库数据库:主要用于存放数据集和目录主要用于存放数据集和目录.创建用户永久库创建用户永久库例：在例：在SASSAS中先用程序方式建立一名为中先用程序方式建立一名为COURSECOURSE的用户库；再用菜单方式建的用户库；再用菜单方式建立一名为立一名为COURSE1COURSE1的用户库。的用户库。操作操作第13页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据集数据集:主要用于存放数据文件。它一般由主要用于存放数据文件。它一般由描述部分描述部分和和数据部分数据部

9、分组成，也有一些数据集包含组成，也有一些数据集包含索引部分索引部分。数据部分数据部分是存放数据的地方，它组织成一张是存放数据的地方，它组织成一张二维表二维表，其列，其列表示变量，行表示观测。表示变量，行表示观测。描述部分说明了数据集的一般信息，如数据集的名字和成员描述部分说明了数据集的一般信息，如数据集的名字和成员类型；建立的日期和时间；观测数以及变量个数、属性等。类型；建立的日期和时间；观测数以及变量个数、属性等。第14页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据集数据集:数据部分数据部分是存放数据的地方，它组织成一张是存放数据的地方，它组织成一张二维表二维

10、表，其列，其列表示变量，行表示观测。表示变量，行表示观测。观测观测变量变量第15页/共129页2.SAS2.SAS数据库与数据集数据库与数据集SASSAS数据集的引用数据集的引用:库名数据集名园点连接符园点连接符第16页/共129页Databodyfat;Databodyfat;inPUtsex$fatpct;inPUtsex$fatpct;fat=fatpct/100;fat=fatpct/100;cards;cards;M13.3F22M19F26M20F16M8F12M18F21.7M13.3F22M19F26M20F16M8F12M18F21.7M22F23.2M20F21M31F28

11、M21F30M12F23M22F23.2M20F21M31F28M21F30M12F23M16M12M24M16M12M24RUNRUN；PROCmeansdata=bodyfat;PROCmeansdata=bodyfat;varfatpct;varfatpct;run;run;DATADATA步步程程序序步步1 1、一个简单、一个简单SASSAS程序程序第17页/共129页该程序的结构包括两部分：该程序的结构包括两部分：以关键词以关键词DATADATA开始的开始的DATADATA步；主要创建步；主要创建数据集数据集,如此例中的如此例中的bodyfatbodyfat，同时也可以产，同时也可以

12、产生一些描述性统计和报表输出；生一些描述性统计和报表输出；以关键词以关键词PROCPROC开头开头PROCPROC步。步。PROCPROC步步则是系统预先写好的执行特殊任务的则是系统预先写好的执行特殊任务的SASSAS子子程序程序,如此例中的如此例中的meansmeans过程过程，它主要是承，它主要是承担数据的分析，同时也可产生相关的数据集。担数据的分析，同时也可产生相关的数据集。1 1、一个简单、一个简单SASSAS程序程序第18页/共129页从该程序可以看出从该程序可以看出SASSAS程序的特点：程序的特点：vvSASSAS程序是由程序是由SASSAS关键词、关键词、SASSAS名字、特殊

13、字符（如例子中的名字、特殊字符（如例子中的$号）和运算符组成的字符串，并以分号（；）结尾的号）和运算符组成的字符串，并以分号（；）结尾的SASSAS语句组成，语句组成，它要求它要求SASSAS系统执行一种操作或给系统执行一种操作或给SASSAS系统提供信息。系统提供信息。vv SASSAS程序不区分大小写，语句位置任意，一条语句的书写既可起程序不区分大小写，语句位置任意，一条语句的书写既可起始于任何行、列，也可跨越多行，同时也可以将多条语句写在一行。始于任何行、列，也可跨越多行，同时也可以将多条语句写在一行。不过一个好的习惯是；不过一个好的习惯是；DATADATA步和步和PROCPROC步顶格

14、书写，在每个步骤内步顶格书写，在每个步骤内的语句以缩进的方式书写，这样阅读和理解起来方便。的语句以缩进的方式书写，这样阅读和理解起来方便。1 1、一个简单、一个简单SASSAS程序程序操作操作第19页/共129页SASSAS变量的命名：第一个字符必须是字母或下划线变量的命名：第一个字符必须是字母或下划线(_);(_);后面可跟后面可跟1 1至至3232个字符（数字、字母或下划线）个字符（数字、字母或下划线）变量的属性：变量的属性：类型类型(type)(type)：数值型、字符型：数值型、字符型(定义时变量后带定义时变量后带$)$)长度长度(length)(length)：存储该变量所用的字节数

15、（缺省为存储该变量所用的字节数（缺省为8 8）输入格式输入格式(informat)/(informat)/输出格式输出格式(format)(format)：标签标签(label)(label)：对变量的作用给予描述。：对变量的作用给予描述。2 2、SASSAS变量变量第20页/共129页输入输入/输出格式：输出格式：类型：类型：字符型：用于字符变量的现成格式字符型：用于字符变量的现成格式数值型：用于一般数值型变量的现成格式数值型：用于一般数值型变量的现成格式日期时间型：用于日期和时间的现成格式日期时间型：用于日期和时间的现成格式用户定义：用于用户自定义格式用户定义：用于用户自定义格式2 2、S

16、ASSAS变量变量第21页/共129页输入输入/输出格式：输出格式：一般形式：一般形式：格式名格式名.2 2、SASSAS变量变量字符字符型变型变量量宽度宽度小数小数点后点后数值数值位数位数第22页/共129页输入输入/输出格式：输出格式：字符格式：字符格式：$w.$w.$CHARw.$CHARw.2 2、SASSAS变量变量其含义是读入或显示字符变量的前其含义是读入或显示字符变量的前w w个字节。其缺省值个字节。其缺省值均为变量的长度。当变量的长度未定时，均为变量的长度。当变量的长度未定时，$CHARw.$CHARw.格式格式宽度的缺省值为宽度的缺省值为8 8；对；对$w.$w.作为输出格式

17、时缺省宽度为作为输出格式时缺省宽度为1 1，作为输入格式时，必须提供作为输入格式时，必须提供w w的值。的值。二者作为输出格式时，均保留字符串前的空格；二者作为输出格式时，均保留字符串前的空格；但作为输入格式时，但作为输入格式时，$CHARw.$CHARw.保留字符串前的空格，而保留字符串前的空格，而$w.$w.则删除字符前的空格。则删除字符前的空格。第23页/共129页输入输入/输出格式：输出格式：数值输入格式：数值输入格式：w.dw.d数值输出格式：数值输出格式：BESTw.BESTw.2 2、SASSAS变量变量BESTw.BESTw.是缺省输出格式。是缺省输出格式。WW的缺省值为的缺省

18、值为1212，最大可取，最大可取3232。含义为将数值在规定的宽度下用最为精确有效的方式显示。含义为将数值在规定的宽度下用最为精确有效的方式显示。若用定点方式表示可提供更多的有效数字就使用定点方式；若用定点方式表示可提供更多的有效数字就使用定点方式；否则自动使用浮点数。否则自动使用浮点数。如：如：数据集中数值数据集中数值输出格式输出格式输出表示输出表示1257000best6.1.26E61257000best6.1.26E61257000best3.1E61257000best3.1E6第24页/共129页输入输入/输出格式：输出格式：数值输入格式：数值输入格式：COMMAw.dCOMMAw

19、.d2 2、SASSAS变量变量如：如：原始数据原始数据输入格式输入格式记入数值记入数值$1,2345,6789comma12.21234567.89$1,2345,6789comma12.21234567.89$1234,5678.9comma12.212345678.9$1234,5678.9comma12.212345678.9($1,2345,6789)comma14.2-1234567.89($1,2345,6789)comma14.2-1234567.89第25页/共129页输入输入/输出格式：输出格式：数值输出格式：数值输出格式：COMMAw.d&DOLLARw.dCOMMAw.

20、d&DOLLARw.d2 2、SASSAS变量变量如：如：原始数据原始数据输出格式输出格式输出表示输出表示123456.7comma12.2123,456.70123456.7comma12.2123,456.70-123456.7comma12.2-123,456.70-123456.7comma12.2-123,456.70123456.7dollar12.2$123,456.70123456.7dollar12.2$123,456.70-123456.7dollar12.2-123,456.70-123456.7dollar12.2-123,456.70第26页/共129页输入输入/输出

21、格式输出格式(日期日期)：日期、时间的处理：日期、时间的处理：将日期和时间值均作为将日期和时间值均作为数值变量数值变量存储存储以以19601960年年1 1月月1 1日日为为0 0点，其它日期值为该日期值与点，其它日期值为该日期值与19601960年年1 1月月1 1日日相差的天数。相差的天数。时间以时间以19601960年年1 1月月1 1日日0 0时时0 0分分0 0秒秒为起点（为起点（0 0点），其余时间为与其点），其余时间为与其的差值。的差值。2 2、SASSAS变量变量第27页/共129页输入输入/输出格式输出格式(日期日期)：I.I.输入格式输入格式:YYMMDDw.:YYMMDD

22、w.输出格式输出格式:YYMMDDxw.:YYMMDDxw.源数据日期源数据日期 输入格式输入格式 SASSAS日期值日期值 输出格式输出格式 输出表示输出表示99-10-15 yymmdd8.14532 yymmddd10.1999-10-1599-10-15 yymmdd8.14532 yymmddd10.1999-10-1599/10/15 ddmmyy8.14532 yymmdds10.1999/10/1599/10/15 ddmmyy8.14532 yymmdds10.1999/10/15151099 ddmmyy6.14532 mmddyyp10.10.15.1999151099

23、ddmmyy6.14532 mmddyyp10.10.15.199910-15-99 mmddyy8.14532 ddmmyyn8.1510199910-15-99 mmddyy8.14532 ddmmyyn8.151019992 2、SASSAS变量变量其中输出格式中的其中输出格式中的x x的取值为的取值为:x:BCDNPSx:BCDNPS含义含义空格空格:-:-无无././第28页/共129页输入输入/输出格式输出格式(日期日期)：II.II.输入输出格式输入输出格式:DATEw.:DATEw.源数据日期源数据日期 输入格式输入格式 SASSAS日期值日期值 输出格式输出格式 输出表示输出

24、表示15OCT99 date7.14532 date9.15OCT199915OCT99 date7.14532 date9.15OCT19992 2、SASSAS变量变量第29页/共129页输入输入/输出格式输出格式(时间时间)：III.III.输入格式输入格式:TIMEw.:TIMEw.和输出格式和输出格式:TIMEw.d:TIMEw.d hh:mm:ss.ss(hh:mm:ss.ss(时时:分分:秒秒.百分秒百分秒)源数据日期源数据日期 输入格式输入格式 SASSAS日期值日期值 输出格式输出格式 输出表示输出表示13:40:35.37 time11.49235.37 time11.2

25、13:40:35.3713:40:35.37 time11.49235.37 time11.2 13:40:35.3713 40 35.37 time11.49235.37 Time2.1313 40 35.37 time11.49235.37 Time2.1313/40/35.37 time11.49235.37 Time5.13:4013/40/35.37 time11.49235.37 Time5.13:4013-40-35.37 time11.49235.37 Time8.13:40:3513-40-35.37 time11.49235.37 Time8.13:40:351 40 35

26、.37 PM time14.49235.37 Timeampm13.2 1:40:35.37 PM1 40 35.37 PM time14.49235.37 Timeampm13.2 1:40:35.37 PM2 2、SASSAS变量变量第30页/共129页输入输入/输出格式输出格式(日期时间日期时间)：IV.IV.输入格式输入格式:DATETIMEw.:DATETIMEw.ddmmmyy hh:mm:ss.ss ddmmmyy hh:mm:ss.ss 或或 ddmmmyyyy hh:mm:ss.ssddmmmyyyy hh:mm:ss.ss 输出格式输出格式:DATETIMEw.d:DATE

27、TIMEw.d ddmmmyy:hh:mm:ss.ss ddmmmyy:hh:mm:ss.ss 或或 ddmmmyyyy:hh:mm:ss.ssddmmmyyyy:hh:mm:ss.ss 2 2、SASSAS变量变量第31页/共129页变量列表：变量列表：x1x2x3xnx1x2x3xn简写为简写为x1-xnx1-xnABCDABCD简写为简写为A-DA-DVXTDUPJPPVXTDUPJPP简写为简写为V-JPPV-JPP2 2、SASSAS变量变量第32页/共129页表达式是由运算符和运算对象形成的指令集表达式是由运算符和运算对象形成的指令集3 3、SASSAS表达式表达式1.SAS1.S

28、AS常数常数数值常数数值常数:110-0.351.8E30.5E-10:110-0.351.8E30.5E-10字符常数字符常数:Smith23JohnSmith:Smith23JohnSmith日期日期,时间和日期时间常数时间和日期时间常数:13JUL2001d:13JUL2001d22:23t13JUL2001:22:23:25DT22:23t13JUL2001:22:23:25DT2.SAS2.SAS运算符运算符:算术运算符算术运算符:*(:*(乘方乘方)*/+-)*/+-比较运算符比较运算符:=(eq)=(ne)(gt)=(ge)(gt)=(ge)(取小取小)()(取大取大)|()|(

29、连接连接)第33页/共129页表达式是由运算符和运算对象形成的指令集表达式是由运算符和运算对象形成的指令集3 3、SASSAS表达式表达式3.3.运算次序运算次序:1)1)括号内的表达式先算括号内的表达式先算2)2)优先级优先级(见下面见下面)3)3)相同优选级别的算符相同优选级别的算符,按顺序按顺序,先左后右先左后右第34页/共129页表达式是由运算符和运算对象形成的指令集表达式是由运算符和运算对象形成的指令集3 3、SASSAS表达式表达式优先级优先级第一级第一级*乘方乘方+仅作前缀时仅作前缀时(取正取正)-仅作前缀时仅作前缀时(取负取负)非非最小最小最大最大第二级第二级*乘乘除除第三级第

30、三级+加加-减减第四级第四级|连接连接第五级第五级,=,=,第六级第六级&与与第七级第七级|第35页/共129页1.SAS1.SAS函数的调用函数的调用:functionname(argument,)-functionname(argument,)-括号不能少括号不能少自变量的表示形式有两种自变量的表示形式有两种:如如x,y,zx,y,z的和的和sum(x,y,z)sum(x,y,z)或或sum(ofxyz)sum(ofx-z)sum(ofxyz)sum(ofx-z)A=sum(x,y,z)A=sum(x,y,z)B=min(sum(ofxyz),1000)B=min(sum(ofxyz),1

31、000)4 4、SASSAS函数函数第36页/共129页2.SAS2.SAS函数的返回值函数的返回值:返回值为返回值为:数值型数值型字符型字符型其返回值的类型一般可通过自变量的类型加以判断其返回值的类型一般可通过自变量的类型加以判断(但有个别有例但有个别有例外外,如如PUTPUT函数函数,总是输出字符型总是输出字符型)4 4、SASSAS函数函数3.SAS3.SAS函数的分类函数的分类:SASSAS函数分为十七类：算术函数、数组函数、截取函数、数学函数、函数分为十七类：算术函数、数组函数、截取函数、数学函数、三角和双曲函数、概率函数、分位数函数、非中心函数、样本统计函数、三角和双曲函数、概率函

32、数、分位数函数、非中心函数、样本统计函数、随机数函数、财政金融函数、逐位逻辑操作函数、数字函数、字符函数、随机数函数、财政金融函数、逐位逻辑操作函数、数字函数、字符函数、日期和时间函数、洲和日期和时间函数、洲和ZipZip码换算函数和特殊函数。码换算函数和特殊函数。第37页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数1.1.概率函数概率函数:SASSAS系统提供了系统提供了1010个常见分布的概率函数，其定义为：个常见分布的概率函数，其定义为：在下面的函数均表示求这些分布函数的值。在下面的函数均表示求这些分布函数的值。其中其中DFDF表示自由度，

33、表示自由度，NCNC表示非中心参数值，表示非中心参数值，NDFNDF表示分表示分子自由度，子自由度，DDFDDF表示分母自由度，其余参数为相关分布的表示分母自由度，其余参数为相关分布的参数项。参数项。第38页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数1.1.概率函数概率函数:标准正态分布函数标准正态分布函数PROBNORMPROBNORM（X X）计算计算PUX(PU0PROBGAM(X,a)a0贝塔分布函数贝塔分布函数PROBBETA(X,a,b)XPROBBETA(X,a,b)X(0,1),a,b0(0,1),a,b0F F分布函数分布函数

34、PROBF(X,NDF,DDF,NC)PROBF(X,NDF,DDF,NC)t t分布函数分布函数PROBT(X,DF,NC)PROBT(X,DF,NC)二项分布的概率分布函数二项分布的概率分布函数PROBBNML(P,N,M)PROBBNML(P,N,M)波松分布波松分布POISSON(,N)POISSON(,N)负二项分布负二项分布PROBNEGB(P,N,M)PROBNEGB(P,N,M)超几何分布超几何分布PROBHYPR(NN,K,N,X,OR)PROBHYPR(NN,K,N,X,OR)表示表示NNNN件产品件产品,其中有其中有K K件件不合格品不合格品,随机地从随机地从NNNN件中

35、抽取件中抽取N N件件,这这N N件中不合格品的个数服从该分布件中不合格品的个数服从该分布.该函数给出该函数给出N N件品中不合格品的个数件品中不合格品的个数=X=X的概率的概率.参数参数OROR是不匀率是不匀率.K=0.K=0是是有效的值有效的值(即即NNNN件产品中没有不合格品件产品中没有不合格品),),由此导出由此导出X=0,X=0,不管其它参数取什么不管其它参数取什么值值,该函数获得值该函数获得值1.1.如果省略如果省略OR,OR,就认为其值为就认为其值为1.1.如如:X=PROBHYPR(10,5,3,2,1.5):X=PROBHYPR(10,5,3,2,1.5)OR=1.5OR=1

36、.5表示抽到不合格品的概率是抽到合格品概率的表示抽到不合格品的概率是抽到合格品概率的1.51.5倍倍.第39页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数1.1.概率函数概率函数:举例举例:DATADD;DATADD;p=PROBCHI(5.6,4);/*p=PROBCHI(5.6,4);/*求自由度为求自由度为4 4的的 2 2分布不超过分布不超过5.65.6的概率值的概率值p*/p*/T=1-2*(1-PROBT(2.4,8);/*T=1-2*(1-PROBT(2.4,8);/*求自由度为求自由度为8 8的的P|t|2.4P|t|2.4*/*

37、/RUN;RUN;在在SASSAS系统中操作系统中操作第40页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数2.2.分位数函数分位数函数:2 2分布的分位数分布的分位数CINVCINV（p p，dfdf，ncnc）nc0nc0）BetaBeta分布的分位数分布的分位数BETAINV(p,a,b)BETAINV(p,a,b)F F分布的分位数分布的分位数FINV(p,ndf,ddf,nc)FINV(p,ndf,ddf,nc)T T分布的分位数分布的分位数TINV(p,df,nc)TINV(p,df,nc)正态分布的分位数正态分布的分位数PROBIT(

38、p)PROBIT(p)标准正态分布的分位数标准正态分布的分位数,其结果在其结果在-5-5和和5 5之间之间.GammaGamma分布的分位数分布的分位数GAMINV(p,a)GAMINV(p,a)第41页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数2.2.分位数函数分位数函数:举例举例:DATAA;DATAA;Q1=FINV(0.025,2,10);/*Q1=FINV(0.025,2,10);/*求自由度为求自由度为(2,10)(2,10)的的F F分布的分布的2.5%2.5%分位数分位数Q1Q1*/*/Q2=FINV(0.975,2,10);/

39、*Q2=FINV(0.975,2,10);/*求自由度为求自由度为(2,10)(2,10)的的F F分布的分布的97.5%97.5%分位数分位数Q2Q2*/*/Q3=FINV(0.95,2,10,3.2);/*Q3=FINV(0.95,2,10,3.2);/*非中心值为非中心值为3.23.2的的F F分布的分布的95%95%的分位数的分位数,其自由度为其自由度为(2,10)(2,10)*/*/RUN;RUN;演示演示第42页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数3.3.随机数函数随机数函数:均匀分布随机数均匀分布随机数UNIFORMUNIF

40、ORM（seedseed）其乘子为其乘子为1680716807，模为，模为2 23131的乘同余发生器和一个的乘同余发生器和一个6464位数的搅乱表形成的组合发生器，位数的搅乱表形成的组合发生器，seedseed必须是常数，它或是必须是常数，它或是0 0或或5 5位，位，6 6位，位，7 7位的奇数。位的奇数。RANUNIRANUNI（seedseed）其乘子为其乘子为397204094397204094，模为，模为2 23131-1-1的素数发生器，的素数发生器，seedseed必须是小于模必须是小于模2 23131-1-1任何常数。任何常数。标准正态分布随机数标准正态分布随机数NORMAL

41、NORMAL（seedseed）RANNORRANNOR（seedseed）一般正态分布一般正态分布N N（,2 2）可由如下变换得到：）可由如下变换得到：第43页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数3.3.随机数函数随机数函数:指数分布随机数指数分布随机数 RANEXP(seed)RANEXP(seed)产生产生=1=1的指数分布的随机数的指数分布的随机数.Y=Y=RANEXP(seed)/RANEXP(seed)/,则产生则产生一般指数分布一般指数分布随机数随机数若若Y=Y=-*LOG(RANEXP(seed),-*LOG(RANEX

42、P(seed),则则Y Y为具有位置参数为具有位置参数 和尺度参数为和尺度参数为 的的极值分布极值分布随机变量随机变量若若Y=Y=FLOOR(-RANEXP(seed)/LOG(p)FLOOR(-RANEXP(seed)/LOG(p),则则Y Y为具有参数为具有参数p p的的几何分布几何分布变量变量第44页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数3.3.随机数函数随机数函数:伽马分布随机数伽马分布随机数 RANGAM(seed,)RANGAM(seed,).它用组合法产生参数为非整数它用组合法产生参数为非整数 的的GAMMAGAMMA随机数随

43、机数.非整数非整数 可表示为整数部分和小数部分的和可表示为整数部分和小数部分的和:=+(-).:=+(-).用逆变换方法产生参数为用逆变换方法产生参数为整数整数的的GAMMAGAMMA随机数随机数X1,X1,用舍选法产生参数为小数用舍选法产生参数为小数(-)(-)的的GAMMAGAMMA随机随机数数X2,X2,那么那么X1+X2X1+X2即为参数为即为参数为 的的GAMMAGAMMA随机数随机数设设X=RANGAM(seed,),X=RANGAM(seed,),则则Y=X/Y=X/为形状参数为为形状参数为,尺度参数为尺度参数为 的的GAMMAGAMMA分分布随机数布随机数;若若2*2*为整数为

44、整数,且且Y=2*X,Y=2*X,则则Y Y是具有自由度为是具有自由度为2*2*的卡方分布随机数的卡方分布随机数;若若=n(=n(正整数正整数)且且Y=X/,Y=X/,那么那么Y Y是是Erlang(Erlang(厄兰厄兰)分布随机数分布随机数.它是它是n n个独立的均值个独立的均值为为 的指数分布变量的和的指数分布变量的和;若若Y1=RANGAM(seed,),Y2=RANGAM(seed,),Y1=RANGAM(seed,),Y2=RANGAM(seed,),则则Y=Y1/(Y1+Y2)Y=Y1/(Y1+Y2)是参数为是参数为 和和 的贝塔分布随机数的贝塔分布随机数第45页/共129页4.

45、4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数3.3.随机数函数随机数函数:三角分布随机数三角分布随机数RANTRIRANTRI（seedseed，h h）0h10h1其密度函数为其密度函数为:的三角分布随机数的三角分布随机数.对一般区间对一般区间a,ba,b上的三角形分布随机数上的三角形分布随机数Y,Y,可由可由0,10,1区间上三角形分布的随机数的区间上三角形分布的随机数的线性变换得到线性变换得到.即即:Y=(b-a)*RANTRI(seed,h)+ah=(c-a)/(b-a)c:Y=(b-a)*RANTRI(seed,h)+ah=(c-a)/(b-a)c a

46、,ba,b第46页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数3.3.随机数函数随机数函数:柯西分布随机数柯西分布随机数RANCAU(seed)RANCAU(seed)用舍选法产生用舍选法产生Y=+*RANCAU(seed),Y=+*RANCAU(seed),则则Y Y是位置参数为是位置参数为,尺度参数为尺度参数为 的一般柯西分布随机数的一般柯西分布随机数 二项分布随机数二项分布随机数RANBIN(seed,n,p)RANBIN(seed,n,p)波松分布随机数波松分布随机数RANPOI(seed,)RANPOI(seed,)离散分布随机数离散分

47、布随机数 RANTBL(seed,p1,p2,pRANTBL(seed,p1,p2,pn n)此函数产生一个由此函数产生一个由p1p1到到p pn n定义的概率定义的概率密度函数的变量密度函数的变量.它返回数值它返回数值1,2,3,n,1,2,3,n,且具以下概率分布且具以下概率分布:123n123np1p2p3pnp1p2p3pn第47页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数3.3.随机数函数随机数函数:vv关于种子数关于种子数(seed):(seed):当当SEED0SEED0时时,则取系统当前的日期时间值作为初始值则取系统当前的日期时

48、间值作为初始值vv随机数函数每次执行时虽然当前的种子不断变化随机数函数每次执行时虽然当前的种子不断变化,但自变量但自变量SEEDSEED的值保持不变的值保持不变.因此不可能控制初始化之后的随机数因此不可能控制初始化之后的随机数.若想控制随机数流若想控制随机数流,可使用可使用CALLCALL子程序子程序 Callranbin(seed,n,p,x)Callranbin(seed,n,p,x)Callrancau(seed,x)Callrancau(seed,x)Callranexp(seed,x)Callranexp(seed,x)Callrannor(seed,x)Callrannor(see

49、d,x)Callranpoi(seed,x)Callranpoi(seed,x)Callrantbl(see,p1,p2,pn,x)Callrantbl(see,p1,p2,pn,x)Callrantri(seed,h,x)Callrantri(seed,h,x)Callranuni(seed,x)Callranuni(seed,x)Callrangam(seed,a,x)Callrangam(seed,a,x)第48页/共129页4.4.与统计有关的与统计有关的SASSAS函数函数4 4、SASSAS函数函数4.4.样本统计函数样本统计函数:均值均值 MEANMEAN（OF X1-XnOF

50、X1-Xn）或或 MEAN(X,Y,Z,)(MEAN(X,Y,Z,)(注注:缺失值不含在内缺失值不含在内)最大值最大值 MAX(OF X1-Xn)MAX(OF X1-Xn)或或 MAX(X,Y,Z,)(MAX(X,Y,Z,)(注注:缺失值不含在内缺失值不含在内)最小值最小值 MIN(OF X1-Xn)MIN(OF X1-Xn)或或 MIN(X,Y,Z,)(MIN(X,Y,Z,)(注注:缺失值不含在内缺失值不含在内)非缺失值的个数非缺失值的个数 N(OF X1-Xn)N(OF X1-Xn)或或 N(X,Y,Z,)N(X,Y,Z,)缺失数据的个数缺失数据的个数 NMISS(OF X1-Xn)NMI