【知识】SAS数据分析完整笔记(3)

2018-03-09 14:08:35 浏览数 (1)

SAS学习笔记(3):SAS一般高级语言

本篇SAS读书笔记主要介绍SAS一般高级语言,主要内容包括赋值语句、输出语句、分支机构、循环结构、数组以及函数等六个部分。

1 赋值语句

在SAS中用赋值语句计算一个值并存放到变量中。

格式为 变量名=表达式:

例如

2 输出语句

SAS数据步的输出一般是数据集,用赋值语句计算的结果会自动写入数据集。SAS也提供了一个PUT语句,可以象其它语言程序的PRINT,WRITE,printf等语句一样立即显示输出结果。PUT语句在关键字后面列出要输出的各项,每一项可以是变量名或字符串,不能为数值常量或表达式,各项之间用空格分开。PUT语句的输出结果显示在LOG窗口。例如:

代码语言:javascript复制
data ;
x=1; y=sin(x);
put 'sin function valueof' x ' is ' y;
run;

结果sin functionvalue of 1 is 0.8414709848

PUT语句的输出项还可以指定具体列位置

put x 10-20 .6; X数值显示在第10-20列,保留6位小数

put x 20.6 y 20.8; X显示在1-20列,6位小数。Y显示在21-40列,8位小数

如果希望PUT语句的输出不产生换行,使下一个PUT的结果可以显示在同一行,只要在PUT语句结尾处加一个@符,

如 put x @;

PUT语句的输出结果缺省情况下被送到运行记录窗口。在PUT语句之前用FILE语句可以改变PUT语句的输出目的地。比如,在PUT语句之前用file print;

可以把PUT语句的输出转向到输出窗口。

在FILE语句中指定一个包含文件名的字符串可以把PUT语句的输出转向到此文件中。比如file 'temp.out';把后续的PUT语句输出转向到当前工作目录下的文件"tmp.out"中,生成输出文件tmp.out。

3 分支结构

格式一: IF 条件 THEN 语句;

Else 语句;

有时我们在条件成立时需要进行的操作无法用一个语句完成,这时可以使用SAS提供的复合语句功能:只要把若干个语句用"DO;"语句和"END;"语句包围起来,就可以把它们看作是一个语句,就可以用在需要指定一个语句的地方。

if x>0 then do;

put 'x是正数';

x=2*x;

put x=;

end;

格式二:

SAS的IF结构允许嵌套,但SAS不提供IF-ELSEIF-ELSE的多分支结构。SAS的SELECT结构提供了更为灵活的多分支结构,可以实现比其它语言的IF -ELSEIF -ELSE结构更强的功能。SELECT结构有两种基本用法,第一种为:

SELECT(选择表达式);

WHEN(值列表)语句;

WHEN(值列表)语句;

OTHERWISE 语句:

END;

另一种形式:

SELECT;

WHEN(条件)语句;

WHEN(条件)语句;

OTHERWISE语句:

END;

4 循环结构

1. 计数DO循环

DO 计数变量=起始值TO结束值BY步长:

循环体语句……:

END;

在循环体中可以用LEAVE语句跳出循环,相当于C语言的break语句。

在循环体内用CONTINUE语句可以立即结束本轮循环并转入下一轮循环的判断与执行。

(2)当型循环

DO WHILE 循环继续条件:

循环体语句……:

END;

(3)直到型

DOUNTIL 循环退出条件:

循环体语句……:

END;

事实上,SAS的循环语句比上面所述还要灵活得多它在DO语句中可以指定一个循环列表,比如

循环变量i取5,7,1 1,14循环体被执行,当i取17时i的平方为289故循环体不被执行,循环结束。注意WHILE条件只作用于用逗号隔开的最后一项。

5 数组

SAS可以把一组同为数值型或同为字符型的变量合在一起,使用同一个名字称呼,用下标来区分。这与通常的程序设计语言中的数组略有区别,通常的程序设计语言中数组元素没有对应的变量名,而SAS数组每个元素都有自己的变量名。

一、数值型数组

定义:

ARRAY 数组名(维数说明)数组元素名列表(初始值表);

例如:array tests(3) math chinese english(0, 0, 0);

初始值表可以省略,这时其初始值为相应数组元素的值(如果其数组元素还没有值则初值为缺失值)

数组元素名列表可以省略,array tests(3)

