《生物统计附试验设计》实习习题

1、1、应用克矽平治疗矽肺患者10人，治疗前后血红蛋白含量（g%）如下 。问该药是否会引起血红蛋白含量的变化？（自身配对t检验） 病人号 治疗前 治疗后 1 11.3 14.02 15.0 13.83 15.0 14.04 13.5 13.55 12.8 13.56 10.0 12.07 11.0 14.78 12.0 11.49 13.0 13.810 12.3 12.0data a;input a x y;d=x-y;cards;1 11.3 14.0 2 15.0 13.8 3 15.0 14.0 4 13.5 13.5 5 12.8 13.5 6 10.0 12.0 7 11.0 14.7

2、 8 12.0 11.4 9 13.0 13.810 12.3 12.0proc means n mean std stderr t prt;var d;run; N 均值 标准偏差 标准误差 t 值 Pr > |t| - 10 -0.6800000 1.6457352 0.5204272 -1.31 0.22372、11只60日龄的雄鼠在x射线照射前后之体重数据见下表（单位：g）：检验雄鼠在照射射线前后体重差异是否显著？（自身配对t检验）编 号1234567891011照射前25.724.421.125.226.423.821.522.923.125.129.5照射后22.523.22

3、0.623.425.420.420.621.922.623.524.3 data e;input a x y;d=x-y;cards;1 25.7 22.5 2 24.4 23.2 3 21.1 20.6 4 25.2 23.45 26.4 25.4 6 23.8 20.4 7 21.5 20.6 8 22.9 21.99 23.1 22.6 10 25.1 23.5 11 29.5 24.3proc means n mean std stderr t prt;var d;run; N 均值 标准偏差 标准误差 t 值 Pr > |t| - 11 1.8454545 1.4807860

4、0.4464738 4.13 0.00203、某生物药品厂研制出一批新的鸡瘟疫苗，为检验其免疫力，用200只鸡进行试验，其中注射100只（经注射后患病的10只，不患病的90只），对照组（注射原疫苗组）100只（经注射后患病的15只，不患病的85只），试问新旧疫苗的免疫力是否有差异。（次数资料：卡方） data w;input a b x;cards;1 1 10 1 2 902 1 15 2 2 85proc freq;table a*b/chisq;weight x;run; 统计量 自由度 值 概率 - 卡方 1 1.1429 0.2850 似然比卡方 1 1.1497 0.2836 连续

5、校正卡方 1 0.7314 0.3924 Mantel-Haenszel 卡方 1 1.1371 0.2863 Phi 系数 -0.0756 列联系数 0.0754 Cramer 的 V -0.07564、随机抽测了10只兔的直肠温度，其数据为：38.7、39.0、38.9、39.6、39.1、39.8、38.5、39.7、39.2、38.4（），已知该品种兔直肠温度的总体平均数=39.5（），试检验该样本平均温度与是否存在显著差异？（单个因素t检验）data r;input x;y=x-39.5;cards;38.7 39.0 38.9 39.6 39.1 39.8 38.5 39.7 39

6、.2 38.4proc means n mean std stderr t prt;var y;run;N 均值 标准偏差 标准误差 Pr > |t|10 -0.4100000 0.4909175 0.1552417 0.02695、甲、乙、丙三个奶牛场高产奶牛、低产奶牛、低产奶牛头数统计如下，试问三个奶牛场高、中、低产奶牛的构成比是否有差异。（半定量（等级）资料：只有卡方）场 地高产奶牛中产奶牛低产奶牛甲321518乙282610丙113816 data q;input a b x;cards;1 1 32 1 2 15 1 3 182 1 28 2 2 26 2 3 103 1 11

7、 3 2 38 3 3 16proc freq;table a*b/chisq;weight x;run; a * b 表的统计量统计量 自由度 值 概率-卡方 4 22.8289 0.0001似然比卡方 4 24.6805 <.0001Mantel-Haenszel 卡方 1 4.8178 0.0282Phi 系数 0.3430列联系数 0.3245Cramer 的 V 0.2426样本大小 = 1946、从7名健康人体中抽取血样，每个受试者的血样分成4份，然后随机地将4份血样分别放置10、20、30、40 min测定其血糖浓度，资料如下表。试检验不同受试者和放置不同时间的血糖浓度有无

