SAS程序及函数详解

SAS程序及函数详解
SAS程序及函数详解

第一章 SAS编程操作预备知识

一、SAS系统简介

SAS是一个庞大的系统,它目前的版本可以在多种操作系统中运行。当前在国内被广泛使用的最新版本是8.2版,功能很强大,我深有体会。

据说9.0版已在国外面世,已经有一些有关它的抢先报道在网络上广为流传,说它如何如何美妙,令人不禁充满期待。

SAS8.2的完整版本包含以下数十个模块。

BASE,GRAPH,ETS,FSP,AF,OR,IML,SHARE,QC,STAT,INSIGHT,ANALYST,ASSIST, CONNECT,CPE,LAB,EIS,WAREHOUSE,PC File Formats,GIS,SPECTRAVIEW,SHARE*NET, R/3,OnlineTutor: SAS Programming,MDDB Server,IT Service Vision Client, IntrNet Compute Services,Enterprise Reporter,MDDB Server common products,Enterprise Miner,AppDev Studio,Integration Technologies 等。

所谓模块,我的理解是将功能相近的程序、代码等集中起来组成相对独立的部分,就称之为模块,类似于办公软件系统office中包含的word、excel、access 等。各模块具有相对独立的功能范围,我们常用的模块有base,graph,stat,insight,assist,analyst模块等,分别执行基本数据处理、绘图、统计分析、数据探索、可视化数据处理等功能。其余模块我用得很少,知道得也很少,所以也就不多说了。

SAS系统的长处,体现于它的编程操作功能的无比强大。SAS一直以来也是注重于其编程语言的发展,对于可视化方式的菜单操作投入较少,其较早的版本仅有很少的菜单操作功能,使用起来也是非常的别扭。这很可能就是在windows 人机交互式操作系统占统治地位的今天,SAS较少被人问津的原因之一。

到了最新的几个版本,SAS也对可视化操作方法投入了一定的关注。从8.0版以后,出现了几个功能强大的可视化操作的模块,如insight模块和analyst 模块等,其菜单操作的方便程度以及人机界面的亲和性绝不亚于SPSS等著名的可视化统计分析软件。

然而要想完全发挥SAS系统强大的功能,充分利用其提供的丰富资源,掌握SAS的编程操作是必要的,也只有这样才能体现出SAS在各个方面的杰出才能。

二、SAS系统基本操作及基本概念

哪位要是连软件的安装和打开都要我啰嗦的话,我劝您还是买一套洪恩的《开天辟地》好好热热身先。

SAS 8.2的界面中间是三个并排(或层叠)的窗口,那个叫做Program Editor 的窗口(窗口标签为Editor)就是用来输入SAS语句的,编程操作的所有内容都是在该窗口内完成的,各位还是要跟它先多熟悉一下。

(一)数据集(dataset)和库

统计学的操作都是针对数据的,SAS中容纳数据的文件称为数据集,数据集又包含在不同的库(暂且理解为数据库吧)中。SAS中的库分为永久性和临时性两种。顾名思义,存在于永久库中的数据集是永久存在的(只要你不去删除它),临时库中的数据集则在你退出SAS后自动被删除。至于SAS中库的概念,最简单的理解就是一个目录,一个存放数据集的目录。

数据集的结构完全等同于我们一般所理解的数据表,由字段和记录所构成,在统计学中我们习惯将字段称为变量,在后面的内容中字段和变量我们就理解为同一种东西吧!建立数据集的方法很多,编程操作中有专门的数据读入方法来建立数据集,但需要将数据现场录入,费时费力。如果数据量大,我劝各位还是先以其它方法将数据集建好,否则程序语句的绝大部分会浪费在数据的输入上。

What are其它方法?各位是不是去参考一下别的书籍或资料。要不您是不是可以等一等,我准备若干年后出一本SAS操作大全。

(二)SAS程序概述

和其它计算机语言一样,SAS语言(称为SCL语言,SAS Component Language)也有其专有的词汇(即关键字)和语法。关键字、名字、特殊字符和运算符等按照语法规则排列组成SAS语句,而执行完整功能的若干个SAS语句就构成了SAS 程序。

