概率统计及SAS应用程序

《概率统计及SAS应用程序》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、ＳＡＳ应用讲义7/22/20211StatisticalAnalysisSystem简称为SAS,可用来分析数据和编写报告.它是美国SAS研究所的产品,在国际上被誉为标准软件,在我国深受医学、农林、财经、社会科学、行政管理等众多领域的专业工作者的好评。有关ＳＡＳ的最新信息，可以查看http://www.sas.com.cn。ＳＡＳ采用积木式模块结构，其中的ＳＡＳ/ＳＴＡＴ模块是目前功能最强的多元统计分析程序集,可以做回归分析、聚类分析、判别分析、主成分分析、因子分析、典型相关分析(下学期介绍)以及各种试

2、验设计的方差分析和协方差分析。本讲义围绕SAS的应用,讲述以下六部分内容：（１）ＳＡＳ应用基础；（２）ＳＡＳ常用语句；（３）ＳＡＳ服务过程；（４）描述性统计程式；（５）方差分析程式；（６）回归分析程式；7/22/20212ＳＡＳ的显示管理系统启动计算机,点击SAS图标后,即可进入SAS的显示管理系统.在View中有四个主要的窗口(其他的先不考虑)：(1)编辑窗口(programeditor):编辑程式和数据文件;(2)日志窗口(log):记录运行情况,显示error信息;(3)输出窗口(output):

3、输出运行的结果;(4)图形窗口(graph):输出图形.点击View菜单中的Programeditor、Log、Output、Graph命令可以进入编辑、日志、输出及图形窗口.按功能键F5、F6、F7也可以进入编辑、日志及输出窗口.退出SAS有两种方法：(1)点击File菜单中的Exit命令；(2)点击窗口右上角的×。7/22/20213概率统计及SAS应用教材中的程序应用SAS计算二项分布的概率，请注意SAS中probbnml(p,n,k)=P(Xk)=因此，当n=5，k=3，p=0.2时，应用S

4、AS直接计算P{X=3}的程序为：dataprobnml;p=probbnml(0.2,5,3)-probbnml(0.2,5,2);procprint;run;输出的结果为：0.0512。7/22/20214当n=5,k=4,p=0.8时,应用SAS直接计算P(X=4)+P(x=5)的程序为：dataex;p=1-probbnml(0.8,5,3);procprint;run;输出的结果为：0.73728。应用SAS直接计算例1.3.1中所求概率的P{8≤X≤12}的程序为：dataex;p=prob

5、bnml(0.5,20,12)-probbnml(0.5,20,7);procprint;run;输出的结果为：0.7368240356。7/22/20215应用SAS中的probnorm(x)近似计算二项分布的概率时,请注意probnorm(x)=因此，应用SAS近似计算P{8≤X≤12}的程序为：dataex;p=probnorm(1.12)-probnorm(-1.12);procprint;run;输出的结果为：0.73729.其中1.12=(12+0.5-10)/sqrt(5)7/22/202

6、162.在SAS中有probnorm(x)函数,用此函数可以求P{X≤x}.当x=1.645，1.96，2.576时,不查标准正态分布的分布函数的函数值表，应用SAS直接计算P{X≤x}的程序为dataex;dox=1.645,1.96,2.576;(给x依次赋值，增加赋值后可全部列出的函数值表)p=probnorm(x);putxp;(计算并输出x对应的概率)end;run;输出的结果如下(在Log窗口中显示)：1.6450.95001509451.960.97500210492.5760.99500

7、246777/22/20217用下列程序更好:dataex;inputx@@;p=probnorm(x);list;cards;1.6451.962.576;procprint;run;输出的结果如下(在Log窗口中显示)：1.6450.95001509451.960.97500210492.5760.99500246777/22/20218以下是用SAS程序绘制的二维正态分布分布密度函数的示意图。所用的SAS程序为：dataex;dox=-3to3by0.25;doy=-3to3by0.25;p=ex

8、p(-((x*x+y*y)*5/4+x*y*3/2)/2)/2/3.1416;output;end;end;procg3d;ploty*x=p;run;7/22/202197/22/2021103．应用SAS计算标准正态分布的分位数在SAS中有probit(p)函数，用此函数可以求p分位数.SAS程序为dataex;dop=0.025,0.05,0.1,0.9,0.95,0.975;u=probit(p);putup@;end;run;输出