8、显著差异。(两因素交叉分组单个观测值方差分析)不同受试者、放置不同时间血滤液的血糖浓度(mg/100ml)受试者编号放置时间(min)1020304019595898329594888431061059790498979590510298978861121121019471051039788data t;input a$ b$ x;cards;a1 b1 95 a1 b2 95 a1 b3 89 a1 b4 83a2 b1 95 a2 b2 94 a2 b3 88 a2 b4 84a3 b1 106 a3 b2 105 a3 b3 97 a3 b4 90a4 b1 98 a4 b2 97 a4

9、 b3 95 a4 b4 90a5 b1 102 a5 b2 98 a5 b3 97 a5 b4 88a6 b1 112 a6 b2 112 a6 b3 101 a6 b4 94a7 b1 105 a7 b2 103 a7 b3 97 a7 b4 88proc anova;class a b;model x=a b;means a b/isd duncan snk;means a b/isd duncan alpha=0.01;run; R-Square Coeff Var Root MSE x Mean 0.949781 2.144480 2.066359 96.35714 Source D

10、F Anova SS Mean Square F Value Pr > Fa 6 629.4285714 104.9047619 24.57 <.0001b 3 824.1428571 274.7142857 64.34 <.00017、研究甲、乙两药对某病的治疗效果，甲药治疗病畜70例，治愈53例；乙药治疗75例，治愈62例，问两药的治愈率是否有显著差异？ （次数资料：卡方）data y;input a b x;cards;1 1 53 1 2 172 1 62 2 2 13proc freq;table a*b/chisq;weight x;run;a * b 表的统计量

11、统计量 自由度 值 概率 - 卡方 1 1.0665 0.3017似然比卡方 1 1.0675 0.3015连续校正卡方 1 0.6849 0.4079Mantel-Haenszel 卡方 1 1.0592 0.3034Phi 系数 -0.0858列联系数 0.0855Cramer 的 V -0.0858Variable Method Variances DF t Value 分析变量：d均值 N 标准偏差 标准误差 t 值 Pr > |t|-0.7919048 2 0.0491608 0.0347619 22.78 0.02798、某医院对9例慢性苯中毒患者用中草药抗苯一号和耦合剂治疗

12、，治疗后测定白细胞总数（109/L)，资料如下。问两种药是否对患者的白细胞总数有影响？（计数资料：配对t检验）病人号 抗苯一号 耦合剂1 6.0 4.22 4.8 5.43 5.0 6.34 3.4 3.85 7.0 4.46 3.8 4.07 6.0 5.98 3.5 8.09 4.3 5.0data q;input a x y;d=x-y;cards;1 6.0 4.2 2 4.8 5.4 3 5.0 6.3 4 3.4 3.85 7.0 4.4 6 3.8 4.0 7 6.0 5.9 8 3.5 8.09 4.3 5.0proc means n mean std stderr t prt

13、;var d;run;、MEANS 过程分析变量：dN 均值 标准偏差 标准误差 t 值 Pr > |t|-9 -0.3555556 1.9956898 0.6652299 -0.53 0.60759、上午和下午各测8只家兔的血糖值，结果如下。试比较其血糖值上、下午有无变化。(配对t检验)兔号 上午血糖(mg%) 下午血糖(mg%)1 100 1052 110 1153 120 1284 98 1045 92 1026 138 1347 100 1108 125 130data d;input a x y;d=x-y;cards;1 100 105 2 110 115 3 120 128

14、 4 98 1045 92 102 6 138 134 7 100 110 8 125 130proc means n mean std stderr t prt;var d;run;分析变量：d均值 N 标准偏差 标准误差 t 值 Pr > |t|-5.6250000 8 4.4380659 1.5690932 -3.58 0.008910、在同样饲养管理条件下，三个品种猪的增重如下表，试对三个品种增重差异是否显著进行检验。(单因素各处理重复数相等的方差分析)品种增 重xij(kg)A116121818131115101718A210131191614815138A3118136715

15、9121011data y;input a$ x;cards;a1 16 a1 12 a1 18 a1 18 a1 13 a1 11 a1 15 a1 10 a1 17 a1 18a2 10 a2 13 a2 11 a2 9 a2 16 a2 14 a2 8 a2 15 a2 13 a2 8a3 11 a3 8 a3 13 a3 6 a3 7 a3 15 a3 9 a3 12 a3 10 a3 11proc glm;class a;model x=a;means a/lsd duncan snk;means a/lsd duncan alpha=0.01;run;Sum ofSource DF