SAS程序包括多个步骤和一些控制语句,一般情况下均包括数据步和过程步,一个或多个、数据步或过程步,它们之间任何形式的组合均可成为一段SAS 程序,只要能完成一个完整的功能。通常情况下SAS程序还包括一些全程语句,用以控制贯穿整个SAS程序的某些选项、变量或程序运行的环境。

SAS程序的语句一般以关键字开始,以一个分号结束,一条语句可占多行(SAS 每看到一个分号,就将其以前、上一个分号以后的所有东东当作一条语句来处理,而不管他们处在多少个不同的行中)。SAS语句对字母的大小写不敏感,你可以根据个人习惯决定字母的大写或小写。

1. 库名(库标记)的定义

为了保存宝贵的数据和方便操作起见,我习惯于指定自己的库名及其路径(目录),因为SAS系统中已有的永久库(SASUSER)无论库名还是其对应的路径都太过繁琐,使用太不方便。程序中用到的数据,都可以永久的保存于该路径下,保证以后可以重复使用。指定库名的语句为全程语句,其格式如下:

Libname 库名‘路径’;

例如我们指定的库名为“a”,路径为:“e:\data\”,SAS语句如下:

libname a ‘e:\data\’;

2. 数据步

SAS的数据步以data语句开始,用于创建和处理数据集。Data语句以关键字“data”开始,格式如下:

data 数据集名;

例如:data a.case; 将创建在库a中名为case的SAS数据集,语句执行后你可在与库a对应的目录下看到刚刚建立的数据集文件case。

Data语句所指定的数据集,一般都是以“库名.数据集名”的格式出现的,也可以单独的“数据集名”出现,此时的数据集系统默认为是临时库中的数据集,退出系统后将会被删除。

data语句有两个重要的功能,标志数据步的开始和命名将要创建的SAS数据集。

除data语句外,数据步一般情况下还包括infile语句、input语句以及datalines语句等。在不同的数据输入方式下对于它们的使用方式也不一样。

SAS程序有两种常见的数据输入方式,即从外部文件读入和直接输入两种方式。

(1)外部文件读入方式

数据若已经包含在某个外部文件(文本文件或数据文件)中,可用此方法输入数据到数据集文件中。

在以上介绍的data语句后,写入以下语句:

infile ‘外部文件的所在位置及名称’选项;

input 变量名1变量名2 …变量名n;

infile语句用于从外部文件读入数据,必须出现在input语句之前。它的功能是指定一个包含原始数据的外部文件。

input语句用于向系统表明如何读入每一条数据记录。它的主要功能有:读入由语句指定的数据列,为相应的数据域定义变量名,确定变量的读入模式。

例如:

libname a ‘e:\data\’;

data a.student;

infile ‘e:\data\student.txt’;

input name height weight;

以上程序将目录“e:\data\”下的文本文件“student.txt”中的数据输入数据集student中,该数据集存放于目录“e:\data\”下。

(2)直接输入方式

数据量较少或操作者意志力坚强的情况下采用此种输入方式,在data语句之后写入如下语句:

input变量名1变量名2 …变量名n;

datalines;(在以前的版本下为cards,新版本下两者可

通用)

…………(数据行)

…………(数据行)

…………(数据行)

;

datalines语句用于直接输入数据,标志着数据块的开始。

注意:这里的数据行中数据之间以空格分隔,当然也可以其它东东如逗号等来分隔,这里大家先以空格来分隔好了。因为不同的分隔方式下input语句要采取相应的控制选项,这些我们以后再讨论,这里我们还是省省力气吧。另外数据行输完后不能像其它语句那样直接在后面加上分号,而要另起一行输入分号,这样SAS才认为这是在输入原始数据而不是在搞别的什么。

例如:

libname a ‘e:\data\’;

data a.student;

input name $ height weight;(name后面的$符号表示变量name为字符型变量)

datalines;

Linda 171 51

Mary 168 50

Selinna 169 49

;

以上程序将直接建立数据集文件student,该数据集文件存放于目录

“e:\data\”下。

3. 过程步

SAS程序的过程步表示一个处理过程,如排序、T检验、方差分析等等。过程步以关键字proc开始,后面紧跟着过程名,用以区分不同的程序步,并以关键字run结束。

一般的格式如下:

proc 过程名选项列表;

……(其它语句);

……(其它语句);

run;

SAS程序中涉及的过程多达数百种,实现统计功能时常用的过程也有数十种之多,现将最为常用的过程名称及其所能实现的功能列入下表(表1.1),以便各位提前热热身。

表1.1 常用的过程名称及其功能

4. 几个常用的重要过程

在进入一般统计学功能实现的内容之前,有关数据预处理和执行重要公共功能的过程大家有必要预先掌握,这里选出几个常用的和重要的过程进行讨论。

(1)对SAS文件进行操作的datasets过程

datasets过程是对数据文件进行管理操作的工具,利用它我们可以实现以下功能:

将SAS文件从一个库中拷入另一个库中;

对SAS文件进行重命名;

修复损坏的SAS文件;

删除SAS文件;

列出某一SAS库中所有的SAS文件;

列出一个SAS数据集的属性,如最后修改时间、数据是否压缩、数据是否索引等;

对SAS文件进行设置密码的操作;

向SAS数据集添加记录;

对SAS数据集的属性以及数据集内变量的属性进行修改;

创建或删除SAS数据集的索引;

创建并管理SAS数据集的核查文件;

创建或删除SAS数据集的完整性规则。

datasets过程的一般格式如下:

proc datasets <选项列表>;

age 当前文件名相关文件名列表

append base=数据集名

audit 文件名<(操作密码)>; initiate;<其它代码;>

change 旧文件名1=新文件名1 <…旧文件名n=新文件名n> <选项列表>;

contents <其它选项>;

copy out=库标记 <其它选项>;

exclude 文件名 <其它选项>;(该语句只能在copy语句后出现,不能和select语句同

时出现)

select 文件名 <其它选项>;(该语句只能在copy语句后出现,不能和exclude语句同

时出现)

delete 文件名 <其它选项>;

exchange文件名1=交换文件名1 <…文件名n=交换文件名n> <选项列表>;

modify 文件名 <选项列表>;

repair文件名 <选项列表>;

save文件名 <选项列表>;

run;

是不是太繁琐了,我都有些不耐烦了。不过我还是得写,已经写了嘛!不忍心浪费掉,最起码可以用来充充数,扩大一下篇幅。至于各位,嫌烦的可以略过此部分,想用的时候再来温习也不迟。下面我就不厌其烦地向各位介绍一下各选项及各条语句的含义和用法。

proc datasets语句后各选项的含义及用法见下表(表1.2)。

表1.2 proc datasets语句后各选项的含义及用法

age语句用于批量地重命名文件,按照当前文件和相关文件的排列顺序,依次将后一个文件名重命名给前一个文件,结果是最后一个文件被删除,当前文件名被废弃。

对此语句我所了解的就这么多,而且还是通过试验得出的结论,但总感觉还是理解的不对,哪位若有高见还请不吝赐教。

append语句执行向数据集添加记录的功能,选项“base=数据集名”用以指定要添加记录的数据集,“data=数据集名”则指定所要添加的记录所在的数据集,此选项若省略则默认为当前数据集(最近一次操作的数据集)。

audit语句用于对文件的核查,生成核查文件并对其进行管理;change语句以新文件名替换旧文件名;contents语句用于显示指定数据集或当前数据集的各种属性;copy语句用于将当前库中相应的文件拷贝到指定的库中,选项“out=库标记”用来指定文件要拷贝到的目标库;delete语句用于删除指定的文件;exchange语句的功能是将等号前后两个文件的文件名进行互换;modify语句用于修改文件各方面的属性;repair语句用于对指定的文件(受到过某种损坏)进行修复,使其恢复到可以使用的状态;save语句的功能是将其指定的文件保留,当前库中的其他所有文件则被删除。

(2)对数据文件中记录进行排序的sort过程

sort过程的功能是对指定数据集中的记录按照指定的变量进行排序。由于诸多过程有对数据集记录进行排序的需要,比如过程步中存在by语句(用以将数据文件分割为若干部分)的情况下,就需要按照by语句后的变量对数据先行排序,所以sort过程非常有用,而且非常常用。

sort过程的一般格式如下:

proc sort 选项列表;

