特别说明:本节【SAS Says】基础篇:开发数据,用的是数说君学习《The little SAS book》时的中文笔记,我们认为这是打基础的最好选择。
转载请在文章开头注明微信号:shushuojun,谢谢!
复习:
前面五节
- 【SAS Says】基础篇:SAS软件入门(上)
- 【SAS Says】基础篇:SAS软件入门(下)
- 【SAS Says】基础篇:读取数据(上)
- 【SAS Says】基础篇:读取数据(中)
- 【SAS Says】基础篇:读取数据(下)
在微信号“shushuojun”中回复“SAS”查看。
前面我们介绍了各种用SAS读取数据的知识,现在数说君考你两个问题:
(1)有一个TXT数据文档,如何让SAS只读取第3到第5行的数据?
(2)我们一般读取的数据都是逗号分隔的,现在有一个数据文件是制表符分隔的,那么要怎么读取,或者说要用到什么命令?
如果答不上来,回复“SAS”查看前面的内容。本节讲的是如何开发数据——如何使用SAS的函数、如何用if-then语句、如何处理SAS烦人的日期格式、如何使用retain语句让SAS保存前一次变量的值、如何快捷的列出SAS的变量名等等。
本节目录:
开发数据
3.1 创建并重新定义变量
3.2 使用SAS函数
3.3 使用IF-THEN语句
3.4 用IF-THEN语句将观测值分组
3.5 构造子集
3.6 处理SAS的日期数据
3.7 可选择的数据格式和函数
3.8 使用retain和sum语句
3.9 用数组简化程序
3.10 列出变量名的快捷方式
开发数据
3.1 创建并重新定义变量
可以通过分配语句来创建并重新定义变量,基本形式为:
Variable=expression
Variable是变量名,expression可以是常量、其他变量、或者数学表达式。分配语句的基本类型有:
Expression是数学表达式时,需要遵循运算法则,先算指数、再算乘除、之后是加减。可以用括号改变运算等级。
例子 如下是一个农产品估重数据,每位农民要求对他们的番茄、南瓜、豌豆、葡萄进行估重:
下面代码从garden.dat原始文件中读取数据,并进行修改:
这个程序包含了5句分配语句,第一个将14赋值给zone,第二个使type等于一个字符串常量……打印出的结果中,既包括旧变量,又包括新变量:
由于观测值susan的peas变量出现了缺失值,因此这个观测值的total和pertom变量也出现了缺失值。日志窗口的说明如下:
NOTE:Missing valueswere generated as a result of performing an operation onmissing values.
3.2 使用SAS函数
SAS有400多个函数,主要涵盖如下领域:
函数基本形式:function-name(argument,argument,…),所有的函数都需要括号,即使没有参数。下面的程序计算使用MDY函数,MDY函数需要三个参数:月、日、年。
Birthday=MDY(MonthBorn,DayBorn,YearBorn);
函数可以嵌套,即一个函数可以是另一个函数的参数。比如:NewValue=INT(LOG(10));
例子 有一个南瓜雕刻比赛的数据,pumpkin.dat数据包含了参赛者的名字、年龄、雕刻的南瓜类型、报名日期、五位裁判给出的分数。
下面的代码读取了数据、创建了两个新变量、转换了一个大小写:
AvgScore 使用均值函数创建的变量,计算参数的均值,这与直接相加再除以5不同的地方在于,当参数中出现缺失值时,直接相加再除的方法返回缺失值,而均值函数计算非缺失参数的均值。
DayEntered变量使用DAY函数,返回日期在一个月里的天数。
Type用大写转换函数将原来的字母转换成大写字母。
结果是:
3.3 使用IF-THEN语句
条件语句IF-THEN的基本形式为:IF 条件 THEN 执行;
比如:IF Model='Mustang' THEN Make='Ford';
条件语句中的一些基本比较符号:
还有IN比较符,比如这句中IF Model IN('Corvette','Camaro') THEN Make='Chevrolet';代表当Model为Corvette或Camaro的时候,将Chevrolet赋给Make。
一个条件只能有一个执行,如果要多个执行,则需要DO和END关键字。
可以用AND和OR来定义多个条件:IFModel='Mustang' AND Year<1975 THEN Status='classic';
例子 如下的数据包含了模型的名字、年份、制造商和颜色:
下面的代码从cars.dat的原始文件中读取数据,使用IF-THEN语句填满缺失值,并创建一个新变量Status
输出结果如下:
3.4 用IF-THEN语句将观测值分组
IF THEN/ELSE的一般形式为:
IF condition THEN action;
ELSE IF condition THEN action;
ELSE IF condition THEN action;
用else语句与直接用多个IF-THEN语句比起来,有两个优势,第一是更有效率,电脑将占用更少的时间;第二是else可以确保你的两个condition之间互斥。
有时候最后一个ELSE只有action,没有IF-THEN:
例子有一个住房改善的数据,home.dat,包括了姓名、改善工作、改善成本:
下面的代码读取数据,并新建了一个CostGroup的变量。根据Cost的值将数据分成high、medium、low和missing三类:
输出结果是:
3.5 构造子集
IF语句可以构造子集,取数据集中的部分数据。
基本形式为: IF expression; 比如: IF Sex='f';
如果IF条件中的数据是真,则数据步将继续执行。
还可以使用DELETE语句,来删除哪些不要的数据: IFexpression THEN DELETE;
这两句话是等价的: IF Sex='f'; IF Sex='m'THEN DELETE;
例子 有关于莎士比亚歌剧的清单,Shakespeare.dat,包含歌剧名、首次表演年份、类型:
下面的代码读取数据,并且用IF语句构造一个只包含喜剧(comedies)的子集:
输出结果如下:
观察日志有时能很好的保证我们截取了我们要的数据:
在这个例子中,用DELETE等价的语句为:
IF Type='tragedy' OR Type='romance' OR Type='history' THEN DELETE;
3.6 处理SAS的日期数据
日期数据的处理很棘手,有的月份有31天、有的30天、有的28天。SAS简化的日期数据,将所有的日期转化成一个以1960年1月1日为起点的数。比如:
SAS处理日期数据的三个工具为:读取数据的informats,使用数据的函数(functions),打印数据的formats
Informats 读取日期数据需要用formatted input。比如,如何告诉SAS用MMDDYY10. imforat读取名为BirthDate的变量:
INPUT BirthDate MMDDYY10.;
设定默认的百年 07/04/76这样的数据可能是1976,也可能是2076、1776。因此需要YEARCUTOFF=来指定一个一百年的第一年,默认的是1920年。下面的语句就是告诉SAS将一个两位年份的日期解释为1960年到2049年之间:
OPTIONS YEARCUTOFF=1950;
SAS表达式中的日期一旦被以SAS日期格式读取之后,可以将此数据想其他数值数据一样用在表达式中。比如像为图书馆的书设定21天的还书日期,只需要在结束日期上加上21:
OPTIONS YEARCUTOFF=1950;
通过在表达式中加入引号和字母D,可以将一个日期当做常数来使用,如下的代码创建了一个EarthDay05的变量,并且等于April22,2005
EarthDay05='22APR2005'D;
函数 SAS日期函数使得操作大大简便,比图today()返回今天的日期。
语句DaysOverDue=TODAY()-DateDue;可以计算一本书应归还还剩余的期限。
Fomats 打印日期数据时,还需要将数值换成日期,下面的FORMAT语句告诉SAS用WEEKDATE17.格式打印变量BirthDate。
FORMAT BirthDate WEEKDATE17.;
例子 图书馆有借书卡数据,Dates.dat,包含持卡人姓名、出生日期、卡办理日期。
下面的代码读取数据,计算变量使用期限(expiredate),使用期限为3年;变量expirequarter计算使用期限的四分之一,使用函数QTR()。接着用IF语句来判断一个卡是否为新卡,在2003年1月1日之后办理的,为新卡:
输出结果为:
注意BirthDate没有用日期格式。
3.7 可选择的Date Informats,Functions和Formats
下面是例子:
3.8 使用retain和sum语句
当开始数据步的每一个观测值迭代时,SAS会先将所有变量值设为缺失,再通过input和分配语句改变。Ratain和sum语句可以改变这种方式,
Retain语句 retain语句可以让SAS保存前一次变量的值。它可以出现在数据步的任何位置,基本形式为:
RETAIN variable-list;
也可以指定一个初始值,而不是用缺失值或前一次的值代替初始值
RETAIN variable-list initial-value;
Sum语句 SUM语句用于你想将一个表达式的值累加到一个变量上去时,基本形式为:
variable expression;
这个语句将表达式的值赋给变量,同时将变量的值保留到下一次迭代。这个变量必须是数值型,且初始值为0。因此,语句等价于如下形式:
RETAIN variable 0;
variable=SUM(variable,expression);
例子 有一个关于本赛季棒球比赛的数据,games.dat,包含比赛日期、参赛队伍、hit数据、runs数据
现在需要增加两个变量,一个反应本赛季的总runs数,一个反应一场比赛中最大的runs数。下面的代码用sum语句实现总run数,用retain和max函数实现最大runs数:
变量maxruns取前面迭代的maxruns和runs中最大值;变量runstodate将每一场比赛的runs都加到自己身上。结果如下:
3.9 用数组简化程序
对于太多变量要处理的程序,数组将大大简化程序。
SAS中,数组是一组变量,变量可以是已存在的,也可以是新创建的。
数组在数据步中用ARRAY来定义,基本形式为:
ARRAY name(n) $ variable-list;
Name是数组名,n是变量数,()也可以用[]和{}代替。如果变量是字符串,则需要$,且变量是新创建的字符串时,$是必须的。变量名依照顺序排列,如数组:
ARRAY store(4) Macys Penneys Sears Target;
则store(1)是Macys,store(2)是Penneys,store(3)是Sear,store(4)是Target。
数组本身不储存在数据集中,只有在数据步中才被定义。命名规则与变量一样(不超过32字节,以字母、下划线开头,只能包含字母、数字、下划线)
例子 广播电台wbrk做了一份关于歌曲的听众调查,对10首歌进行打分,分值在1-5,如果没听过则填9。数据文件wbrk.dat包括了被访者姓名、年龄、以及十首歌的打分。
下面的代码将所有打分为9的改为缺失值:
十首歌被放入song的数组中。输出结果如下:
注意这里数组没有被保存到数据集中,而i被保存了。
3.10 列出变量名的快捷方式
如果想把100个变量放入数组,并不需要一个一个变量名的输入,有快捷方式可以列出变量名。
Number range lists 开始于同一个单词,结尾于连续的数字的,可以使用Numberrange list。比如:
Name range lists这种列表是依据变量在数据集中的排列顺序来的,比如,创建如下数据步:
则变量的排列顺序就为:Y A C H R
那么可以依照这个顺序用“put 第一个变量--最后一个变量”来简化:
如果不能确定数据集中变量的顺序,可以用proc contents的postion选项来查看。下面的代码列出了永久数据集distance的变量顺序:
例子 广播电台WBRK想要修改前面的代码(将9改为缺失值),使用mean函数计算平均分数:
程序中,当原始变量(domk-ttr)值为9时,song变量值为缺失值,否则就把原始变量的值赋给song变量。另外avgscore计算平均值: