SAS学习第12章：习题练习

例1：11只60日龄的雄鼠在x射线照射前后之体重数据见下表（单位：g）：检验雄鼠在照射射线前后体重差异是否显著？

编  号	1	2	3	4	5	6	7	8	9	10	11
照射前	25.7	24.4	21.1	25.2	26.4	23.8	21.5	22.9	23.1	25.1	29.5
照射后	22.5	23.2	20.6	23.4	25.4	20.4	20.6	21.9	22.6	23.5	24.3

data a1;
input id x1 x2@@;
d=x1-x2;
cards;
1 25.7 22.5
2 24.4 23.2
3 21.1 20.6
4 25.2 23.4
5 26.4 25.4
6 23.8 20.4
7 21.5 20.6
8 22.9 21.9
9 23.1 22.6
10 25.1 23.5
11 29.5 24.3
;
proc means mean stderr t prt;
var d;
run;

例2：某生物药品厂研制出一批新的鸡瘟疫苗，为检验其免疫力，用200只鸡进行试验，其中注射100只（经注射后患病的10只，不患病的90只），对照组（注射原疫苗组）100只（经注射后患病的15只，不患病的85只），试问新旧疫苗的免疫力是否有差异。

data a2;
input r c f@@;
cards;
1 1 10
1 2 90
2 1 15
2 2 85
;
proc freq;
table r*c/chisq;
weight f;
run;

例3：随机抽测了10只兔的直肠温度，其数据为：38.7、39.0、38.9、39.6、39.1、39.8、38.5、39.7、39.2、38.4（℃），已知该品种兔直肠温度的总体平均数=39.5（℃），试检验该样本平均温度与是否存在显著差异？

data a3;
input x@@;
y=x-39.5;
cards;
38.7 39.0 38.9 39.6 39.1 39.8 38.5 39.7 39.2 38.4
;
proc means n mean stderr t prt;
var y;
run;

例4：

甲、乙、丙三个奶牛场高产奶牛、低产奶牛、低产奶牛头数统计如下，试问三个奶牛场高、中、低产奶牛的构成比是否有差异。

场  地	高产奶牛	中产奶牛	低产奶牛
甲	32	15	18
乙	28	26	10
丙	11	38	16

data a4;
do r=1 to 3;
do f=1 to 3;
input y@@;
output;
end;
end;
cards;
32 15 18 28 26 10 11 38 16
;
PROC freq;
table r*f/chisq;
weight y;
run;

例5：从7名健康人体中抽取血样，每个受试者的血样分成4份，然后随机地将4份血样分别放置10、20、30、40 min测定其血糖浓度，资料如下表。试检验不同受试者和放置不同时间的血糖浓度有无显著差异。

不同受试者、放置不同时间血滤液的血糖浓度(mg/100ml)

受试者编号	放置时间(min)
	10	20	30	40
1	95	95	89	83
2	95	94	88	84
3	106	105	97	90
4	98	97	95	90
5	102	98	97	88
6	112	112	101	94
7	105	103	97	88

data a5;
input A$ B$ y@@;
cards;
A1 B1 95 A1 B2 95 A1 B3 89 A1 B4 83
A2 B1 95 A2 B2 94 A2 B3 88 A2 B4 84
A3 B1 106 A3 B2 105 A3 B3 97 A3 B4 90
A4 B1 98 A4 B2 97 A4 B3 95 A4 B4 90
A5 B1 102 A5 B2 98 A5 B3 97 A5 B4 88
A6 B1 112 A6 B2 112 A6 B3 101 A6 B4 94
A7 B1 105 A7 B2 103 A7 B3 97 A7 B4 88
;
proc anova;
class A B;
model y=A B;
means A B/snk;
means A B/snk alpha=0.01;
run;

例6：研究甲、乙两药对某病的治疗效果，甲药治疗病畜70例，治愈53例；乙药治疗75例，治愈62例，问两药的治愈率是否有显著差异？ 

data a6;
input r c f@@;
cards;
1 1 53
1 2 17
2 1 62
2 2 13
;
proc freq;
table r*c/chisq;
weight f;
run;

例7：上午和下午各测8只家兔的血糖值，结果如下。试比较其血糖值上、下午有无变化。

