重庆医科大学
卫生统计学统计软件包SAS
上机实习题
重庆医科大学卫生统计学教研室 彭斌 编写
2010年12月
SAS上机练习题(一) 本SAS习题集由彭斌编写 2010.11.16
卫生统计学统计软件包上机实习题
SAS上机练习题(一)
1、SAS常用的窗口有哪三个?请在三个基本窗口之间切换并记住这些命令或功能键。
2、请在PGM窗口中输入如下几行程序,提交系统执行,并查看OUTPUT窗和LOG窗中内容,注意不同颜色的含义;并根据日志窗中的信息修改完善程序。 DATS EX01; 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、试根据某班12名学生3门功课成绩表完成后面的问题:
表1 某班12名学生3门功课成绩表
学号 083 084 085 087 084 085 086 089 118 083 084 099 问题:
1) 建立数据集;
2) 打印至少有1门功课不及格同学的信息;(提示,使用if语句)
生化 68 74 73 75 85 78 80 67 70 70 62 82 物理 71 61 75 71 85 79 76 73 54 66 73 70 病理 65 68 46 68 87 75 79 71 75 84 65 79
1
SAS上机练习题(一) 本SAS习题集由彭斌编写 2010.11.16
5、(选做)根据下列数据建立数据集
表2 销售数据
开始时间 2005/04/28 2005 09 18
20040508
终止时间 25MAY2009 05OCT2009 30JUN2009
费用 $123,345,000 $33,234,500 $432,334,500
提示:(格式化输入;数据之间以空格分隔,数据对齐;注意格式后面的长度应以前一个位置结束开始计算,如果读入错误,可试着调整格式的宽度;显示日期需要使用输出格式) 开始时间,输入格式yymmdd10. 终止时间,输入格式date10. 费用,输入格式dollar12.
例如:input x1 yymmdd10. X2 date10. X3 dollar12.
6、手机号码一编码规则一般是:YYY-XXXX-ZZZZ,其YYY为号段;XXXX一般为所在地区编码;ZZZZ为对应的个人识别编号。下面有一组电话号码(来源于网络,末位以X替换),请用程序完成下列要求: (1)分别列出属于联通、移动、电信的号码; (2)分别提取地区编号及个人识别编号。
(提示:列输入方式或者字符串操作substr()函数) 说明:
移动:134-139、150、151、152、157、158、159、188 联通:130、131、132、155、156 电信:133、153、180、189
手机号码 1508320464X 1510291126X 1339644156X 1532464828X 1365721647X 1523105754X 1592624347X 1331237668X 1370048578X 1507244457X 1321246707X 1368464734X 1394913242X 手机号码 1313654836X 1351425709X 1303612231X 1383824855X 1522322763X 1357851051X 1508311759X 1327313520X 1556443719X 1804346016X 1513441713X 1308279203X 1897876466X
7、下表是某班学生几门功课的成绩,其中形势课是考查课,其它均为考试课。请完成以下处理并且保存程序、结果和日志:(注意函数在DATA STEP中使用的位置) (1)用Mean()函数求出每位同学的平均分,将其保存在变量中;
(2)用sum()函数或者表达式求出每位同学的总分,将其保存在变量中;
(3)如果每门考试课权重是0.75,每门考查课权重为0.25,试求每位同学的加权平均分; (4)将数据集保存为永久数据集,保存位置为“D:\\mydata\\”; (5)调用上面的永久数据集,并将男、女生的成绩分别保存在数据集Dmale和Dfemale中。
2
SAS上机练习题(一) 本SAS习题集由彭斌编写 2010.11.16
表3 某班同学几门功课的成绩
学号 0083 0414 0564 0774 0873 0874 0875 0876 0877 0878 0880 0881 0883 0884 0885 0886 0887 0888 0889 0890 0891 0892 性别 高数 (0=女,1=男) 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 73 90 97 40 68 74 73 79 75 76 81 76 77 85 78 80 85 77 67 75 69 66 生理 73 79 87 60 65 68 46 79 68 60 73 81 60 87 75 79 82 75 71 70 76 71 人解 64 71 89 61 60 56 65 74 55 64 60 84 60 84 62 77 81 74 71 75 76 60 数理统计 形势(考查) 74 85 91 65 76 60 66 89 60 71 75 89 66 85 66 73 84 67 69 69 79 60 75 78 80 75 75 75 87 87 70 80 82 78 75 90 80 83 80 80 90 80 90 78
8、下面是3个大类疾病的ICD-10编码及对应的疾病名。请完成以下任务: (1)建立数据集;
(2)提取每种疾病的大类编码;
(3)分别将3个大类的疾病存入3个数据集。
(提示:ICD10编码中小数点前面的三位表示大类;length语句定义字符变量长度;字符串取子串函数)
表4 4类疾病的ICD10编码及对应疾病名
ICD10 A01.001 A01.002 A01.101 A01.201 A01.301 A01.401 A02.001 A02.002 A02.004 A02.006 A02.007
Disease 伤寒
伤寒杆菌性败血症 甲型副伤寒 乙型副伤寒 丙型副伤寒 副伤寒
B群沙门氏菌肠炎 C群沙门氏菌肠炎 沙门氏菌性肠炎 沙门氏菌胃肠炎 鼠伤寒沙门氏菌性肠炎
3
SAS上机练习题(一) 本SAS习题集由彭斌编写 2010.11.16
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、现有两个文件,内容如下: 婴儿沙门氏菌肠炎 沙门氏菌败血症 沙门氏菌感染 鼠伤寒沙门氏菌感染 沙门氏菌属食物中毒 什密氏志贺菌痢疾 弗氏志贺菌痢疾 鲍氏志贺菌痢疾 宋内氏志贺菌痢疾 菌痢混合感染 细菌性痢疾
慢性细菌性痢疾急性发作 慢性细菌性痢疾 中毒性痢疾
A文件 学号 统计 体育 年龄 1903 87 83 23 1902 56 96 22 1901 93 75 19 1904 77 84 24 1908 88 55 18 B文件 学号 性别 班级 1901 男 1 1902 女 1 1903 男 3 1904 男 2 1905 女 2 1908 女 3 要求:(1)根据A、B文件分别建立数据集A、B; (2)打印至少有一门功课不及格的同学的年龄、性别和班级。
10、某研究员欲分析急性染毒对肝脏功能的影响,将40只小鼠分为两组,雄雌各半,试验组进行急性染毒试验,染毒后2小时测定血液中的ALT,整理的结果见下表。
ALT(丙氨酸转氨酶) NO 1 2 3 4 5 6 7 8 9 38.4 30.6 25.5 29.9 29.0 33.6 39.4 42.5 46.7 Male Control Treatment 57.8 62.4 57.5 49.9 55.6 47.3 50.1 48.7 57.9 60.4 Control 39.7 34.6 40.5 44.3 41.2 25.6 30.6 37.6 39.5 25.1 Female Treatment 55.7 63.8 59.9 49.7 48.8 51.3 53.8 42.6 61.7 65.3 10 49.3 要求:根据上面的数据建立恰当的SAS数据集
4
SAS上机练习题(一) 本SAS习题集由彭斌编写 2010.11.16
11、某职业病防治所对30名矿工分别测定血清铜蓝蛋白含量(μmol/L),资料如下。问各期血清铜蓝蛋白含量的测定结果有无差别?
疾病分期 0期
I期 II期
测定结果
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
要求:根据上面的数据建立恰当的SAS数据集
12、24名志愿者随机分成两组,每组12人,接受降胆固醇试验,甲组为特殊饮食组,乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇(mmol/L),数据见下表:
甲组 受试者 1 2
3 4 5 6 7 8 9 10 11
试验前 6.11 6.81 6.48 7.59 6.42 6.94 9.17 7.33 6.94 7.67 8.15
试验后 6.00 6.83 4.49 7.28 6.30 6.64 8.42 7.00 6.58 7.22 6.57
受试者 1 2 3 4 5 6 7 8 9 10 11 12
乙组 试验前 6.90 6.40 6.48 7.00 6.53 6.70 9.10 7.31 6.96 6.81 8.16 6.98
试验后 6.93 6.35 6.41 7.10 6.41 6.68 9.05 6.83 6.91 6.73 7.65 6.52
12 6.60 6.17
要求:根据上面的数据建立恰当的SAS数据集
(彭斌,2010-10-10)
5
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16
SAS上机练习题(二)
1、某农村地区1999年14岁女孩的身高资料列于表1。 142.3 137.9 142.4 148.8 140.5 150.8 141.8 144.0
148.8 143.9 141.5 147.8 143.4 147.9 148.1 129.4
表1 某农村地区1999年14岁女孩身高资料(cm)
142.7 144.4 144.7 145.1 143.3 154.2 152.3 141.2 141.9 146.7 152.9 136.9 145.4 142.8
139.3 147.9 132.7 147.5 146.5 134.6 132.1
145.8 125.8 149.7 147.7 130.6 130.5 141.8
142.2 139.9 154.0 162.6 142.5 145.2 143.3
137.9 148.9 158.2 141.6 149.0 146.2 143.8
141.2 154.9 138.2 143.6 145.4 146.4 134.7 142.9 146.7 138.9 142.9
150.6 145.7 149.8 144.0 139.5 142.4 147.1 141.4 138.5 143.5 126.8
142.7 142.7 140.8 151.1 150.6 148.9 137.1 140.9 134.7 139.6 139.9
156.6 151.3 139.6 140.1 138.9 144.5 141.4 137.4 138.5 139.2 134.4
142.5 146.6 135.5 146.8 156.3 150.0 147.3 146.6 134.5 135.1 141.9 142.1 138.1 134.9 148.8 150.3 140.7 143.5 140.2 143.6 138.7 133.1 145.9 139.2 137.4 142.3 160.9 137.7 问题:(1)计算均数、中位数;计算均数的95%可信区间; (2)计算四分位间距、标准差、变异系数;计算标准误; (3)请进行正态性检验; (4)观察频数分布情况;
2、(t检验、配对t检验)24名志愿者随机分成两组,每组12人,接受降胆固醇试验,甲组为特殊饮食组,乙组为药物治疗组。受试者试验前后各测量一次血清胆固醇(mmol/L),数据见下表:
甲组
受试者 1 2 3 4 5 6 7 8 9 10 11 12
试验前 6.11 6.81 6.48 7.59 6.42 6.94 9.17 7.33 6.94 7.67 8.15 6.60
试验后 6.00 6.83 4.49 7.28 6.30 6.64 8.42 7.00 6.58 7.22 6.57 6.17
受试者 1 2 3 4 5 6 7 8 9 10 11 12
6.90 6.40 6.48 7.00 6.53 6.70 9.10 7.31 6.96 6.81 8.16 6.98
乙组 试验前
试验后 6.93 6.35 6.41 7.10 6.41 6.68 9.05 6.83 6.91 6.73 7.65 6.52
问题:
(1)试建立能够满足以下统计分析的SAS数据集。
(2)试判断两组受试者试验前血清胆固醇水平是否相等? (3)分别判断两种降胆固醇措施是否有效?
(4)试判断两种降胆固醇措施的效果是否相同?
6
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16
3、(秩和检验)某职业病防治所对30名矿工分别测定血清铜蓝蛋白含量(μmol/L),资料如下。问各期血清铜蓝蛋白含量的测定结果有无差别?
疾病分期 测定结果 0期
I期 II期
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
请用秩和检验完成分析。 4、(频数表等)以下是对某年级的问卷调查结果,问题表内容如下: id sex date height weight money 机器号 性别:1=male,2=female 记忆中影响最大的日子是哪一天? 身高是多少 体重是多少 昨天吃饭花多少钱(元)? question1 question2 46级改革之后学习英语的积极性:a=增加,b=不变,c=减弱 结核病的早期症状有(可以多选): a=咳嗽,b=痰中带血,c=盗汗,d=贫血,e=胸痛,f=午后潮热,g=胃肠不适
请进行以下分析:
1) 此次调查中男性的比例为多少;
2) 对46级改革之后学习英语的积极性的情况进行描述,三种选择占的比例分别是多少? 3) 男生与女性的对46级改革的反应有无差异;
4) 计算身体质量指数(Body Mass Index,BMI?体重(Kg)身高2(m2)),比较男女生的BMI有无差异? 5) 关于BMI与人体的肥胖程度通常可以根据下面的标准来判断,请列出每一类的人数及比例;
偏瘦: 正常体重: 超重: 轻度肥胖: BMIBMIBMIBMI指数 指数 指数 指数 < 18 = 18 - 25 = 25 - 30 > 30 中度肥胖: BMI指数 > 35 重度肥胖: BMI指数 > 40
机器号 id
cnw1l01 cnw1l02 cnw1l03 cnw1l04 cnw1l05 cnw1l06 cnw1r01 cnw1r02 cnw1r03 cnw1r05 cnw1r06 cnw2l02 cnw2l03
性别 sex 1 1 2 2 2 2 2 2 1 2 1 2 2
身高 height 179 175 157 163 161 150 158 158 168 150 162 160 164
体重 weight 70 70 47 48 52 44 56 56 64 44 53 53 58
开支 money 5.7 7.5 4.5 5 5 5.5 4.4 6.2 6.5 5.5 7 8 6
英语改革 question1 a a a c b b b b a b b c c
TB症状 question2 ABDEF ABE ABE ABDF ABF ABD ADFG ABDE ABE ABD ABDE ABEF ABE
7
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16 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 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 8
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16 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
cnw8l03
2
160
60
7
A
ABDEF
9
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16
补充练习题
(摘自《医学统计学》,颜虹主编,人民卫生出版社,第一版)
P132:(t检验)第2,4题
2. 某医生研究脑缺氧对脑组织中生化指标的影响,将出生状况相近的乳猪按出生体重配成7对共2组,随机接受两种处理,一组设为对照组,一组设为脑缺氧模型组,试验结果如表。试比较两种处理的脑组织钙泵的含量有无差别。 乳猪号 1 2 3 4 5 6 7
对照组 0.3550 0.2000 0.3130 0.3630 0.3544 0.3450 0.3050
实验组 0.2755 0.2545 0.1800 0.1800 0.3113 0.2955 0.2870
4. 随机抽取20只小鼠分配到A,B两个不同饲料组,每组10只,在喂养一定时间后,测得鼠肝中铁的含量(μg/g),数据如下。试问不同饲料对鼠肝中铁的含量有无影响。 A组 B组
P149:(方差分析)第1,3,4题
1. 将12份同一个体的血液标本随机等分成3组,分别用甲,乙,丙3种抗凝剂处理后测红细胞沉降率。试验结果如下,问3种抗凝剂作用的差异是否有统计学意义? 抗凝剂 甲 乙 丙
红细胞沉降率(mm/h) 17 10 11
16 11 9
16 12 8
15 12 9
3. 比较下列试验结果(如下表)雌激素3个剂量组之间的差异有无统计学意义。 种系 甲 乙 丙 丁
10
3.59 2.23
0.96 1.14
3.89 2.63
1.23 1.00
1.61 1.35
2.94 2.01
1.96 1.64
3.68 1.13
1.54 1.01
2.59 1.70
雌激素注射剂量 0.2 106 42 70 42
0.4 116 68 111 63
0.8 145 115 133 87
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16
4.分析下列拉丁方设计的试验结果(下表),其中字母A,B,C,D表示胰岛素注射的4种剂量(U/kg),分别为A:0.32,B:0.47,C:0.62,D:0.77。问4种剂量下的测量结果的差异有无统计学意义? 实验日期 1 2 3 4
1 C B A D
家兔编号 32.7 26.2 14.0 33.2
2 A D C B
11.2 31.8 14.0 16.5
3 B C D A
23.2 28.9 27.5 21.2
4 D A B C
48.1 18.7 25.6 40.2
P168:(卡方检验)第1,3题
1.某医师用两种疗法治疗脑血管梗死,结果如下表,是比较两种疗法的疗效是否相同? 疗法 甲疗法 乙疗法 合计
有效 25 29 54
无效 6 3 9
合计 31 32 63
有效率(%) 80.65 90.63 85.71
3.有50份痰液标本,每份分别接种在甲乙两种培养基中,观察结核杆菌的生产情况,结果如下表,是比较两种培养基的效果。 甲培养基 + - 合计
P184:(秩和检验)第2,3,6题
2. 某研究中心观察溶脲脲原体(UU)感染对家兔精子质量的影响,分别测定十只兔子感染UU前后的精子密度(10E9/L),数据见下表,试分析溶脲脲原体是否影响家兔的精子密度。 编号 感染前 感染后
1 336 258
2 371 291
3 386 300
4 364 285
5 377 298
6 292 303
7 288 312
8 304 260
9 333 339
10 302 290
乙培养基 + 23(a) 7(c) 30
- 12(b) 8(d) 20
合计 35 15 50
11
SAS上机练习题(二) 本SAS习题集由彭斌编写 2010.12.16
3. 用中草药治疗不同类型的小儿肺炎,其疗效分为4个等级,结果见表,试比较该药物治疗不同类型的小儿肺炎的疗效有无差别? 疗效 控制 显效 有效 无效 合计
病毒性肺炎 65 18 30 13 126
细菌性肺炎 42 6 23 11 82
合计 107 24 53 24 208
6. 4种疾病患者痰液内嗜酸性粒细胞的检查结果如下表,问4种疾病患者痰液内嗜酸性粒细胞的等级分布有无差别?
白细胞等级 支气管扩张 肺水肿 - + ++ +++ 合计
(彭斌,2010-10-10)
1 3 10 7 21
4 6 6 2 18
肺癌 7 9 5 4 25
病毒性呼吸道感染 5 7 5 1 18
12
SAS上机练习题(三) 本SAS习题集由彭斌编写 2010.12.16
SAS上机练习题(三)
1、(方差分析)某研究人员研究3种刺激因素(高葡萄糖、高胰岛素和过氧化氢)刺激大鼠肾小球系膜细胞对单核细胞趋化蛋白表达的影响,在四种条件下对大鼠肾小球系膜细胞进行培养,检测结果如表1。
表1 不同刺激因素作用下蛋白表达结果
分组 对照组
高葡萄糖组 高胰岛素组 过氧化氢组
0.3547 1.2247 1.1453 1.3582
表达结果(光密度值) 0.3652 1.2354 1.1526 1.3657
0.3842 1.2136 1.1375 1.3518
0.3598 1.2286 1.1547 1.3259
0.3742 1.2415 1.1294 1.3726
0.3925 1.2275 1.1659 1.3318
问题1:计算每一组的平均光密度,观察各个组间均数的差异大小;
问题2:进行方差分析,将结果整理成方差分析表,观察变异、均方及F值之间的关系; 问题3:分别采用SNK法和Dunnett法进行多重比较; 问题4:对上述结果进行解释。
2、(方差分析)研究枸杞多糖对脂肪肝的预防作用,按窝别作为区组标志,每一区组3只大鼠,随机分配到三个组:生理盐水组、酒精组、酒精+枸杞多糖组,三组灌胃5周,检测肝脏中谷胱甘肽(GSH)的含量(mg/gprot),结果见表2。
表2 三组小鼠肝脏中谷胱甘肽(GSH)的含量(mg/gprot)
区组 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
酒精组 30.48 31.25 33.28 34.61 28.35 29.17 27.34 30.58 34.25 27.31 28.09 30.45 33.25 34.04 34.25
LBP预防组 65.08 63.04 67.59 68.58 64.12 66.55 66.89 67.15 68.05 65.48 64.38 65.04 66.84 67.56 67.46
生理盐水 79.15 75.46 79.32 75.98 76.55 80.34 84.35 88.14 87.35 72.15 74.61 86.33 94.35 92.05 96.42
问题1:按照随机区组设计进行分析,处理组间、区组间效应是否有差异?
问题2:忽略区组,按照完全随机设计进行分析,处理组间效应是否有差异?
问题3:对上述两个方差分析表进行比较,找出相同项与不同项?观察处理组间比较的F值变化,会得出怎样的结论?
问题4:根据分析结果对该研究下一个结论。
13
SAS上机练习题(三) 本SAS习题集由彭斌编写 2010.12.16
3、(方差分析)研究人员研究某种物质的毒性,将40只小鼠分为两组,雄雌各半,试验组给予研究物质,2小时后测定血液中碱性磷酸酶的含量,结果如下。
表3 40只小鼠给予不同物质后血液中碱性磷酸酶的含量
性别 分组
碱性磷酸酶
374.5 432.7 401.3 399.4 367.5 421.6 489.0 432.5 421.0 420.4 341.3 322.5 365.4 321.6 401.9 460.5 420.1 394.4 389.6 420.5
雄性 对照组 367.9 408.6 375.6 354.9 421.7 试验组 423.8 446.9 432.5 478.1 437.5 雌性 对照组 378.1 345.2 390.6 399.0 421.1 试验组 420.4 473.2 450.3 405.5 427.4 问题1:研究者考虑了几个影响因素? 问题2:性别对碱性磷酸酶有无影响? 问题3:该物质对碱性磷酸酶有无影响?
4、(回归分析)以下是制作标准曲线时对应的浓度和吸光度值,试求出标准曲线。并求出A=1.15时的浓度。
表4 不同浓度对应的吸光度值 浓度 0.97 1.27 1.57 1.88
2.18
5、(相关分析、回归分析)27名糖尿病人的血清总胆固醇、甘油三脂、空腹胰岛素、糖化血红蛋白、空腹血糖的测量值列于表中,试分析血糖与其它几项指标关系的关系。
表5 27名糖尿病人的血糖及有关变量的测量结果
总胆固醇
序号i 1 2 3 4 5 6 7 8 9 10 11 12 13 14
A值 0.3 0.56 0.93 1.35 1.51
甘油三脂 (mmol/L) X2 1.90 1.64 3.56 1.07 2.32 0.64 8.50 3.00 2.11 0.63 1.97 1.97 1.93 1.18
(mmol/L)
X1 5.68 3.79 6.02 4.85 4.60 6.05 4.90 7.08 3.85 4.65 4.59 4.29 7.97 6.19
胰岛素 (?U/ml) X3 4.53 7.32 6.95 5.88 4.05 1.42 12.60 6.75 16.28 6.59 3.61 6.61 7.57 1.42
糖化血 红蛋白(%)
X4 8.2 6.9 10.8 8.3 7.5 13.6 8.5 11.5 7.9 7.1 8.7 7.8 9.9 6.9
血 糖 (mmol/L)
Y 11.2 8.8 12.3 11.6 13.4 18.3 11.1 12.1 9.6 8.4 9.3 10.6 8.4 9.6
14
SAS上机练习题(三) 本SAS习题集由彭斌编写 2010.12.16
15 16 17 18 19 20 21 22 23 24 25 26 27
6.13 5.71 6.40 6.06 5.09 6.13 5.78 5.43 6.50 7.98 11.54 5.84 3.84
2.06 1.78 2.40 3.67 1.03 1.71 3.36 1.13 6.21 7.92 10.89 0.92 1.20
10.35 8.53 4.53 12.79 2.53 5.28 2.96 4.31 3.47 3.37 1.20 8.61 6.45
10.5 8.0 10.3 7.1 8.9 9.9 8.0 11.3 12.3 9.8 10.5 6.4 9.6
10.9 10.1 14.8 9.1 10.8 10.2 13.6 14.9 16.0 13.2 20.0 13.3 10.4
问题1:计算各指标之间的pearson相关系数及spearman相关系数; 问题2:以血糖为应变量Y其它为自变量X,进行回归分析。
6、(综合分析)以下是对某年级的问卷调查结果,问题表内容如下: id sex date height weight money question1 机器号 性别:1=male,2=female 记忆中影响最大的日子是哪一天? 身高是多少 体重是多少 昨天吃饭花多少钱(元)? 46级改革之后学习英语的积极性:a=增加,b=不变,c=减弱 question2 结核病的早期症状有(可以多选): a=咳嗽,b=痰中带血,c=盗汗,d=贫血,e=胸痛,f=午后潮热,g=胃肠不适
请进行以下分析:
1) 试对男女生的开支进行分析;
2) 析被调查者对“结核病的早期症状”的选择情况,即列出每种选项的频数及比例;(提示:可分别对每
个选项单独进行分析,使用index()函数判断是否选择某个选项,具体用法请查帮助。) 3) 分析身高与体重的相关性;建立身高预测体重的预测模型。
(原始数据请参见“SAS练习题(二)第4题”Page 7)
(彭斌 2010-10-10 整理)
15
SAS上机练习题(四) 本SAS习题集由彭斌编写 2010.12.16
SAS上机练习题(四)
1. 用三种不同治疗方法治疗慢性支气管炎的疗效如表1所示,试比较三种治疗方法治疗慢性支气管炎的疗效。
表1 三种不同治疗方法治疗慢性支气管炎的疗效
组别 A药 B药 C药 合计
要求:
(1)用卡方检验分析; (2)用logistic回归分析。
2、欲研究胎儿颅内出血的危险因素,某研究人员收集了下列相关资料,请进行logistic回归分析。
妊高症等 高危因素 1=有,0=无
1 1 1 1 0 0 0 1 1 0 1 0 0 0 0 0 1 1 1 1 0 1 1 0 1 0 1 1
胎盘异常 1=有,0=无
0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 1 0 0 1 0 0 0 0 0 0
脐带绕颈 1=有,0=无
0 0 0 0 0 1 1 0 1 0 0 0 0 0 1 0 0 0 1 1 1 1 0 0 1 0 0 0
羊水异常 1=有,0=无
0 1 1 1 1 1 1 0 1 0 0 0 0 0 0 1 0 0 1 1 0 0 1 1 1 0 0 0
产时窒息 1=有,0=无
0 1 0 0 0 1 0 1 1 1 0 0 0 0 0 1 1 1 0 1 1 0 0 0 0 0 0 0
酸中毒 1=有,0=无
0 1 1 1 0 0 1 0 1 1 0 0 0 1 1 1 0 0 0 1 0 0 1 1 1 1 0 0
颅内出血 1=有,0=无
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0
16
有效 35 20 7 62
无效 5 10 25 40
合计 40 30 32 102
有效率% 87.50 66.67 21.88 60.78
注:数据来源《卫生统计学》(方积乾主编,第六版)
SAS上机练习题(四) 本SAS习题集由彭斌编写 2010.12.16
0 0 0 0 1 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0
1 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 1 0 0
1 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 1 0 0
0 1 1 0 0 0 1 0 0 1 1 1 0 0 0 0 1 0 1 0 0 0 0
1 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
17