by 变量名1 <变量名2 … 变量名n>;

run;

proc sort语句后各选项含义及其用法见下表(表1.3)。

表1.3 proc sort语句后各选项含义及其用法

by语句即用以指定排序所要依据的变量,变量可为数值型也可为字符型,其后可指定多个变量,sort过程在按照靠前的变量进行排序的情况下再按照靠后的变量进行排序。by语句中每个变量前可用descending/ascending选项来指定按照其排序的方式(降序或升序),默认状态为升序。

(3)将数据文件输出显示的print过程

print过程的功能是将SAS数据集的记录以一定的方式显示到输出设备(显示屏),可以显示其全部的变量或部分变量。利用此过程,你可以创建从简单列表到可进行数据汇总的各种报告的各种不同的表单。

print过程的一般格式如下:

proc print 选项列表;

by 变量名1 <变量名2 … 变量名n 其

他选项>;

pageby 变量名;

sumby 变量名;

id 变量名;

sum变量名;

var变量名;

run;

proc print语句后可跟的选项含义及其用法见下表(表1.4)

表1.4 proc print语句后选项含义及其用法

by语句在所有过程中的用法都相同,即将数据集分割为若干小数据集分别进行处理。pageby语句用来控制换页时变量的显示方式,对于其后所指定的变量,相同的值不会显示在不同的页中,该变量某一值的记录在一页的剩余部分显示不下时,则从该值的第一条记录开始换行显示。sumby语句的作用和pageby语句相似,只不过是将换页的动作换为求和,对指定变量的每一值计算var变量的总计值。id语句的作用是用指定的变量值代替记录编号对每一条记录进行标识。sum语句用于指定报告中要进行求和操作的变量,var语句用于指定要在报告中显示的变量。

以上过程作用较为普遍,使用频率较高,有必要预先了解,以便于后面所讨论内容的顺利进行。

为节省篇幅,这里不进行实例演示。内容过于枯燥,可能的错误也难免,还请各位多多包涵。

SAS程序操作的大概情况就草草的这样介绍一下吧,说得太多的话我怕各位没有耐心看下去。更多的内容我想还是留在实际的例子中来介绍,这样大家可以好好的切身体会一下,然后就会印象深刻,实际运用起来也就得心应手了。

第二章计量资料的统计描述

一、计量资料的统计描述

计量资料统计描述的内容主要包括平均指标和变异指标的计算、资料分布形态(或特征)的图形表现等。

(一)实现描述性统计功能的几个过程

描述性统计指标的计算可以用四个不同的过程来实现,它们分别是means过程、summary过程、univariate过程以及tabulate过程。它们在功能范围和具体的操作方法上存在一定的差别,下面我们大概了解一下它们的异同点。

相同点:他们均可计算出均数、标准差、方差、标准误、总和、加权值的总和、最大值、最小值、全距、校正的和未校正的离差平方和、变异系数、样本分布位置的t检验统计量、遗漏数据和有效数据个数等,均可应用by语句将样本分割为若干个更小的样本,以便分别进行分析。

不同点:(1)means过程、summary过程、univariate过程可以计算样本的偏度(skewness)和峰度(kurtosis),而tabulate过程不计算这些统计量;(2)univariate过程可以计算出样本的众数(mode),其它三个过程不计算众数;(3)summary过程执行后不会自动给出分析的结果,须引用output语句和print过程来显示分析结果,而其它三个过程则会自动显示分析的结果;(4)univariate过程具有统计制图的功能,其它三个过程则没有;(5)tabulate

过程不产生输出资料文件(存储各种输出数据的文件),其它三个均产生输出资料文件。

以上是它们的主要异同点,其它更为具体的异同点需要在实际应用中去体会。掌握了各种过程的异同点,就可以根据具体需要选择最佳的过程进行工作。

统计制图的过程均可以实现对样本分布特征的图形表示,一般情况下可以使用的有chart过程、plot过程、gchart过程和gplot过程。大家有没有发现前两个和后两个只有一个字母‘g’(代表graph)的差别,其实它们之间(只差

一个字母g的过程之间)的统计描述功能是相同的,区别仅在于绘制出的图形的复杂和美观程度。

chart过程和plot过程绘制的图形类似于我们用文本字符堆积起来的图形,只能概括地反映出资料分布的大体形状,实际上这两个过程绘制的图形并不能称之为图形,因为他根本就没有涉及一般意义上图形的任何一种元素(如颜色、分辨率等)。而gchart过程和gplot过程给出的是真正意义上的图形,可以用很多的语句和选项来控制图形的各方面的性质和特征。

chart和gchart与plot和gplot的区别则体现在不同的作图功能,前两个过程可以绘制出的图形主要有条形图(包括横条和竖条)、圆图、环形图和星形图等,后两个过程通常用一个记录中的两个变量值表示点的坐标来绘制图形,如散点图和线图等。

(二)描述性统计过程的一般格式

1. means过程的一般格式

proc means 选项列表;

by 变量名称(分组变量);

class 变量名称(分组变量);

freq变量名称(数值变量,用以表示相应记录出现的频

数)

weight变量名称(数值变量,用以表示相应记录的权重

系数)

var 变量名称(待分析的数值变量);

run;

Proc means 语句后的选项主要用来指定所要计算的统计量,默认情况下,Means过程会给出频数、均数、标准差、最大值和最小值等,其余统计量的计算均需要在选项中指定。class语句所指定的分组变量用来进行分组,而by语句所指定的分组变量是用来将数据分为若干个更小的样本,以便SAS分别在各小样本内进行各自独立的处理。freq语句和weight语句分别引导代表记录出现频数和权重系数的数值变量。var语句引导所要进行分析的所有变量的列表,SAS将对var语句所引导的所有变量分别进行描述性统计分析。

2. summary过程的一般格式

proc summary 选项列表;

by 变量名称(分组变量);

class 变量名称(分组变量);

freq变量名称(数值变量,用以表示相应记录出现的频

数)

weight变量名称(数值变量,用以表示相应记录的权重

系数)

output <统计量关键字=自定义变量名>

var 变量名称(待分析的数值变量);

run;

summary过程的格式和means过程可以说是完全相同的,各条语句和选项的含义也是相同的,包括在means过程中未列出的output语句也可以应用于means 过程,只是此语句在summary过程应用较多(这样才能将分析结果显示出来),所以才将其列入一般格式中。output语句用来对分析结果输出为数据文件进行控制,其后的选项可有可无,若无则SAS按照默认方式进行。“out=数据集名”用来定义输出数据文件的文件名称,文件名的格式和数据步中数据文件名相同。“统计量关键字=自定义变量名”用来自定义输出数据文件中各种统计量的变量名称,前者是系统定义的(和proc语句后选项中的统计量关键字完全相同),必须正确无误,后者可自行定义。默认状态下输出统计量只有频数、均数、标准差、最大值和最小值,在默认状态不能满足需要时这一选项则是必需的。

3. univariate过程的一般格式

proc univariate 选项列表;

by 变量名称(分组变量);

class 变量名称(分组变量);

freq变量名称(数值变量,用以表示相应记录出现的频

数)

weight变量名称(数值变量,用以表示相应记录的权重

系数)

histogram 变量名称/选项列表

output <统计量关键字=自定义变量名>

var 变量名称(待分析的数值变量);

run;

univariate过程和以上两个过程的格式非常相似,相同的语句和选项其含义也相同,所不同的是某些统计量只能在univariate过程中计算(如众数),以及univariate过程中所具有的绘图功能。histogram语句即用来指示SAS对其后所指定的变量绘制直方图,其后的选项用来指示SAS添加不同类型的拟合图形(如正态分布的分布密度曲线)。

4. tabulate过程的一般格式

proc tabulate 选项列表;

by 变量名称(分组变量);

class 变量名称(分组变量);

freq变量名称(数值变量,用以表示相应记录出现的频数)

weight变量名称(数值变量,用以表示相应记录的权重系数) table <<页变量表达式>,<行变量表达式>,<列变量表达式>>

选项>

var 变量名称(待分析的数值变量,统计量列入相应的表单元格);

run;

tabulate过程和上述几个过程的格式也基本相似,相同的语句和选项也代表相同的含义。最大的不同也是tabulate过程中最为重要的是table语句,他用来定义表格的具体格式以及表格中所要包括的统计量。

