SAS上机练习题(全部,含参考答案)

1、重庆医科大学-卫生统计学统计软件包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题的程序、结果及日

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

3、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提示：(格式化输入；数据之间以空格分隔，数据对齐；注意格式后面的长度应以前一个位置结束开始计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出

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 print;run;proc print;format x1 yymmdd10. x2 d

5、ate9. 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手机号码手机号码1508320464X1313654836X15102911

6、26X1351425709X1339644156X1303612231X1532464828X1383824855X1365721647X1522322763X1523105754X1357851051X1592624347X1508311759X1331237668X1327313520X1370048578X1556443719X1507244457X1804346016X1321246707X1513441713X1368464734X1308279203X1394913242X1897876466X参考程序：data a;input y 1-3 x $ 4-7 z $ 8-11;if

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

11、+xs*0.25)/(0.75*4+0.25);cards;008317373647475041419079718578056419787899180077414060616575087316865607675087417468566075087517346656687087617979748987087717568556070087817660647180088018173607582088117681848978088317760606675088408587848590088507875626680088608079777383088708582818480088807775746780

12、088906771716990089007570756980089106976767990089206671606078089317084758085089416265656981089519178616678089615160424478089816773677385089917972786880090017668638490090108170708279;data dmale dfemale;set pb.score;if xb=1 then output dmale;if xb=0 then output dfemale;run;8、下面是3个大类疾病的ICD-10编码及对应的疾病名。请

13、完成以下任务：（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.00

14、8婴儿沙门氏菌肠炎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.;ca

15、t=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 鼠伤寒沙

16、门氏菌感染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&

17、#39; 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

18、 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

19、=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.614

20、2.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.813

21、2.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）计算四分位间距、标准差、变异系数；计算标准

22、误；（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

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

24、.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） 比较

25、实验组与对照组的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=&#

26、39;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名志愿者随机分成两组

27、，每组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.17

28、126.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.41

29、66.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='乙组'r

30、un;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

31、='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

32、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;els

33、e 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=a

34、4 ;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 d

35、ata=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

36、.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题，等级资料，配对设

37、计，用Friedman检验data a3;do id=1 to 10; do grp= 1 to 2; input x$ ; if x="治愈" then x1=3; if x="有效" then x1=2; if x="无效" then x1=1; output; end;end;cards;治愈 有效有效 无效治愈 有效治愈 治愈有效 有效治愈 有效治愈 无效治愈 有效有效 无效治愈 治愈;proc freq data=a3;table id*grp*x1/cmh scores=rank noprint;run; (彭斌，2010

38、-4-18)SAS上机练习题（三）1、（方差分析）某研究人员研究3种刺激因素（高葡萄糖、高胰岛素和过氧化氢）刺激大鼠肾小球系膜细胞对单核细胞趋化蛋白表达的影响，在四种条件下对大鼠肾小球系膜细胞进行培养，检测结果如表1。表1 不同刺激因素作用下蛋白表达结果分组表达结果（光密度值）对照组0.35470.36520.38420.35980.37420.3925高葡萄糖组1.22471.23541.21361.22861.24151.2275高胰岛素组1.14531.15261.13751.15471.12941.1659过氧化氢组1.35821.36571.35181.32591.37261.331

39、8问题1：计算每一组的平均光密度，观察各个组间均数的差异大小；问题2：进行方差分析，将结果整理成方差分析表，观察变异、均方及F值之间的关系；问题3：分别采用SNK法和Dunnett法进行多重比较；问题4：对上述结果进行解释。参考程序：data a;do grp=1 to 4; do i=1 to 6; input x; output; end; end;cards;0.35470.36520.38420.35980.37420.39251.22471.23541.21361.22861.24151.22751.14531.15261.13751.15471.12941.16591.35821.

40、36571.35181.32591.37261.3318;或者data a;length grp $10;input grp $ ; do i=1 to 6; input x; output; end;cards;对照组 0.35470.36520.38420.35980.37420.3925高葡萄糖组1.22471.23541.21361.22861.24151.2275高胰岛素组1.14531.15261.13751.15471.12941.1659过氧化氢组1.35821.36571.35181.32591.37261.3318;proc glm;class grp;model x=gr

41、p;means grp/snk dunnett;means grp;run;quit;2、（方差分析）研究枸杞多糖对脂肪肝的预防作用，按窝别作为区组标志，每一区组3只大鼠，随机分配到三个组：生理盐水组、酒精组、酒精+枸杞多糖组，三组灌胃5周，检测肝脏中谷胱甘肽（GSH）的含量（mg/gprot），结果见表2。表2 三组小鼠肝脏中谷胱甘肽（GSH）的含量（mg/gprot）区组酒精组LBP预防组生理盐水130.4865.0879.15231.2563.0475.46333.2867.5979.32434.6168.5875.98528.3564.1276.55629.1766.5580.3472

42、7.3466.8984.35830.5867.1588.14934.2568.0587.351027.3165.4872.151128.0964.3874.611230.4565.0486.331333.2566.8494.351434.0467.5692.051534.2567.4696.42问题1：按照随机区组设计进行分析，处理组间、区组间效应是否有差异？问题2：忽略区组，按照完全随机设计进行分析，处理组间效应是否有差异？问题3：对上述两个方差分析表进行比较，找出相同项与不同项？观察处理组间比较的F值变化，会得出怎样的结论？问题4：根据分析结果对该研究下一个结论。参考程序：data a;i

43、nput blk ;do grp='酒精组 ','LBP预防组','生理盐水' input x; output;end;cards;130.4865.0879.15231.2563.0475.46333.2867.5979.32434.6168.5875.98528.3564.1276.55629.1766.5580.34727.3466.8984.35830.5867.1588.14934.2568.0587.351027.3165.4872.151128.0964.3874.611230.4565.0486.331333.2566.8494.

44、351434.0467.5692.051534.2567.4696.42;proc glm;class grp blk;model x= grp blk;means grp/snk;run;quit;proc glm;class grp;model x= grp;means grp/snk;run;quit;3、（方差分析）研究人员研究某种物质的毒性，将40只小鼠分为两组，雄雌各半，试验组给予研究物质，2小时后测定血液中碱性磷酸酶的含量，结果如下。表3 40只小鼠给予不同物质后血液中碱性磷酸酶的含量性别分组碱性磷酸酶雄性对照组367.9408.6375.6354.9421.7374.5432.

45、7401.3399.4367.5试验组423.8446.9432.5478.1437.5421.6489.0432.5421.0420.4雌性对照组378.1345.2390.6399.0421.1341.3322.5365.4321.6401.9试验组420.4473.2450.3405.5427.4460.5420.1394.4389.6420.5问题1：研究者考虑了几个影响因素？问题2：性别对碱性磷酸酶有无影响？问题3：该物质对碱性磷酸酶有无影响？参考程序：data a;do sex='M','F' do grp='C','T

46、9; do i= 1 to 10; input x; output; end; end;end;cards;367.9408.6375.6354.9421.7374.5432.7401.3399.4367.5423.8446.9432.5478.1437.5421.6489.0432.5421.0420.4378.1345.2390.6399.0421.1341.3322.5365.4321.6401.9420.4473.2450.3405.5427.4460.5420.1394.4389.6420.5;proc glm;class sex grp;model x= sex grp;run;quit;4、（回归分析）以下是制作标准曲线时对应的浓度和吸光度值，试求出标准曲线。并求出A=1.15时的浓度。表4 不同浓度对应的吸光度值浓度A值0.970.31.270.561.570.931.