Đang tải... (xem toàn văn)
Điều tra, khảo sát Phỏng vấn, thu thập thông tin Tổng hợp thông tin, báo cáo Kiểm tra, đối chiếu, rà soát thông tin và kết thống kê Xây dựng báo cáo thống kê
LATIN1 - SECTION ENTRY B C A D NSUAT 1504 1209 1522 1358 MEANS 1398 (3) OVERALL: MEANS STD ERR 5% LSD C.V RES DF % EFFCY 1398 62.72 216.9 90 Từ kết chạy ñược, ta quan tâm phân tích ba phần: (1) Là bảng phân tích phương sai: trường hợp này, xác suất chấp nhận Ho cột cuối bảng phân tích phương sai 0,039 nhỏ 0,05 ta bác bỏ Ho kết luận công thức (giống) khác dẫn tới kết thí nghiệm (năng suất) khác cách có ý nghĩa thống kê với ñộ tin cậy 95% (2) Là suất trung bình cơng thức (giống) (3) Là Kết tổng hợp gồm có trung bình tồn thí nghiệm, giá trị sai khác nhỏ có ý nghĩa mức xác suất nhỏ 5% sai số thí nghiệm C.V % Dùng kết suất trung bình phần (2) ñể so sánh hiệu số cặp trung bình với so với giá trị sai khác nhỏ (5% LSD) phần (3) ñể kết luận đơi trung bình so sánh có khác thực hay khơng Ví dụ so trung bình giống A với trung bình giống B ta hiệu chúng 18 nhỏ 216.9 suất hai giống khơng khác hay nói cách xác khác chưa có ý nghĩa Nhưng so sánh suất giống A với giống ñối chứng C, ta ñược hiệu chúng 313 lớn 216.9 có nghĩa suất hai giống khác cách có ý nghĩa thống kê với mức tin cậy 95% Còn sai số thí nghiệm trường hợp C.V % chấp nhận thí nghiệm xác (đối với thí nghiệm so sánh giống số tác giả đề nghị cần đạt độ xác cao hơn: C.V % nên 5) tuỳ trường hợp cụ thể người nghiên cứu có kết luận cuối Cách biểu diễn kết cụ thể báo cáo trình bày phần sau 4.1.4 Phân tích phương sai kết thí nghiệm bố trí kiểu Lattice cân đối Trong phân tích cần xác ñịnh thành phần biến ñộng phân tích phương sai bao gồm: nhắc lai, cơng thức chưa ñiều chỉnh, khối ñược ñiều chỉnh, sai số khối, cơng thức điều chỉnh, sai số hiệu Ta có bảng phân tích phương sai sau, bảng 4.7 Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….56 Bảng 4.7 ANOVA cho thiết kế thí nghiệm kiểu Lattice cân ñối Source of variation Reps Trea (unadj) Block (adj) Intrablock error Trea.(adj) Effective error Total Degree of freedom K = t1/2 K2 - K2 - ( k-1)(K2 - 1) K2 -1 ( k-1)(K2 - 1) K2( k+1) – Sum of square Mean square F Ví dụ: có thí nghiệm so sánh hiệu 16 cơng thức bón phân khác thiết kế Lattice cân đối Như thí nghiệm phải có lần nhắc lại (r = 5), số khối không ñầy ñủ k Kết thí nghiệm thu được trình bày theo sơ đồ thiết kế bảng 4.8 ñây: Bảng 4.8 Số dảnh m2 16 cơng thức bón phân, thiết kế kiểu Lattice cân ñối x Khối số Nlai (1) 147 (5) 127 (9) 147 (13) 155 Tổng Nlại Khối số (1) 155 10 (5) 182 11 (9) 137 12 (13) 185 Tổng Nlại (2) 152 (6) 155 (10) 100 (14) 195 (3) 167 (7) 162 (11) 192 (15) 192 T.khố Khối i số (4) 616 150 (8) 616 172 (12) 616 177 (16) 747 205 Nlai (6) 162 (2) 130 (14) 185 (10) 122 (11) 177 (15) 177 (3) 152 (7) 182 (16) 152 (12) 165 (8) 152 (4) 192 Nlai (1) 140 (10) 97 (7) 155 (16) 182 2595 Tổng Nlại T.khố Khối i số 646 13 (1) 220 654 14 (13) 205 626 15 (5) 165 681 16 (9) 155 2607 Tổng Nlại (5) 165 (2) 155 (15) 182 (8) 207 (9) 182 (14) 192 (3) 192 (12) 232 T.khối (13) 152 (6) 142 (11) 192 (4) 162 (7) 175 (11) 180 (3) 200 (15) 185 586 721 783 2729 T.khối Nlai (14) 202 (2) 152 (10) 150 (6) 177 639 (12) 205 (8) 187 (16) 160 (4) 172 802 724 675 689 2890 Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….57 Khối số 17 (1) 147 18 (9) 180 19 (13) 172 20 (5) 177 Tổng Nlại Nlai (10) 112 (2) 205 (6) 212 (14) 220 (15) 177 (7) 190 (3) 197 (11) 205 T.khố i (8) 583 147 (16) 742 167 (12) 773 192 (4) 827 225 2925 ðể tiến hành phân tích phương sai, thực bước sau: • Vào số liệu IRRISTAT, kết ghi lại File có sys thể hình 4.18 Hình 4.18 • Các bước phân tích IRRISTAT Làm tương tự trường hợp thiết kế kiểu ô vng latin có hình 4.19 sau: Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….58 Hình 4.19 Từ hình 4.19, chọn SODANH từ hộp Data File Variable đưa vào hộp Analysis Variate tự động vào hộp Analysis Variate w/ANOVA, ta có hình 4.20 sau Hình 4.20 Từ hình 4.20, kích chuột vào menu SSA Model, tiếp tục chọn Lattice Blocks hộp Type of Design, ñưa NLAI vào hộp Replicate, KHOI vào hộp Block, CTHUC vào hộp Treatment để có hình 4.21 sau Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….59 Hình 4.21 Từ hình 4.21, chọn option ñể hộp Heading mở, ñánh dòng chữ “Thiet ke Lattice can doi” vào hộp Heading, hộp Sorting Variate hộp Variate for Percentage Checks chọn SODANH để có hình 4.22 sau Hình 4.22 Từ hình 4.22, kích chuột vào hộp OK để chạy chương trình cho kết đây: • Kết chạy mơ hình phân tích ANOVA FOR SINGLE VARIATES - LATTICE BLOCK FILE LATTICE 7/ 1/** 9:39 PAGE Thiet ke Lattice can doi VARIATE V004 SODANH SOURCE D.F S.S M.S REPS 5946.05 1486.51 Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….60 TREATMENTS (UNADJ) 15 26994.4 1799.62 BLOCKS ADJ FOR TRTS 15 11381.8 758.789 INTRA-BLOCK ERROR 45 14533.3 322.962 TOTAL 79 58855.6 VARIANCE COMPONENT FOR BLOCKS 136.2 TREATMENTS (ADJ) 15 20988.0 1399.20 RANDOMIZED BLOCK ERROR 60 25915.2 431.919 EFFECTIVE PLOT ERROR 45 16620.2 369.337 F FOR TREATMENTS(ADJ)/INTRA-BLOCK ERROR 4.33 FPROB= 0.000 F FOR TREATMENTS(UNADJ)/RANDOMIZED BLOCK ERROR 4.17 FPROB= 0.000 EFFICIENCY 116.9% (Có phần mơ tả tất dư thừa tất công thức qua lần nhắc lại vào vị trí ta khơng quan tâm nên bỏ ) L.S.D (5%) 24.48 COEFFICIENT OF VARIATION 11.18 MEANS FOR EACH VARIETY - LATTICE BLOCK FILE LATTICE 7/ 1/** 9:39 PAGE Thiet ke Lattice can doi LATICEKQ ENTRY - SECTION 14 12 11 15 SODANH 197.2 190.5 188.2 185.7 183.9 13 176.9 175.7 173.8 169.5 168.4 16 167.8 165.8 163.0 162.9 161.0 10 118.8 MEANS 171.8 OVERALL: MEANS STD ERR 5% LSD C.V RES DF % EFFCY 171.8 8.595 24.48 11 45 117 • Phân tích kết Trong bảng phân tích phương sai, IRRISTAT tính tất nguồn biến ñộng cần thiết bảng ANOVA tổng quát ñã nêu trên, ñộ tự do, tổng bình phương phương sai cho nguồn sau: Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….61 Dòng số nhắc lại Dịng cơng thức chưa điều chỉnh Dịng khối điều chỉnh theo cơng thức Dịng 10 sai số khối Dịng 11 tổng tồn Dòng 12 hợp phần phương sai cho khối Dịng 13 cho cơng thức điều chỉnh Dòng 14 cho sai số khối ngẫu nhiên Dòng 15 cho sai số hiệu thí nghiệm Dịng 16 trị số F thực nghiệm xác suất chấp nhận H0 cho công thức ñã ñược ñiều chỉnh so với sai số khối - Dòng 17 tương tự dòng 16 cho cơng thức chưa điều chỉnh so với sai số khối ngẫu nhiên - Dòng 18 hiệu cách dùng Lattice làm tăng độ xác thí nghiệm lên 17% - Dịng 19 gí trị LSD0.05 sai số thí nghiệm - Các kết phía biểu thị giống cho ví dụ nêu với cách thiết kế khác Các quan tâm cho phân tích ví dụ thuộc dịng 16, 17, 18 19 cho thấy hai trường hợp kết điều chỉnh hay khơng ñiều chỉnh so với sai số khối khối ngẫu nhiên cho kết luận cơng thức khác dẫn đến kết khác (vì xác suất chấp nhận H0 cho hai trường hợp ñều nhỏ) Hiệu áp dụng thiết kế kiểu Lattice cân đối làm tăng độ xác thí nghiệm lên 17% Giá trị LSD0.05 24.48 Sai số thí nghiệm đạt 11,18% Nếu so trung bình cơng thức 14 với 12 có hiệu sai xấp xỉ nhỏ giá trị LSD0.05 chứng tỏ chúng khác khơng đáng tin cậy, so cơng thức 14 với 10 có kết luận ngược lại hiệu chúng lớn giá trị LSD0.05 nghĩa chúng khác ñáng tin cậy mức ñộ tin 95% Ta so sánh tất đơi trung bình với biểu diễn kết dạng số a, b, c, - Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……….62 4.2 Thí nghiệm hai nhân tố 4.2.1 Phân tích kết thí nghiệm hai nhân tố thiết kế kiểu hoàn toàn ngẫu nhiên (CRD) Trong trường hợp này, mức nhân tố ñược phối hợp thành tổ hợp công thức khác tổ hợp coi cơng thức đơn lẻ thí nghiệm nhân tố cho việc bố trí thí nghiệm theo kiểu hồn tồn ngẫu nhiên Ví dụ: Có thí nghiệm nghiên cứu ảnh hưởng việc bón phối hợp a mức ñạm (ao, a1) b mức lân (bo, b1, b2) khác tới suất lúa (tạ/ha), thí nghiệm bố trí kiểu hồn tồn ngẫu nhiên với bốn lần nhắc lại ô xi măng Năng suất thí nghiệm ghi lại bảng 4.9 sau ñây: Bảng 4.9 Năng suất lúa tổ hợp công thức lần nhắc lại Năng suất (tạ/ha) ðạm (A) Lân (B) Nlại Nlại Nlại Nlại bo 24 26 23 27 ao b1 28 30 30 27 28 31 32 28 b2 bo 31 34 34 31 a1 b1 47 45 47 46 b2 59 51 65 60 Trong phân tích phương sai cho trường hợp này, ta cần phân tích để thấy ñược vai trò ảnh hưởng ñạm, lân phối hợp ñạm lân khác ñến suất lúa Bảng phân tích phương sai hình thành bảng 4.10 sau: Bảng 4.10 Bảng phân tích phương sai Source of Degree of Sum of Mean square variation freedom square Factor A a-1 Factor B b-1 AxB (a-1)(b-1) Error (r-1)ab Total rab-1 Các bước tiến hành cho phân tích phương sai • Vào số liệu IRRISTAT Save file với SYS để có kết hình 4.23 sau: F Hình 4.23 Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 63 • Các bước phân tích Làm liên tục bốn bước đầu để có hình tương tự hình 4.2, 4.3, 4.4, 4.5 trường hợp phân tích ANOVA thiết kế kiểu CRD để có hình 4.24 sau Hình 4.24 Từ hình 4.24, chọn biến NSUAT đưa vào hộp Analysis Variates, ba biến cịn lại đưa vào hộp Factor Tiếp tục ñưa biến DAM$, LAN$ DAM$ x LAN$ từ hộp Facor vào hộp ANOVA Model Specification để hình 4.25 Hình 4.25 Từ hình 4.25, vào menu Options để hộp Heading mở đánh dịng chữ “ANOVA cho thi nghiem hai nhan to thiet ke kieu hoan toan ngau nhien” xong kích chuột vào hộp OK để chạy kết ñây BALANCED ANOVA FOR VARIATE NSUAT FILE FSAIHNTO 6/ 1/** 23: PAGE ANOVA cho thi nghiem thiet ke kieu hoan toan ngau nhien VARIATE V004 NSUAT Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 64 LN ER SOURCE OF VARIATION DF SUMS OF MEAN SQUARES SQUARES F RATIO PROB LN ============================================================================ = DAM$ 2016.67 2016.67 226.88 0.000 LAN$ 962.333 481.167 54.13 0.000 DAM$*LAN$ 462.333 231.167 26.01 0.000 * RESIDUAL 18 159.999 8.88885 -* TOTAL (CORRECTED) 23 3601.33 156.580 TABLE OF MEANS FOR FACTORIAL EFFECTS FILE FSAIHNTO 6/ 1/** 23: PAGE ANOVA cho thi nghiem thiet ke kieu hoan toan ngau nhien MEANS FOR EFFECT DAM$ DAM$ NOS 12 12 ao a1 NSUAT 27.5000 45.8333 SE(N= 12) 0.860661 5%LSD 18DF 2.55715 MEANS FOR EFFECT LAN$ LAN$ NOS 8 bo b1 b2 NSUAT 28.7500 37.0000 44.2500 SE(N= 8) 1.05409 5%LSD 18DF 3.13186 MEANS FOR EFFECT DAM$*LAN$ DAM$ ao ao ao a1 a1 a1 LAN$ bo b1 b2 bo b1 b2 NOS 4 4 4 NSUAT 25.0000 27.7500 29.7500 32.5000 46.2500 58.7500 SE(N= 4) 1.49071 5%LSD 18DF 4.42911 - Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 65 ANALYSIS OF VARIANCE SUMMARY TABLE FILE FSAIHNTO 6/ 1/** 23: PAGE ANOVA cho thi nghiem thiet ke kieu hoan toan ngau nhien F-PROBABLIITY VALUES FOR EACH EFFECT IN THE MODEL SECTION - VARIATE |DAM$ |LAN$ |$ GRAND MEAN STANDARD DEVIATION C OF V |DAM$*LAN| (N= 24) SD/MEAN | | | NO BASED ON BASED ON OBS TOTAL SS RESID SS 12.513 2.9814 % | | | | | | | | NSUAT 0.0000 24 36.667 8.1 0.0000 0.0000 • Phân tích kết Nhìn vào bảng phân tích phương sai ta thấy tương tự trường hợp thí nghiệm nhân tố bố trí kiểu hồn tồn ngẫu nhiên, bảng khơng có thành phần biến động bình phương nhắc lại trường hợp coi đơn vị thí nghiệm hồn tồn đồng Vì thí nghiệm hai nhân tố nên bảng ANOVA có ba nguồn biến ñộng nhân tố nghiên cứu gây nên: biến ñộng ñạm, biến ñộng lân biến ñộng tương tác ñạm lân Kết cho thấy đạm khác nhau, lân bón khác phối hợp ñạm với lân khác ñều dẫn tới suất khác có ý nghĩa thống kê cao (xác suất chấp nhận H0 bé = 0,000) Phần bảng suất trung bình theo mức nhân tố A giá trị LSD0.05 tương ứng chúng, từ ñây dễ dàng so sánh chênh lệch suất hai mức A với LSD0.05 ñể kết luận Cụ thể cho ví dụ này, chênh lệch suất hai mức nhân tố A lớn so với LSD0.05 nên đến kết luận mức nhân tố A khác ñã dẫn ñến suất khác cách có ý nghĩa Tương tự mức lân bón khác cho suất khác có ý nghĩa thống kê, cụ thể mức lân b2 cho suất cao sau đến mức b1 thấp b0 Bảng suất tác ñộng kết hợp mức nhân tố A B, cách so sánh tương tự ta thấy kết hợp mức a1b2 ñã cho suất cao nhất, sau đến a1b1 thấp suất kết hợp hai mức a0b0 Nhìn vào dòng cuối cùng, cột C of V thấy CV% 8,1%, với thí nghiệm xi măng độ xác khơng cao, theo lý thiuyết chung phải 5% tốt 4.2.2 Phân tích phương sai kết thí nghiệm hai nhân tố thiết kế kiểu khối ngẫu nhiên đầy đủ (RCB) Phân tích phương sai cho trường hợp ñược xây dựng ñể xác ñịnh hiệu nhân tố A, B tương tác chúng tới kết thí nghiệm ñồng thời cần xác ñịnh biến ñộng khối ñể loại trừ khỏi sai số thí nghiệm Giả Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 66 thiết thí nghiệm nhân tố A có a mức, nhân tố B có b mức với r lần nhắc lại, bảng phân tích phương sai cấu tạo bảng 4.11 sau: Bảng 4.11 ANOVA cho thí nghiệm hai nhân tố thiết kế kiểu khối ngẫu nhiên ñầy ñủ Source of Degree of Sum of Mean square F variation freedom square Reps r-1 Factor A a-1 Factor B b-1 AxB (a-1)(b-1) Error (r-1)(ab-1) Total rab-1 Ví dụ: Dùng file số liệu có sẵn IRRISTAT ví dụ (thí nghiệm hai nhân tố thiết kế kiểu hồn tồn ngẫu nhiên) để tiến hành phân tích ANOVA cho trường hợp Cách tiến hành cụ thể sau: • Thực bước phân tích tương tự trường hợp phân tích ANOVA cho kiểu thiết kế CRD đến bước có hình 4.24 ví dụ ñây Tiếp tục ñưa biến NSUAT vào hộp Analysis Variate, ba biến cịn lại đưa vào hộp Factor Sau đưa biến NLAI, DAM$, LAN$ từ hộp Factor DAM$ x LAN$ xuống hộp ANOVA Model Specification ñể có hình 4.26 sau: Hình 4.26 Từ hình 4.26, kích chuột vào menu Options để hộp Heading mở xong ñánh dòng chữ “ANOVA cho thi nghiem hai nhan to thiet ke kieu khoi ngau hien day du” vao hộp sau kích chuột vào hộp OK để chạy phân tích cho kết BALANCED ANOVA FOR VARIATE NSUAT FILE FSAIHNTO 6/ 1/** 22:35 PAGE ANOVA cho thi nghiem hai nhan to thiet ke kieu khoi ngau nhien day du VARIATE V004 NSUAT Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 67 ... kết thí nghiệm đồng thời cần xác ñịnh biến ñộng khối ñể loại trừ khỏi sai số thí nghiệm Giả Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT? ??…… 66 thiết thí. .. học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT? ??…… .62 4.2 Thí nghiệm hai nhân tố 4.2.1 Phân tích kết thí nghiệm hai nhân tố thiết kế kiểu hồn tồn ngẫu nhiên (CRD) Trong... phương sai • Vào số liệu IRRISTAT Save file với SYS để có kết hình 4.23 sau: F Hình 4.23 Trường ðại học Nơng nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT? ??…… 63 • Các bước