5. gchart过程的一般格式

proc gchart 选项列表;

图形关键词变量名称/选项列表

run;

此过程格式简单,复杂的地方在于图形关键字(每个图形关键字对应一种图形类型)所引导的语句,这里是控制图形类型及图形要素的地方,涉及到众多的关键字和选项。gchart过程可以使用的图形关键字及其所绘制的图形类型见下表(表2.1)。

表2.1 gchart过程可以使用的图形关键字及其所绘制的图形类型

图形关键字后的变量名用以指定进行图形描述时的分组变量,可以是数值型的(此时以各组的组中值为分组的标志),也可以是字符型的。其后的选项比较重要的有:(1)type=统计量关键字,表示以图形对变量(sumvar所指定的变量)的哪一种统计量进行描述,比如频数(freq)、均数(mean)、总计(sum)、频数百分比(pctn)等;(2)subgroup=变量名(分组变量),指定要进行分组(各组段内再分组)的变量;(3)sumvar=变量名(数值变量),指定要进行统计计算的变量,也就是“type=统计量关键字”选项中统计量的计算所依据的变

量。其它的选项较少用到或系统默认值即可基本满足要求,这里还是少罗嗦,以后用到再说。

6. gplot过程的一般格式

proc gplot 选项列表;

bubble 散点图表达式

bubble2 散点图表达式

plot散点图表达式

plot2散点图表达式

run;

从gplot过程的一般格式中我们就可看出,此过程只能绘制两种类型的图形,bubble语句指示SAS绘制泡状散点图,plot语句指示SAS绘制点状散点图。bubble2语句和plot2语句指示SAS在同一区域内(bubble2和bubble在同一区域,plot2和plot在同一区域)绘制第二个图形,两者的横坐标相同(同一变量),纵坐标分别位于左右两侧(可以是同一变量,也可以是两个不同的变量)。

散点图表达式的一般形式为:

(1)bubble和bubble2语句:纵坐标变量名*横坐标变量名=泡尺寸变量名(变量值以泡的大小表示),三者均应为数值变量;

(2)plot和plot2语句:纵坐标变量名*横坐标变量名<=n/分类变量名>,此处等号及其后的部分可以省略,此时SAS以默认的散点类型绘制散点图;若等号后为n(n为正整数,是散点类型的编号),SAS则以指定的编号对应的散点类型绘制散点图;若等号后为分类变量名(可为字符型或数值型,为数值型时作为离散型变量处理,每一个值将被当作一个类别),此变量的具体值(或与每个具体值对应的图形)将被作为散点用来绘制散点图。

chart过程和plot过程的一般格式及各选项使用方法分别与gchart过程和gplot过程是基本相同的,不同之处仅在于后两者中涉及到有关三维和图形元素(颜色等)的语句和选项在前两者中是无效的。例如vbar3d语句在chart过程中无效,bubble语句在plot过程中无效。其余的语句和选项使用方法完全相同,所以在掌握了gchart过程和gplot过程后,chart过程和plot过程你会不学自通。

(三)描述性统计关键字及其含义

SAS中可计算的描述性统计量多达二十余种,大部分可在以上介绍的前四个过程中计算,个别统计量在某些过程中不能计算,大家需要注意,要不然系统显示错误信息时还不知道是怎么回事。

我经常遇到这种情况,系统提示错误(此类提示信息显示在log窗口中)时总是摸不着头脑,费半天劲才能搞明白。没办法,摸着石头过河嘛!不过这样也并非一无是处,最起码可以积累很多使用经验。

下表(表2.2)列出SAS中可以计算的所有描述性统计量关键字及其含义,供大家使用时参考。

表2.2 SAS中可以计算的描述性统计量关键字及其含义

(四)应用实例演示

1. 大样本资料的统计描述

对于大样本资料的描述性统计,医学统计学教材上一般要求采用编制频数表的方法进行各种指标的计算。在用软件进行工作时大可不必这么辛苦,尤其是在玩SAS时(据说这家伙在处理成百上千万的数据时都绝不在话下)。但要用直方图表现资料分布特征时,频数表的编制却是非常必要的。下面我们就用一个实例来演示大样本资料统计描述的SAS编程实现方法。

例2-1 下表(表2.3)为长沙市某大学160名正常成年女子1999年体检资料中的血清甘油三酯的测量结果(mmol/L),请对其进行统计描述。

表2.3 160名正常成年女子的血清甘油三酯(mmol/L)

按题目要求我们来计算该指标的均数、中位数、标准差和若干个重要的百分位数以及四分位数,并以直方图来表现此资料的分布特征。

在用SAS进行处理前,我们最好先将数据转换成SAS数据集的格式,编程时直接引用就可以了。在这里我就将我一贯的做法展现给大家,欢迎各位参考,更希望大家积极尝试,找出更为便捷的方法。

谁要有了更好的方法,给我也指点指点,大家共同提高嘛!吝啬鬼是成不了大事的。

先将原始数据以Excel电子表格的形式录入,用SAS的import功能转换成SAS数据集文件,或以Analyst模块直接读入Excel文件,自动生成SAS数据集文件。

噢,差点忘了,我还没给大家介绍这两部分的内容。而且这些内容也是说来话长,我们留待以后再说吧,大家知道有此方法就行了。感兴趣的可以自己去尝试,或以后掌握了这些内容后再来操练也行。

更容易(不一定简便)的方法是将数据以纯文本(数据间以空格或逗号隔开,一个变量占一列或多列均可)的形式录入,在SAS程序中以“infile”语句来引用该文件。

当然你要不怕麻烦,也可以将原始数据全部列在Data步的datalines(或cards)语句中。在这里我就不这么做了,我嫌麻烦,而且易出错。

我们先来计算所要的描述性统计量。

计量资料的描述性统计功能通常用Summary过程、Means过程和univariate 过程来实现。tabulate过程也可以实现描述性统计功能,但它和以上三个过程来比就显得复杂了一些,其实tabulate过程的过人之处在于对含有较多分类变量的资料的汇总性分析,以后大家用到的时候会慢慢体会到这一点。

这里我们分别以三个过程来完成对该例题的描述性统计分析,使各位能深刻体会这些过程的用法。好了,我们下面就开始操作。

将上表中的数据以文本文件的形式录入,文件名为case2_1,存入目录“e:\data\”下。然后进入SAS操作环境,在Program Editor窗口中编写相应程序。

(1)以means过程进行分析

具体操作程序如下

libname a 'e:\data\';

data a.case2_1;

infile 'e:\data\data2_1.txt';

input x@@;

proc means n mean median p1 p5 p95 p99 q1 q3 max min;

var x;

run;

以上程序建立了一个永久性数据集,以后可重新利用。

默认情况下,Means过程会给出频数、均数、标准差、最大值和最小值等,若默认情况满足需要,则Proc Means后的所有选项(代表所要求计算的统计量)皆可省略,否则需一一写出,如以上程序。

程序编写完成后,点击菜单Run>>Submit或直接点击快捷工具栏内的那个向右跑的小人,将程序提交给SAS执行。以上程序提交执行后的结果如下。

程序执行结果给出了所有指定的统计量,列表中各项目的含义不用我再啰嗦了吧!我想大家都能看得懂。

很是简单吧!的确,我们在做一般的统计分析时编写的程序比较简单,SAS 中相当部分的语句和选项是很少用到的,但要想充分利用SAS的强大功能,尽可能的掌握更多的内容是非常必要的。在这里我不想和那些没完没了的语句和选项过多地纠缠,等以后的例子中涉及到相应的内容时我们再来仔细研究。

(2)以summary过程进行分析

操作程序如下

libname a 'e:\data\';

data a.case2_1;

infile 'e:\data\data2_1.txt';

input x@@;

proc summary;

var x;

output out=temp1 mean=m median=md p1=pt1 p5= pt5 p95=pt95 p99=pt99 q1=qt1 q3=at3;

run;

proc print data=temp1;

run;

大家要注意在summary过程中对output语句的使用,他将需要显示的统计量以自定义变量的形式存储在自命名的输出数据文件中(但还是不能自动在结果中显示)。最后使用print过程将输出数据文件显示出来。

将程序提交执行,结果如下。

相关主题
相关文档
最新文档