最新医用SAS统计分析(一)精品课件.ppt

《最新医用SAS统计分析(一)精品课件.ppt》由会员分享，可在线阅读，更多相关《最新医用SAS统计分析(一)精品课件.ppt（58页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、医用医用SAS统计分析统计分析(一一)SAS软件简介vSAS发展概述 SAS(Statistical Analysis System)是一个管理数据、分析数据和编写打印各种形式报告的组合软件系统，是国际上非常流行的统计分析软件之一。1985年，美国SAS研究所推出了可以在微机上运行的SAS/PC版本，此后，又不断出新的版本，功能不断增强，越来越多的人在用它进行统计分析的同时，也把它用作数据管理软件。在英美等国，能熟练使用SAS进行统计分析是许多公司和科研机构选材的条件之一。美国FDA新药审批程序中，新药试验结果的统计分析规定只能用SAS进行 ，目前已到SAS9.0版本。qSAS运行菜单v“Fi

2、le” 调用、储存、打印文件v “Edit” 用于编辑文件v“Local” 用于运行本机文件 其中：Submit是运行程序，Recall是调回前面用过的程序。v“Globals” 可供选择窗口v “Options” 选择SAS的功能v“Help” 可随时提供帮助咨询q菜单下方的工具条 Submit（运行），New（新文件），Open（打开文件），Save（储存文件），Print（打印），View（预览），Cut（裁剪），Copy（复制），Paste（粘贴），Undo（恢复），DOS Prompt（DOS方式），Browse（连接英特网）， Libraries（库）, SAS/ Assist（应

3、用SAS/ASSIST模块）， Help（提供帮助）q SAS文件系统q*.sd2 SAS数据集(6.12版), 8.0以上版为 *.sas7bdatq*.sas EDITOR视窗输出SAS程序文件 q*.log LOG视窗输出文件 q*.lst OUTPUT视窗口输出文件qSAS中的常用变量 v数值型变量 变量名有18个字符组成，以英文字母开头，其余可以是英文字母、数字或下划线，不能有空格。 如AGE， X2， X1_1，等都是合法的， 1X ， XY-1，X 1，等都是不合法的。 变量值过大或过小的数可用科学记数法，如：1.785E-19 即为1.78510-19 。 v字符型变量 字符型

4、变量在变量名后加“$”号表示,如NAME$ ， SEX$等。一个字符变量是由1-200个字符组成的字符串，如：Zhang Hua , M等都是一个字符变量。qSAS运算符与常用函数v算术运算符 （乘方）（乘） /（ 除） +（ 加） （减）v比较运算符 = 或 EQ （等于） = 或 NE（不等于） 或 GT（大于） = 或 GE（大于等于） 或 LT（小于） 25 & B30 ； A25 | B30 ；NOT(A3)v 常用函数 ABS(x) x的绝对值； SQRT(x) x的平方根； EXP(x) e(2.71828)的x次方； LOG(x) x的自然对数（以e为底）； LOG2(x) x

5、的以2为底的对数； LOG10(x) x的普通对数（以10为底）qSAS程序结构：程序结构：由三部分组成；SAS语句以一个关键词开始, 以分号(;)结束。v数据步 以data语句开始，用于创建和处理数据集。Data语句的功能：语句的功能：(1)标志数据步的开始；标志数据步的开始；(2)命名将要创建的命名将要创建的SAS数据集。数据集。 data语句格式：语句格式： data 数据集名；数据集名；Input语句的功能：语句的功能：(1)读入由语句指定的数据列；读入由语句指定的数据列；(2)为相应数据域定义变量名；为相应数据域定义变量名；(3)确定变量的读入模式。确定变量的读入模式。Input语句

