SAS上机练习题全部含参考答案.docx

1、SAS上机练习题全部含参考答案重庆医科大学-卫生统计学统计软件包SAS上机练习题（一）1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中内容，注意不同颜色的含义；并根据日志窗中的信息修改完善程序。 DATS EX0101; INPUTT NAME $ AGE SEX; CARDS; XIAOMIN 19 1 LIDONG 20 1 NANA 18 2 ; PROD PRONT DATS=EX1; RUN; PROC PRINT DATA=EX1; VAR NAME AGE; RUN;

2、 3、将第2题的程序、结果及日志保存到磁盘。4、试根据如下例1的程序完成后面的问题：表1 某班16名学生3门功课成绩表如下问题：1） 建立数据集；2） 打印至少有1门功课不及格同学的信息；(提示，使用if语句)参考程序：data a;input id sh wl bl;cards;083 68 71 65084 74 61 68085 73 75 46086 79 80 79087 75 71 68084 85 85 87085 78 79 75086 80 76 79087 85 80 82088 77 71 75089 67 73 71080 75 81 70118 70 54 75083

3、 70 66 84084 62 73 65099 82 70 79;run;data b;set a;if sh60 or wl60 or bl60 then output;run;proc print data=b;var id sh wl bl;run;5、根据下列数据建立数据集表2销售数据开始时间终止时间费用2005/04/2825MAY2009$123,345,0002005 09 1805OCT2009$33,234,5002007/08/1222SEP2009$345,60030JUN2009$432,334,500提示：(格式化输入；数据之间以空格分隔，数据对齐；注意格式后面的长

4、度应以前一个位置结束开始计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式) 开始时间，输入格式yymmdd10. 终止时间，输入格式date10. 费用，输入格式dollar12. 参考程序：data a;input x1 yymmdd10. x2 date10. x3 dollar13.;cards;2005/04/28 25MAY2009 $123,345,0002005 09 18 05OCT2009 $33,234,500 2007/08/12 22SEP2009 $345,600 20040508 30JUN2009 $432,334,500;run;proc pr

5、int;run;proc print;format x1 yymmdd10. x2 date9. x3 dollar13.;run;6、手机号码一编码规则一般是：YYY-XXXX-ZZZZ，其YYY为号段；XXXX一般为所在地区编码；ZZZZ为对应的个人识别编号。下面有一组电话号码（来源于网络，末位以X替换），请用程序完成下列要求：（1）分别列出属于联通、移动、电信的号码；（2）分别提取地区编号及个人识别编号。（提示：列输入方式或者字符串操作）说明：移动：134、139、150、151、152、157、158、159、188联通：130、131、132、155、156电信：133、153、18

6、0、189手机号码手机号码1508320464X13X1510291126X1351425709X1339644156X1303612231X1532464828X1383824855X1365721647X1522322763X1523105754X1357851051X1592624347X1508311759X1331237668X1327313520X1370048578X1556443719X1507244457X1804346016X1321246707X1513441713X1368464734X1308279203X1394913242X1897876466X参考程序：data

7、 a;input y 1-3 x $ 4-7 z $ 8-11;if 134=y=139 or y in (150,151,152,157,158,159,188) then ygrp=移动;if y in (130,131,132,155,156) then ygrp=联通;if y in (133,153,180,189) then ygrp=电信;cards;1508320464X1510291126X1339644156X。1513441713X1308279203X1897876466X;proc print;run;7、下表是某班学生几门功课的成绩，其中形势课是考查课，其它均为考试

8、课。请完成以下处理并且保存程序、结果和日志：（注意函数在DATA STEP中使用的位置）（1）用Mean()函数求出每位同学的平均分，将其保存在变量中；（2）用sum()函数或者表达式求出每位同学的总分，将其保存在变量中；（3）如果每门考试课权重是0.75，每门考查课权重为0.25，试求每位同学的加权平均分；（4）将数据集保存为永久数据集，保存位置为“D:mydata”；（5）调用上面的永久数据集，并将男、女生的成绩分别保存在数据集Dmale和Dfemale中。表3 某班同学几门功课的成绩学号性别(0=女，1=男)高数生理人解数理统计形势(考查)008317373647475041419079

9、718578056419787899180077414060616575087316865607675087417468566075087517346656687087617979748987087717568556070087817660647180088018173607582088117681848978088317760606675088408587848590088507875626680088608079777383088708582818480088807775746780088906771716990089007570756980089106976767990089206671