───────────────────────

兔号         上午血糖(mg%)       下午血糖(mg%)

───────────────────────

1                100               105

2                110               115

3                120               128

4                98                104

5                92                102

6                138               134

7                100               110

8                125               130

───────────────────────

data a7;
input x1 x2@@;
d=x1-x2;
cards;
100 105
110 115
120 128
98 104
92 102
138 134
100 110
125 130
;
proc means mean stderr t prt;
var d;
run;

例8：在同样饲养管理条件下，三个品种猪的增重如下表，试对三个品种增重差异是否显著进行检验。

品种	增       重xij(kg)
A1	16	12	18	18	13	11	15	10	17	18
A2	10	13	11	9	16	14	8	15	13	8
A3	11	8	13	6	7	15	9	12	10	11

data a8;
do a=1 to 3;
do i=1 to 10;
input y@@;
output;
end;
end;
drop i;
cards;
16 12 18 18 13 11 15 10 17 18
10 13 11 9 16 14 8 15 13 8
11 8 13 6 7 15 9 12 10 11
;
PROC ANOVA;
class A;
MODEL y=A;
means A/lsd duncan SnK;
means A/lsd duncan SnK alpha=0.01;
run;

例9：为了比较4种饲料(A)和猪的3个品种(B)，从每个品种随机抽取4头猪(共12头)分别喂以4种不同饲料。随机配置，分栏饲养、位置随机排列。从60日龄起到90日龄的时期内分别测出每头猪的日增重(g),数据如下，试检验饲料及品种间的差异显著性。

4种饲料3个品种猪60～90日龄日增重

	A1	A2	A3	A4
B1	505	545	590	530
B2	490	515	535	505
B3	445	515	510	495

data a9;
input B$ A$ y@@;
cards;
B1 A1 505 B1 A2 545 B1 A3 590 B1 A4 530
B2 A1 490 B2 A2 515 B2 A3 535 B2 A4 505
B3 A1 445 B3 A2 515 B3 A3 510 B3 A4 495
;
PROC ANOVA;
class B A;
model y=B A;
means A B/snk;
means A B/snk alpha=0.01;
run;

例10：为了从3种不同原料和3种不同温度中选择使酒精产量最高的水平组合，设计了试验，每一水平组合重复3次，结果如下表，试进行分析。

用不同原料及不同温度发酵的酒精产量

原料	温   度B
	B1(30℃)				B2(35℃)				B3(40℃)
A1	41	49	23		11	12	25		6	22	26
A2	47	59	50		43	38	33		8	22	18
A3	48	35	53		55	38	47		30	33	26

data a10;
do A=1 to 3;
do B=1 to 3;
do n=1 to 3;
input y@@;
output;
end;
end;
END;
drop n;
cards;
41 49 23 11 12 25 6 22 26
47 59 50 43 38 33 8 22 18
48 35 53 55 38 47 30 33 26
;
proc anova;
class A B;
model y=A B A*B;
means A B A*B/SNK;
means A B A*B/SNK alpha=0.01;
run;

例11：

对陕西三个秦川牛保种基地县进行秦川牛肉用性能外形调查，划分为优良中下4个等级，试问三个地区秦川牛肉用性能各级构成比差异是否显著。

地区	优	良	中	下
甲	10	10	60	10
乙	10	5	20	10
丙	5	5	23	6

data a11;
do r=1 to 3;
do c=1 to 4;
input y@@;
output;
end;
end;
cards;
10 10 60 10
10 5 20 10
5 5 23 6
;
proc freq;
table r*c/chisq;
weight y;
run;

例12：10头育肥猪的饲料消耗（x）和增重（y）资料如下表（单位：kg），试对增重与饲料消耗进行直线回归分析。

x	191	167	194	158	200	179	178	174	170	175
y	33	11	42	24	38	44	38	37	30	35

data a12;
input x y@@;
cards;
191 33 167 11 194 42 158 24 200 38 179 44 178 38 174 37 170 30 175 35
;
proc reg corr;
model y=x/clm cli;
run;
plot y*x="+";
run;

例13：