6、格式：语句格式：input 变量名变量名 变量类型变量类型；Cards语句：语句：标志数据块的开始。格式为：标志数据块的开始。格式为：cards;data esr ; /*建立临时数据集esr.sd2,自动放在saswork子目录下*/Input x ; /* 指明要输入的变量 */Cards ; /*标志数据区开始*/ 3 9 8 6 5 5 7 3 10 8 10 4 ;v过程步 以PROC语句开始，对数据进行分析和处理时，需要调用SAS过程，SAS过程是一系列应用程序。语句格式：PROC SAS过程名选项；常用的选项为：data=数据集名 如：proc print ; /*调用现成的SA

7、S过程，进行统计分析*/v提交程序步 Run ; /*在程序最后,指示过程步或数据步结束,可以提交程序*/ Log窗输出 Output窗输出窗输出q程序运行程序运行-在编辑窗输入程序后在编辑窗输入程序后点击图标按功能键“F8”在命令行（窗）键入“submit”在local下拉菜单中选“submit ”q窗口内容保存窗口内容保存-激活窗口后激活窗口后在命令行（窗）键入 file 路径文件名点击图标在file下拉菜单中选“save ”或“save as”只有编辑窗口(edit 窗)才能打开文件q常用的快捷键常用的快捷键F1 Help命令，进入帮助视窗。F4 recall命令,用于edit窗,将前面

8、提交的程序 调回，可多次使用。F5 切换到edit视窗。F6 切换到log视窗。F7 切换到output视窗。F8 submit命令。F9 Keys命令，进入功能键视窗。Ctrl-E 清除当前窗口的内容。 edit窗被清除的内容不能用F4调回。SAS数据集建立数据集建立-用用input和和cards语句直接输入语句直接输入 data child; /*建立临时数据集child.sd2,自动放在saswork子目录下*/ input id x1 $ x2 x3 x4 x5 x6; /* 指明要输入的变量 , $为字符型变量*/ cards; /*标志数据区开始,数据之间以一个或几个空格分隔*/

9、1 m 32 95.5 14.0 53.5 49.64 2 m 35 92.0 13.0 52.0 41.61 3 m 33 89.0 12.5 53.5 35.81 254 m 176 168.0 53.5 82.0 100.14 255 f 30 91.0 11.0 48.0 35.39 256 f 33 91.0 11.5 47.0 44.98 521 f 178 163.0 51.0 79.0 87.42 ; /* 标志数据区结束，分号必须单独一行 */ run;SAS数据集建立数据集建立-数据读入数据读入data a1;input id x1$ x2 x3 x4;cards;1 m

10、32 95.52 m 35 92.0 13.03 m 33 89.0 12.5;proc print;run;结果:OBS ID X1 X2 X3 X41 1 m 32 95.5 2.02 3 m 33 89.0 12.5data a1;input id x1$ x2 x3;cards;1 m 32 95.52 m 35 92.0 13.03 m 33 89.0 12.5;proc print;run;结果:OBS ID X1 X2 X3 1 1 m 32 95.5 2 2 m 35 92.0 3 3 m 33 89.0q不分行符的使用 data esr; input x ； /*每读入一条记

11、录后,数据指针保持原位不换行, 继续读下一条记录*/ cards; 3 9 8 6 5 5 7 3 10 8 10 4 ; proc means; run; 结果： Analysis Variable : X N Mean Std Dev Minimum Maximum - 12 6.5000000 2.5405797 3.0000000 10.0000000 -SAS数据集建立数据集建立q缺失值的输入- 以 . 表示，缺失值不进入分析 data child; input id x1 $ x2 x3 x4 x5 x6; cards; 1 m 32 95.5 14.0 53.5 49.64 2

12、m 35 . 13.0 52.0 41.61 3 m 33 89.0 12.5 53.5 35.81 254 m 176 168.0 53.5 82.0 100.14 255 f 30 91.0 11.0 48.0 35.39 256 f 33 91.0 11.5 47.0 44.98 521 f 178 163.0 51.0 79.0 87.42 ;proc means;run;SAS数据集建立数据集建立结果Variable N Mean Std Dev Minimum Maximum - ID 7 184.5714286 194.9639771 1.0000000 521.0000000

13、X2 7 73.8571429 70.4779601 30.0000000 178.0000000 X3 6 116.2500000 38.2410120 89.0000000 168.0000000 X4 7 23.7857143 19.4825930 11.0000000 53.5000000 X5 7 59.2857143 14.7361719 47.0000000 82.0000000 X6 7 56.4271429 26.2540446 35.3900000 100.1400000 -SAS数据集建立数据集建立-永久数据集和临时数据集永久数据集和临时数据集qSAS数据集名称的表示方式