10、606078089317084758085089416265656981089519178616678089615160424478089816773677385089917972786880090017668638490090108170708279参考程序：libname pb D:mydata;data pb.score;input id$ xb gs sl rj sltj xs;mscore=mean(of gs sl rj sltj xs);tscore=sum(of gs sl rj sltj xs);wmscore= (gs+sl+rj+sltj)*0.75+xs*0.25)/(

11、0.75*4+0.25);cards;0083 1 73 73 64 74 750414 1 90 79 71 85 780564 1 97 87 89 91 800774 1 40 60 61 65 750873 1 68 65 60 76 750874 1 74 68 56 60 750875 1 73 46 65 66 870876 1 79 79 74 89 870877 1 75 68 55 60 700878 1 76 60 64 71 800880 1 81 73 60 75 820881 1 76 81 84 89 780883 1 77 60 60 66 750884 0 8

12、5 87 84 85 900885 0 78 75 62 66 800886 0 80 79 77 73 830887 0 85 82 81 84 800888 0 77 75 74 67 800889 0 67 71 71 69 900890 0 75 70 75 69 800891 0 69 76 76 79 900892 0 66 71 60 60 780893 1 70 84 75 80 850894 1 62 65 65 69 810895 1 91 78 61 66 780896 1 51 60 42 44 780898 1 67 73 67 73 850899 1 79 72 7

13、8 68 800900 1 76 68 63 84 900901 0 81 70 70 82 79;data dmale dfemale;set pb.score;if xb=1 then output dmale;if xb=0 then output dfemale;run;8、下面是3个大类疾病的ICD-10编码及对应的疾病名。请完成以下任务：（1）建立数据集；（2）提取每种疾病的大类编码；（3）分别将3个大类的疾病存入3个数据集。（提示：ICD10编码中小数点前面的三位表示大类；length语句定义字符变量长度；字符串取子串函数）表44类疾病的ICD10编码及对应疾病名ICD10Dis

14、easeA01.001伤寒A01.002伤寒杆菌性败血症A01.003+伤寒性脑膜炎A01.101甲型副伤寒A01.201乙型副伤寒A01.301丙型副伤寒A01.401副伤寒A02.001B群沙门氏菌肠炎A02.002C群沙门氏菌肠炎A02.004沙门氏菌性肠炎A02.006沙门氏菌胃肠炎A02.007鼠伤寒沙门氏菌性肠炎A02.008婴儿沙门氏菌肠炎A02.101沙门氏菌败血症A02.201+沙门氏菌性肺炎A02.202+沙门氏菌性关节炎A02.203+沙门氏菌性脑膜炎A02.901沙门氏菌感染A02.902鼠伤寒沙门氏菌感染A02.903沙门氏菌属食物中毒A03.001什密氏志贺菌痢疾A

15、03.101弗氏志贺菌痢疾A03.201鲍氏志贺菌痢疾A03.301宋内氏志贺菌痢疾A03.802菌痢混合感染A03.901细菌性痢疾A03.902慢性细菌性痢疾急性发作A03.903慢性细菌性痢疾A03.905慢性迁延型细菌性痢疾A03.906中毒性痢疾参考程序：data a;input ICD10 $ Disease:$24.;cat=substr(ICD10,1,3);cards;A01.001 伤寒A01.002 伤寒杆菌性败血症A01.003+ 伤寒性脑膜炎A01.101 甲型副伤寒A01.201 乙型副伤寒A01.301 丙型副伤寒A01.401 副伤寒A02.001 B群沙门氏菌

16、肠炎A02.002 C群沙门氏菌肠炎A02.004 沙门氏菌性肠炎A02.006 沙门氏菌胃肠炎A02.007 鼠伤寒沙门氏菌性肠炎A02.008 婴儿沙门氏菌肠炎A02.101 沙门氏菌败血症A02.201+ 沙门氏菌性肺炎A02.202+ 沙门氏菌性关节炎A02.203+ 沙门氏菌性脑膜炎A02.901 沙门氏菌感染A02.902 鼠伤寒沙门氏菌感染A02.903 沙门氏菌属食物中毒A03.001 什密氏志贺菌痢疾A03.101 弗氏志贺菌痢疾A03.201 鲍氏志贺菌痢疾A03.301 宋内氏志贺菌痢疾A03.802 菌痢混合感染A03.901 细菌性痢疾A03.902 慢性细菌性痢疾急

