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

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.

2、0 14.7 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 检验）x编 号 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

3、 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.3data 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 均值 标准偏差 标准误差

4、 t 值 Pr |t|-11 1.8454545 1.4807860 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;统计量 自由度

5、 值 概率-卡方 1 1.1429 0.2850似然比卡方 1 1.1497 0.2836连续校正卡方 1 0.7314 0.3924Mantel-Haenszel 卡方 1 1.1371 0.2863Phi 系数 -0.0756列联系数 0.0754Cramer 的 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（） ，试检验该样本平均温度与 是否存在显著差异？（单个因00素 t 检验）data r;input x;y=x-

6、39.5;cards;38.7 39.0 38.9 39.6 39.1 39.8 38.5 39.7 39.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、甲、乙、丙三个奶牛场高产奶牛、低产奶牛、低产奶牛头数统计如下，试问三个奶牛场高、中、低产奶牛的 构成比构成比 是否有差异。 （半定量（等级）资料：只有卡方）场 地 高产奶牛 中产奶牛 低产奶牛甲 32 15 18乙 28 26 10丙 11 38 16data q

7、;input a b x;cards;1 1 32 1 2 15 1 3 182 1 28 2 2 26 2 3 103 1 11 3 2 38 3 3 16proc freq;table a*b/chisq;weight x;run;a * b 表的统计 量统计量 自由度 值 概率-卡方 4 22.8289 0.0001似然比卡方 4 24.6805 Fa 6 629.4285714 104.9047619 24.57 |t|-0.7919048 2 0.0491608 0.0347619 22.78 0.02798、某医院对9例慢性苯中毒患者用中草药抗苯一号和耦合剂治疗，治疗后测定白细胞总

8、数（10 9/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;var d;ru

9、n;、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 4 98 1045

10、 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)A1 16 12 18 18 13 11 15 10 17 18A2 10 13 11 9 16 14 8 15 13 8A3 1

11、1 8 13 6 7 15 9 12 10 11data 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;ru

12、n;Sum ofSource DF 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.00

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

14、a2 b2 515 a2 b3 515a3 b1 590 a3 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 2

15、668.750000 13.21 0.0063b 3 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 次，结果如下表，试进行分析。 （两因素交互分组有重复观测值方差分析）用不同原料及不同温度发酵的酒精产量原料 温 度

16、BB1(30 ) B2(35 ) B3(40 )A1 41 49 23 11 12 25 6 22 26A2 47 59 50 43 38 33 8 22 18A3 48 35 53 55 38 47 30 33 26data 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

17、48 a3 b1 35 a3 b1 53 a3 b2 55 a3 b2 38 a3 b2 47 a3 b3 30 a3 b3 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.814815

18、Corrected Total 26 5955.185185R-Square Coeff Var Root MSE 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 |t|Intercept 1 -47.80836 37.54743 -1.27 0.2387 x 1 0.45357 0.20973 2.16 0.062515、某猪场从

19、10 窝大白猪的仔猪中，每窝抽出性别相同、体重接近的仔猪 2 头，将每窝两头仔猪随机地 分配 到两个饲料组，进行饲料对比试验，试验时间 30 天，增重结果见下表。试检验两种饲料喂饲的仔猪平均增重差异是否显著？（同源配对 t 检验）窝号 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.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

20、 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、试对下列资料进行直线相关和回归分析。x 36 30 26 23 26 30 20 19 20 16y 0.89 0.80 0.74 0.80 0.85 0.68 0.73 0.68 0.80 0.

21、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.0010x 1 0.00949 0.00406 2.34 0.047717、 为了研究不同种类饲料对奶牛产奶量的

22、影响，设置了 A、B、C 、D 、E 5 种饲料，用5 头奶牛进行试验，试验根据泌乳阶段分为 5 期，每期 4 周，试验结果列于下表，试对其进行方差分析。 （3 因素交互分组单个观测值方差分析 拉丁方实验设计）饲料种类对奶牛产奶量影响的试验结果 (单位：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)

23、 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 c3 390 a4 b4 c5 28

24、0 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 Fa 4 3224.00000 806.00000 1.32 0.3193b 4 2144.00000 536.00000 0.87 0.5071 c 4 50504.00000 12626.00000 20.61 .000118、用青木香治疗高血压病前后的舒张压(mmHg)资料如下，比较治疗前后舒张压有无显著差异。（自身配对 t 检验）