14、数据集名称的表示方式:两级文件名两级文件名v定义库关联名定义库关联名： libname 命令 libname 库关联名库关联名 盘符盘符:路径路径; 如： libname dw c:teach； 路径必须实际存在路径必须实际存在。v库关联名库关联名. 数据集名数据集名： 如 dw.esrq临时数据集临时数据集：存放在sassaswork子目录中的数据集，一旦退出sas，会被自动删除。它的库关联名库关联名默认为work，均省略不写。前面所生成的数据集都是临时数据集，退出sas后都会被删除。q永久数据集永久数据集：通过两级文件名的形式，将生成的数据集存放到saswork以外的路径中（库关联名所指向

15、的路径库关联名所指向的路径），即可以生成永久数据集。即使退出SAS也还存在，可供以后调用。例： libname dw c:teach； data dw.esr; input x1 x2 x3 ; (略) q从已建立的SAS数据集中读入数据建立新的SAS数据集-set命令libname L c:teach;data L.child ;set chlid;proc print;run; data a;set L.child;proc print;run; data b;set a;run;注：Proc 步的操作对象总是最新生成的SAS数据集。只有proc步时，也可通过在proc语句后加data=指

16、定要分析的数据集。如:proc means data=L.chlid;run;SAS数据集建立数据集建立-从其它数据集转入从其它数据集转入q利用file下拉菜单中的import功能 ，直接将其他类型的数据文件，转为SAS数据集。可转入：v文本型数据文件（*.txt) 可用Notebook(记事本）输入数据，首行为变量名，制表键(Tab)分隔。vdbase II 和 III数据库文件（*.dbf）vexcel数据库 SAS8.0以上可转入Excel 97及2000的数据库，用Excel输入数据，首行为变量名。 有的数据可通过粘贴的方法贴到edit窗后，加入input和cards等语句，生成sas

17、数据集。q用infile语句建立SAS数据集 语句格式 infile 盘符：文件夹文件名.txt (或dat) ; 首先用任何文字处理软件来建立一个数据文件，但一定用纯文本格式保存。该纯本数据文件是不包含变量名的，必须知道那一列数据是哪一个变量值，调用使用input x1 x2 指明。例 新建纯文本文件“xy”，路径为 e:sasxy.txt 。数据为： 56 34 57 35 58 36 59 37 60 38 61 39 62 40 63 41 64 43 61 28data Li_1;infile e:sasxy.txt;input x y ;proc print;run;OBS X Y

18、 1 56 34 2 57 35 3 58 36 4 59 37 5 60 38 6 61 39 7 62 40 8 63 41 9 64 43 10 61 28SAS数据集的整理数据集的整理-变量操作变量操作q变量重命名：变量重命名： rename 旧变量名=新变量名； 如 rename id=x0; q变量赋值或修改变量值：变量赋值或修改变量值： 变量名=表达式；如 y=log10(x);q累加语句 累加变量+表达式 累加语句能把语句中表达式的运算结果或变量值累加到累加变量上，并作为记录送到数据集中。 累加变量必须是数值型，表达式是符合SAS语言规则的表达式。累加变量的初始值由系统自动设置

19、为零。例data abc;input a b c ;sum=a+b+c;sumz+sum;cards;21 22 26 5 7 4 12 3 67 16 23 19 2 3;proc print;run; OBS A B C SUM SUMZ 1 21 22 26 69 69 2 5 7 4 16 85 3 12 3 6 21 106 4 7 16 23 46 152 5 19 2 3 24 176q用keep或drop语句整理变量名表 keep ；（指定要保留在数据集的变量） drop ；（指定要从数据集中删除的变量）q结合if语句进行条件赋值 if then ； else ；例 data

20、keep; input x1 x2 x3 ; x4=log10(x1); keep x1 x2 x4; rename x4=lgx1; Cards; 20 34 32 24 50 93 34 43 56 42 56 88 ; proc print; run;注：keep语句等价于 drop x3 ;OBS X1 X2 LGX1 1 20 34 1.30103 2 24 50 1.38021 3 34 43 1.53148 4 42 56 1.62325data drop;set child;if id200 then group=1; else group=2;drop x3 x4 x5 x6

21、;rename id=x0; proc print;run;注：drop语句等价于 keep id x1 x2 group;OBS X0 X1 X2 GROUP 1 1 m 32 1 2 2 m 35 1 3 3 m 33 1 4 254 m 176 2 5 255 f 30 2 6 256 f 33 2 7 521 f 178 2例SAS数据集的整理数据集的整理-记录操作记录操作q数据子集化数据子集化 v记录输出记录输出语句格式：if ； “将符合表达式的记录输出到新建的数据集。”语句格式：if then output ； else output；Output 的使用：一般SAS在data步

22、语句全部执行完后才将数据送入新数据集。但遇到output语句时，会马上将当前的数据送入新数据集；此时，如果output语句后还有其他data步语句，则不再被执行。 例1 data a;input x ; cards;2 3 5 6 8 9 10;data b;set a;if x3; /*只将x3的变量值由数据集a 赋值给数据集b*/proc print;run;v记录删除记录删除 if then delete; “ 将符合表达式的变量从数据集中删除。” OBS X 1 5 2 6 3 8 4 9 5 10data b1; set child; if x1=m; proc print; Run

23、; OBS ID X1 X2 X3 X4 X5 X6 1 1 m 32 95.5 14.0 53.5 49.64 2 2 m 35 92.0 13.0 52.0 41.61 3 3 m 33 89.0 12.5 53.5 35.81 4 254 m 176 168.0 53.5 82.0 100.14data b1 b2;set child; if x1=m then output b1 ;else output b2;run；例2例例3data a;input id x1 x2 x3 ;cards;1 0 12 12 1 55 23 0 15 24 1 14 35 0 42 16 0 25

24、27 0 35 48 0 27 3; data c1;set a;if x1=1 then delete ;proc print;run;结果：OBS ID X1 X2 X3 1 1 0 12 1 2 3 0 15 2 3 5 0 42 1 4 6 0 25 2 5 7 0 35 4 6 8 0 27 3 例例4data c2;set a;if x1=0 then output ;if x3=2 then x3=2; （执行了output语句，此语句不再执行。）proc print;run; 结果： OBS ID X1 X2 X3 1 1 0 12 1 2 3 0 15 2 3 5 0 42

25、1 4 6 0 25 2 5 7 0 35 4 6 8 0 27 3 例例5data c3;set a; if x3=2 then x3=2;if x1=0 then output ;proc print;run;结果：OBS ID X1 X2 X3 1 1 0 12 1 2 3 0 15 2 3 5 0 42 1 4 6 0 25 2 5 7 0 35 2 6 8 0 27 2 SAS数据集的整理数据集的整理-记录操作记录操作q数据排序数据排序 对数据集数据进行排序。 语句格式：proc sort 选择项； by descending 变量名； sort过程按by变量名的值进行排序。 选择项

26、：1.如果对非当前数据集排序，须用data指出要排序的数据集名。如：proc sort data=L.px1; 2.排序后产生一新数据集，数据集名可用“out=“指出，如：proc sort out=L.px2; 如缺省，则原数据集将被输出的数据集替换。 descending 选择为降序，不选择为升序。如： by descending x1; 例例1. Libname Ld:sassas1;data L.class1;input name $ kc1 kc2 kc3 kc4;sum=(kc1+kc2+kc3+kc4);cards;zhao 83 76 68 88qian 76 77 80 69

27、sun 75 82 78 78li 80 78 85 81;proc print;proc sort;by descending sum;proc print;Proc sort out=L.class2;by sum;proc print;run;OBS NAME KC1 KC2 KC3 KC4 SUM 1 zhao 83 76 68 88 315 2 qian 76 77 80 69 302 3 sun 75 82 78 78 313 4 li 80 78 85 81 324OBS NAME KC1 KC2 KC3 KC4 SUM 1 li 80 78 85 81 324 2 zhao 8

28、3 76 68 88 315 3 sun 75 82 78 78 313 4 qian 76 77 80 69 302OBS NAME KC1 KC2 KC3 KC4 SUM 1 qian 76 77 80 69 302 2 zhao 83 76 68 88 315 3 sun 75 84 78 78 315 4 li 80 78 85 81 324q数据集的纵向合并 将两个或多个数据集的变量值纵向合 并。 语句格式：set ; 例如：有两个数据集例如：有两个数据集A和和B，将两者纵向合并成新数据集，将两者纵向合并成新数据集AB的程序如下：的程序如下： data AB; set A b; ru

29、n;SAS数据集的整理数据集的整理-记录操作记录操作 若A、B两数据集的变量完全相同，则新数据集AB的变量与两个数据集变量一致；若A、B量数据集的变量不完全一致，则新数据集AB除含有两个数据集共同的变量外，还包括各自的变量，对于原数据集中没有的变量，其值为缺失值。例如，数据集A、B的变量值分别如下： 数据集 A中的数据 NO AGE SEX HEIGHT 01 16 m 178 02 16 m 170 03 17 m 175 04 18 m 169 数据集B中的数据NO AGE SEX HEIGHT WEIGHT 16 f 160 55 17 f 165 6011 17 f 159 50新数据

30、集AB的数据OBS NO AGE SEX HEIGHT WEIGHT 1 01 16 M 178 . 2 02 16 M 170 . 3 03 17 M 175 . 4 04 18 M 175 . 5 11 16 F 160 55 6 12 17 F 165 60 7 13 17 F 159 50q数据集的横向合并 将两个或多个数据集的变量值横向合 并。 v语句格式 merge ; by 识别字段;使用by 识别字段时，要先按识别字段排序。 例data a;input id x1 x2 ;cards;1 11 122 21 224 41 425 51 52;data b;input id x3

31、;cards;1 133 334 43data ab;merge a b ;by id ;proc print;run;OBS ID X1 X2 X3 1 1 11 12 13 2 2 21 22 . 3 3 . . 33 4 4 41 42 43 5 5 51 52 .例data a; input id x1 x2 ;cards;1 11 122 21 224 41 424 411 4225 51 52;data b; input id x3;cards;1 133 334 435 535 533;data c;merge a b;by id;proc print;run;OBS ID X1

32、 X2 X3 1 1 11 12 13 2 2 21 22 . 3 3 . . 33 4 4 41 42 43 5 4 411 422 43 6 5 51 52 53 7 5 51 52 533 q在PROC过程的几个常用语句vWHERE语句 PROC过程步中定义表达式，在分析之前用该表达式选择符合条件的变量值作为分析的对象。语句格式 where ；例 data a; input id x1 x2 ; cards; 1 11 12 2 21 22 3 41 42 4 45 44 5 51 52 ;proc print;where x120;run; OBS ID X1 X2 2 2 21 22

33、 3 3 41 42 4 4 45 44 5 5 51 52v SUM语句 语句格式 SUM 变量表 对SUM语句中的数值变量求合并打印输出。例data abc;input a b c ;cards;21 22 26 5 7 4 12 3 67 16 23 19 2 3;proc print;sum a b c;run; OBS A B C 1 21 22 26 2 5 7 4 3 12 3 6 4 7 16 23 5 19 2 3 = = = 64 50 62v标题title语句 在PROC过程步中指定打印到输出文件上的标题内容。 语句格式 title n 标题内容 n为标题行，n的范围从1

34、到10，缺省值为1。标题内容是要打印的标题字符串，字符不超过132个。 一旦为一行规定了一个标题，则在后继的输出中都使用该标题，直至取消该标题或为该行定义其它的标题为止。若要取消现存的标题，使用下面语句： title n ； 如 title; title2; 等。例data a;input x ;cards;4 5 6 7 8 9 3 4 5;data b;input x ;cards;45 34 56 33 44 54 32 36;data ab;set a b;proc means data=a;title sample a;proc means data=b;title sample b

35、;proc means data=ab;title sample ab;proc print data=a;title;run; sample a Analysis Variable : X N Mean Std Dev Minimum Maximum - 9 5.6666667 2.0000000 3.0000000 9.0000000 - sample b Analysis Variable : X N Mean Std Dev Minimum Maximum - 8 41.7500000 9.5131488 32.0000000 56.0000000 - sample ab Analysis Variable : X N Mean Std Dev Minimum Maximum - 17 22.6470588 19.6530569 3.0000000 56.0000000 - OBS X 1 4 2 5 3 6 4 7 5 8 6 9 7 3 8 4 9 5