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

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 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.220.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 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 连续校正卡方 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.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 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 .0001b 3 824.1428571 274.7142857 64.34 |t|-0.7919048 2 0.0491608 0.0347619 22.78 0.02798、某医院对9例慢性苯中毒患者用中草药抗苯一号和耦合剂治疗，治疗后测定白细胞总数（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;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 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)A116121818131115101718A210131191614815138A31181367159121011data 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 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、为了比较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 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 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)A141492311122562226A247595043383382218A3483553553847303326data 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 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 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、某猪场从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.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、试对下列资料进行直线相关和回归分析。x36302623263020192016y0.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.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(390)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 c3 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 2144.00000 536.00000 0.87 0.5071 c 4 50504.00000 12626.00000 20.61 |t|-10 11.8000000 14.2890790 4.5186035 2.61 0.028219、将20名某病患者随机分成两组，分别用甲、乙两药治疗，测得治疗前后的血沉（mm/h）如下表：甲乙两药治疗前后的血沉（mm/h）甲 药 病人号 1 2 3 4 5 6 7