|
很好的一篇文章,特适合入门者。
/ p1 d6 M! J$ u从本质上讲,SAS是一种完善的第四代计算机语言。因此要真正掌握它,我们仍然要抛开其华丽的外表,从学习它的核心――SAS程序开始。 现在,让我们将SAS看成一个计算能力极强的统计学白痴(之所以这样说,是因为它计算能力虽然极强,却只能帮你计算而不能提出自己的实验设计方案或研究方向来),而你有一个非常小的关于数据分析的问题要请它帮忙。自然你要开口提出请求,无论措辞是委婉动听还是直截了当,你的大实话无非是“喂,老兄,我有这样一些数据,我想做这样一种统计分析,您能帮我吗?”。他迅速检查您的要求,在认为合情、合理并且合法后,就开始进行计算,并且在很短的时间内给出计算结果。 好的,计算机语言就是我们和计算机对话时所用的语言。和以上人类对话的例子相对照,“我有这样一些数据”对应SAS程序中的数据步;而“我想做这样一些统计处理”对应SAS程序中的程序步;最后的那一句哀求“您能帮我吗”则对应了Submit命令。在默认情况下,运行结果或者出错信息将自动给出(当然你也可以强制不让它输出)。此外,还有一些系统环境控制语句,如Libname、title等。SAS程序就是由一个或多个数据步和/或程序步加上一些乱七八糟的环境控制语句组成。在本章中,我们将首先学习SAS对数据的管理方式,然后会了解到什么是数据步和程序步,最后则介绍一下SAS语言中结构化语句的语法。 §3.1 SAS系统对数据的管理 在SAS系统中只有SAS数据集才能被SAS过程直接调用,SAS数据集的结构和DBF数据库完全相同,因此无须多讲。而SAS数据集存储在被称为SAS数据库的文件集中,在PC系统中,SAS数据库与某一个文件夹相对应,我们要为每一个数据库指定一个库标记(库名)来识别该库,使用Libname命令可以指定库标记。它的一般格式如下: Libname 库标记 ’文件夹位置’ 选项; 例如要指定目录“C:\USER”为库标记A,可以在视窗中提交如下语句: libname a ’c:\user’; 数据库可分为永久库和临时库两种。临时库只有1个,名为WORK,它在每次启动SAS系统后自动生成,关闭SAS时库中的数据集被自动删除;永久库可有多个,用户可以使用Libname语句指定永久库的库标记,永久库中的所有文件都将被保留。但库标记仍是临时的,每次启动SAS系统后都要重新指定。不过为了方便用户,SAS在每次启动时都会自动指定两个库标记: SASUSER:指明为永久库,即库中的数据集被保存起来,以便下次启动系统时使用。
3 i; G( ]$ h2 \, V; uWORK:指明为临时库,对于开发和检查新程序非常有用,但每次SAS运行结束后WORK库中的所有文件将被删除。在程序中引用该库中的数据集可以省略库标记,即它被认为是缺省的数据库。
. Z) y4 s. w- H0 `8 B SAS数据库对应文件夹,如SASUSER对应c:\sas\sasuser,WORK对应c:\sas\saswork, SAS数据集则对应文件,每个数据集实际上是在硬盘的相应文件夹内产生一个文件名为数据集名,而扩展名为.sd2的文件,因此如果你成功的建立了一个数据集,理论上你就可以反复使用它直到硬盘报废。 每一个数据集都有一个两级名,第一级是库标记,第二级是数据集名,中间用“.”隔开,在程序中通过指定两级名来识别文件。文件两级名的一般形式如下: 库标记.数据集名 如在a库中的数据集abc可以这样来引用:a.abc 。 §3.2 SAS程序的数据步 3.2.1 预备知识 SAS程序的语句从上一个语句结束处开始,以一个分号结束,可占多行。 3.2.2 数据步 以DATA语句开始,用于创建和处理数据集。其中最常用的语句有: DATA语句 它的主要功能是:
( s2 z( p- B# B6 E( R5 C l标志数据步的开始。 a6 ~' u0 y, Y" A+ V
命名将要创建的SAS数据集。5 Z9 N. s( @2 w* m
DATA语句的一般形式如下: DATA 数据集名; INFILE语句 用于从外部文件读入数据,必须出现在INPUT语句之前。它的主要功能是:( S8 m4 I( R# y) V5 G* P8 s, M4 ~
确定一个包含原始数据的外部文本文件。
% C; J8 \8 |/ m! S' H) u) gINFILE语句的格式如下: INFILE ’外部文件的所在位置及名称’ 选项; CARDS语句 用于直接输入数据,标志着数据块的开始。格式如下:. ?; \, N" O% J5 N2 u' r( ]
CARDS;
3 e9 ~3 Z/ E" _5 S- G& G9 C数据块2 Q/ H4 o3 D3 i. L
;
6 R2 U |* K4 b6 r* iINPUT语句 用于向系统表明如何读入每一条记录。它的主要功能有:( O6 a% ]- G6 p0 U! J6 {
读入由语句指定的数据列。, H3 n3 x) c' l, l7 s6 H
为相应的数据域定义变量名。- h8 H) u+ E4 r2 @9 p
确定变量的读入模式。
# j* ` a% \. E6 \INPUT语句的格式如下: INPUT 变量名 [变量类型 起止列数]...; 方括号表示其中的内容为可选,如果不输入,系统会以默认值代替。 3.2.3 数据的两种输入方式 上面我们学到了INFILE语句和CARDS语句,它们分别对应了两种数据输入的方式: 【直接输入方式】 实际上我们在第一章已经用到了直接输入方式,仍以那个程序为例,其中的数据步如下: data temp;+ a9 ?; i2 C6 b3 N. x2 c o7 x. I
命名将要建立的数据集为work.temp
9 t9 O* c ~! v2 \ x- K- Finput x y@@;
' f P; u$ w R8 E a 要输入的变量为x、y,并且连续输入 cards;- y7 p6 }0 [$ h) O2 H6 i2 p' R8 }5 Q
直接输入数据,数据块开始 34 56 78 90 35 67 89 10 23 65 77 456 u- k2 G/ S4 B) p) k
数据块) ?" s/ x. ~/ M2 Q" h
;
4 ^1 g7 [0 L- l4 ?" l/ d( N0 H 数据块结束 变量y后面的@@表示数据可以在一行里连续读入,SAS默认按列来分隔变量,可是这里只有两个变量,输成两列数据太长,因此加上两个@@,SAS见到这个符号,在按变量名依次读取完数据后,不是跳到下一行,而是继续在该行读数据,直至本行结束或到达分号为止。 如果你对这段叙述无法理解,请去掉两个@@,重新运行一遍程序,看看LOG视窗中的提示都有些什么变化。
) i+ {5 c& [ C3 |5 N5 r$ e【外部文件读入方式】 如果刚才的数据已经事先输好,在硬盘上的“C:\USER”文件夹内存为temp.dat文件,该文件内容如下(纯文本): 34 56 78 90 35 67 89 10 23 65 77 45 则我们可以输入程序如下: data temp; 命名将要建立的数据集为work.temp
) {$ _4 E% b! p! A0 C5 X1 ZInfile ’c:\user\temp.dat’; 指定外部数据文本文件名% C& J6 f8 O7 m9 k! M2 v# v' n
input x y@@; 要输入的变量为x、y,并且连续输入 Run; 数据步结束,以上语句可以执行了 数据步以DATA语句开始,那么在哪里结束?在遇到run语句或另一个数据步/程序步时就结束了。 为什么程序的最后都要加一句run?这个问题还不太好解释,每一个程序应至少在最后有一个run语句,表明前面的所有语句可以提交运行了,如果没有这一句,SAS会以为你后面还有语句要输入,从而将一些非立即执行的语句保留在编译缓冲区中,等待后续命令发出后一起执行,这有时会把事情弄得很糟。出于减少麻烦的需要,建议大家养成在每一个数据步或程序步后都写上一句run的习惯。 Submit命令和run语句有什么区别?这是初学者最爱问的一个问题,也同样是一个难以回答的问题。简单地讲,run是程序语句,可在程序中多次出现,它表示前面的程序段已经全部写完,可以作为一个或几个整体提交运行了,而Submit则是SAS命令,只在最后程序运行的时候发出,表示将程序正式提交运行。如果你还不明白,那么我还是来给你打个比方,run语句是你对SAS说“我想让您做的就是这些了”,而Submit命令就是对SAS说“现在开始为我做这些事,好吗”,明白了吗? 【读入其他格式的数据文件】 除了以上的两种通过数据步创建数据集的方法,SAS还提供了一些其他的方式可以用来读入其他格式的数据文件。6.11版本以上的SAS可以利用FILE菜单上的import命令将其他格式的数据文件导入SAS系统,创建SAS自己的数据集。可以导入的数据文件格式有:dBase数据库,EXCEL工作表,LOTUS的数据库,纯文本的数据文件等。 导入的操作完全是对话式的,界面友好,简便实用。以下简单叙述导入的步骤,假如例1.1的数据输成一个dBASEIII数据库文件temp.dbf, 已经存放在"c:\user"下,要导入成数据集work.li1_1。选择FILE菜单上的import,弹出一个对话框,按照向导的提示进行下去。 1. 选择导入的数据格式,从下拉式菜单上选择DBF格式,单击NEXT按钮。 2. 给出数据文件的位置和文件名,在对话框中键入C:\USER\temp.DBF,或点BROWSE直接从上面选择文件,选好后单击NEXT按钮。 3. 选择导入的目的地,即指定要创建的数据集的名字和存放的数据库名,先在左面的对话框选择数据库名WORK(临时库),在右面的对话框键入数据集的名字li1_1,此名可任意起,少于8个字符,选择完后,单击FINISH按钮,就完成了此次操作。 这时已经建好了一个数据集,名为WORK.li1_1,与我们前面建立的数据集完全一致。 §3.3 SAS程序的程序步 通俗的讲,SAS的程序步(有的书中也称过程或过程步)就是已经编好了的用于数据整理和统计的计算机程序,你只需要调用它们就是了。程序步总是用一个PROC语句开始,后面紧跟着程序步名,用以区分不同的程序步。以下是一些常用的程序步的名称及功能。 程序步名
& h' [& X/ b' F' [7 |功 能 SORT 将指定的数据集按指定变量排序. h. L8 V# m2 a2 _8 v
PRINT 将数据集中的数据列表输出6 D5 p2 N: ~5 p/ P$ @5 R9 V9 A7 W/ X
MEANS 对指定的数值变量进行简单的统计描述0 L* }. m0 X5 w; v6 K+ F4 D7 B3 e; ~
FREQ 对指定的分类变量进行简单的统计描述 W( W, F) J* A1 ~1 u
TTEST 对指定的变量做t检验' g1 n9 Q% B( F1 \0 e
ANOVA 对指定的变量做方差分析5 Z/ D) v Q9 K7 u6 `4 ^
NPAR1WAY 对指定的变量做非参数检验- |5 u* [" i* t9 l7 r
REG 对指定的变量做回归分析& _8 n# s. g" c
CORR 对指定的变量做相关分析4 R0 L; ], {5 T1 \# Y- k/ q
CHART 绘出低分辨率的统计图 PROC就是程序(procedure)的缩写,而程序步的名字大都是其功能相对应的单词或词组的缩写。 可见SAS的许多功能就是通过程序步来体现的,可如此多的种类也使我们难以总结出一个程序步的通式来。不过仍可大致给出如下结构: PROC 过程名 [DATA=数据集名] [选项];, o q+ \; n" Z+ f7 j
该过程的专用语句描述;/ Z1 t* \3 K& S) ~
[VAR 变量序列;]5 e4 E( n8 f" ^5 x9 [3 k$ g8 I) X
[WHERE 条件表达式...;]3 X$ F) E# k" R) u9 ]# w0 _- Q8 c# U8 I3 C
[BY 变量序列;]$ a! o& N4 ^9 m7 O, Y/ Z( @2 V
Run; 方括号里的语句均可以省略,在这时该过程按最通常的情况来处理,即: 处理最新建立的SAS数据集。
$ w6 ]% O; ~* t7 z4 Q; ?% z# i处理所有的变量(或对一个计算过程来说处理全部数值变量)。& O3 A6 {6 e/ _) ~! o. T; N
一次处理整个数据集而不是某个子集。5 P* K& r, g" N; k8 \5 v
我们在第一章中用到的PRINT程序步就是采用的这种默认方式。 但有时我们的要求超出了默认方式所提供的范围,这时就要动用方括号里的秘密武器了。 DATA参数 指明所需处理的数据集名,请注意在这里它是一个参数而不是语句。
. b* Z' s& C' z) u! tVAR语句 如果只想分析某一个或几个特定的变量,则可用VAR语句指定它们。例如只想显示变量x的列表,则PRINT过程如下:
1 H/ O0 L; X6 y* T$ p1 N% ^Proc print ;
9 @$ Q" j& R) c+ y ?8 O Var x;* q8 I" R) y. k" ]) n
Run; 如果想显示x和y两个变量,则将VAR语句改为如下形式即可: var x y; 当然,也可以将VAR语句删除,结果相同。 Where语句 如果你想处理的不是整个数据集而只是其中符合某种条件的子集,那么WHERE语句将会非常有用,如上例中我们只想显示大于50的x变量的值,则在PRINT过程中加入where语句如下:9 o( N) n& D+ H5 u9 L/ H
Where x>50; 如果条件变为x、y两个变量的值都要大于50,则where语句改为: where x>50 and y>50; SAS语言中常用的逻辑表达符号有and(和)、or(或)、xor(异或)、ne(不等)等。 By语句 如果你需要分组处理数据,例如要按性别分组输出统计结果,你有两种选择,一是用不同的where语句将同一个程序步反复写几遍;另一种显然更酷的方法就是采用by语句。如上例我们想按不同的y值输出x值,则加入by语句如下:
3 | q! E+ v# gBy y; 使用BY语句要求数据集已经按BY语句中指定的变量排序。如果没有排序,则程序无法正确运行。可以用SORT过程来排序,语法结构如下: PROC SORT DATA=数据集名;' _. w1 ^/ ], D4 I
BY 变量名列;5 ]9 ~( [- P) A, i* F( W4 a
RUN; §3.4 结构化语句简介 每一种结构化语言编写的程序都由顺序、分支、循环三种结构构成, SAS语言也不例外。在这里简要介绍一下分支和循环语句的语法。这些语句均可直接在数据步和程序步中使用,适当地使用它们可以大大简化我们的工作。 3.4.1 分支(条件)语句 【语法格式】 语法格式如下: IF 条件 THEN
, U! M# x# Q* e5 B 程序块;
) ?* [0 x$ B0 C# V0 t8 `1 }* o- ?/ |ELSE 程序块; 可见其语法和FOXBASE语言十分相似,只是前面多了THEN,结尾没有ENDIF,可以将两者相比较来理解。其中程序块如果只有一句,则可直接写出,否则应以DO开头,以END结束。 【应用实例】 例3.1 在产生数据集temp的同时为其增加变量class,当x>50时class=1,否则class=2。 解:程序如下: data temp;
- m5 n! h% W/ C 数据步开始,定义要建立的数据集为WORK库的TEMP
' J- L7 L7 w! {6 W input x y@@;
: p3 T5 Q, ^- A( v1 ~1 Z+ b3 A 要输入的变量为X和Y,并且采用数据连续读入方式
# V. F, s* Z* u if x>50 then class=1;) ?& y2 F# ^6 t& `
建立新变量CLASS,如果X>50,则CLASS=1
# `# Q. G7 v6 p# Y j9 J else class=2;: g w+ I- U/ ^6 n
否则,CLASS=2
$ z- @5 C3 U' b& u, Z( zcards;+ k1 R' r) G8 z& y$ k9 `' U" T
数据块开始2 L" o( ?' R) X0 Z" T3 K
34 56 78 90 35 67 89 10 23 65 77 45
8 U) z- k( h, {* C2 ] 数据块
- o) a; r, q2 Y C( X5 L;3 Z# M0 u0 e/ f. J X
数据块结束* o: ]7 V/ e: j) G2 m q" C9 |
proc print; 列表输出数据集中的数据,检查有无错误$ q3 _% }" N" E' z d9 O1 j
run; 程序结束,开始运行以上程序 以上程序在书写时采用了缩进格式,使程序的结构更清楚。当然大家可以左对齐写完所有的语句,但这样书写的程序在较长时难以阅读。 为什么可以这样写?在这里有必要解释一下数据步的执行过程,input语句按变量顺序将数据读入内存缓冲区,直至每一个变量都有值相对应。然后数据步继续向下执行,从而其它语句可以修改内存缓冲区,如修改值、增加变量等,在执行到数据步结束后,程序又回到input语句处继续执行,直到数据读完为止。现在再看看这个程序,大家可以理解它了吧! 3.4.2 循环语句 【语法格式】 语法格式如下: DO 起始条件 TO 终止条件;
8 e; Z- x: e; |% j3 G2 S 程序块;: \5 k1 X, i- L
END; 该语句主要用于建立数据集。 【应用实例】 例3.2 在产生数据集temp的同时为其增加变量class,取值依次为1、2。 解:程序如下: Data temp; 数据步开始,定义要建立的数据集为WORK库的TEMP
R; {- K' w% b* ? do class = 1 to 2; 循环开始,循环控制变量为CLASS,取值从1到2- [/ q% m, X/ D3 H
input x y@@; 要输入的变量为X和Y,并且采用数据连续读入方式
6 L! y6 Z; W- W# G1 J: ^4 ] output; 用OUTPUT语句将循环控制变量写入数据集中
6 [* T) s, X( K& ?1 X end; 循环结束
. k) O7 ?0 h, h! \cards; 数据块开始: a, h* |* b6 @ d: M
34 56 78 90 35 67 89 10 23 65 77 45
& D. Y; {6 o1 R) L6 g5 |+ Z. N 数据块- A% C% J6 P/ \1 R6 g
; 数据块结束
' h2 K1 K+ c0 {) j" E3 [2 v( drun; 程序结束,开始运行以上程序 注意在数据步中,我们不需要用语句改变循环变量的大小,系统会自动改变。 现在,SAS语言的基本知识算是介绍完了。在理论上,你现在可以坐在计算机前,独立编写程序以做出你的统计作业或者是向SAS公司发射核导弹(狗头)。! `" J Q& u/ ~ l/ Y7 f
$ @' e) y1 |4 j |