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

1、SAS上机练习题全部含参考答案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;

2、VAR NAME AGE; RUN; 2004050830JUN2009$432,334,500提示：(格式化输入；数据之间以空格分隔，数据对齐；注意格式后面的长度应以前一个位置结束开始计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式) 开始时间，输入格式yymmdd10. 终止时间，输入格式date10. 费用，输入格式dollar12. 参考程序：data a;input x1 yymmdd10. x2 date10. x3 dollar13.;cards;2005/04/28 25MAY2009 $123,345,0002005 09 18 05OCT2009 $33,

3、234,500 2007/08/12 22SEP2009 $345,600 20040508 30JUN2009 $432,334,500;run;proc print;run;proc print;format x1 yymmdd10. x2 date9. x3 dollar13.;run;6、手机号码一编码规则一般是：YYY-XXXX-ZZZZ，其YYY为号段；XXXX一般为所在地区编码；ZZZZ为对应的个人识别编号。下面有一组电话号码（来源于网络，末位以X替换），请用程序完成下列要求：（1）分别列出属于联通、移动、电信的号码；（2）分别提取地区编号及个人识别编号。（提示：列输入方式或者字

4、符串操作）说明：移动：134-139、150、151、152、157、158、159、188联通：130、131、132、155、156电信：133、153、180、189手机号码手机号码1508320464X1313654836X1510291126X1351425709X1339644156X1303612231X1532464828X1383824855X1365721647X1522322763X1523105754X1357851051X1592624347X1508311759X1331237668X1327313520X1370048578X1556443719X15072444

5、57X1804346016X1321246707X1513441713X1368464734X1308279203X1394913242X1897876466X参考程序：data 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;1508320464X1510291126X13396

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

7、别保存在数据集Dmale和Dfemale中。表3 某班同学几门功课的成绩学号性别(0=女，1=男)高数生理人解数理统计形势(考查)0083173736474750414190797185780564197878991800774140606165750873168656076750874174685660750875173466566870876179797489870877175685560700878176606471800880181736075820881176818489780883177606066750884085878485900885078756266800886080797

8、77383088708582818480088807775746780088906771716990089007570756980089106976767990089206671606078089317084758085089416265656981089519178616678089615160424478089816773677385089917972786880090017668638490090108170708279参考程序：libname pb D:mydata;data pb.score;input id$ xb gs sl rj sltj xs;mscore=mean(of g

9、s sl rj sltj xs);tscore=sum(of gs sl rj sltj xs);wmscore= (gs+sl+rj+sltj)*0.75+xs*0.25)/(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 70087

10、8 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 85 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

11、 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 78 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个大类的

12、疾病存入3个数据集。（提示：ICD10编码中小数点前面的三位表示大类；length语句定义字符变量长度；字符串取子串函数）表44类疾病的ICD10编码及对应疾病名ICD10DiseaseA01.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.

13、202+沙门氏菌性关节炎A02.203+沙门氏菌性脑膜炎A02.901沙门氏菌感染A02.902鼠伤寒沙门氏菌感染A02.903沙门氏菌属食物中毒A03.001什密氏志贺菌痢疾A03.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.00

14、2 伤寒杆菌性败血症A01.003+ 伤寒性脑膜炎A01.101 甲型副伤寒A01.201 乙型副伤寒A01.301 丙型副伤寒A01.401 副伤寒A02.001 B群沙门氏菌肠炎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 什密氏志贺菌痢疾

15、A03.101 弗氏志贺菌痢疾A03.201 鲍氏志贺菌痢疾A03.301 宋内氏志贺菌痢疾A03.802 菌痢混合感染A03.901 细菌性痢疾A03.902 慢性细菌性痢疾急性发作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

16、 56 96 2203 93 75 1904 77 84 2405 88 55 18B文件学号 性别 班级01 男 102 女 103 男 304 男 2 05 女 208 女 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 d

17、ata=fileA;by no;proc sort data=fileB;by no;data fileAB; merge fileA fileB; by no; if min(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;

18、keep i;run;proc print;run;(彭斌，2010-3-23)SAS上机练习题（二）1、某农村地区1999年14岁女孩的身高资料列于表1。表1 谋农村地区1999年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.61

19、48.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 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

20、 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.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%可信区间；

21、（2）计算四分位间距、标准差、变异系数；计算标准误；（3）请进行正态性检验；（4）观察频数分布情况；参考程序：data a;input x;cards;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

22、 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 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.

23、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.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 stder

24、r;var x;proc univariate data=a normal;histogram x;var x;run;2、某研究员欲分析急性染毒对肝脏功能的影响，将40只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后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

25、.8633.647.325.651.3739.450.130.653.8842.548.737.642.6946.757.939.561.71049.360.425.165.3参考程序：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

26、 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.

27、3056.536.4166.946.6466.706.6879.178.4279.109.0587.337.0087.316.8396.946.5896.966.91107.677.22106.816.73118.156.57118.167.65126.606.17126.986.52问题：（1）试建立能够满足以下统计分析的SAS数据集。（2）试判断两组受试者试验前血清胆固醇水平是否相等？（3）分别判断两种降胆固醇措施是否有效？（4）试判断两种降胆固醇措施的效果是否相同？参与程序：data a3;do grp=甲组,乙组; input id before after ; cha=before-after; output;end;cards;1 6.11 6.00 1 6.90 6.932 6.81 6.83 2 6.40 6.353 6.48 4.49 3 6.48 6.414 7.59 7.28 4 7.00 7.105 6