SAS上机练习题(全部,含参考答案)
重庆医科大学--卫生统计学统计软件包
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
提示：(
计算，如果读入错误，可试着调整格式的宽度；显示日期需要使用输出格式) 开始时间，输入格式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
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
参考程序：
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）将数据集保存为永久数据集，保存位置为“D:\mydata\”；
（5）调用上面的永久数据集，并将男、女生的成绩分别保存在数据集Dmale和Dfemale中。
参考程序：
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
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=1then output dmale;
if xb=0then 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 丙型副伤寒
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 沙门氏菌性肠炎
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、现有两个文件，内容如下：
要求：打印至少有一门功课不及格的同学的年龄、性别和班级。
参考程序：
data fileA;
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 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=100to200;
flag=0;
do j=2to i-1;
if mod(i,j)=0then flag=1;
end;
if flag=0then output;
end;
keep i;
run;
proc print;
run;
(彭斌，2010-3-23)
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 means data=a n mean median clm qrange std cv stderr;
var x;
proc univariate data=a normal;
histogram x;
var x;
2、某研究员欲分析急性染毒对肝脏功能的影响，将40只小鼠分为两组，雄雌各半，试验组进行急性染毒试验，染毒后2小时测定血液中的ALT，整理的结果见下表。
（1）比较实验组与对照组的ALT有无差异；
参考程序：
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;
3、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）试判断两种降胆固醇措施的效果是否相同？
参与程序：
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.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
;
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期';
else if _n_<=20then 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;
教材上的练习题
(卫生统计学，方积乾主编，第六版)
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_<=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;
第四题
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，2×2表，卡方检验
data a1;
do row=1to2;
do col=1to2;
input f@@;
output;
end;
end;
cards;
64 21
51 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 2
1 2 10
2 1 14
2 2 14
;
run;
proc freq data=a2 ;
tables row*col/chisq;
weight f;
run;
例9.4，3×2表，卡方检验
data a3;
do row=1to3;
do col=1to2;
input f@@;
output;
end;
end;
cards;
35 5
20 10
7 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.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=1to4;
do col=1to3;
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;
p193 7题，等级资料，配对设计，用Friedman检验
data a3;
do id=1to10;
do grp= 1to2;
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-4-18)
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=1to4;
do i=1to6;
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.1453 1.1526 1.1375 1.1547 1.1294 1.1659
1.3582 1.3657 1.3518 1.3259 1.3726 1.3318
;
或者
data a;
length grp $10;
input grp $ @@;
do i=1to6;
input x@@;
output;
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.1453 1.1526 1.1375 1.1547 1.1294 1.1659
过氧化氢组 1.3582 1.3657 1.3518 1.3259 1.3726 1.3318
;
SAS上机练习题(二)参考答案
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上机练习题(全部,含参考答案)
重庆医科大学--卫生统计学统计软件包 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上机实习-8页文档资料
中科院研究生院统计分析与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（进料
统计实验与SAS上机简易过程步
数据统计分析一般可遵循以下思路: (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上机练习题(全部-含参考答案).docx
重庆医科大学“卫生统计学统计软件包 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上机练习题(全部_含参考问题详解)
医科大学--卫生统计学统计软件包 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. 参考程序：
SAS上机实习考试
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课程上机练习一
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上机练习题(全部-含参考答案)
————————————————————————————————作者：————————————————————————————————日期：
重庆医科大学--卫生统计学统计软件包 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.
SAS上机实验
实验内容 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 ; 实验结果：
0九研究生SAS软件上机试卷(1)
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）
最新SAS上机实习
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
new_SAS课程上机练习四
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。试检验针刺不同穴位的镇痛效果有无差别？