17、性发作A03.903 慢性细菌性痢疾A03.905 慢性迁延型细菌性痢疾A03.906 中毒性痢疾;data cat01 cat02 cat03;set a;if cat=A01 then output cat01;if cat=A02 then output cat02;if cat=A03 then output cat03;run;9、现有两个文件，内容如下：A文件学号 统计 体育 年龄01 87 83 2302 56 96 2203 93 75 1904 77 84 2405 88 55 18B文件学号 性别 班级01 男 102 女 103 男 304 男 2 05 女 208 女

18、3要求：打印至少有一门功课不及格的同学的年龄、性别和班级。参考程序：data fileA;input no$ tj ty age;cards;01 87 83 2302 56 96 2203 93 75 1904 77 84 2405 88 55 18;data fileB;input no$ sex$ bj;cards;01 男 102 女 103 男 304 男 2 05 女 208 女 3;proc sort data=fileA;by no;proc sort data=fileB;by no;data fileAB; merge fileA fileB; by no; if min(

19、tj,ty) . ;run;proc print;run;10、显示100到200以内的素数。（提示：素数即只能被1和本身整除的数，如2，3，5，7等等；循环语句，条件语句，取余数函数mod(a,b)）参考程序：data a;do i=100 to 200; flag=0; do j=2 to i-1; if mod(i,j)=0 then flag=1; end; if flag=0 then output;end;keep i;run;proc print;run;(彭斌，2010-3-23)SAS上机练习题（二）1、某农村地区1999年14岁女孩的身高资料列于表1。表1 谋农村地区199

20、9年14岁女孩身高资料（cm）142.3 148.8 142.7 144.4 144.7 145.1 143.3 154.2 152.3 142.7 156.6137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 150.6 142.7 151.3142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 145.7 140.8 139.6148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 149.8 151.1 140.1140.5 143.4 152.9 147.5 1

21、47.7 162.6 141.6 143.6 144.0 150.6 138.9150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 139.5 148.9 144.5141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 142.4 137.1 141.4144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 147.1 140.9 137.4142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 141.4 134.7 138.5

22、146.6 134.5 135.1 141.9 142.1 138.1 134.9 146.7 138.5 139.6 139.2148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 143.5 139.9 134.4133.1 145.9 139.2 137.4 142.3 160.9 137.7 142.9 126.8问题：（1）计算均数、中位数；计算均数的95%可信区间；（2）计算四分位间距、标准差、变异系数；计算标准误；（3）请进行正态性检验；（4）观察频数分布情况；参考程序：data a;input x;cards;142.3 148.8 1

23、42.7 144.4 144.7 145.1 143.3 154.2 152.3 142.7 156.6137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 150.6 142.7 151.3142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 145.7 140.8 139.6148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 149.8 151.1 140.1140.5 143.4 152.9 147.5 147.7 162.6 141.6 143.6 144.0

24、 150.6 138.9150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 139.5 148.9 144.5141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 142.4 137.1 141.4144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 147.1 140.9 137.4142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 141.4 134.7 138.5146.6 134.5 135.1 141.9 142.

25、1 138.1 134.9 146.7 138.5 139.6 139.2148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 143.5 139.9 134.4133.1 145.9 139.2 137.4 142.3 160.9 137.7 142.9 126.8;proc means data=a n mean median clm qrange std cv stderr;var x;proc univariate data=a normal;histogram x;var x;run;2、某研究员欲分析急性染毒对肝脏功能的影响，将40只小鼠分

26、为两组，雄雌各半，试验组进行急性染毒试验，染毒后2小时测定血液中的ALT，整理的结果见下表。 （1） 比较实验组与对照组的ALT有无差异；ALT(丙氨酸转氨酶)NOmale control TreaetmentFemalecontrol Treaetment138.457.839.755.7230.662.434.663.8325.557.540.559.9429.949.944.349.7529.055.641.248.8633.647.325.651.3739.450.130.653.8842.548.737.642.6946.757.939.561.71049.360.425.165.3

27、参考程序：data a2;do sex=M,F; do grp=C,T; input alt; output; end;end;cards;38.4 57.8 39.7 55.730.6 62.4 34.6 63.825.5 57.5 40.5 59.929.9 49.9 44.3 49.729.0 55.6 41.2 48.833.6 47.3 25.6 51.339.4 50.1 30.6 53.842.5 48.7 37.6 42.646.7 57.9 39.5 61.749.3 60.4 25.1 65.3;proc anova data=a2;class sex grp;model alt=sex grp;means grp;run;quit;3、24名志愿者随机分成两组，每组12人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L），数据见下表：甲组乙组受试者试验前试验后受试者试验前试验后16.116.0016.906.9326.816.8326.406.3536.484.4936.486.4147.597.2847.007.1056.426.3