16、 Squares Mean Square F Value Pr > FModel 2 110.0666667 55.0333333 6.42 0.0052Error 27 231.3000000 8.5666667Corrected Total 29 341.3666667R-Square Coeff Var Root MSE x Mean0.322429 23.92551 2.926887 12.23333Source DF Anova SS Mean Square F Value Pr > Fa 2 110.0666667 55.0333333 6.42 0.005211、为了

17、比较4种饲料(A)和猪的3个品种(B)，从每个品种随机抽取4头猪(共12头)分别喂以4种不同饲料。随机配置，分栏饲养、位置随机排列。从60日龄起到90日龄的时期内分别测出每头猪的日增重(g),数据如下，试检验饲料及品种间的差异显著性。（两因素交叉分组单个观测值方差分析）4种饲料3个品种猪6090日龄日增重A1A2A3A4B1505545590530B2490515535505B3445515510495data o;input a$ b$ x;cards;a1 b1 505 a1 b2 490 a1 b3 445a2 b1 545 a2 b2 515 a2 b3 515a3 b1 590 a3

18、 b2 535 a3 b3 510a4 b1 530 a4 b2 505 a4 b3 495proc anova;class a b;model x=a b;means a b/lsd duncan snk;means a b/lsd duncan alpha=0.01;run;R-Square Coeff Var Root MSE x Mean0.908835 2.760311 14.21560 515.0000Source DF Anova SS Mean Square F Value Pr > Fa 2 5337.500000 2668.750000 13.21 0.0063b 3

19、 6750.000000 2250.000000 11.13 0.0073Source DF Anova SS Mean Square F Value Pr > Fa 3 6225.000000 2075.000000 12.49 0.0054b 2 6078.333333 3039.166667 18.30 0.002812、为了从3种不同原料和3种不同温度中选择使酒精产量最高的水平组合，设计了试验，每一水平组合重复3次，结果如下表，试进行分析。（两因素交互分组有重复观测值方差分析）用不同原料及不同温度发酵的酒精产量原料温 度BB1(30)B2(35)B3(40)A1414923111

20、22562226A247595043383382218A3483553553847303326data g;input a$ b$ x;cards;a1 b1 41 a1 b1 49 a1 b1 23 a1 b2 11 a1 b2 12 a1 b2 25 a1 b3 6 a1 b3 22 a1 b3 26a2 b1 47 a2 b1 59 a2 b1 50 a2 b2 43 a2 b2 38 a2 b2 33 a2 b3 8 a2 b3 22 a2 b3 18a3 b1 48 a3 b1 35 a3 b1 53 a3 b2 55 a3 b2 38 a3 b2 47 a3 b3 30 a3 b3

21、 33 a3 b3 26proc anova;class a b;model x=a b a*b;means a b a*b/lsd duncan snk;means a b a*b/lsd duncan alpha=0.01;run;Sum ofSource DF Squares Mean Square F Value Pr > FModel 8 4680.518519 585.064815 8.26 0.0001Error 18 1274.666667 70.814815Corrected Total 26 5955.185185R-Square Coeff Var Root MSE

22、 x Mean0.785957 25.30169 8.415154 33.25926Source DF Anova SS Mean Square F Value Pr > Fa 2 1308.074074 654.037037 9.24 0.0017b 2 2545.407407 1272.703704 17.97 <.0001a*b 4 827.037037 206.759259 2.92 0.050413、 对陕西三个秦川牛保种基地县进行秦川牛肉用性能外形调查，划分为优良中下4个等级，试问三个地区秦川牛肉用性能各级构成比差异是否显著。（半定量资料：卡方）地区优良中下甲10106

23、010乙1052010丙55236data v;input a b x;cards;1 1 10 1 2 10 1 3 60 1 4 102 1 10 2 2 5 2 3 20 2 4 103 1 5 3 2 5 3 3 23 3 4 6proc freq;table a*b/chisq;weight x;run;a * b 表的统计量统计量 自由度 值 概率-卡方 6 7.7268 0.2588似然比卡方 6 7.5059 0.2766Mantel-Haenszel 卡方 1 0.0405 0.8406Phi 系数 0.2107列联系数 0.2062Cramer 的 V 0.1490样本大小

