Đề cương thiết kế thí nghiệm chăn nuôi

PHẦN I: ƯỚC LƯỢNG VÀ KIỂM ĐỊNH GIÁ TRỊ TRUNG BÌNH: A ƯỚC LƯỢNG KIỂM ĐỊNH GIÁ TRỊ TRUNG BÌNH: Ví dụ M-1.3: ( Kiểm định số C có phải thuộc quần thể hay khơng?) Các bước tiến hành: Bước 1: Tóm tắt trình bày liệu: Bước 2: Giả thiết H0 đối thiết H1: Bước 3: Kiểm tra điều kiện: Bước 4: Tính xác suất P: Bước 5: So sánh P với α rút kết luận: BÀI LÀM Bước 1: Tóm tắt trình bày liệu: Vào Stat/Basic Statistics/ Display Descriptive Statistics… Variable N(n) P 36 N* Mean( X 599.19 ) SE Mean(SE) 3.11 StDev(S) 18.66 Variance 348.05 CoefVar(Cv) 3.11 S 18,66 Minimum 559.00 Maximum 636.00 Cv(%) 3,11 Tăng trọng trung bình giống lợn Landrace: Chỉ tiêu Tăng trọng Đơn vị tính gram/ngày n 36 X ± SE 599,19 ± 3,11 Bước 2: Giả thiết H0 đối thiết H1: Ho: Tăng trọng trung bình tồn đàn lợn trại 607 gram/ngày µ = 607 gram/ngày H1: Tăng trọng trung bình tồn đàn lợn trại khác 607 gram/ngày µ ≠ 607 gram/ngày Bước 3: Kiểm tra điều kiện: Vào Stat/Basic Statistics/ Normality Test… Kiểm tra phải tuân theo phân phối chuẩn: P-Value = 0,997 > 0,05 Số liệu có phân phối chuẩn Điều kiện: P-Value ≥ 0,05 tuân theo phân phối chuẩn P-Value < 0,05 khơng tn theo phân phối chuẩn Bước 4: Tính xác suất P: Trường hợp biết độ lệch chuẩn σ ( Standard Deviation) µ ( Mean cần kiểm định ) = 607 Ví dụ M-1.3 có σ = 21,75 gram Sử dụng phép thử Z ( 1Z): Stat/Basic Statistics/ 1-Sample Z… Nhập giá trị σ = 21,75 gram vào Standard Deviation Nhập giá trị µ = 607 gram vào ô Hypothesized mean One-Sample Z: P Test of mu = 607 vs not = 607 The assumed standard deviation = 21.75 Variable P N 36 Mean 599.19 StDev 18.66 SE Mean 3.63 95% CI (592.09, 606.30) Z -2.15 P = 0.031 Bước 5: So sánh P với α rút kết luận: P = 0.031 < 0.05 nên bác bỏ giả thiết H0 chấp nhận giả thiết H1 P 0.031 Kết luận: Tăng trọng trung bình tồn đàn lợn trại khác 607 gram/ngày ( µ ≠ 607 gram/ngày Khoảng ước lượng: (592.09, 606.30) Điều kiện: P > α ( =0,05) chấp nhận H0 P < α ( =0,05) bác bỏ H0, chấp nhận H1 Trường hợp chưa biết độ lệch chuẩn σ ( Standard Deviation) biết µ ( Mean cần kiểm định ) = 607 Bước 4: Tính xác suất P: Sử dụng phép thử T (1T): Stat/Basic Statistics/ 1-Sample t… Nhập giá trị µ = 607 gram vào Hypothesized mean One-Sample T: P Test of mu = 607 vs not = 607 Variable P N 36 Mean 599.19 StDev 18.66 SE Mean 3.11 95% CI (592.88, 605.51) T -2.51 P 0.017 P = 0,017 Bước 5: So sánh P với α rút kết luận: P = 0,017 < 0,05 nên bác bỏ giả thiết H0 chấp nhận giả thiết H1 Kết luận: Tăng trọng trung bình tồn đàn lợn trại khác 607 gram/ngày ( µ ≠ 607 gram/ngày) Khoảng ước lượng: (592.88, 605.51) B SO SÁNH GIÁ TRỊ TRUNG BÌNH CỦA BIẾN CHUẨN: I SO SÁNH GIÁ TRỊ TRUNG BÌNH CỦA BIẾN CHUẨN KHI LẤY MẪU ĐỘC LẬP: khơng có tương quan với Phương sai đồng nhất: VÍ DỤ: M-1.4: So sánh µ biết phương sai (σ ) đồng Các bước tiến hành: Bước 1: Tóm tắt trình bày liệu: Bước 2: Giả thiết H0 đối thiết H1: Bước 3: Kiểm tra điều kiện: Bước 4: Tính xác suất P: Bước 5: So sánh P với α rút kết luận: BÀI LÀM Bước 1: Tóm tắt trình bày liệu: Vào Stat/Basic Statistics/ Display Descriptive Statistics… Descriptive Statistics: KL Variable KL GIONG N 12 15 N* 0 Variable KL GIONG Mean 196.17 153.70 Maximum 221.10 181.80 Khối lượng hai bị giống: Cơng thức TN Đơn vị tính Giống bị Kg Giống bị Kg SE Mean 3.06 3.18 n 12 15 Bước 2: Giả thiết H0 đối thiết H1: StDev 10.62 12.30 Variance 112.70 151.32 CoefVar 5.41 8.00 X ± SE 196,17 ± 3,06 153,70 ± 3,18 Minimum 180.30 135.30 S 10,62 12,30 Cv(%) 5,41 8,00 Ho: Khối lượng hai giống bò µGiống = µGiống H1: Khối lượng hai giống bị khơng µGiống ≠ µGiống Bước 3: Kiểm tra điều kiện: Vào Stat/Basic Statistics/ Normality Test… a Kiểm tra phải tuân theo phân phối chuẩn: Giống 1: P-Value =0,530 > 0,05 Số liệu có phân phối chuẩn Giống 2: P-Value =0.407 > 0,05 Số liệu có phân phối chuẩn b Kiểm tra phương sai đồng nhất: Vào: Stat/ Basic Statistics/ Variances… Test and CI for Two Variances: GIONG1, GIONG2 Method Null hypothesis Alternative hypothesis Significance level Sigma(GIONG1) / Sigma(GIONG2) = Sigma(GIONG1) / Sigma(GIONG2) not = Alpha = 0.05 Statistics Variable GIONG1 GIONG2 N 12 15 StDev 10.616 12.301 Variance 112.700 151.324 Ratio of standard deviations = 0.863 Ratio of variances = 0.745 95% Confidence Intervals Distribution of Data Normal Continuous CI for Variance Ratio (0.241, 2.502) (0.152, 2.749) CI for StDev Ratio (0.491, 1.582) (0.390, 1.658) Tests Method F Test (normal) Levene's Test (any continuous) DF1 11 DF2 14 25 Test Statistic 0.74 0.46 P-Value 0.631 0.503 P-Value F Test: P-Value = 0,631 > 0,05 nên chấp nhận H0, hai phương sai đồng Bước 4: Tính xác suất P: Sử dụng phép thử T( 2t): Stat/ Basic Statistics/ 2-Sample t… Chú ý: - Nếu hai phương sai đồng (V) vào Assume equal variances - Nếu hai phương sai đồng khơng tích (V) vào ô Assume equal variances Two-sample T for GIONG1 vs GIONG2 GIONG1 GIONG2 N 12 15 Mean 196.2 153.7 StDev 10.6 12.3 SE Mean 3.1 3.2 Difference = mu (GIONG1) - mu (GIONG2) Estimate for difference: 42.48 95% CI for difference: (33.23, 51.72) T-Test of difference = (vs not =): T-Value = 9.46 Both use Pooled StDev = 11.5901 P-Value = 0.000 DF = 25 P-Value = 0,000 < 0,05( α) Bước 5: So sánh P với α rút kết luận: P-Value = 0,000 < 0,05( α) nên bác bỏ giả thiết H0, chấp nhận giả thiết H1 Kết luận: Khối lượng hai giống bị khơng (µGiống ≠ µGiống 2) ( P 0,05 Số liệu có phân phối chuẩn Giống 2: P-Value =0.407 > 0,05 Số liệu có phân phối chuẩn b Kiểm tra phương sai đồng nhất: Vào: Stat/ Basic Statistics/ Variances… Test and CI for Two Variances: GIONG1, GIONG2 Variance 112.70 151.32 CoefVar 5.41 8.00 X ± SE 196,17 ± 3,06 153,70 ± 3,18 Minimum 180.30 135.30 S 10,62 12,30 Cv(%) 5,41 8,00 Method Null hypothesis Alternative hypothesis Significance level Sigma(GIONG1) / Sigma(GIONG2) = Sigma(GIONG1) / Sigma(GIONG2) not = Alpha = 0.05 Statistics Variable GIONG1 GIONG2 N 12 15 StDev 10.616 12.301 Variance 112.700 151.324 Ratio of standard deviations = 0.863 Ratio of variances = 0.745 95% Confidence Intervals Distribution of Data Normal Continuous CI for Variance Ratio (0.241, 2.502) (0.152, 2.749) CI for StDev Ratio (0.491, 1.582) (0.390, 1.658) Tests Method F Test (normal) Levene's Test (any continuous) DF1 11 DF2 14 25 Test Statistic 0.74 0.46 P-Value 0.631 0.503 P-Value F Test: P-Value = 0,631 > 0,05 nên H0 chấp nhận, hai phương sai đồng Bước 4: Tính xác suất P: Sử dụng phép thử T( 2t): Stat/ Basic Statistics/ 2-Sample t… Chú ý: - Nếu hai phương sai không đồng khơng tích (V) vào Assume equal variances Two-Sample T-Test and CI: GIONG1, GIONG2 Two-sample T for GIONG1 vs GIONG2 GIONG1 GIONG2 N 12 15 Mean 196.2 153.7 StDev 10.6 12.3 SE Mean 3.1 3.2 Difference = mu (GIONG1) - mu (GIONG2) Estimate for difference: 42.48 95% CI for difference: (33.37, 51.58) T-Test of difference = (vs not =): T-Value = 9.62 P-Value = 0.000 DF = 24 P-Value = 0,000 < 0,05( α) Bước 5: So sánh P với α rút kết luận: P-Value = 0,000 < 0,05( α) nên bác bỏ giả thiết H0, chấp nhận giả thiết H1 Kết luận: Khối lượng hai giống bị khơng (µGiống ≠ µGiống 2) ( P 0,05 Số liệu có phân phối chuẩn B4: Tìm P-value kết luận Sử dụng phép thử T cặp đôi( t-t): Vào Stat/ Basic Statistics/ Paired t… Paired T-Test and CI: A, B Paired T for A - B A B Difference N 10 10 10 Mean 42.80 38.20 4.600 StDev 3.82 4.13 1.955 SE Mean 1.21 1.31 0.618 95% CI for mean difference: (3.201, 5.999) T-Test of mean difference = (vs not = 0): T-Value = 7.44 P-Value = 0.000 P-Value = 0,000 < 0,05 nên bác bỏ giả thiết H0, chấp nhận H1 Kết luận: Tăng trọng trung bình hai cách chăm sóc khác (µcách A ≠ µcách B) PHẦN II: PHÂN TÍCH PHƯƠNG SAI I THÍ NGHIỆM YẾU TỐ HỒN TỒN NGẪU NHIÊN: Bước 1: Tóm tắt trình bày liệu Bước 2: Giả thiết H0 đối thiết H1 Bước 3: Kiểm tra điều kiện Bước 4: Tính xác suất P Bước 5: So sánh P với α rút kết luận BÀI LÀM: Bước 1: Tóm tắt trình bày liệu Vào Stat/Basic Statistics/ Display Descriptive Statistics… Cv(%) 8,93 10,81 Descriptive Statistics: KL Variable KL CT A B C D E N 4 4 N* 0 0 Mean 0.8875 0.4250 0.7625 0.9375 0.9500 Variable KL CT A B C D E SE Mean 0.0239 0.0104 0.0473 0.0239 0.0204 StDev 0.0479 0.0208 0.0946 0.0479 0.0408 Variance 0.0023 0.0004 0.0090 0.0023 0.0017 CoefVar 5.39 4.90 12.41 5.11 4.30 Minimum 0.8500 0.4000 0.7000 0.9000 0.9000 Maximum 0.9500 0.4500 0.9000 1.0000 1.0000 Cơng thức TN Đơn vị tính n S Cv(%) X ± SE Công thức nuôi A Kg 0,8875 ± 0,0239 0,0479 5,39 Công thức nuôi B Kg 0,4250 ± 0,0104 0.0208 4,90 Công thức nuôi C Kg 0,7625 ± 0,0473 0.0946 12,41 Công thức nuôi D Kg 0,9375 ± 0,0239 0,0479 5,11 Công thức nuôi E Kg 0,9500 ± 0,0204 0,0408 4,30 Bước 2: Giả thiết H0 đối thiết H1 Ho: Tăng trọng cá công thức nuôi µA =µ B= µC =µ D=µE H1: Tăng trọng cá công thức nuôi khác µA ≠µ B≠ µC ≠µ D≠µE Bước 3: Kiểm tra điều kiện a Phân phối chuẩn: - Tìm cột RESI( phần dư) bằng: Vào Stat/ ANOVA/ One-Ways… tích vào cột Store residuals - Vào Stat/Basic Statistics/ Normality Test… kiểm tra phân phối chuẩn cột RESI P-Value = 0,159 > 0,05 nên phân phối chuẩn b Phương sai đồng nhất: Vào Stat/ ANOVA/ Test for Equal Variances… P-Value Levene’s Test: P-Value = 0,539 > 0,05 nên chấp nhận H0, phương sai đồng Bước 4: Tính xác suất P Vào Stat/ ANOVA/ One-Ways… Test for Equal Variances: KL versus CT One-way ANOVA: KL versus CT Source CT Error Total DF 15 19 SS 0.76325 0.04692 0.81017 S = 0.05593 Level A B C D E N 4 4 MS 0.19081 0.00313 R-Sq = 94.21% Mean 0.8875 0.4250 0.7625 0.9375 0.9500 StDev 0.0479 0.0208 0.0946 0.0479 0.0408 Pooled StDev = 0.0559 F 60.99 P 0.000 R-Sq(adj) = 92.66% Individual 95% CIs For Mean Based on Pooled StDev + -+ -+ -+ ( * ) ( * ) ( * ) ( * ) (-* ) + -+ -+ -+ 0.40 0.60 0.80 1.00 Nguồn biến động CT Sai số Tổng số Bậc tự ( DF) Tổng bình phương(SS) Trung bình bình F phương(MS) P 0.76325 0.19081 0.000 15 0.04692 0.00313 19 0.81017 60.99 P = 0,000 Bước 5: So sánh P với α rút kết luận P= 0,000 < 0,05 nên bác bỏ H0 chấp nhận H1 Kết luận: Tăng trọng cá công thức nuôi khác nhau( µA ≠µ B≠ µC ≠µ D≠µE ) So sánh cặp: bác bỏ H0 chấp nhận H1: Minitab 16: Vào Stat/ ANOVA/ One-Ways… tích vào Comparison… sau tích vào Turkey’s, family error rate Grouping Information Using Tukey Method CT E D A C B N 4 4 Mean 0.95000 0.93750 0.88750 0.76250 0.42500 Grouping A A A B C Means that not share a letter are significantly different Tukey 95% Simultaneous Confidence Intervals All Pairwise Comparisons among Levels of CT Individual confidence level = 99.25% CT = A subtracted from: CT B C D E Lower -0.58471 -0.24721 -0.07221 -0.05971 Center -0.46250 -0.12500 0.05000 0.06250 Upper -0.34029 -0.00279 0.17221 0.18471 + -+ -+ -+( -* ) ( * -) ( * -) ( -* ) + -+ -+ -+-0.35 0.00 0.35 0.70 CT = B subtracted from: CT C D E Lower 0.21529 0.39029 0.40279 Center 0.33750 0.51250 0.52500 Upper 0.45971 0.63471 0.64721 + -+ -+ -+( -* ) ( -* ) ( * ) + -+ -+ -+-0.35 0.00 0.35 0.70 CT = C subtracted from: CT D E Lower 0.05279 0.06529 Center 0.17500 0.18750 Upper 0.29721 0.30971 + -+ -+ -+( * ) ( * -) + -+ -+ -+-0.35 0.00 0.35 0.70 CT = D subtracted from: CT Lower Center Upper + -+ -+ -+- E -0.10971 0.01250 0.13471 ( * -) + -+ -+ -+-0.35 0.00 0.35 0.70 - So sánh cơng thức: + Xắp xếp giá trị trung bình theo thứ tự giảm dần: E = 0.95000 | a D = 0.93750 | | a A = 0.88750 | | | a C = 0.76250 | b B = 0.42500 | c + Bằng nối nét liền, khác dừng lại sang phải xuống + Đường nối dài lấy, đường nối ngắn xóa bỏ + Bằng dấu Lower Upper khác dấu, Khác dấu Kết : Cơng thức: Trung bình Cơng thức A: 0,8875a Cơng thức B: 0,425c Công thức C: 0,7625b Công thức D: 0,9375a Cơng thức E: 0,95a II THÍ NGHIỆM YẾU TỐ KHỐI NGẪU NHIÊN ĐẦY ĐỦ: Ví dụ: M-1.7 Bước 1: Tóm tắt trình bày liệu Bước 2: Giả thiết H0 đối thiết H1 Bước 3: Kiểm tra điều kiện Bước 4: Tính xác suất P Bước 5: So sánh P với α rút kết luận BÀI LÀM: Bước 1: Tóm tắt trình bày liệu Vào Stat/Basic Statistics/ Display Descriptive Statistics… Descriptive Statistics: TB Variable TB THUOC A B C D N 5 5 N* 0 0 Mean 6.420 5.720 6.060 5.660 Variable TB THUOC A B C D SE Mean 0.271 0.297 0.340 0.273 StDev 0.606 0.665 0.760 0.611 Variance 0.367 0.442 0.578 0.373 CoefVar 9.44 11.62 12.55 10.79 Maximum 7.100 6.700 7.100 6.700 Số lượng tế bào lympho chuột sử dụng loại thuốc : Công thức TN Đơn vị tính n X ± SE -3 Thuốc A X1000 tb mm máu 6.420 ± 0.271 -3 Thuốc B X1000 tb mm máu 5.720 ±0.297 Thuốc C X1000 tb mm-3 máu 6.060 ± 0.340 -3 Thuốc D X1000 tb mm máu 5.660 ± 0.273 Bước 2: Giả thiết H0 đối thiết H1 Ho: Ảnh hưởng thuốc đến tế bào lympho µA =µ B= µC =µ D Minimum 5.600 5.100 5.000 5.200 S Cv(%) 0.606 9.44 0.665 11.62 0.760 12.55 0.611 10.79 H1: Ảnh hưởng thuốc đến tế bào lympho khác µA ≠µ B≠ µC ≠µ D Bước 3: Kiểm tra điều kiện a Phân phối chuẩn: - Tìm cột RESI( phần dư) bằng: Vào Stat/ ANOVA/ Two-Ways… tích vào cột Store residuals - Vào Stat/Basic Statistics/ Normality Test… kiểm tra phân phối chuẩn cột RESI P-Value = 0,104 > 0,05 nên số liệu có phân phối chuẩn b Phương sai đồng nhất: Vào Stat/ ANOVA/ Test for Equal Variances… P-Value Levene’s Test: P-Value = 0,977 > 0,05 nên chấp nhận H0, phương sai đồng Bước 4: Tính xác suất P Vào Stat/ ANOVA/ Two-Ways… Two-way ANOVA: TB versus THUOC, LUA Source THUOC LUA Error Total DF 12 19 SS 1.8455 6.4030 0.6370 8.8855 S = 0.2304 MS 0.61517 1.60075 0.05308 R-Sq = 92.83% F 11.59 30.16 P 0.001 0.000 R-Sq(adj) = 88.65% Nguồn biến Bậc tự ( DF) Tổng bình động phương(SS) 1.8455 THUOC 6.4030 LUA 12 0.6370 Sai số 19 8.8855 Tổng số Cách 2: Sử dụng (GLM) Vào Stat/ ANOVA/ General Linear Model… Trung bình bình F phương(MS) P 0.61517 11.59 0.001 1.60075 30.16 0.000 0.05308 Test for Equal Variances: TB versus THUOC General Linear Model: TB versus THUOC, LUA Factor THUOC LUA Type fixed fixed Levels Values A, B, C, D 1, 2, 3, 4, Analysis of Variance for TB, using Adjusted SS for Tests Source THUOC LUA Error Total DF 12 19 S = 0.230398 Seq SS 1.8455 6.4030 0.6370 8.8855 Adj SS 1.8455 6.4030 0.6370 R-Sq = 92.83% Adj MS 0.6152 1.6008 0.0531 F 11.59 30.16 P 0.001 0.000 R-Sq(adj) = 88.65% P = 0,001 Bước 5: So sánh P với α rút kết luận P= 0,001 < 0,05 nên bác bỏ H0 chấp nhận H1 Kết luận: Ảnh hưởng thuốc đến tế bào lympho khác (µA ≠µ B≠ µC ≠µ D) So sánh cặp: bác bỏ H0 chấp nhận H1: Minitab 16: Vào Stat/ ANOVA/ One-Ways… tích vào Comparison… sau tích vào Turkey’s, family error rate General Linear Model: TB versus THUOC, LUA Factor THUOC LUA Type fixed fixed Levels Values A, B, C, D 1, 2, 3, 4, Analysis of Variance for TB, using Adjusted SS for Tests Source THUOC LUA Error Total DF 12 19 Seq SS 1.8455 6.4030 0.6370 8.8855 S = 0.230398 Adj SS 1.8455 6.4030 0.6370 Adj MS 0.6152 1.6008 0.0531 R-Sq = 92.83% F 11.59 30.16 P 0.001 0.000 R-Sq(adj) = 88.65% Grouping Information Using Tukey Method and 95.0% Confidence THUOC A C B D N 5 5 Mean 6.4 6.1 5.7 5.7 Grouping A A B B B Means that not share a letter are significantly different Tukey 95.0% Simultaneous Confidence Intervals Response Variable TB All Pairwise Comparisons among Levels of THUOC THUOC = A subtracted from: THUOC B C D Lower -1.133 -0.793 -1.193 THUOC = B THUOC C D THUOC D Upper -0.2672 0.0728 -0.3272 + -+ -+ -+ -( * -) ( * ) ( * -) + -+ -+ -+ 1.20 -0.60 0.00 0.60 subtracted from: Lower -0.0928 -0.4928 THUOC = C Center -0.7000 -0.3600 -0.7600 Center 0.34000 -0.06000 Upper 0.7728 0.3728 + -+ -+ -+ -( -* ) ( * ) + -+ -+ -+ 1.20 -0.60 0.00 0.60 subtracted from: Lower -0.8328 Center -0.4000 Upper 0.03276 + -+ -+ -+ -( * -) + -+ -+ -+ 1.20 -0.60 0.00 0.60 - So sánh công thức: + Xắp xếp giá trị trung bình theo thứ tự giảm dần: + Bằng nối nét liền, khác dừng lại sang phải xuống + Đường nối dài lấy, đường nối ngắn xóa bỏ + Bằng dấu Lower Upper khác dấu, Khác dấu Kết : Công thức: Trung bình Cơng thức A: 6,4a Cơng thức B: 5,7b Công thức C: 6,1 ab Công thức D: 5,7b III THÍ NGHIỆM KIỂU Ơ VNG LA TINH: Ví dụ: M-1.9a Bước 1: Tóm tắt trình bày liệu Bước 2: Giả thiết H0 đối thiết H1 Bước 3: Kiểm tra điều kiện Bước 4: Tính xác suất P Bước 5: So sánh P với α rút kết luận BÀI LÀM: Bước 1: Tóm tắt trình bày liệu: Vào Stat/Basic Statistics/ Display Descriptive Statistics… Descriptive Statistics: KLCO Variable KLCO TA A B C D N 4 4 N* 0 0 Mean 11.500 11.075 10.925 9.225 Variable KLCO TA A B C D SE Mean 0.445 0.373 0.256 0.304 StDev 0.891 0.746 0.512 0.608 Variance 0.793 0.556 0.263 0.369 CoefVar 7.75 6.73 4.69 6.59 Minimum 10.800 10.000 10.200 8.500 Maximum 12.800 11.700 11.400 9.900 Khối lượng cỏ khơ trung bình 10 ngày thí nghiệm là: Cơng thức TN Thức ăn A Thức ăn B Thức ăn C Thức ăn D Đơn vị TN Kg n X ± SE S Cv(%) Bước 2: Giả thiết H0 đối thiết H1 Ho: Ảnh hưởng bốn loại thức ăn bổ sung đến thu nhận cỏ khơ bê µA =µ B= µC =µ D H1: Ảnh hưởng bốn loại thức ăn bổ sung đến thu nhận cỏ khô bê µA ≠µ B≠ µC ≠µ D Bước 3: Kiểm tra điều kiện a Phân phối chuẩn: - Tìm cột RESI( phần dư) bằng: Vào Stat/ ANOVA/ General Linear Model…dí vào Storage… tích vào cột Residuals - Vào Stat/Basic Statistics/ Normality Test… kiểm tra phân phối chuẩn cột RESI P-Value = 0,760 > 0,05 nên số liệu có phân phối chuẩn b Phương sai đồng nhất: Vào Stat/ ANOVA/ Test for Equal Variances… P-Value Levene’s Test: P-Value = 0,936 > 0,05 nên chấp nhận H0, phương sai đồng Bước 4: Tính xác suất P Stat/ ANOVA/ General Linear Model… General Linear Model: KLCO versus GD, BE, TA Factor GD BE TA Type fixed fixed fixed Levels 4 Values 1, 2, 3, 1, 2, 3, A, B, C, D Analysis of Variance for KLCO, using Adjusted SS for Tests Source GD BE TA Error Total DF 3 15 S = 0.380515 Seq SS 1.4819 3.5919 12.0219 0.8688 17.9644 Adj SS 1.4819 3.5919 12.0219 0.8688 Adj MS 0.4940 1.1973 4.0073 0.1448 R-Sq = 95.16% F 3.41 8.27 27.68 P 0.094 0.015 0.001 R-Sq(adj) = 87.91% Nguồn biến động Bậc tự ( DF) Tổng bình phương(SS) Trung bình bình F phương(MS) P GD 3.41 0.094 BE 8.27 0.015 TA 27.68 0.001 Sai số Tổng số 15 P = 0,001 Bước 5: So sánh P với α rút kết luận P= 0,001 < 0,05 nên bác bỏ H0 chấp nhận H1 Kết luận: Ảnh hưởng bốn loại thức ăn bổ sung đến thu nhận cỏ khơ bê (µA ≠µ B≠ µC ≠µ D) So sánh cặp: bác bỏ H0 chấp nhận H1: Minitab 16: Vào Stat/ ANOVA/ General Linear Model… dí vào Comparison… sau tích vào Turkey Unusual Observations for KLCO Obs 11 KLCO 12.8000 Fit 12.2875 SE Fit 0.3008 Residual 0.5125 St Resid 2.20 R R denotes an observation with a large standardized residual Grouping Information Using Tukey Method and 95.0% Confidence TA A B C D N 4 4 Mean 11.5 11.1 10.9 9.2 Grouping A A A B Means that not share a letter are significantly different Tukey 95.0% Simultaneous Confidence Intervals Response Variable KLCO All Pairwise Comparisons among Levels of TA TA = A subtracted from: TA B C D Lower -1.357 -1.507 -3.207 TA = B Center -0.425 -0.575 -2.275 Upper 0.507 0.357 -1.343 subtracted from: -+ -+ -+ ( * -) ( -* -) ( -* -) -+ -+ -+ -2.4 -1.2 0.0 TA C D Lower -1.082 -2.782 TA = C TA D Center -0.150 -1.850 Upper 0.7823 -0.9177 -+ -+ -+ ( -* -) ( -* ) -+ -+ -+ -2.4 -1.2 0.0 subtracted from: Lower -2.632 Center -1.700 Upper -0.7677 -+ -+ -+ ( -* -) -+ -+ -+ -2.4 -1.2 0.0 Tukey Simultaneous Tests Response Variable KLCO All Pairwise Comparisons among Levels of TA TA = A subtracted from: TA B C D Difference of Means -0.425 -0.575 -2.275 TA = B TA C D TA D T-Value -1.580 -2.137 -8.455 Adjusted P-Value 0.4536 0.2427 0.0006 T-Value -0.557 -6.876 Adjusted P-Value 0.9411 0.0019 T-Value -6.318 Adjusted P-Value 0.0030 subtracted from: Difference of Means -0.150 -1.850 TA = C SE of Difference 0.2691 0.2691 0.2691 SE of Difference 0.2691 0.2691 subtracted from: Difference of Means -1.700 SE of Difference 0.2691 - So sánh cơng thức: + Xắp xếp giá trị trung bình theo thứ tự giảm dần: + Bằng nối nét liền, khác dừng lại sang phải xuống + Đường nối dài lấy, đường nối ngắn xóa bỏ + Bằng dấu Lower Upper khác dấu, Khác dấu Kết : Cơng thức: Trung bình Cơng thức A: 11,5a Công thức B: 11,1b Công thức C: 10,9ab Công thức D: 9,2b PHẦN III: TƯƠNG QUAN HỒI QUY VÀ BẢN TƯƠNG QUAN I HỆ SỐ TƯƠNG QUAN: Ví dụ M-1.11: Vào Stat/ Basis Statistics/ Correlation… Correlations: KL, DL, DN DL DN KL 0.897 0.000 DL 0.905 0.000 0.648 0.001 Cell Contents: Pearson correlation P-Value Nhận xét: - Hệ số tương quan Khối lượng đường kính lớn 0,897 - Hệ số tương quan Khối lượng đường kính bé 0,905 - Hệ số tương quan đường kính lớn đường kính bé 0,648 - Xác suất hệ số tương quan < 0,05 Kết luận: Mối quan hệ tiêu khác Chú ý: H0: ρ = khơng tồn quan hệ tuyến tính H1: ρ ≠ có quan hệ tuyến tính Với hệ số tương quan: R R = ko có quan hệ tuyến tính R ≠ có quan hệ tuyến tính R > tương quan thuận R < tương quan nghịch đảo R = tương quan dương lý tưởng R = -1 tương quan âm lý tưởng | R| < 0,4 có quan hệ, mức quan hệ yếu 0,4 ≤ | R| ≤ 0,8 quan hệ trung bình | R| > 0,8 quan hệ chặt II PHƯƠNG TRÌNH HỒI QUY TUYẾN TÍNH: a = P-Value > α chấp nhận H0  b = a ≠ P-Value < α chấp nhận H1  tồn a, b nên hệ số a, b có ý nghĩa thiết kế b ≠ Vào Stat/ Regression/ Regression… Hồi quy đơn biến: y = ax + b Ví dụ: KL = a0 DL + b1 KL = a0 DN + b1 a Phần Response đánh KL Phần Predictor đánh DL Regression Analysis: KL versus DL The regression equation is KL = - 53.7 + 2.04 DL Predictor Constant DL S = 2.69651 Coef -53.67 2.0379 SE Coef 12.78 0.2250 R-Sq = 80.4% T -4.20 9.06 P 0.000 0.000 R-Sq(adj) = 79.4% Analysis of Variance Source Regression Residual Error Total DF 20 21 SS 596.60 145.42 742.02 MS 596.60 7.27 F 82.05 P 0.000 Unusual Observations Obs DL 57.1 58.2 KL 67.900 59.000 Fit 62.629 64.871 SE Fit 0.579 0.658 Residual 5.271 -5.871 St Resid 2.00R -2.25R R denotes an observation with a large standardized residual - Phương trình hồi quy ước tính khối lượng (Y) thơng qua đường kính lớn (X) là: Y = - 53.7 + 2.04 X - Với P = 0,000 nên hệ số phương trình hồi quy khác khơng - Hệ số xác định phương trình R2 =( R-Sq )= 80.4% hiệu chỉnh R2 =( R-Sq(adj) = 79.4% b Phần Response đánh KL Phần Predictor đánh DN Regression Analysis: KL versus DN The regression equation is KL = - 116 + 4.01 DN Predictor Constant DN Coef -115.54 4.0110 S = 2.58776 SE Coef 18.64 0.4209 R-Sq = 82.0% T -6.20 9.53 P 0.000 0.000 R-Sq(adj) = 81.0% Analysis of Variance Source Regression Residual Error Total DF 20 21 SS 608.09 133.93 742.02 MS 608.09 6.70 F 90.81 P 0.000 Unusual Observations Obs DN 44.9 KL 59.000 Fit 64.397 SE Fit 0.607 Residual -5.397 St Resid -2.15R R denotes an observation with a large standardized residual - Phương trình hồi quy ước tính khối lượng (Y) thơng qua đường kính nhỏ (X) là: Y = - 116 + 4.01 - Với P = 0,000 nên hệ số phương trình hồi quy khác khơng - Hệ số xác định phương trình R2 =( R-Sq )= 82.0% hiệu chỉnh R2 =( R-Sq(adj) = 81.0% Hồi quy đa biến: y = a + b1x1 + b2x2 + … Ví dụ: KL = a + b1DL + b2 DN Vào Stat/ Regression/ Regression… Phần Response đánh KL Phần Predictor đánh DL DN Regression Analysis: KL versus DL, DN The regression equation is KL = - 117 + 1.21 DL + 2.48 DN X Predictor Constant DL DN Coef -116.555 1.21473 2.4764 S = 0.759757 SE Coef 5.472 0.08323 0.1623 T -21.30 14.60 15.26 R-Sq = 98.5% P 0.000 0.000 0.000 R-Sq(adj) = 98.4% Analysis of Variance Source Regression Residual Error Total Source DL DN DF 1 DF 19 21 SS 731.05 10.97 742.02 MS 365.53 0.58 F 633.24 P 0.000 Seq SS 596.60 134.46 - Phương trình hồi quy ước tính khối lượng (Y) thơng qua đường kính lớn (X1) đường kính nhỏ (X2) là: Y = - 117 + 1.21 X1 + 2.48 X2 - Với P = 0,000 nên hệ số phương trình hồi quy khác khơng - Hệ số xác định phương trình R2 =( R-Sq )= 98.5% hiệu chỉnh R2 = ( R-Sq(adj) = 98.4% III BẢNG TƯƠNG LIÊN ( BIẾN ĐỊNH TÍNH): Ví dụ M-1.12 Bước 1: Tóm tắt trình bày liệu Bước 2: Giả thiết H0 đối thiết H1 Bước 3: Kiểm tra điều kiện Bước 4: Tính xác suất P Bước 5: So sánh P với α rút kết luận Bài làm: Bước 1: Tóm tắt trình bày liệu Vào Stat/ Tables/ Cross Tabulation and Chi-Square… Nhập vào: For rows: GIONG For colums: KETQUA Frequencies are in: TANSUAT Dí vào Chi-Square… sau tích vào Expected cell counts Tabulated statistics: GIONG, KETQUA Using frequencies in TANSUAT Rows: GIONG Columns: KETQUA - + All H 400 80.00 421.4 100 20.00 78.6 500 100.00 500.0 J 190 95.00 168.6 10 5.00 31.4 200 100.00 200.0 All 590 84.29 590.0 110 15.71 110.0 700 100.00 700.0 Cell Contents: Count % of Row Expected count Giống bò n Holstein 500 Jersey 200 Viêm nội mạc tử cung Có 100 (20%) 190 (95%) Khơng 400 (80%) 10 (5%) Bước 2: Giả thiết H0 đối thiết H1 H0: Tỷ lệ viêm nội mạc tử cung hai giống bò Holstein Jersey πHolstein= πJersey H1: Tỷ lệ viêm nội mạc tử cung hai giống bò Holstein Jersey khác πHolstein ≠ πJersey Bước 3: Kiểm tra điều kiện - Tần suất lý thuyết (tần suất ước tính) > 5: sử dụng phép thử Khi bình phương (Chi-square) - Trường hợp có giá trị tần suất ước tính lý thuyết ≤ sử dụng phép thử xác Fisher với bảng tương liên 2x2 (Fisher’s exact for 2x2 tables) Bước 4: Tính xác suất P a Với tần suất ước tính > dùng Stat/ Table/ Cross Tabulation and Chi-Square… Dí vào Chi-Square… tích vào Chi-Square analysis Tabulated statistics: GIONG, KETQUA Using frequencies in TANSUAT Rows: GIONG Columns: KETQUA - + All H 400 80.00 421.4 100 20.00 78.6 500 100.00 500.0 J 190 95.00 168.6 10 5.00 31.4 200 100.00 200.0 All 590 84.29 590.0 110 15.71 110.0 700 100.00 700.0 Cell Contents: Count % of Row Expected count Pearson Chi-Square = 24.268, DF = 1, P-Value = 0.000 Likelihood Ratio Chi-Square = 29.054, DF = 1, P-Value = 0.000 P-Value lấy Pearson Chi-Square với P-Value = 0,000 b Với tần suất ước tính ≤ dùng Stat/ Table/ Cross Tabulation and Chi-Square… Dí vào ô Other Stats… tích vào ô Fisher’s exact for 2x2 tables Tabulated statistics: GIONG, KETQUA Using frequencies in TANSUAT Rows: GIONG Columns: KETQUA - + All H 400 80.00 421.4 100 20.00 78.6 500 100.00 500.0 J 190 95.00 168.6 10 5.00 31.4 200 100.00 200.0 All 590 84.29 590.0 110 15.71 110.0 700 100.00 700.0 Cell Contents: Count % of Row Expected count Fisher's exact test: P-Value = 0.0000001 P-Value = 0.0000001 Bước 5: So sánh P với α rút kết luận a P-Value = 0,000 < 0,05 nên bác bỏ H0 chấp nhận H1.Tỷ lệ viêm nội mạc tử cung hai giống bò Holstein Jersey khác (πHolstein ≠ πJersey) Kết luận: Có mối liên hệ bệnh viêm nội mạc tử cung giống bò b P-Value = 0.0000001 < 0,05 nên bác bỏ H0 chấp nhận H1.Tỷ lệ viêm nội mạc tử cung hai giống bò Holstein Jersey khác (πHolstein ≠ πJersey) Kết luận: Có mối liên hệ bệnh viêm nội mạc tử cung giống bò -the end ... 605.51) T -2.51 P 0.017 P = 0,017 Bước 5: So sánh P với α rút kết luận: P = 0,017 < 0,05 nên bác bỏ giả thiết H0 chấp nhận giả thiết H1 Kết luận: Tăng trọng trung bình tồn đàn lợn trại khác 607 gram/ngày... P-Value = 0,000 < 0,05( α) Bước 5: So sánh P với α rút kết luận: P-Value = 0,000 < 0,05( α) nên bác bỏ giả thiết H0, chấp nhận giả thiết H1 Kết luận: Khối lượng hai giống bị khơng (µGiống ≠ µGiống... P-Value = 0,000 < 0,05( α) Bước 5: So sánh P với α rút kết luận: P-Value = 0,000 < 0,05( α) nên bác bỏ giả thiết H0, chấp nhận giả thiết H1 Kết luận: Khối lượng hai giống bị khơng (µGiống ≠ µGiống