重庆医科大学--卫生统计学统计软件包
SAS上机练习题（一）
1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。
2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中内容，注意不同颜色的含义；并根据日志窗中的信息修改完善程序。
3、将第2题的程序、结果及日志保存到磁盘。
4、试根据某班12名学生3门功课成绩表完成后面的问题：
问题：
1）建立数据集；
2）打印至少有1门功课不及格同学的信息；(提示，使用if语句)
5、（选做）根据下列数据建立数据集
表2销售数据
开始时间终止时间费用
2005/04/28 25MAY2009 $123,345,000
2005 09 18 05OCT2009 $33,234,500
提示：(
计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式)
开始时间，输入格式yymmdd10.
终止时间，输入格式date10.
费用，输入格式dollar12.
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
7、下表是某班学生几门功课的成绩，其中形势课是考查课，其它均为考试课。请完成以下处理并且保存程序、结果和日志：（注意函数在DATA STEP中使用的位置）
（1）用Mean()函数求出每位同学的平均分，将其保存在变量中；
（2）用sum()函数或者表达式求出每位同学的总分，将其保存在变量中；
（3）如果每门考试课权重是0.75，每门考查课权重为0.25，试求每位同学的加权平均分；
（4）将数据集保存为永久数据集，保存位置为“D:\mydata\”；
表3某班同学几门功课的成绩
8、下面是3个大类疾病的ICD-10编码及对应的疾病名。请完成以下任务：
（1）建立数据集；
（2）提取每种疾病的大类编码；
（3）分别将3个大类的疾病存入3个数据集。
（提示：ICD10编码中小数点前面的三位表示大类；length语句定义字符变量长度；字符串取子串函数）
表44类疾病的ICD10编码及对应疾病名
ICD10 Disease
A01.001 伤寒
A01.002 伤寒杆菌性败血症
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.901 沙门氏菌感染
A02.902 鼠伤寒沙门氏菌感染
A02.903 沙门氏菌属食物中毒
A03.001 什密氏志贺菌痢疾
A03.101 弗氏志贺菌痢疾
A03.201 鲍氏志贺菌痢疾
A03.301 宋内氏志贺菌痢疾
A03.802 菌痢混合感染
A03.901 细菌性痢疾
A03.902 慢性细菌性痢疾急性发作
A03.903 慢性细菌性痢疾
A03.906 中毒性痢疾
9、现有两个文件，内容如下：
要求：(1)根据A、B
(2)打印至少有一门功课不及格的同学的年龄、性别和班级。
10、某研究员欲分析急性染毒对肝脏功能的影响，将40只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后2小时测定血液中的ALT，整理的结果见下表。
要求：根据上面的数据建立恰当的SAS数据集
11、某职业病防治所对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
要求：根据上面的数据建立恰当的SAS数据集
12、24名志愿者随机分成两组，每组12人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L），数据见下表：
甲组乙组
受试者试验前试验后受试者试验前试验后
1 6.11 6.00 1 6.90 6.93
2 6.81 6.8
3 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
(彭斌，2012-11-6)
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）观察频数分布情况；
2、（t检验、配对t检验）24名志愿者随机分成两组，每组12人，接受降胆固醇试验，甲组为特殊饮食组，乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇（mmol/L），数据见下表：
甲组乙组受试者试验前试验后受试者试验前试验后
1 6.11 6.00 1 6.90 6.93
2 6.81 6.8
3 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）试判断两种降胆固醇措施的效果是否相同？
3、（秩和检验）某职业病防治所对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
请用秩和检验完成分析。
请进行以下分析：
1） 此次调查中男性的比例为多少；
2） 对46级改革之后学习英语的积极性的情况进行描述，三种选择占的比例分别是多少？ 3） 男生与女性的对46级改革的反应有无差异； 4） 计算身体质量指数（Body Mass Index ，)
()
(2m Kg BMI 身高体重
），比较男女生的BMI 有无差
异？
5） 关于BMI
某年级的问卷调查数据
机器号 性别 身高 体重 开支 英语改革 TB 症状 id sex height weight money question1
question2 cnw1l01 1 179 70 5.7 a ABDEF cnw1l02 1 175 70 7.5 a ABE cnw1l03 2 157 47 4.5 a ABE cnw1l04 2 163 48 5 c ABDF cnw1l05
2
161
52
5
b
ABF
cnw1r02 2 158 56 6.2 b ABDE cnw1r03 1 168 64 6.5 a ABE cnw1r05 2 150 44 5.5 b ABD cnw1r06 1 162 53 7 b ABDE cnw2l02 2 160 53 8 c ABEF cnw2l03 2 164 58 6 c ABE cnw2l04 1 167 60 7 c ABCDEF cnw2l05 2 164 54 5 c AB cnw2r02 2 151 44 6 b B cnw2r03 2 150 51 7.8 b ABDF cnw2r04 2 160 55 5 b ABE cnw2r05 2 160 47 3.9 b ABEF cnw2r06 2 156 48 8.5 c ABCDEFG cnw3l01 2 162 48 8 c ABDF cnw3l02 1 163 60 7.2 a ABE cnw3l03 1 161 50 7.3 b ADFG cnw3r01 1 169 57 5.3 a ABDF cnw3r02 2 162.5 46.5 7 a ABCDEF cnw3r03 2 155 53 7.5 b ABCDEF cnw3r05 2 163 50 4.4 b ABED cnw3r06 1 163 65 8 b ABCDEF cnw4l01 1 168 60 5 c ABDG cnw4l02 1 178 68 8 a AC cnw4l03 2 158 45 6 c ADEF cnw4l04 1 173 68 8.5 b ABDEF cnw4r01 1 175 53 5 b ABDEF cnw4r02 1 178 69 6 b A cnw4r03 2 158 50 8 b ABDEF cnw4r04 2 165 53 4.5 b ABDE cnw4r05 2 157 46 5.3 b ABE cnw4r06 2 155 48 6.7 c ABD cnw5l01 2 160 49 4.2 b ABDEF cnw5l02 2 160 49 7 a ABE cnw5l03 2 170 57 8.5 c AB cnw5r01 2 160 49 4.2 b ABDE cnw5r02 1 168 65 9.2 a ABDEF cnw5r02 2 170 57 8.5 c AD cnw6l02 2 160 50 4 b ABE cnw6l03 1 181 63 3.5 b ADEF cnw6l04 1 175 53 20 a AB
cnw6r02 2 166 54 4.5 c ABDF cnw6r03 2 155 43 5.4 c ABEF cnw6r04 2 153 45 5.4 c EG cnw6r05 2 160 48 5 a AB cnw6r06 2 164 46 7.4 b ABEF cnw7l02 2 160 46 8 b A cnw7l05 1 174 70 7 b ADF cnw7r01 2 155 55 8 c AB cnw7r02 1 168 60 10 c ABEF cnw7r03 1 181 60 10 c ABE cnw7r05 1 166 60 7 c ABDE cnw7r06 2 164 46 7.4 b ABEF cnw8l01 1 160 50 9 c ADF cnw8r01 1 174 70 8 a ADEG cnw8r02 2 172 53 5 c ABE cnw1r01 2 154 49 4.2 C ABDEF cnw1r02 1 169 64 10.4 B ABDE cnw1r04 1 169 67 15 A ABDEFG cnw1r05 2 153 50 10 B ABCD cnw1r06 2 160 50 3.6 A ABE cnw2r02 2 156 46 4 A ADG cnw2r03 1 173 62 7.5 B ABFG cnw2r04 1 168 57 8.9 B ABDF cnw2r05 2 155 45 6 B ABD cnw2r06 1 170 60 7 C ABD cnw3r02 1 162 56 6.3 C ABD cnw3r03 2 161 50 7 B ABE cnw3r04 2 154 48 5 C ABE cnw3r05 2 155 46 6 B ABDF cnw3r06 2 155 51 5 A ABDEG cnw4r02 2 153 50 3.5 C BE cnw4r03 2 150 48 5 A ABDE cnw4r04 2 152.5 46 5.2 A ADE cnw4r05 1 185 70 7.5 A ABDG cnw4r06 1 170 58 6.5 B AD cnw5r01 1 173 62 6 B ABDE cnw5r02 1 168 50 6.1 A ABDE cnw5r03 2 160 50 7 C ABCEF cnw5r04 2 155 42 6 B ABDF cnw5r05 1 168 55 7 B ADF
cnw6r01 2 158 50 8.1 C ADF cnw6r02 1 175 61 10.5 B ABDE cnw6r03 2 155 44 11 C ABDE cnw6r04 1 162 60 10 C ADF cnw6r05 1 173 64 15 C ADF cnw7r01 1 175 86 11 C ABE cnw7r02 1 172 62 12 B ABEG cnw7r03 1 172 65 6.3 B ABEG cnw7r05 2 151 48 7 B ABDEFG cnw8r02 1 167 50 9 B ABDEF cnw8r03 2 160 50 5 B ABDEFG cnw8r04 2 170 65 10 B ABDEFG cnw1l01 1 165 47 7.3 B AB cnw1l02 2 152 44 4.2 B AB cnw1l03 2 162 49.5 4.5 B ABE cnw1l04 2 157 52 3.2 B ABE cnw1l05 2 150 45 3.8 A ABE cnw2l02 2 170 55 4.9 B ADF cnw2l03 2 155 47 5.9 B ABDF cnw2l04 2 163 48 5.8 A ADFG cnw2l05 2 150 47.5 8 A ABE cnw3l01 2 153 46 4 B ABDFG cnw3l02 1 165 54 35 A ADF cnw4l01 2 160 52 12.4 A ABDF cnw4l02 1 170 65 5 A D cnw4l04 1 165 51 5 B AF cnw4l05 2 155 42 8 B ABCDE cnw5l01 2 160 45 7 C ABDF cnw5l02 1 179 63 10 C ABDF cnw5l04 1 170 67 9 C ADF cnw6l01 1 170 58 9.6 B ABCDEF cnw6l02 2 160 54 5.5 C ABCDEF cnw6l03 2 155 48 4.4 A ABCDEF cnw6l04 1 170 65 10 C ABCDEF cnw6l05 2 165 47 7 C ABDF cnw7l01 2 150 48 5.8 C ABE cnw7l02 1 160 50 7 A ABCDEF cnw7l03 1 175 66 10 A ABCDEF cnw7l04 2 163 51 8 B ABD cnw7l05 2 165 57 4.9 A ABD cnw8l01 1 160 60 10 B ABE cnw8l02 2 154 50 4.3 A BE cnw8l03 2 160 60 7 A ABDEF
教材上的练习题
(医学统计学，颜虹主编，人民卫生出版社，第一版)
请根据第二版教材相应题目练习。
P132：（t检验）第1，2，4题
P149：（方差分析）第1，3，4题
P168：（卡方检验）第1，3题
P184：（秩和检验）第2，3，6题
(彭斌，2012-11-6)
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：对上述结果进行解释。
2、（方差分析）研究枸杞多糖对脂肪肝的预防作用，按窝别作为区组标志，每一区组3只大鼠，随机分配到三个组：生理盐水组、酒精组、酒精+枸杞多糖组，三组灌胃5周，检测肝脏中谷胱甘肽（GSH）的含量（mg/gprot），结果见表2。
表2 三组小鼠肝脏中谷胱甘肽（GSH）的含量（mg/gprot）
区组酒精组LBP预防组生理盐水
130.4865.0879.15
231.2563.0475.46
333.2867.5979.32
434.6168.5875.98
528.3564.1276.55
629.1766.5580.34
727.3466.8984.35
830.5867.1588.14
934.2568.0587.35
1027.3165.4872.15
1128.0964.3874.61
1230.4565.0486.33
1333.2566.8494.35
1434.0467.5692.05
1534.2567.4696.42
问题1：按照随机区组设计进行分析，处理组间、区组间效应是否有差异？
问题2：忽略区组，按照完全随机设计进行分析，处理组间效应是否有差异？
问题3：对上述两个方差分析表进行比较，找出相同项与不同项？观察处理组间比较的F值变化，会得出怎样的结论？
问题4：根据分析结果对该研究下一个结论。
3、（方差分析）研究人员研究某种物质的毒性，将40只小鼠分为两组，雄雌各半，试验组给予研究物质，2小时后测定血液中碱性磷酸酶的含量，结果如下。
表3 40只小鼠给予不同物质后血液中碱性磷酸酶的含量
性别分组碱性磷酸酶
雄性对照组367.9 408.6 375.6 354.9 421.7 374.5 432.7 401.3 399.4 367.5 试验组423.8 446.9 432.5 478.1 437.5 421.6 489.0 432.5 421.0 420.4 雌性对照组378.1 345.2 390.6 399.0 421.1 341.3 322.5 365.4 321.6 401.9 试验组420.4 473.2 450.3 405.5 427.4 460.5 420.1 394.4 389.6 420.5
问题1：研究者考虑了几个影响因素？
问题2：性别对碱性磷酸酶有无影响？
问题3：该物质对碱性磷酸酶有无影响？
4、（回归分析）以下是制作标准曲线时对应的浓度和吸光度值，试求出标准曲线。并求出A=1.15时的浓度。
表4 不同浓度对应的吸光度值
浓度A值
0.97 0.3
1.27 0.56
1.57 0.93
1.88 1.35
2.18 1.51
5、(相关分析、回归分析)27名糖尿病人的血清总胆固醇、甘油三脂、空腹胰岛素、糖化血红蛋白、空腹血糖的测量值列于表中，试分析血糖与其它几项指标关系的关系。
表5 27名糖尿病人的血糖及有关变量的测量结果
序号i
总胆固醇甘油三脂胰岛素糖化血血糖（mmol/L）(mmol/L) ( U/ml) 红蛋白(%) (mmol/L) X1 X2 X3 X4 Y
1 5.68 1.90 4.53 8.
2 11.2
2 3.79 1.64 7.32 6.9 8.8
3 6.02 3.56 6.95 10.8 12.3
4 4.8
5 1.07 5.88 8.3 11.6
5 4.60 2.32 4.05 7.5 13.4
6 6.05 0.64 1.42 13.6 18.3
7 4.90 8.50 12.60 8.5 11.1
8 7.08 3.00 6.75 11.5 12.1
9 3.85 2.11 16.28 7.9 9.6
10 4.65 0.63 6.59 7.1 8.4
11 4.59 1.97 3.61 8.7 9.3
12 4.29 1.97 6.61 7.8 10.6
13 7.97 1.93 7.57 9.9 8.4
14 6.19 1.18 1.42 6.9 9.6
15 6.13 2.06 10.35 10.5 10.9
16 5.71 1.78 8.53 8.0 10.1
17 6.40 2.40 4.53 10.3 14.8
18 6.06 3.67 12.79 7.1 9.1
19 5.09 1.03 2.53 8.9 10.8
20 6.13 1.71 5.28 9.9 10.2
21 5.78 3.36 2.96 8.0 13.6
22 5.43 1.13 4.31 11.3 14.9
23 6.50 6.21 3.47 12.3 16.0
24 7.98 7.92 3.37 9.8 13.2
25 11.54 10.89 1.20 10.5 20.0
26 5.84 0.92 8.61 6.4 13.3
27 3.84 1.20 6.45 9.6 10.4
问题1：计算各指标之间的pearson相关系数及spearman相关系数；
问题2：以血糖为应变量Y其它为自变量X，进行回归分析。
请进行以下分析：
1）试对男女生的开支进行分析；
2）析被调查者对“结核病的早期症状”的选择情况，即列出每种选项的频数及比例；（提示：可分别对每个选项单独进行分析，使用index()函数判断是否选择某个选项，具体用法请查帮助。）
3）分析身高与体重的相关性；建立身高预测体重的预测模型。
机器号性别身高体重开支英语改革TB症状
id sex height weight money question1 question2
cnw1l01 1 179 70 5.7 a ABDEF
cnw1l02 1 175 70 7.5 a ABE
cnw1l03 2 157 47 4.5 a ABE
cnw1l04 2 163 48 5 c ABDF
cnw1l05 2 161 52 5 b ABF
cnw1l06 2 150 44 5.5 b ABD
cnw1r01 2 158 56 4.4 b ADFG
cnw1r06 1 162 53 7 b ABDE cnw2l02 2 160 53 8 c ABEF cnw2l03 2 164 58 6 c ABE cnw2l04 1 167 60 7 c ABCDEF cnw2l05 2 164 54 5 c AB cnw2r02 2 151 44 6 b B cnw2r03 2 150 51 7.8 b ABDF cnw2r04 2 160 55 5 b ABE cnw2r05 2 160 47 3.9 b ABEF cnw2r06 2 156 48 8.5 c ABCDEFG cnw3l01 2 162 48 8 c ABDF cnw3l02 1 163 60 7.2 a ABE cnw3l03 1 161 50 7.3 b ADFG cnw3r01 1 169 57 5.3 a ABDF cnw3r02 2 162.5 46.5 7 a ABCDEF cnw3r03 2 155 53 7.5 b ABCDEF cnw3r05 2 163 50 4.4 b ABED cnw3r06 1 163 65 8 b ABCDEF cnw4l01 1 168 60 5 c ABDG cnw4l02 1 178 68 8 a AC cnw4l03 2 158 45 6 c ADEF cnw4l04 1 173 68 8.5 b ABDEF cnw4r01 1 175 53 5 b ABDEF cnw4r02 1 178 69 6 b A cnw4r03 2 158 50 8 b ABDEF cnw4r04 2 165 53 4.5 b ABDE cnw4r05 2 157 46 5.3 b ABE cnw4r06 2 155 48 6.7 c ABD cnw5l01 2 160 49 4.2 b ABDEF cnw5l02 2 160 49 7 a ABE cnw5l03 2 170 57 8.5 c AB cnw5r01 2 160 49 4.2 b ABDE cnw5r02 1 168 65 9.2 a ABDEF cnw5r02 2 170 57 8.5 c AD cnw6l02 2 160 50 4 b ABE cnw6l03 1 181 63 3.5 b ADEF cnw6l04 1 175 53 20 a AB cnw6l05 1 174 60 6.5 b AD cnw6l06 1 162 52 4.7 b ABE cnw6r01 2 154 54 10 a ABDF cnw6r02 2 166 54 4.5 c ABDF
cnw6r05 2 160 48 5 a AB cnw6r06 2 164 46 7.4 b ABEF cnw7l02 2 160 46 8 b A cnw7l05 1 174 70 7 b ADF cnw7r01 2 155 55 8 c AB cnw7r02 1 168 60 10 c ABEF cnw7r03 1 181 60 10 c ABE cnw7r05 1 166 60 7 c ABDE cnw7r06 2 164 46 7.4 b ABEF cnw8l01 1 160 50 9 c ADF cnw8r01 1 174 70 8 a ADEG cnw8r02 2 172 53 5 c ABE cnw1r01 2 154 49 4.2 C ABDEF cnw1r02 1 169 64 10.4 B ABDE cnw1r04 1 169 67 15 A ABDEFG cnw1r05 2 153 50 10 B ABCD cnw1r06 2 160 50 3.6 A ABE cnw2r02 2 156 46 4 A ADG cnw2r03 1 173 62 7.5 B ABFG cnw2r04 1 168 57 8.9 B ABDF cnw2r05 2 155 45 6 B ABD cnw2r06 1 170 60 7 C ABD cnw3r02 1 162 56 6.3 C ABD cnw3r03 2 161 50 7 B ABE cnw3r04 2 154 48 5 C ABE cnw3r05 2 155 46 6 B ABDF cnw3r06 2 155 51 5 A ABDEG cnw4r02 2 153 50 3.5 C BE cnw4r03 2 150 48 5 A ABDE cnw4r04 2 152.5 46 5.2 A ADE cnw4r05 1 185 70 7.5 A ABDG cnw4r06 1 170 58 6.5 B AD cnw5r01 1 173 62 6 B ABDE cnw5r02 1 168 50 6.1 A ABDE cnw5r03 2 160 50 7 C ABCEF cnw5r04 2 155 42 6 B ABDF cnw5r05 1 168 55 7 B ADF cnw5r06 1 173 61 8 B ABD cnw6r01 2 158 50 8.1 C ADF cnw6r02 1 175 61 10.5 B ABDE cnw6r03 2 155 44 11 C ABDE
cnw6r04 1 162 60 10 C ADF
cnw6r05 1 173 64 15 C ADF
cnw7r01 1 175 86 11 C ABE
cnw7r02 1 172 62 12 B ABEG
cnw7r03 1 172 65 6.3 B ABEG
cnw7r05 2 151 48 7 B ABDEFG
cnw8r02 1 167 50 9 B ABDEF
cnw8r03 2 160 50 5 B ABDEFG
cnw8r04 2 170 65 10 B ABDEFG
cnw1l01 1 165 47 7.3 B AB
cnw1l02 2 152 44 4.2 B AB
cnw1l03 2 162 49.5 4.5 B ABE
cnw1l04 2 157 52 3.2 B ABE
cnw1l05 2 150 45 3.8 A ABE
cnw2l02 2 170 55 4.9 B ADF
cnw2l03 2 155 47 5.9 B ABDF
cnw2l04 2 163 48 5.8 A ADFG
cnw2l05 2 150 47.5 8 A ABE
cnw3l01 2 153 46 4 B ABDFG
cnw3l02 1 165 54 35 A ADF
cnw4l01 2 160 52 12.4 A ABDF
cnw4l02 1 170 65 5 A D
cnw4l04 1 165 51 5 B AF
cnw4l05 2 155 42 8 B ABCDE
cnw5l01 2 160 45 7 C ABDF
cnw5l02 1 179 63 10 C ABDF
cnw5l04 1 170 67 9 C ADF
cnw6l01 1 170 58 9.6 B ABCDEF
cnw6l02 2 160 54 5.5 C ABCDEF
cnw6l03 2 155 48 4.4 A ABCDEF
cnw6l04 1 170 65 10 C ABCDEF
cnw6l05 2 165 47 7 C ABDF
cnw7l01 2 150 48 5.8 C ABE
cnw7l02 1 160 50 7 A ABCDEF
cnw7l03 1 175 66 10 A ABCDEF
cnw7l04 2 163 51 8 B ABD
cnw7l05 2 165 57 4.9 A ABD
cnw8l01 1 160 60 10 B ABE
cnw8l02 2 154 50 4.3 A BE
（彭斌2012-11-6 整理）
6$6? ?д乬??? ?1乬 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 means data=a n mean median clm qrange std cv stderr maxdec=2; var x; proc univariate data=a normal; histogram x; var x; run; ??? The MEANS Procedure Analysis Variable : x N Mean Median Lower 95% CL for Mean Upper 95% CL for Mean Quartile Range Std Dev Coeff of Variation Std Error ? Ё? 95% ??? 95% ??? ? ? ? ? ? 130 143.22142.75142.10144.337.80 6.43 4.490.56
重庆医科大学--卫生统计学统计软件包 SAS上机练习题（一） 1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中内容，注意不同颜色的含义；并根据日志窗中的信息修改完善程序。 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 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 sh<60 or wl<60 or bl<60then output; run; proc print data=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. 参考程序： 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
中科院研究生院统计分析与SAS实现第1次上机实习题 一、定量资料上机实习题 要求： （1）先判断定量资料所对应的实验设计类型； （2）假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算； （3）摘录主要计算结果并合理解释，给出统计学结论和专业结论。 【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较。 表1 未成年大白鼠注射不同剂量雌激素后的子宫重量 子宫重量(g) 窝别 剂量(μg/100g):0.2 0.4 0.8 合计 1 106 116 145 367 2 42 68 115 225 3 70 111 133 314 4 42 63 87 192 合计260 358 480 1098 定量资料的随机区组设计（区组因素：窝别；实验因素：剂量）【SAS程序】： DATA PGM15G; DO A=1TO4; /*A为窝别*/ DO B=1TO3; /*B为雌激素剂量*/ INPUT X @@; /*X为子宫重量*/ OUTPUT; END; END; CARDS; 106 116 145 42 68 115 70 111 133 42 63 87 RUN; ods html; /*将结果输出成网页格式，SAS9.0以后版本可用*/ PROC GLM DATA=PGM15G; CLASS A B; MODEL X=A B / SS3; MEANS A B; /*给出因素A、B各水平下的均值和标准差*/ MEANS B / SNK; /*对因素B（即剂量）各水平下的均值进行两两比较*/ RUN; ODS HTML CLOSE; 【练习2】一位工程师研究由钻头压力产生的冲力。考察了A（钻孔速度）和B（进料
数据统计分析一般可遵循以下思路: (1)先确定研究目的，根据研究目的选择方法。不同研究目的采用的统计方法不同，常见的研究目的主要有三类：①差异性研究，即比较组间均数、率等的差异，可用的方法有t检验、方差分析、χ2检验、非参数检验等。②相关性分析，即分析两个或多个变量之间的关系，可用的方法有相关分析。③影响性分析，即分析某一结局发生的影响因素，可用的方法有线性回归、logistic回归、Cox 回归等。 (2)明确数据类型，根据数据类型进一步确定方法：①定量资料可用的方法有t检验、方差分析、非参数检验、线性相关、线性回归等。②分类资料可用的方法有χ2检验、对数线性模型、logistic回归等。下图简要列出了不同研究目的、不同数据类型常用的统计分析方法。 (3)选定统计方法后，需要利用统计软件具体实现统计分析过程。SAS中，不同的统计方法对应不同的命令，只要方法选定，便可通过对应的命令辅之以相应的选项实现统计结果的输出。 (4)统计结果的输出并非数据分析的完成。一般统计软件都会输出很多结果，需要从中选择自己需要的部分，并做出统计学结论。但统计学结论不同于专业结论，最终还需要结合实际做出合理的专业结论。 第一部分：统计描述
1.定量资料的统计描述指标及SAS实现； （1）数据分布检验：PROC UNIVARIATE ①基本格式： ②语句格式示例： 1.PROC UNIVARIATE normal;/*normal选项表示进行正态性检验*/ 2.CLASS group;/*指定group为分组变量*/ 3.VAR weight;/*指定分析变量为weight*/ 4.RUN; ③结果：正态性检验(tests for normality)结果，常用的是Shapiro-Wilk 检验和Kolmogorov-Smirnov检验。当例数小于2000时，采用Shapiro-Wilk检验W值为标准；当例数大于2000时，SAS中不显示Shapiro-Wilk检验结果，采用Kolmogorov-Smirnov检验D值为判断标准。正态性检验的P≤0.05提示不服从正态分布，P＞0.05提示服从正态分布。 注：若服从正态分布，进行PROC MEANS过程步；若不服从则计算百分位数，转（3） （2）数据描述（符合正态分布的数据）：PROC MEANS ①基本格式： 关键字（可以无视）：不写任何关键字时默认输出n，mean，std，max，min； n：有效数据记录数（有效样本量） median：中位数 mean：均数 qrange：四分位数间距 std：标准差 var：方差 clm：95%可信区间 max、min：最大、最小值 ②语句格式示例： 1.PROC MEANS n mean std median qrange clm;/*关调用proc means过程， 要求输出的指标有例数、均值、标准差、中位数、四分位数间距、95% 可信区间*/ 2.CLASS group;/*指定group为分组变量*/ 3.VAR weight;/*指定分析变量为weight*/ 4.Run; ③结果以“均数±标准差”表示 （3）偏正态分布的统计描述： ①基本思想：计算中位数和百分位数，并且用“中位数（Q1～Q3）”表示 ②语句格式示例： 1.proc univariate data=aa; 2.var x; 3.output out=c pctlpre=P pctlpts=0 to 100 by 2.5;/*计算0到100
重庆医科大学“卫生统计学统计软件包 SAS上机练习题（一） 1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中内容，注意不同颜色的含义；并根据口志窗屮的信息修改完善程序。 DATS EXO 101; INPUTT NAME $ AGE SEX; CARDS; XIAOMIN 19 1 LIDONG 20 1 NANA 18 2 PROD PRONT DATS=EX1; RUN; PROC PRINT DATA=EX 1; VAR NAME AGE; RUN; 3、将第2题的程序、结果及tl志保存到磁盘。 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 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 sh<60 or wl<60 or bl<60 then output; run; proc print data=b; var id sh wl bl; run; 5、根据下列数据建立数据集 表2销售数据 提示：（格式化输入；数据之间以空格分隔，数据对齐；注意格式后面的长度应以前一个位置结朿开始计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式） 开始时间，输入格式yymmddlO. 终止时间，输入格式血tel0. 费用，输入格式dollar 12.
医科大学--卫生统计学统计软件包 SAS上机练习题（一） 1、SAS常用的窗口有哪三个？请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序，提交系统执行，并查看OUTPUT窗和LOG窗中容，注意不同颜色的含义；并根据日志窗中的信息修改完善程序。 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 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 sh<60 or wl<60 or bl<60then output; run; proc print data=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. 参考程序：
1、为比较两种方法对乳酸饮料中脂肪含量测定结果是否不同，某人随机抽取了10份乳酸饮料制品，分别用脂肪酸水解法和哥特里－罗紫法测定其结果如表3-3第(1)~(3)栏。问两法测定结果是否不同？ 表两种方法对乳酸饮料中脂肪含量的测定结果(%) 编号(1) 哥特里－罗紫法 (2) 脂肪酸水解法 (3) 1 0.840 0.580 2 0.591 0.509 3 0.67 4 0.500 4 0.632 0.316 5 0.687 0.337 6 0.978 0.517 7 0.750 0.454 8 0.730 0.512 9 1.200 0.997 10 0.870 0.506 2、为研究国产四类新药阿卡波糖胶囊的降血糖效果，某医院用20名II型糖尿病病人进行同期随机对照试验。试验者将这些病人随机等分到试验组(用阿卡波糖胶囊)和对照组(用拜唐苹胶囊)，分别测得试验开始前和8周后的空腹血糖，算得空腹血糖下降值见表3-4，能否认为该国产四类新药阿卡波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果不同？ 表3-4 试验组和对照组空腹血糖下降值(mmol/L) 试验组 (n1=10) -0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50 对照组 (n2=10) 3.70 6.50 5.00 5.20 0.80 0.20 0.60 3.40 6.60 -1.10 3. 两组肿瘤患者，单纯放疗组（A）13 例，口服平消胶囊+放疗组（B）12 例，接受放疗前后，血清Sil-2R水平（U/ml）如表所示。试评价平消胶囊对接受放疗患者血清Sil-2R水平的影响？ 两组肿瘤患者的血清Sil-2R水平
SAS课程上机练习一 下面的数据是一次对20岁以上人群进行心血管病随访研究的结果，数据中各变量的含义如下： V1—随访号V2—年龄V3—开始随访时的求诊机构V4/V5—开始随访时的收缩压/舒张压V6/V7—体重/身高V8—开始随访时的胆固醇V9—社会经济地位（1-高2-中3-底）V10—原有心血管病（0-其他心脏病1-冠心病2-冠心与高心3-高心4-高心及风心5-风心6-可疑心脏病7-高血压8-正常）V11—随访结束时的求诊机构V12/V13—随访结束时的收缩压/舒张压V14—随访结束时的胆固醇V15—随访结束时的体重V16—随访结束时的最后诊断（0-未诊断1-3心肌梗塞4-7心绞痛8-9其他）V17—死亡年份（63-68）0-未死。 一、数据集练习操作 1、用data步建立SAS永久数据集（下次练习可以调用此数据集） 2、将年龄（V2）按20—30—40—50—及60岁以上分组 3、计算体重指数={体重/身高2}*100 4、将随访结束时最后诊断的结果（V16）分为四组：(1)心肌梗塞(2) 心绞痛(3)其他心脏 病(4)未诊断 5、按体重指数的大小( 0.30)分为二类：1-超重0-未超重 6、将死亡年份（V17）按是否死亡分为二类（1为死亡，0为未死亡）。 7、将原有心血管病（V10）分为患有心血管病与不患有心血管病二类（1有0无） 8、将数据集按是否死亡分为二个数据集。 9、将社会经济地位（V9）设置为哑变量。 10、计算开始随访时的收缩压与/舒张压（V4/V5）之差。 二.。计算分析 1．年龄分布的特征。 2．开始时收缩压频数分布的特征。 3·超重病人开始时胆固醇（V8）的分布的特征。 3·比较不同体重的人开始时的收缩压。 4·比较不同社会经济地位与患心血管病的差别。 5·比较心肌梗塞的病人开始时的收缩压与随访结束时的收缩压的差别。 6·比较不同年龄、不同体重(指数)间开始时的胆固醇的差别。 7·试分析年龄与血压（V4）之间的关系。 9、比较不同社会经济地位开始随访是收缩压的差别。 三、结果存盘 1.要求将你编辑好的SAS程序保存。 2.要求将所有的计算结果(OUTPUT)保存.。 3.文件名的格式为：你的姓名+你的学号 例如：王永，学号为20918003，则你存盘的文件名为：王永208918003
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; V AR 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 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 sh<60 or wl<60 or bl<60then output; run; proc print data=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.
实验内容 1试用产生标准正态分布的随机数normal(seed)产生参数为10的卡方分布随机数100个。 2根据数据集： （1）创建一个仅包含地区、销售的产品类型、销售数量和销售额的数据集。 （2）分别创建一个仅包含产品类型a100和产品类型a200和SAS数据集。 （3）选择一个人口在50000以上的部分子集。 3 以下数据来自7位同学的高考语文、数学和英语成绩，试用编程的方法计算出平均成绩在75以上的男同学的人数。 实验步骤： 实验1： data a (drop=i) ; do i=1to100by1; z=normal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2+no rmal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2; output; end; 结果如下： 实验2： （1） data biao1;
input region$ product$ quantity price; cards; es a100 150 3750 so a100 410 10250 es a100 350 8750 so a100 710 17750 es a100 750 18750 so a100 760 19000 es a100 150 3000 so a100 410 8200 es a100 350 7000 so a100 710 14200 es a100 750 15000 so a100 760 152000 es a200 165 4125 so a200 425 10425 es a200 365 9125 ne a100 200 5000 we a100 180 4500 ne a100 600 15000 we a100 780 19500 ne a100 800 20000 we a100 880 22000 ne a100 200 4000 we a100 180 3600 ne a100 600 12000 we a100 780 15600 ne a100 800 16000 we a100 880 17600 ne a200 215 5375 we a200 195 4875 ne a200 615 15375 ; 实验结果：
09硕/博士研究生课程《医学统计软件应用》上机试卷 (春学期) 一·内容 有一个名为：Hypertension的SAS数据文件，文件中的数据来自于杭州市某社区居民家庭健康档案一份2000例的随机样本，数据中各变量的含义如下： V1 地区编号：①城镇②农村 V2性别①男②女 V3 年龄（岁） V4文化程度 1.文盲或半文盲2小学 3.初中 4.中专或高中 5.大专6大学或以上 V5职业 1.专业技术人员(医生、教师等) 2.行政管理、文职人员 3.轻体力劳动工人（服务行业、钳工、检验工等） 4.重体力劳动工人（搬运工、建筑工等） 5.在校学生 6. 家务 V6你的家人是否有高血压①父母双方均有②父一方有③母一方有④父母双方均没有⑤不详 V7你的家人是否有糖尿病①父母双方均有②父一方有③母一方有④父母双方均没有⑤不详 V8你的家人是否有其他疾病①肿瘤②冠心病③脑卒中④高血脂 V9. 吸烟①是②否③已戒 V10.饮酒①是②不饮③已戒酒 V11..饮茶①是②偶尔③不饮 V12.锻炼你是否经常参加体育运动（包括散步.太极拳.跑步游泳跳舞等运动持续半小时以上为1次） ①每天参加②每周3次以上③每周1-2次 ④每月1-2次⑤不参加 V13饮食习惯①以荤食为主②以素食为主③荤素各半 V14吃.腌、霉食品①经常吃(每周3次以上) ②偶尔吃③不吃 V15食物口味咸淡程度①偏咸②偏淡③一般 V16甜食①.经常吃(每周3次以上) ②偶尔吃③不吃 V17工作压力①大②一般③没压力 V18 5年内家庭重大事件(亲属病故.伤亡.离异)①有②无 V19 身高（厘米）V20—体重（公斤）V21—腰围（厘米）V22—臀围（厘米） V23是否服降压药①是②否 V24—收缩压（mmHg） V25—舒张压（mmHg）
S A S上机实习
中科院研究生院统计分析与SAS实现第1次上机实习题 一、定量资料上机实习题 要求： （1）先判断定量资料所对应的实验设计类型； （2）假定资料满足参数检验的前提条件，请选用相应设计的定量资料的方差分析，并用SAS软件实现统计计算； （3）摘录主要计算结果并合理解释，给出统计学结论和专业结论。 【练习1】取4窝不同种系未成年的大白鼠，每窝3只，随机分配到三个实验组中，分别注射不同剂量雌激素，经过一定时间后处死大白鼠测子宫重量，资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义？若剂量间差别有统计学意义，请作两两比较。 表1 未成年大白鼠注射不同剂量雌激素后的子宫重量 窝别 子宫重量(g) 剂量(μg/100g):0.2 0.4 0.8 合计 1 106 116 145 367 2 42 68 115 225 3 70 111 133 314 4 42 63 87 192 合计260 358 480 1098 定量资料的随机区组设计（区组因素：窝别；实验因素：剂量）【SAS程序】： DATA PGM15G; DO A=1TO4; /*A为窝别*/ DO B=1TO3; /*B为雌激素剂量*/ INPUT X @@; /*X为子宫重量*/ OUTPUT; END; END; CARDS; 106 116 145 42 68 115 70 111 133 42 63 87 ; RUN; ods html; /*将结果输出成网页格式，SAS9.0以后版本可用*/ PROC GLM DATA=PGM15G; CLASS A B; MODEL X=A B / SS3;
上机练习题 1.括弧匹配检验work1.pas 假设表达式中允许包含两种括号：圆括号和方括号，其嵌套的顺序随意，如（［］（））或［（［］［］）］等为正确的匹配，［（]）或（［］（）或( ( ) ) ）均为错误的匹配。 现在的问题是，要求检验一个给定表达式中的括弧是否正确匹配？ 输入一个只包含圆括号和方括号的字符串，判断字符串中的括号是否匹配，匹配就输出“OK”，不匹配就输出“Wrong”。 输入一个字符串： ［（［］［］）］ 输出： OK 【输入】 输入仅一行字符（字符个数小于255） 【输出】 匹配就输出“OK”，不匹配就输出“Wrong”。 【样例】 输入（work1.in） ［（］） 输出（work1.out） Wrong 2.自然数的拆分work2.pas 任何一个大于1的自然数总可以拆分成若干个自然数之和。 1:4=1+1+1+1 2:4=1+1+2 3:4=1+3 4:4=2+2 5:4=4 分析：设拆分出的数s1≤s2≤…≤sk。定义数组s为一个栈，用来存放因子。从1开始搜索因子，求和，若sum ≤n就将因子压栈；若sum ＝n，则输出解，出栈；若sum ＞n，则修改栈顶元素的值，即回溯。 【输入】 输入仅一行，包含一个自然数n。 【输出】 这个自然数拆分结果。拆分出的因子要求满足s1＋s2＋…＋sk＝n 且s1≤s2≤…≤sk。具体要求见样例。 【样例】 输入（work2.in） 4 输出（work2.out）
1:4=1+1+1+1 2:4=1+1+2 3:4=1+3 4:4=2+2 5:4=4 3.糟糕的头发 work3.pas 农夫约翰有N (1 <= N <= 80,000)头奶牛，自从奶牛们独自照看自己的头发开始，它们的发型变得越来越糟糕。约翰决定统计奶牛们所能看到其他奶牛的头顶的数目总和。 第i头奶牛加上她头发的高度为h[i] (1 <= h[i] <= 1,000,000,000)，奶牛们站成一排面向东方（见图1），所以，从第i头奶牛的发尖看过去，她能看见前面所有比她矮的奶牛头顶（这些可能的奶牛依次为i+1,i+2,……）。 例如 = = = = - = 奶牛面向右 --> = = = = - = = = = = = = = = 1 2 3 4 5 6 Cow#1 可以看到的奶牛编号为#2, 3, 4 Cow#2 什么也看不到 Cow#3 可以看到的奶牛编号为#4 Cow#4 什么也看不到 Cow#5 可以看到的奶牛编号为#6 Cow#6 什么也看不到 约翰用c[i] 表示每头奶牛能看到的奶牛数量，请你帮约翰统计所有c[i]（1 <= i <= N）的总和。上面例子的统计结果就是：3 + 0 + 1 + 0 + 1 + 0 = 5 时间限制：0.5 seconds 【输入】 输入的第一行为一个单独的整数n表示奶牛的数目，以下i＋1行中n个整数分别代表第i头羊的高度。 【输出】只需一行，c[1]~c[n]的和。 【样例】 输入（work3.in） 6 10 3 7 4 12 2 输出（work3.out） 5
SAS课程上机练习四 1、用两种化疗方法治疗淋巴系统肿瘤，其疗效如下表。检验两种化疗方法的缓解率 有否差别？ 两种化疗方法的缓解率比较 治疗方法缓解未缓解合计缓解率（％） 单纯化疗组 2 10 12 16．7 复合化疗组14 13 27 51．9 合计16 23 39 2．用两种方法检查110例肺结核患者，结果如下表，问用两种方法检查的结果有无差别？ 110例肺结核患者分别用两种方法检查结果 甲法乙法合计 + — + 50 28 78 —12 20 32 合计62 48 110 3、某市中心血站采集血清标本673例，检出乙肝核心抗体55例，其血型分布如下 表。问不同血型的乙肝核心抗体阳性率有否差别？ 乙肝核心抗体的阳性率与血型关系 血型检测数乙肝核心抗体阳性数阳性率（%） A 221 13 5．88 B 169 16 9．47 O 235 23 9．79 AB 48 3 6．25 合计673 55 8．17 4、为了解年龄与总胆固醇之间的数量关系，随机抽查了20名高血脂病人治疗前的总胆固醇水平（mmol/L）与年龄，测定值见表9.1。试求胆固醇对年龄的直线回归方程和相关系数。
5、某医院用中药结合化疗（中药组）和化疗（对照组）两种疗法治疗18例白血病患者，最初的缓解期（天）如下，试绘制两组的缓解率曲线并作比较。 中药组：10，13+，26，32，89，169+，198+，205+，220，254+ 对照组：4，12，20+，32，45，69，103，254+ （注：有+号的数据为截尾值） 6、某疾控中心对该地2496名55岁以下男性做了一次吸烟和喝咖啡对心肌梗死影响的病例对照研究，调查结果见下表。 吸烟及喝咖啡与心肌梗死关系调查结果 喝咖啡（杯/天） 吸烟（支/天） 0 1-20 21-35 35以上 病例对照病例对照病例对照病例对照 0 66 123 30 52 15 12 36 13 1-2 141 179 59 45 53 22 69 25 3-4 113 106 63 65 55 16 119 30 5以上129 80 102 58 118 44 373 85 请拟合下列模型： （1）模型1将喝咖啡(cofe)和吸烟(smoke)两个变量作为有序分类变量（0，1，2，3，），拟合模型。 （2）模型2将喝咖啡(cofe)和吸烟(smoke)两个变量以0为基准设为哑变量（0，1），拟合模型 （3）模型3将喝咖啡和吸烟两个变量作为交互效应变量，拟合模型。 7、根据练习一的数据，分析影响血压的主要因素有哪些？ 8、某单位研究三种不同穴位的镇痛效果，资料如表4-16。试检验针刺不同穴位的镇痛效果有无差别？