24、 = 17414、10头育肥猪的饲料消耗（x）和增重（y）资料如下表（单位：kg），试对增重与饲料消耗进行直线回归分析。 x191167194158200179178174170175y33114224384438373035data m;input x y;cards;191 33 167 11 194 42 158 24 200 38 179 44 178 38 174 37 170 30 175 35proc reg corr;model y=x;run;Parameter Estimates Parameter Standard Variable DF Estimate Er

25、ror t Value Pr > |t| Intercept 1 -47.80836 37.54743 -1.27 0.2387 x 1 0.45357 0.20973 2.16 0.062515、某猪场从10窝大白猪的仔猪中，每窝抽出性别相同、体重接近的仔猪2头，将每窝两头仔猪随机地分配到两个饲料组，进行饲料对比试验，试验时间30天，增重结果见下表。试检验两种饲料喂饲的仔猪平均增重差异是否显著？（同源配对t检验）窝号12345678910饲料10.011.212.110.511.19.810.812.512.09.9饲料9.510.511.89.512.08.89.711.211.09

26、.0data h;input a x y;d=x-y;cards;1 10 9.5 2 11.2 10.5 3 12.1 11.8 4 10.5 9.5 5 11.1 12.0 6 9.8 8.8 7 10.8 9.7 8 12.5 11.2 9 12.0 11.0 10 9.9 9.0 proc means mean n std stderr t prt;var d;run;分析变量：d均值 N 标准偏差 标准误差 t 值 Pr > |t|-0.6900000 10 0.6314885 0.1996942 3.46 0.007216、试对下列资料进行直线相关和回归分析。x3630262

27、3263020192016y0.890.800.740.800.850.680.730.680.800.58data l;input x y;cards;36 0.89 30 0.80 26 0.74 23 0.80 26 0.8530 0.68 20 0.73 19 0.68 20 0.80 16 0.58proc reg corr;model y=x;run;Parameter EstimatesParameter StandardVariable DF Estimate Error t Value Pr > |t|Intercept 1 0.52150 0.10272 5.08 0

28、.0010x 1 0.00949 0.00406 2.34 0.047717、 为了研究不同种类饲料对奶牛产奶量的影响，设置了A、B、C、D、E 5种饲料，用5头奶牛进行试验，试验根据泌乳阶段分为5期，每期4周，试验结果列于下表，试对其进行方差分析。（3因素交互分组单个观测值方差分析拉丁方实验设计） 饲料种类对奶牛产奶量影响的试验结果 (单位：kg)牛 号时 期一二三四五E(300)A(320)B(390)C(390)D(380)IID(420)C(390)E(280)B(370)A(270)IIIB(350)E(360)D(400)A(260)C(400)IVA(280)D(400)C(39

29、0)E(280)B(370)VC(400)B(380)A(350)D(430)E(320)data x;input a$ b$ c$ x;cards;a1 b1 c5 300 a1 b2 c1 320 a1 b3 c2 390 a1 b4 c3 390 a1 b5 c4 380a2 b1 c4 420 a2 b2 c3 390 a2 b3 c5 280 a2 b4 c2 370 a2 b5 c1 270a3 b1 c2 350 a3 b2 c5 360 a3 b3 c4 400 a3 b4 c1 260 a3 b5 c3 400a4 b1 c1 280 a4 b2 c4 400 a4 b3 c

30、3 390 a4 b4 c5 280 a4 b5 c2 370a5 b1 c3 400 a5 b2 c2 380 a5 b3 c1 350 a5 b4 c4 430 a5 b5 c5 320proc anova;class a b c;model x=a b c;means a b c/lsd duncan snk;means a b c/lsd duncan alpha=0.01;run; The ANOVA ProcedureDF Anova SS Mean Square F Value Pr > F a 4 3224.00000 806.00000 1.32 0.3193 b 4

31、2144.00000 536.00000 0.87 0.5071 c 4 50504.00000 12626.00000 20.61 <.000118、用青木香治疗高血压病前后的舒张压(mmHg)资料如下，比较治疗前后舒张压有无显著差异。（自身配对t检验）病例号 治疗前舒张压 治疗后舒张压 1 110 902 146 1203 110 964 115 1165 90 936 133 1017 116 1208 134 1329 110 10810 150 120data z;input a x y;d=x-y;cards;1 110 90 2 146 120 3 110 96 4 115 116 5 90 936 133 101 7 116 120 8 134 132 9 110 1