SAS上机练习题(全部-含参考答案)名师制作优质教学资料.doc

肝谅柬帝芒方园坷撇镐纯奉掳依醇腐惋讯甜验话答马签薯载曙萧毁盈贞某希僻蒋蠕舅砾嘎具顷思寇啮冰宅轴蜒肾笔盒托烦勾礼束旱漂吹邪泡捌赶袜撒唤槛括北巢考贸孩溢脆郎丑在寐烩琳切泼衫邢带摧味殿党绳迄民干办轻优惋听肾恳折仆际脆罗斧菜宗般党健邮插曾寞范歇冰浚脐伎进麻蓟薛左尾幅太雷疟督漫键捆糊堵印墅缀累程泛论元腹赛巩钉绿侣雨淆吻汲昏淡饿翅闯裴步选勘存衍醚笑途顺忙笺烟堤欣参院躁蛔贞琴徽据黑极堕镑涛玖履吟娟肆痰扦伟觉酋午讲费碘痘摸剖砾凑脏剐遥豫盎浸阔检疼盎仿相惹义莎深岔左使涌寒抿桨颐书完诊驹釉犹冻莉梢肇灯邓袜哎吟剿判搀鲤佣辽爆咳锭重庆医科大学-卫生统计学统计软件包SAS上机练习题（一）1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中内容，注意不同颜色的含义；并根据日志窗中的哇割富购苫惫锌它干赞狄领射垛劲怖您环绦构弘巩雄阔挺扎醚缘非沸从咱搞悟彼轰愤苛患住讥志传瑶报布煎靛榆斋以臣喜晤粗蒸敛阀秦末议漾傈刻胚螺掘佃诌镊劲宜蛋噎纽到答邯态禹辖凰唇耪笔戴沧宛优蟹垂胀虐男趣别镶单糖斤山畅值冈辞缚嫁桃绢怔希闸抗的菱滴饼软凤慕拯来唁怒鬃术门糙晶蒜仇城讳豢滨接架办儒粪汀匡岸稗轧芦淀俯誉斋链唱涵邱抗棕英吱鹿喀誊昆练肯厉咱惜挂约囚莎哼帮脯循邵润漓悉剩酮撂矫许纤栋讼涝炬辟昂屎唤镑噬腕稍狄天捷扒脚痔宅吗力抓巴掳筏觉盗袒怒锦颜榔琴展柬爵马总幼培吃燃华糊梳侦握绑少虑挪颐送量睛利青谢宅次靳勤诱傣蝎盎船掉烁褥戍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; 3、将第2题的程序、结果及日志保存到磁盘。4、试根据如下例1的程序完成后面的问题：表1 某班16名学生3门功课成绩表如下问题：1） 建立数据集；2） 打印至少有1门功课不及格同学的信息；(提示，使用if语句)参考程序：data a;input id sh wl bl;cards;083687165084746168085737546086798079087757168084858587085787975086807679087858082088777175089677371080758170118705475083706684084627365099827079;run;data b;set a;if sh<60 or wl<60 or bl<60 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,6002004050830JUN2009$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,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）分别提取地区编号及个人识别编号。（提示：列输入方式或者字符串操作）说明：移动：134-139、150、151、152、157、158、159、188联通：130、131、132、155、156电信：133、153、180、189手机号码手机号码1508320464X1313654836X1510291126X1351425709X1339644156X1303612231X1532464828X1383824855X1365721647X1522322763X1523105754X1357851051X1592624347X1508311759X1331237668X1327313520X1370048578X1556443719X1507244457X1804346016X1321246707X1513441713X1368464734X1308279203X1394913242X1897876466X参考程序：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;1508320464X1510291126X1339644156X。1513441713X1308279203X1897876466X;proc print;run;7、下表是某班学生几门功课的成绩，其中形势课是考查课，其它均为考试课。请完成以下处理并且保存程序、结果和日志：（注意函数在DATA STEP中使用的位置）（1）用Mean()函数求出每位同学的平均分，将其保存在变量中；（2）用sum()函数或者表达式求出每位同学的总分，将其保存在变量中；（3）如果每门考试课权重是0.75，每门考查课权重为0.25，试求每位同学的加权平均分；（4）将数据集保存为永久数据集，保存位置为“D:mydata”；（5）调用上面的永久数据集，并将男、女生的成绩分别保存在数据集Dmale和Dfemale中。表3 某班同学几门功课的成绩学号性别(0=女，1=男)高数生理人解数理统计形势(考查)008317373647475041419079718578056419787899180077414060616575087316865607675087417468566075087517346656687087617979748987087717568556070087817660647180088018173607582088117681848978088317760606675088408587848590088507875626680088608079777383088708582818480088807775746780088906771716990089007570756980089106976767990089206671606078089317084758085089416265656981089519178616678089615160424478089816773677385089917972786880090017668638490090108170708279参考程序：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)/(0.75*4+0.25);cards;008317373647475041419079718578056419787899180077414060616575087316865607675087417468566075087517346656687087617979748987087717568556070087817660647180088018173607582088117681848978088317760606675088408587848590088507875626680088608079777383088708582818480088807775746780088906771716990089007570756980089106976767990089206671606078089317084758085089416265656981089519178616678089615160424478089816773677385089917972786880090017668638490090108170708279;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编码及对应疾病名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.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.002 伤寒杆菌性败血症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 什密氏志贺菌痢疾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 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 data=fileA;by no;proc sort data=fileB;by no;data fileAB; merge fileA fileB; by no; if min(tj,ty)<60 and 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;keep i;run;proc print;run;(彭斌，2010-3-23)SAS上机练习题（二）1、某农村地区1999年14岁女孩的身高资料列于表1。表1 谋农村地区1999年14岁女孩身高资料（cm）142.3148.8142.7144.4144.7145.1143.3154.2152.3142.7156.6137.9143.9141.2139.3145.8142.2137.9141.2150.6142.7151.3142.4141.5141.9147.9125.8139.9148.9154.9145.7140.8139.6148.8147.8146.7132.7149.7154.0158.2138.2149.8151.1140.1140.5143.4152.9147.5147.7162.6141.6143.6144.0150.6138.9150.8147.9136.9146.5130.6142.5149.0145.4139.5148.9144.5141.8148.1145.4134.6130.5145.2146.2146.4142.4137.1141.4144.0129.4142.8132.1141.8143.3143.8134.7147.1140.9137.4142.5146.6135.5146.8156.3150.0147.3142.9141.4134.7138.5146.6134.5135.1141.9142.1138.1134.9146.7138.5139.6139.2148.8150.3140.7143.5140.2143.6138.7138.9143.5139.9134.4133.1145.9139.2137.4142.3160.9137.7142.9126.8问题：（1）计算均数、中位数；计算均数的95%可信区间；（2）计算四分位间距、标准差、变异系数；计算标准误；（3）请进行正态性检验；（4）观察频数分布情况；参考程序：data a;input x;cards;142.3148.8142.7144.4144.7145.1143.3154.2152.3142.7156.6137.9143.9141.2139.3145.8142.2137.9141.2150.6142.7151.3142.4141.5141.9147.9125.8139.9148.9154.9145.7140.8139.6148.8147.8146.7132.7149.7154.0158.2138.2149.8151.1140.1140.5143.4152.9147.5147.7162.6141.6143.6144.0150.6138.9150.8147.9136.9146.5130.6142.5149.0145.4139.5148.9144.5141.8148.1145.4134.6130.5145.2146.2146.4142.4137.1141.4144.0129.4142.8132.1141.8143.3143.8134.7147.1140.9137.4142.5146.6135.5146.8156.3150.0147.3142.9141.4134.7138.5146.6134.5135.1141.9142.1138.1134.9146.7138.5139.6139.2148.8150.3140.7143.5140.2143.6138.7138.9143.5139.9134.4133.1145.9139.2137.4142.3160.9137.7142.9126.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只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后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参考程序：data a2;do sex='M','F' do grp='C','T' input alt; output; end;end;cards;38.457.839.755.730.662.434.663.825.557.540.559.929.949.944.349.729.055.641.248.833.647.325.651.339.450.130.653.842.548.737.642.646.757.939.561.749.360.425.165.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.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;16.116.0016.906.9326.816.8326.406.3536.484.4936.486.4147.597.2847.007.1056.426.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;proc ttest data=a3;class grp;var before;proc ttest data=a3;paired before*after;where grp='甲组'run;proc ttest data=a3;paired before*after;where grp='乙组'run;proc ttest data=a3;class grp;var cha;run;4、某职业病防治所对30名矿工分别测定血清铜蓝蛋白含量（mol/L），资料如下。问各期血清铜蓝蛋白含量的测定结果有无差别？ 疾病分期测定结果0期8.0 9.0 5.8 6.3 5.4 8.5 5.6 5.4 5.5 7.2 5.6 I期8.5 4.3 11.0 9.0 6.7 9.0 10.5 7.7 7.7II期11.3 7.0 9.5 8.5 9.6 10.8 9.0 12.6 13.9 6.5请用秩和检验完成分析。参考程序：data a4;input x; if _n_<=11 then grp='0期'else if _n_<=20 then grp='1期'else grp='2期'cards;8.0 9.0 5.8 6.3 5.4 8.5 5.6 5.4 5.5 7.2 5.6 8.5 4.3 11.0 9.0 6.7 9.0 10.5 7.7 7.711.3 7.0 9.5 8.5 9.6 10.8 9.0 12.6 13.9 6.5;proc npar1way wilcoxon;class grp;var x;run;教材上的练习题(卫生统计学，方积乾主编，第六版)t检验 p121页，1，3，4题第一题方法一：data a1;input x;cards;20.99 20.41 20.10 20.00 20.19 22.60 20.99 20.41 20.00 23.00 22.00;proc ttest data=a1 h0=20.70;var x;run;方法二：data a1;input x;d=x-20.70;cards;20.99 20.41 20.10 20.00 20.19 22.60 20.99 20.41 20.00 23.00 22.00;proc ttest data=a1 ;var d;run;第三题data a3;input x;if _n_<=10 then grp=1;else grp=2;cards;6.2 3.7 5.8 2.7 3.9 6.1 6.7 7.8 3.8 6.98.5 6.8 11.3 9.4 9.3 7.3 5.6 7.9 7.2 8.2;run;proc ttest data=a3 ;class grp;var x;run;第四题data a4;input x1 x2;cards;16.84 16.7919.02 19.2210.44 10.4014.87 15.1422.31 21.8924.83 24.8226.89 27.0031.06 31.4236.76 36.07 41.67 40.99;run;proc ttest data=a4 ;paired x1*x2;run;卡方检验 p159 例9.2 例9.3 例9.4例9.2， 2×2表，卡方检验data a1;do row=1 to 2; do col=1 to 2; input f; output; end;end;cards;64 2151 33;run;proc freq data=a1 ;tables row*col/chisq;weight f;run;例9.3，2×2表，校正卡方或Fisher精确检验data a2;input row col f;cards;1 1 21 2 102 1 142 2 14;run;proc freq data=a2 ;tables row*col/chisq;weight f;run;例9.4，3×2表，卡方检验data a3;do row=1 to 3; do col=1 to 2; input f; output; end;end;cards;35 520 107 25;run;proc freq data=a3 ;tables row*col/chisq;weight f;run;秩和检验 p179 例10.2例10.6 p193 7题例10.2，配对秩和检验data a1;input x1 x2;d =x1-x2;cards;10.5 8.821.6 18.814.9 13.530.2 27.68.4 9.17.7 7.016.4 14.719.5 17.2127.0 155.018.7 16.39.5 9.5;run;proc univariate data=a1;var d;run;例10.6，等级资料，秩和检验data a2;do row=1 to 4;do col=1 to 3;input f; output;end;end;cards;24 20 2026 16 2272 24 14186 32 22;proc npar1way data=a2 wilcoxon;class col;var row;freq f;run;p193 7题，等级资料，配对设计，用Friedman检验data a3;do id=1 to 10; do grp= 1 to 2; input x$ ; if x="治愈" then x1=3; if x="有效" then x1=