也可以在说明维数时用"下标下界:下标上界"来说明一个其它的下标下界,如 ARRAY sales (95:97) yr95-yr97;这时sales(95)为yr95,sales(96)为yr96,sales(97)为yr97。

一维数组的维数说明还可以是一个星号,比如:

array tests(*) math chinese english (0, 0, 0);

二维数组定义,只要在维数说明中指定用逗号分开的两个下标界说明,例如:array table(2,2) x11 x12 x21 x22;注意,二维数组元素按行排列。

二、字符数组

定义:ARRAY 数组名(维数说明)$ 元素长度说明 数组元素名列表(初始值表);

三、临时数组

SAS也提供了与其它程序设计语言相同的数组,即数组元素只由数组名和序号决定,没有对应的变量名。这种数组叫做临时数组。

ARRAY数组名(维数说明)_TEMPORARY_ (初始值表)

6 函数

表示形式:

SUM(xl,x2,x3) SUM(OF xl x2 x3)或SUM(OF xl-x3)

1 数学函数

·ABS(x)求x的绝对值。

·MAX(xl,x2,…,xn)求所有自变量中的最大一个。

·MIN(xl,x2,…,xn)求所有自变量中的最小一个。

·MOD(x,y)求x除以y的余数。

·SQRT(x)求x的平方根。 ·ROUND(x,eps)求x按照eps指定的精度四舍五入后的结果,比如ROUND (5654.5654,0.01)结果为5654.57,ROUND(5654.5654,10)结果为5650。

·CEEL(x)求大于等于x的最小整数。当x为整数时就是x本身,否则为x右边最近的整数。

·FLOOR(x)求小于等于x的最大整数。当x为整数时就是x本身,否则为x左边最近的整数。

·INT(x)求x扔掉小数部分后的结果。

·FUZZ(x)当x与其四舍五入整数值相差小于IE-12时取四舍五入。

·LOG(x)求x的自然对数。

·LOG10(x)求x的常用对数。

·EXP(x) 指数函数

·SIN(x),COS(x),TAN(x)求x的正弦、余弦、正切函数。

·ARSIN(y),ARCOS(y),ATAN(y)

·SINH(x),COSH(x),TANH(x)双曲正弦、余弦、正切。

此外还有符号函数SIGN,一阶导数函数DIGAMMA,二阶导数函数TRIGAMMA。

2 数组函数

·DIM(x) 求数组x第一维的元素的个数(注意当下界为1时元素个数与上界相同,否则元素个数不一定与上界相同)

·DIMk(x) 求数组x第k维的元素的个数。如DIM2(x)计算二维数组x第二维长度。

·LBOUND(x)求数组x第一维的下界。

·HBOUND(x)求数组x第一维的上界。

·LBOUNDk(x)求数组x第k维的下界。

·HBOUNDk(x)求数组x第k维的上界。

3 字符函数

·TRIM(s)返回去掉字符串s的尾随空格的结果。

· UPCASE(s)把字符串s中小写字母转换为大写字母后的结果。

·LOWCASE(s)把字符串s中大写字母转换为小写字母后的结果。

·INDEX(s,sl)查找sl在s中出现的位置。找不到时返回0。

·RANK(s) 字符s的ASCII码值。

·BYTE(n)第n个ASCII码值的对应字符。

·REPEAT(s,n)字符表达式s重复n次。

·SUBSTR(s,p,n)从字符串s中的第p个字符开始抽取n个字符长的子串

·TRANWRD(s,sl,s2)从字符串s中把所有字符串sl替换成字符串s2后的结果。

4 日期和时间函数

·MDY(m,d,yr)生成yr年m月d日的SAS日期值

·YEAR(date)由SAS日期值date得到年

·MONTH(date)由SAS日期值date得到月

·DAY(date)由SAS日期值date得到日

·WEEKDAY(date)由SAS日期值date得到星期几

·QTR(date)由SAS日期值date得到季度值

·DHMS(d,h,m,s) 由SAS日期值d、小时h、分钟m、秒s生成SAS日期时间值

5 分布密度函数、分布函数

·分布函数值=CDF('分布',x<,参数表>);

·密度值= PDF('分布',x<,参数表>);

·概率值= PMF('分布',x<,参数表>);

·对数密度值= LOGPDF('分布',x<,参数表>);

·对数概率值= LOGPMF('分布',x<,参数表>);

0 人点赞