《生物統(tǒng)計附試驗設(shè)計》實習習題

1、1、應(yīng)用克矽平治療矽肺患者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射線照射前后之體重數(shù)據(jù)見下表（單位：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只（經(jīng)注射后患病的10只，不患病的90只），對照組（注射原疫苗組）100只（經(jīng)注射后患病的15只，不患病的85只），試問新舊疫苗的免疫力是否有差異。（次數(shù)資料：卡方） 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; 統(tǒng)計量 自由度 值 概率 - 卡方 1 1.1429 0.2850 似然比卡方 1 1.1497 0.2836 連續(xù)

5、校正卡方 1 0.7314 0.3924 Mantel-Haenszel 卡方 1 1.1371 0.2863 Phi 系數(shù) -0.0756 列聯(lián)系數(shù) 0.0754 Cramer 的 V -0.07564、隨機抽測了10只兔的直腸溫度，其數(shù)據(jù)為：38.7、39.0、38.9、39.6、39.1、39.8、38.5、39.7、39.2、38.4（），已知該品種兔直腸溫度的總體平均數(shù)=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、甲、乙、丙三個奶牛場高產(chǎn)奶牛、低產(chǎn)奶牛、低產(chǎn)奶牛頭數(shù)統(tǒng)計如下，試問三個奶牛場高、中、低產(chǎn)奶牛的構(gòu)成比是否有差異。（半定量（等級）資料：只有卡方）場 地高產(chǎn)奶牛中產(chǎn)奶牛低產(chǎn)奶牛甲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 表的統(tǒng)計量統(tǒng)計量 自由度 值 概率-卡方 4 22.8289 0.0001似然比卡方 4 24.6805 <.0001Mantel-Haenszel 卡方 1 4.8178 0.0282Phi 系數(shù) 0.3430列聯(lián)系數(shù) 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例，問兩藥的治愈率是否有顯著差異？ （次數(shù)資料：卡方）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 表的統(tǒng)計量

11、統(tǒng)計量 自由度 值 概率 - 卡方 1 1.0665 0.3017似然比卡方 1 1.0675 0.3015連續(xù)校正卡方 1 0.6849 0.4079Mantel-Haenszel 卡方 1 1.0592 0.3034Phi 系數(shù) -0.0858列聯(lián)系數(shù) 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、某醫(yī)院對9例慢性苯中毒患者用中草藥抗苯一號和耦合劑治療

12、，治療后測定白細胞總數(shù)（109/L)，資料如下。問兩種藥是否對患者的白細胞總數(shù)有影響？（計數(shù)資料：配對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只家兔的血糖值，結(jié)果如下。試比較其血糖值上、下午有無變化。(配對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、在同樣飼養(yǎng)管理條件下，三個品種豬的增重如下表，試對三個品種增重差異是否顯著進行檢驗。(單因素各處理重復數(shù)相等的方差分析)品種增 重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種不同飼料。隨機配置，分欄飼養(yǎng)、位置隨機排列。從60日齡起到90日齡的時期內(nèi)分別測出每頭豬的日增重(g),數(shù)據(jù)如下，試檢驗飼料及品種間的差異顯著性。（兩因素交叉分組單個觀測值方差分析）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種不同溫度中選擇使酒精產(chǎn)量最高的水平組合，設(shè)計了試驗，每一水平組合重復3次，結(jié)果如下表，試進行分析。（兩因素交互分組有重復觀測值方差分析）用不同原料及不同溫度發(fā)酵的酒精產(chǎn)量原料溫 度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、 對陜西三個秦川牛保種基地縣進行秦川牛肉用性能外形調(diào)查，劃分為優(yōu)良中下4個等級，試問三個地區(qū)秦川牛肉用性能各級構(gòu)成比差異是否顯著。（半定量資料：卡方）地區(qū)優(yōu)良中下甲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 表的統(tǒng)計量統(tǒng)計量 自由度 值 概率-卡方 6 7.7268 0.2588似然比卡方 6 7.5059 0.2766Mantel-Haenszel 卡方 1 0.0405 0.8406Phi 系數(shù) 0.2107列聯(lián)系數(shù) 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天，增重結(jié)果見下表。試檢驗兩種飼料喂飼的仔豬平均增重差異是否顯著？（同源配對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、試對下列資料進行直線相關(guān)和回歸分析。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、 為了研究不同種類飼料對奶牛產(chǎn)奶量的影響，設(shè)置了A、B、C、D、E 5種飼料，用5頭奶牛進行試驗，試驗根據(jù)泌乳階段分為5期，每期4周，試驗結(jié)果列于下表，試對其進行方差分析。（3因素交互分組單個觀測值方差分析拉丁方實驗設(shè)計） 飼料種類對奶牛產(chǎn)奶量影響的試驗結(jié)果 (單位：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