某猪场从10窝大白猪的仔猪中，每窝抽出性别相同、体重接近的仔猪2头，将每窝两头仔猪随机地分配到两个饲料组，进行饲料对比试验，试验时间30天，增重结果见下表。试检验两种饲料喂饲的仔猪平均增重差异是否显著？

窝号	1	2	3	4	5	6	7	8	9	10
饲料Ⅰ	10.0	11.2	12.1	10.5	11.1	9.8	10.8	12.5	12.0	9.9
饲料Ⅱ	9.5	10.5	11.8	9.5	12.0	8.8	9.7	11.2	11.0	9.0

data a13;
input x1 x2@@;
d=x1-x2;
cards;
10.0 9.5
11.2 10.5
12.1 11.8
10.5 9.5
11.1 12.0
9.8 8.8
10.8 9.7
12.5 11.2
12.0 11.0
9.9 9.0
;
proc means mean stderr t prt;
var d;
run;

例14：

试对下列资料进行直线相关和回归分析。

x	36	30	26	23	26	30	20	19	20	16
y	0.89	0.80	0.74	0.80	0.85	0.68	0.73	0.68	0.80	0.58

data a14;
input x y@@;
cards;
36 0.89
30 0.80
26 0.74
23 0.80
26 0.85
30 0.68
20 0.73
19 0.68
20 0.80
16 0.58
;
proc reg corr;
model y=x/clm cli;
run;
plot y*x="+";
run;

例15：为了研究不同种类饲料对奶牛产奶量的影响，设置了A、B、C、D、E 5种饲料，用5头奶牛进行试验，试验根据泌乳阶段分为5期，每期4周，试验结果列于下表，试对其进行方差分析。

  饲料种类对奶牛产奶量影响的试验结果         (单位：kg)

牛  号	时     期
	一	二	三	四	五
Ⅰ	E(300)	A(320)	B(390)	C(390)	D(380)
II	D(420)	C(390)	E(280)	B(370)	A(270)
III	B(350)	E(360)	D(400)	A(260)	C(400)
IV	A(280)	D(400)	C(390)	E(280)	B(370)
V	C(400)	B(380)	A(350)	D(430)	E(320)

data cow;
input a b c d@@;
cards;
1 1 5 300 1 2 1 320 1 3 2 390 1 4 3 390 1 5 4 380
2 1 4 420 2 2 3 390 2 3 5 280 2 4 2 370 2 5 1 270
3 1 2 350 3 2 5 360 3 3 4 400 3 4 1 260 3 5 3 400
4 1 1 280 4 2 4 400 4 3 3 390 4 4 5 280 4 5 2 370
5 1 3 400 5 2 2 380 5 3 1 350 5 4 4 430 5 5 5 320
;
proc anova;
class a b c;
model d=c;
means c/lsd duncan snk;
means c/lsd duncan snk alpha=0.01;
run;

例16：

将20名某病患者随机分成两组，分别用甲、乙两药治疗，测得治疗前后的血沉（mm/h）如下表：

甲乙两药治疗前后的血沉（mm/h）

─────────────────────────────────

甲 药 病人号   1   2   3   4   5    6   7   8   9   10

───────────────────────────── 

      治疗前  10  13   6   11  10   7   8   8   5   9

      治疗后   6   9   3   10  10   4   2   5   3   3

─────────────────────────────────

乙药 病人号    1   2   3   4   5    6   7   8   9   10

─────────────────────────────

     治疗前    9  10   9  13   8    6  10  11   10  10

     治疗后    6   3   5   3   3    5   8   2    7   4

─────────────────────────────────

问：（1）甲乙两药是否均有效（使血沉下降）？

data a1601;
input x1 x2@@;
d=x1-x2;
cards;
10 6
13 9
6 3
11 10
10 10
7 4
8 2
8 5
5 3
9 3
;
proc means mean stderr t prt;
var d;
run;

（2）甲乙两药的疗效有无差别？

data a1602;
input x y@@;
cards;
1 4 1 4 1 3 1 1 1 0 1 3 1 6 1 3 1 2 1 6 
2 3 2 7 2 4 2 10 2 5 2 1 2 2 2 9 2 3 2 6 
;
proc ttest;
class x;
var y;
run;