2022年完整word版,SAS上机练习题 .pdf

重庆医科大学 -卫生统计学统计软件包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 ; 083 68 71 65 084 74 61 68 085 73 75 46 086 79 80 79 087 75 71 68 084 85 85 87 085 78 79 75 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 32 页 - - - - - - - - - 086 80 76 79 087 85 80 82 088 77 71 75 089 67 73 71 080 75 81 70 118 70 54 75 083 70 66 84 084 62 73 65 099 82 70 79 ; run; data b; set a; if sh60 or wl 60 or bl60 then output ; run; proc printdata=b; var id sh wl bl; run;5、根据下列数据建立数据集表 2销售数据开始时间终止时间费用2005/04/28 25MAY2009 $123,345,000 2005 09 18 05OCT2009 $33,234,500 2007/08/12 22SEP2009 $345,600 20040508 30JUN2009 $432,334,500 提示： (格式化输入；数据之间以空格分隔，数据对齐；注意格式后面的长度应以前一个位置结束开始计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式) 开始时间，输入格式yymmdd10. 终止时间，输入格式date10. 费用，输入格式dollar12. 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 32 页 - - - - - - - - - 参考程序：data a; input x1 yymmdd10. x2 date10. x3 dollar13.; cards ; 2005/04/28 25MAY2009 $123,345,000 2005 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 手机号码手机号码1508320464X 1313654836X 1510291126X 1351425709X 1339644156X 1303612231X 1532464828X 1383824855X 1365721647X 1522322763X 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 32 页 - - - - - - - - - 1523105754X 1357851051X 1592624347X 1508311759X 1331237668X 1327313520X 1370048578X 1556443719X 1507244457X 1804346016X 1321246707X 1513441713X 1368464734X 1308279203X 1394913242X 1897876466X 参考程序：data a; input y 1-3 x $ 4-7 z $ 8- 11; if134=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 ; 1508320464X 1510291126X 1339644156X 。1513441713X 1308279203X 1897876466X ; proc print ; run;7、下表是某班学生几门功课的成绩，其中形势课是考查课，其它均为考试课。请完成以下处理并且保存程序、结果和日志： （注意函数在DATA STEP 中使用的位置）（1）用 Mean()函数求出每位同学的平均分，将其保存在变量中；（2）用 sum()函数或者表达式求出每位同学的总分，将其保存在变量中；（3）如果每门考试课权重是0.75，每门考查课权重为0.25，试求每位同学的加权平均分；（4）将数据集保存为永久数据集 ，保存位置为“” ；（5）调用上面的永久数据集，并将男、女生的成绩分别保存在数据集Dmale 和 Dfemale 中。表 3某班同学几门功课的成绩学号性别(0=女， 1=男) 高数生理人解数理统计形势 (考查 ) 0083 1 73 73 64 74 75 0414 1 90 79 71 85 78 0564 1 97 87 89 91 80 0774 1 40 60 61 65 75 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 32 页 - - - - - - - - - 0873 1 68 65 60 76 75 0874 1 74 68 56 60 75 0875 1 73 46 65 66 87 0876 1 79 79 74 89 87 0877 1 75 68 55 60 70 0878 1 76 60 64 71 80 0880 1 81 73 60 75 82 0881 1 76 81 84 89 78 0883 1 77 60 60 66 75 0884 0 85 87 84 85 90 0885 0 78 75 62 66 80 0886 0 80 79 77 73 83 0887 0 85 82 81 84 80 0888 0 77 75 74 67 80 0889 0 67 71 71 69 90 0890 0 75 70 75 69 80 0891 0 69 76 76 79 90 0892 0 66 71 60 60 78 0893 1 70 84 75 80 85 0894 1 62 65 65 69 81 0895 1 91 78 61 66 78 0896 1 51 60 42 44 78 0898 1 67 73 67 73 85 0899 1 79 72 78 68 80 0900 1 76 68 63 84 90 0901 0 81 70 70 82 79 参考程序：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 ; 0083 1 73 73 64 74 75 0414 1 90 79 71 85 78 0564 1 97 87 89 91 80 0774 1 40 60 61 65 75 0873 1 68 65 60 76 75 0874 1 74 68 56 60 75 0875 1 73 46 65 66 87 0876 1 79 79 74 89 87 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 32 页 - - - - - - - - - 0877 1 75 68 55 60 70 0878 1 76 60 64 71 80 0880 1 81 73 60 75 82 0881 1 76 81 84 89 78 0883 1 77 60 60 66 75 0884 0 85 87 84 85 90 0885 0 78 75 62 66 80 0886 0 80 79 77 73 83 0887 0 85 82 81 84 80 0888 0 77 75 74 67 80 0889 0 67 71 71 69 90 0890 0 75 70 75 69 80 0891 0 69 76 76 79 90 0892 0 66 71 60 60 78 0893 1 70 84 75 80 85 0894 1 62 65 65 69 81 0895 1 91 78 61 66 78 0896 1 51 60 42 44 78 0898 1 67 73 67 73 85 0899 1 79 72 78 68 80 0900 1 76 68 63 84 90 0901 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 编码及对应疾病名ICD10 Disease A01.001 伤寒A01.002 伤寒杆菌性败血症A01.003+ 伤寒性脑膜炎A01.101 甲型副伤寒A01.201 乙型副伤寒A01.301 丙型副伤寒名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 32 页 - - - - - - - - - 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 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 沙门氏菌性肠炎名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 32 页 - - - - - - - - - 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= A01then output cat01; if cat= A02then output cat02; if cat= A03then output cat03; run; 9、现有两个文件，内容如下：A 文件学号统计体育年龄01 87 83 23 02 56 96 22 03 93 75 19 04 77 84 24 05 88 55 18 B 文件学号性别班级01 男1 02 女1 03 男3 04 男2 05 女2 08 女3 要求：打印至少有一门功课不及格的同学的年龄、性别和班级。参考程序：data fileA; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 32 页 - - - - - - - - - input no$ tj ty age; cards ; 01 87 83 23 02 56 96 22 03 93 75 19 04 77 84 24 05 88 55 18 ; data fileB; input no$ sex$ bj; cards ; 01 男 1 02 女 1 03 男 3 04 男 2 05 女 2 08 女 3 ; proc sortdata =fileA; by no; proc sortdata =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; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 32 页 - - - - - - - - - if flag= 0 then output ; end; keep i; run; proc print ; run;(彭斌， 2010-3-23) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 32 页 - - - - - - - - - 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.6 137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 150.6 142.7 151.3 142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 145.7 140.8 139.6 148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 149.8 151.1 140.1 140.5 143.4 152.9 147.5 147.7 162.6 141.6 143.6 144.0 150.6 138.9 150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 139.5 148.9 144.5 141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 142.4 137.1 141.4 144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 147.1 140.9 137.4 142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 141.4 134.7 138.5 146.6 134.5 135.1 141.9 142.1 138.1 134.9 146.7 138.5 139.6 139.2 148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 143.5 139.9 134.4 133.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 142.7 144.4 144.7 145.1 143.3 154.2 152.3 142.7 156.6 137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 150.6 142.7 151.3 142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 145.7 140.8 139.6 148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 149.8 151.1 140.1 140.5 143.4 152.9 147.5 147.7 162.6 141.6 143.6 144.0 150.6 138.9 150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 139.5 148.9 144.5 141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 142.4 137.1 141.4 144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 147.1 140.9 137.4 142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 141.4 134.7 138.5 146.6 134.5 135.1 141.9 142.1 138.1 134.9 146.7 138.5 139.6 139.2 148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 143.5 139.9 134.4 133.1 145.9 139.2 137.4 142.3 160.9 137.7 142.9 126.8 ; proc meansdata=a n meanmedianclm qrangestd cv stderr; var x; proc univariate data=a normal; histogram x; var x; run ;名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 32 页 - - - - - - - - - 2、某研究员欲分析急性染毒对肝脏功能的影响，将40 只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后2 小时测定血液中的ALT，整理的结果见下表。（1）比较实验组与对照组的ALT有无差异；ALT( 丙氨酸转氨酶) NO male control Treaetment Female control Treaetment 1 38.4 57.8 39.7 55.7 2 30.6 62.4 34.6 63.8 3 25.5 57.5 40.5 59.9 4 29.9 49.9 44.3 49.7 5 29.0 55.6 41.2 48.8 6 33.6 47.3 25.6 51.3 7 39.4 50.1 30.6 53.8 8 42.5 48.7 37.6 42.6 9 46.7 57.9 39.5 61.7 10 49.3 60.4 25.1 65.3 参考程序：data a2; do sex=M,F; do grp=C,T; input alt; output; end; end; cards; 38.4 57.8 39.7 55.7 30.6 62.4 34.6 63.8 25.5 57.5 40.5 59.9 29.9 49.9 44.3 49.7 29.0 55.6 41.2 48.8 33.6 47.3 25.6 51.3 39.4 50.1 30.6 53.8 42.5 48.7 37.6 42.6 46.7 57.9 39.5 61.7 49.3 60.4 25.1 65.3 ; proc anova data=a2; class sex grp; model alt=sex grp; means grp; run ; quit ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 32 页 - - - - - - - - - 3、24 名志愿者随机分成两组，每组12 人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L ） ，数据见下表：甲组乙组受试者试验前试验后受试者试验前试验后1 6.11 6.00 1 6.90 6.93 2 6.81 6.83 2 6.40 6.35 3 6.48 4.49 3 6.48 6.41 4 7.59 7.28 4 7.00 7.10 5 6.42 6.30 5 6.53 6.41 6 6.94 6.64 6 6.70 6.68 7 9.17 8.42 7 9.10 9.05 8 7.33 7.00 8 7.31 6.83 9 6.94 6.58 9 6.96 6.91 10 7.67 7.22 10 6.81 6.73 11 8.15 6.57 11 8.16 7.65 12 6.60 6.17 12 6.98 6.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.93 2 6.81 6.83 2 6.40 6.35 3 6.48 4.49 3 6.48 6.41 4 7.59 7.28 4 7.00 7.10 5 6.42 6.30 5 6.53 6.41 6 6.94 6.64 6 6.70 6.68 7 9.17 8.42 7 9.10 9.05 8 7.33 7.00 8 7.31 6.83 9 6.94 6.58 9 6.96 6.91 10 7.67 7.22 10 6.81 6.73 11 8.15 6.57 11 8.16 7.65 12 6.60 6.17 12 6.98 6.52 ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 32 页 - - - - - - - - - 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.7 II 期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_=11then grp=0期; elseif _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.7 11.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 ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 32 页 - - - - - - - - - 教材上的练习题(卫生统计学，方积乾主编，第六版) 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 ; procttestdata=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 ; procttestdata=a1 ; var d; run ; 第三题data a3; input x; if _n_=10then 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.9 8.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 ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 32 页 - - - - - - - - - 第四题data a4; input x1 x2; cards; 16.84 16.79 19.02 19.22 10.44 10.40 14.87 15.14 22.31 21.89 24.83 24.82 26.89 27.00 31.06 31.42 36.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， 22 表，卡方检验data a1; do row=1 to 2; do col=1 to 2; input f; output; end; end; cards; 64 21 51 33 ; run ; proc freq data=a1 ; tables row*col/ chisq; weight f; run ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 32 页 - - - - - - - - - 例 9.3，22 表，校正卡方或Fisher 精确检验data a2; input row col f; cards; 1 1 2 1 2 10 2 1 14 2 2 14 ; run ; proc freq data=a2 ; tables row*col/ chisq; weight f; run ; 例9.4，32表，卡方检验data a3; do row=1 to 3; do col=1 to 2; input f; output; end; end; cards; 35 5 20 10 7 25 ; run ; proc freq data=a3 ; tables row*col/ chisq; weight f; run ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 32 页 - - - - - - - - - 秩和检验p179 例 10.2例 10.6 p193 7 题例 10.2，配对秩和检验data a1; input x1 x2; d =x1-x2; cards; 10.5 8.8 21.6 18.8 14.9 13.5 30.2 27.6 8.4 9.1 7.7 7.0 16.4 14.7 19.5 17.2 127.0 155.0 18.7 16.3 9.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 20 26 16 22 72 24 14 186 32 22 ; proc npar1way data=a2 wilcoxon; class col; var row; freq f; run ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 32 页 - - - - - - - - - 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=2; if x=无效 then x1=1; output; end; end; cards; 治愈有效有效无效治愈有效治愈治愈有效有效治愈有效治愈无效治愈有效有效无效治愈治愈; proc freq data=a3; table id*grp*x1/ cmhscores =rank noprint; run ; (彭斌， 2010-4-18) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 32 页 - - - - - - - - - SAS 上机练习题（三）1、 （方差分析） 某研究人员研究3 种刺激因素（高葡萄糖、高胰岛素和过氧化氢）刺激大鼠肾小球系膜细胞对单核细胞趋化蛋白表达的影响，在四种条件下对大鼠肾小球系膜细胞进行培养，检测结果如表1。表 1 不同刺激因素作用下蛋白表达结果分组表达结果（光密度值）对照组0.3547 0.3652 0.3842 0.3598 0.3742 0.3925 高葡萄糖组1.2247 1.2354 1.2136 1.2286 1.2415 1.2275 高胰岛素组1.1453 1.1526 1.1375 1.1547 1.1294 1.1659 过氧化氢组1.3582 1.3657 1.3518 1.3259 1.3726 1.3318 问题 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.3547 0.3652 0.3842 0.3598 0.3742 0.3925 1.2247 1.2354 1.2136 1.2286 1.2415 1.2275 1.14