ðây là cuốn sách hướng dẫn thực hành trên máy vi tính ñể phân tích các kết quả thực nghiệm trong nông nghiệp bằng IRRISTAT 4.0 trong Window. Với cách chỉ dẫn ngắn gọn dùng nhiều hình ảnh minh hoạ các bước cụ thể nên rất tiện cho người sử dụng. ðặc biệt là với những người không am hiểu nhiều về máy tính cũng như kiến thức về thống kê cũng có thể sử dụng ñược. Sách rất cần thiết cho sinh viên và các nghiên cứu viên nông nghiệp trong nghiên cứu, thiết kế thí nghiệm và xử lý các kết quả thí nghiệm của họ
Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam BỘ GIÁO DỤC VÀ ðÀO TẠO TRƯỜNG ðẠI HỌC NÔNG NGHIỆP HÀ NỘI Phạm Tiến Dũng THIẾT KẾ THÍ NGHIỆM VÀ XỬ LÝ KẾT QUẢ Bằng phần mềm thống kê IRRISTAT Hà nội 2008 https://sites.google.com/site/lophocphank57vnua/ Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam LỜI TỰA ðây sách hướng dẫn thực hành máy vi tính để phân tích kết thực nghiệm nơng nghiệp IRRISTAT 4.0 Window Với cách dẫn ngắn gọn dùng nhiều hình ảnh minh hoạ bước cụ thể nên tiện cho người sử dụng ðặc biệt với người không am hiểu nhiều máy tính kiến thức thống kê sử dụng ñược Sách cần thiết cho sinh viên nghiên cứu viên nông nghiệp nghiên cứu, thiết kế thí nghiệm xử lý kết thí nghiệm họ Nội dung gồm Chương Chương 1: Giới thiệu khái quát nội dung chức IRRISTAT 4.0 Window Chương 2: Giới thiệu cách quản lý số liệu IRRISTAT 4.0 Chương 3: Thiết kế thí nghiệm Chương 4: Phân tích phương sai kết thí nghiệm Chương 5: Phân tích tương quan hồi quy Chương 1và ñược dịch từ Tutorial Manual “IRRISTAT 4.0 for Windows” Biometric Unit Internationnal Rice Research Institute biên soạn Các chương 3, 4, tác giả ñưa ví dụ cụ thể hướng dẫn cách thực máy tính cho trường hợp thiết kế phân tích kết chạy ðây ví dụ mẫu cho người học làm theo IRRISTAT 4.0 5.0 Trong khn khổ điều kiện cịn nhiều hạn chế, khơng tránh khỏi thiếu sót biên soạn Rất mong nhận ñược nhiều ý kiến ñóng góp quý báu bạn ñọc Tác giả https://sites.google.com/site/lophocphank57vnua/ 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……….i Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam LỜI GIỚI THIỆU IRRISTAT chương trình xử lý số liệu thống kê Viện nghiên cứu lúa Quốc tế ðây chương trình quen thuộc với cán ngành trồng trọt, phục vụ thiết thực cho việc bố trí thí nghiệm xử lý kết nghiên cứu lúa trồng khác Nội dung cách dùng IRRISTAT ver 92-1 ñã ñược giới thiều tập sách tin học cho Cao học ngành sinh học trường ðại học Nơng nghiệp IRRISTAT lúc gồm phần: 1.Bố trí thí nghiệm đồng ruộng với kiểu chính: hồn tồn ngẫu nhiên (RCD), khối ngẫu nhiên (RCB), chia (Split-plot), chia băng (Strip-plot) Mỗi kiểu có 1, nhiều nhân tố (ngồi khối coi nhân tố dặc biệt, thực chất hạn chế điều kiện khơng thể lựa chọn nhiều thí nghiệm đồng nhất) Phần ngày giúp cho việc thết kế thí nghiệm chia khối, chọn kích thước ơ, xếp ngẫu nhiên theo yêu cầu kiểu thí nghiệm 2.Phân tích phương sai kiểu thí nghiệm trên, ngồi cịn có phần Pooled ANOVA Combined ANOVA ñể liên kết kết thí nghiệm nhiều vùng, qua nhiều năm cân đo nhiều đợt Sau phân tích phương sai so sánh trung bình qua phương pháp LSD Duncan, sâu để đánh giá số kết thơng qua phương pháp lập tương phản Phần viết tỉ mỉ sát với nội dung Statistical Procedures for Agricultural Research Kwanchai A.Gomez & Arturo A.Gomez viện nghiên cứu lúa quốc tế 3.Tương quan, hồi quy tuyến tính đa thức Phần có phân tích hồi quy ñơn, hồi quy bội tuyến tính, hồi quy ña thức, sau tính hệ số hồi quy có bảng phân tích phương sai phân tích phần dư 4.Một số tiện ích để quản lý số liệu chép tệp, biến đổi số liệu Nhìn chung số liệu cỡ trung bình, thiết thực cho phân tích số liệu thí nghiệm trồng Tuy nhiên, chương trình cịn nhiều nhược điểm chạy chậm, có nhiều câu hỏi mà người dùng khơng hưỡng dẫn đầy đủ khơng có tài liệu chi tiết khơng dùng được, phần trợ giúp khơng có, vẽ thơ sơ xấu, nội dung hẹp (chỉ bao gồm vấn ñề sở mơn thống kê sinh học) chưa đáp ứng ñược yêu cầu nghiên cứu nông nghiệp Bộ chương trình IRRISTAT ver 4.0 bước nhảy vọt chuyển hồn tồn sang chạy Windows với nhiều nội dung Có thể giới thiệu nét sau: Tồn chương trình hoạt ñộng dạng Menu theo ñúng khuynh hướng chung chương trình chạy Windows Các liệu ñược truy nhập dễ dàng thuận tiện trước nhiều Bộ chương trình giữ lại phần phân tích phương sai cũ tổng quát phân chia thành trường hợp cân ñối phần hồi quy tổng qt với hình vẽ đẹp chương trình thống kê khác, ngồi bổ sung nhiều phần phân tích ảnh hưởng giống mơi trường, phân tích tính ổn định, số vấn đề định lượng phân tích ổ gen phân chia nhóm theo số mơ hình https://sites.google.com/site/lophocphank57vnua/ 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……….ii Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nơng Nghiệp Việt Nam Tồn phần cần thiết cho việc nhiên cứu sâu giống ảnh hưởng môi trường Phần chia nhóm theo số mơ hình mang tính phác hoạ, tim cách ứng dụng số mơ hình thống kê nhiều chiều ñể xem xét gần gũi mốt số giống hay số mơi trường Tuy có nhiều ưu ñiểm ñây bước ñầu chuyển sang chạy Windows nên số thiếu sót, khơng thuận tiện tổng qt SAS, SPSS Về nội dung ñã bổ sung thêm nhiều phần tập trung xung quanh việc đánh giá giống mơi trường Phần trợ giúp nghèo nàn khơng thoả mãn cho người dùng gặp khó khăn Trên số nhận xét chung chương trình IRRISTAT để bạn đọc có số ý niệm sơ lược khái quát IRRISTAT Sách không nhằm giới thiệu tổng quát IRRISTAT mà tập trung vào phần: (1) Giới thiệu cách dùng IRRISTAT ver 4.0 (2) Dùng IRRISTAT ñể thiết kế phân tích số mẫu phân tích phương sai hồi quy hay gặp nghiên cứu trồng Phần giới thiệu ñược soạn tỉ mỉ theo tài liệu hướng dẫn cách dùng IRRISTAT ver 4.0 phận Biometric Viện lúa Quốc tế biên soạn Qua việc giới thiệu loại cửa sổ (cửa sổ chính, cửa sổ quản lý liệu, cửa sổ văn bản, cửa sổ đồ hoạ) hình dung cách làm việc với IRRISTAT Trong cửa sổ có cơng việc (menu) hình xuất nhiều loại cơng cụ nút ñiều khiển, nút tuỳ chọn, hộp danh sách, hộp kiểm soát, hộp soạn thảo… Việc sử dụng menu, nút hộp ñã ñược giới thiệu kỹ lưỡng với hình minh hoạ cụ thể Trong phần II tác giả ñã chọn số mẫu thiết kế thí nghiệm thường gặp nghiên cứu trồng có kèm theo ví dụ kết xử lý IRRISTAT ver 4.0 ðối với thí nghiệm nhân tố có kiểu thiết kế: Thí nghiệm nhân tố thiết kế kiểu hoàn toàn ngẫu nhiên (CRD) Thí nghiệm nhân tố thiết kế kiểu khối ngẫu nhiên đầy đủ (RCB) Thí nghiệm nhân tố bố trí kiểu Ơ vng La tinh (LS) Thí nghiệm nhân tố bố trí kiểu lưới vng (Lattice design) Trừ kiểu lưới vng giới thiệu kỹ cịn kiểu đầu đề cập sơ qua có nhiều tài liệu hướng dẫn dạy giáo trình phương pháp thí nghiệm ðối với thí nghiệm nhân tố có kiểu thiết kế: Thí nghiệm nhân tố kiểu hồn tồn ngẫu nhiên (CRD) Thí nghiệm nhân tố kiểu khối ngẫu nhiên ñầy ñủ (RCB) Thí nghiệm nhân tố kiểu chia lớn nhỏ (Split-plot) Thí nghiệm nhân tố kiểu chia băng (Strip-plot) Cả kiểu ñều ñược giới thiệu kỹ tài liệu khác khơng viết kỹ không nêu rõ lý phải thiết kế kèm theo việc phân tích ưu khuyết kiểu ðối với thí nghiệm nhân tố tác giả giới thiệu kiểu quen thuộc RCB chia ô lớn vừa ô nhỏ (Split-split-plot) Ngoài kiểu thiết kế thí nghiệm nói tác giả chọn vấn ñề ñang ñược nhiều người quan tâm phân tích số liệu thiết kế nhiều địa điểm thiết kế qua nhiều năm Phân tích hồi quy gồm: https://sites.google.com/site/lophocphank57vnua/ 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……….iii Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam - Hồi quy tuyến tính đơn Phần độc giả quen nhanh chóng hiểu kết phân tích - Hồi quy bội tuyến tính Phần giới thiệu kỹ phức tạp có nhiều chi tiết cần phải nói rõ để khơng phạm phải sai lầm sử dụng kết phân tích Hiện sinh viên nghành kỹ thuật nơng nghiệp khối trồng học tin học, xác suất thống kê, phương pháp thí nghiệm với khối lượng kiến thức trình bày khơng nhỏ thiếu tài liệu hướng dẫn tỉ mỉ việc dùng chương trình thống kê để xử lý số liệu Tơi hy vọng sách đáp ứng ñược yêu cầu ñông ñảo bạn ñọc Nhà giáo Ưu tú Nguyễn ðình Hiền https://sites.google.com/site/lophocphank57vnua/ 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……….iv Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Chương GIỚI THIỆU CƠ BẢN VỀ IRRISTAT 4.0 IRRITAT chương trình phần mềm máy tính viết IRRI (Viện ghiên cứu lúa quốc tế Philippine) qua nhiều phiên khác nhau, phiên trước thường chạy mơi trường DOS, sử dụng thuận tiện, phiên 4.0 chạy ñược môi trường Window nên sử dụng thuận tiện IRRISTAT có chức để thiết kế thí nghiệm, quản lý phân tích thống kê bản, phân tích phương sai hồi qui số liệu thực nghiệm thu nghiên cứu khoa học Chương trình chạy hệ thống hoạt ñộng Window 32 bít Khi cài đăt chương trình vào máy, chạy file SETUP EXE ñĩa cài ñặt Khi khởi động chương trình, chạy file IRRISTAT EXE kích chuột vào biểu tượng cài đặt, xem biểu tượng hình 1.1 Hiện phần mềm có phiên 5.0 xong khơng khác nhiều Vậy cho mục đích sách người đọc sử dụng hai phiên sở giới thiệu sách Hình 1.1 Biểu tượng Shortcut IRRISTAT hình IRRISTAT biên soạn chủ yếu để thiết kế thí nghiệm, phân tích số liệu thu từ thí nghiệm ñược thiết kế thực ñúng phương pháp, có nhiều nội dung sử dụng để phân tích số liệu từ nguồn nghiên cứu khác Các mơ hình phân tích chức khác bao gồm: Quản lý số liệu spreadsheet Trang ghi kết Phân tích phương sai https://sites.google.com/site/lophocphank57vnua/ 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……… Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Phân tích hồi qui Phân tích ảnh hưởng giống mơi trường Phân tích di truyền số lượng Phân tích theo nhóm Vẽ đồ thị Ngẫu nhiên hóa sơ đồ thí nghiệm, phương sai trung bình nhiều nhân tố đa thức trực giao 1.1 Các cửa sổ Window 1.1.1 Các menu Window (Main Window) Cửa sổ Window cung cấp cho người sử dụng menu công cụ hình 1.2 Hình 1.2 Các menu Window 1.1.2 Cửa sổ quản lý số liệu (Data Editor) ðể mở ñược cửa sổ data editor, chọn Window cửa sổ xong kích đúp chuột vào data editor, data editor mở hình 1.3 Hình 1.3 Cửa sổ quản lý số liệu (data editor) https://sites.google.com/site/lophocphank57vnua/ 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……… Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nơng Nghiệp Việt Nam Data editor trình bày số liệu dạng hàng cột Mỗi hàng trường hợp cột biến Có thể ghi số liệu vào hàng, cột nhập số liệu từ file trang văn file excel dbase Một lúc mở ñồng thời nhiều file số liệu Số file ñược mở nhiều hay tuỳ thuộc khả nhớ máy tính Ví dụ bảng số liệu data editor sau (Hình 1.4): Hình 1.4 Bảng số liệu ñược lưu trữ Data editor 1.1.3 Text Editor (Thường cửa sổ quản lý file kết quả) ðể có text editor, chọn từ menu Window cửa sổ xong chọn text editor có text editor mở hình 1.5 Hình 1.5 Cửa Text editor https://sites.google.com/site/lophocphank57vnua/ 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……… Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nơng Nghiệp Việt Nam Ví dụ kết tính tốn ghi text editor 1.6: Hình 1.6 Kết phân tích lưu lại Text editor 1.1.4 Cửa sổ ñể vẽ ñồ thị (Graph Window) Mở cửa sổ từ menu Analysis cửa sổ Window hình 1.7 dây Hình 1.7 Cửa sổ ñể vẽ ñồ thị (Graph Window) Ví dụ ñồ thị vẽ Graph Window hình 1.8 sau ñây: https://sites.google.com/site/lophocphank57vnua/ 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……… Lớp Học Phần VNUA - Khoa Nơng Học - Học Viện Nơng Nghiệp Việt Nam Hình 1.8 ðồ thị ñược vẽ Graph Window 1.2 Các menu (menus) Mỗi cửa Window IRRISTAT có menu riêng có lựa chọn thích hợp với chức riêng chúng 1.2.1 Các menu Window - Menu File (hình 1.9): dùng để mở file, dùng data editor text editor ñể mở file (tuỳ thuộc vào kiểu file), dùng ñể save the log file, save-as the log file in the log file từ menu Hình 1.9 Menu File Main Window - Menu Edit (hình 1.10): sử dụng để cut, copy, past để xố Log window https://sites.google.com/site/lophocphank57vnua/ 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……… Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Kích chuột vào OK để chạy mơ hình, có kết sau: BALANCED ANOVA FOR VARIATE YIELD FILE QUANOI 11/ 9/ 14:43 :PAGE VARIATE V004 YIELD LN SOURCE OF VARIATION DF SUMS OF MEAN F RATIO PROB ER SQUARES SQUARES LN ============================================================================= SITE 5.87069 2.93534 4.67 0.060 2 SITE*REP 3.77451 629085 1.46 0.224 NTRT 62.9067 12.5813 29.26 0.000 NTRT*SITE 10 11.2305 1.12305 2.61 0.020 * RESIDUAL 30 12.8985 429949 * TOTAL (CORRECTED) 53 96.6809 1.82417 TABLE OF MEANS FOR FACTORIAL EFFECTS FILE QUANOI 11/ 9/ 14:43 :PAGE MEANS FOR EFFECT SITE SITE NOS 18 18 18 YIELD 5.43567 5.99800 5.21478 SE(N= 18) 0.186947 5%LSD 6DF 0.646680 MEANS FOR EFFECT SITE*REP SITE 1 2 3 REP 3 NOS 6 6 6 6 YIELD 5.63183 5.51183 5.16333 5.78100 5.93250 6.28050 5.66500 5.18383 4.79550 SE(N= 6) 0.267691 5%LSD 30DF 0.773113 MEANS FOR EFFECT NTRT NTRT NOS 9 9 9 YIELD 3.39267 5.04533 5.77333 6.38678 6.51178 6.18700 SE(N= 9) 0.218568 5%LSD 30DF 0.631244 MEANS FOR EFFECT NTRT*SITE NTRT 1 SITE NOS 3 YIELD 2.38467 3.71200 4.08133 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……… https://sites.google.com/site/lophocphank57vnua/ 95 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam 2 3 4 5 6 3 3 3 3 3 3 3 3 3 3 4.34800 5.65433 5.13367 5.78100 6.24900 5.29000 6.66100 6.68100 5.81833 6.79367 7.08900 5.65267 6.64567 6.60267 5.31267 SE(N= 3) 0.378572 5%LSD 30DF 1.09335 ANALYSIS OF VARIANCE SUMMARY TABLE FILE QUANOI 11/ 9/ 14:43 :PAGE F-PROBABLIITY VALUES FOR EACH EFFECT IN THE MODEL SECTION - VARIATE |NTRT*SIT| GRAND MEAN STANDARD (N= SD/MEAN | 54) DEVIATION C OF V |SITE |SITE*REP|NTRT | | |E | | | | 0.0000 | | NO BASED ON BASED ON OBS 54 5.5495 TOTAL SS 1.3506 RESID SS 0.65571 % | | YIELD 0.0205 | 11.8 0.0600 | 0.2238 Nhìn vào dịng 10 phần kết ta kết luận: nơi khác không cho kết khác dùng kết phân tích phía để phân tích thơng thường trường hợp thí nghiệm yếu tố Nếu người nghiên cứu quan tâm chặt chẽ kết luận nơi khác cho kết khác (vì xác suất nhận Ho chưa diển hình dịng 7), cần phân tích riêng cho nơi Trong trường hợp này, người nghiên cứu nên làm lại thí nghiệm để có đánh giá xác 4.5 Phân tích phương sai kết thí nghiệm bố trí nhiều năm Kết thí nghiệm thu ñược bảng sau: Bảng Năng suất giống lúa thí nghiệm thiết kế RCB với lần nhắc lại thời vụ năm liên tiếp Giống Năm Năm số Nlai Nlai Nlai T.số Nlai Nlai Nlai 3036 4177 3884 11097 1981 3198 3726 T.số 8905 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……… 96 https://sites.google.com/site/lophocphank57vnua/ Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam 1369 1554 1899 4822 3751 2391 3714 9856 5311 5091 4839 15241 3868 3134 3487 10489 2559 3980 3853 10392 2729 2786 2598 8113 1291 1705 2130 5126 3222 3554 2452 9228 3452 3548 4640 11640 4250 4134 3339 11723 1812 2914 958 5684 3336 4073 2885 10294 T.theo 64002 68608 Năm Nguồn: Gomez, 1984 Ký hiệu: số năm: y; Nhắc lại: r; Công thức: t Bảng đề cương phân tích phương sai hình thành sau: Source of variation Year(Y) Rep Within year Treatment(T) TxY Pooled Error Total Degree of freedom y –1 = y (r – ) = t –1 = (y-1)(t-1) = y(r-1)(t-1) = 24 yrt-1 = 41 Vào số liệu IRRISTAT dạng 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……… https://sites.google.com/site/lophocphank57vnua/ 97 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Các bước phân tích phương sai IRRISTAT Làm tương tự phân tích ANOVA qua nhiều nơi thay nơi năm thể bảng dây Tiếp tục làm bướcc phân tích qua nhiều nơi có kết sau: BALANCED ANOVA FOR VARIATE NSUAT FILE QUANAM 13/ 9/ 8:41 :PAGE Phan Tich ANOVA qua nam RCB VARIATE V004 NSUAT LN SOURCE OF VARIATION DF SUMS OF MEAN F RATIO PROB ER SQUARES SQUARES LN ============================================================================= NAM 505124 505124 1.36 0.309 2 NAM*NLAI 148246E+07 370616 1.01 0.422 GIONG 191589E+08 319315E+07 8.71 0.000 NAM*GIONG 154960E+08 258267E+07 7.04 0.000 * RESIDUAL 24 879934E+07 366639 * TOTAL (CORRECTED) 41 454418E+08 110834E+07 TABLE OF MEANS FOR FACTORIAL EFFECTS FILE QUANAM 13/ 9/ 8:41 :PAGE Phan Tich ANOVA qua nam RCB MEANS FOR EFFECT NAM - 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……… https://sites.google.com/site/lophocphank57vnua/ 98 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam NAM NOS 21 21 NSUAT 3047.71 3267.05 SE(N= 21) 132.847 5%LSD 4DF 520.732 MEANS FOR EFFECT NAM*NLAI NAM 1 2 NLAI NOS 7 7 7 3 NSUAT 2690.00 3281.29 3171.86 3305.29 3324.29 3171.57 SE(N= 7) 228.860 5%LSD 24DF 667.979 MEANS FOR EFFECT GIONG GIONG NOS 6 6 6 6 NSUAT 3333.67 2446.33 4288.33 3084.17 2392.33 3893.83 2663.00 SE(N= 6) 247.197 5%LSD 24DF 721.500 MEANS FOR EFFECT NAM*GIONG NAM 1 1 1 2 2 2 GIONG NOS 3 3 3 3 3 3 3 7 NSUAT 3699.00 1607.33 5080.33 3464.00 1708.67 3880.00 1894.67 2968.33 3285.33 3496.33 2704.33 3076.00 3907.67 3431.33 SE(N= 3) 349.590 5%LSD 24DF 1020.35 ANALYSIS OF VARIANCE SUMMARY TABLE FILE QUANAM 13/ 9/ 8:41 :PAGE Phan Tich ANOVA qua nam RCB F-PROBABLIITY VALUES FOR EACH EFFECT IN THE MODEL SECTION - VARIATE |NAM*GION| GRAND MEAN STANDARD (N= SD/MEAN | 42) DEVIATION C OF V |NAM |NAM*NLAI|GIONG | | |G | | | | NO BASED ON BASED ON % | | 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……… https://sites.google.com/site/lophocphank57vnua/ 99 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam OBS TOTAL SS RESID SS 1052.8 605.51 | | | | | NSUAT 0.0002 42 3157.4 19.2 0.3086 0.4224 0.0000 Cách phân tích kết tương tự phân tích kết phân tích phương sai qua nơi, coi năm nơi phân tích phương sai qua nơi cách nhận xét tương tự 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……… https://sites.google.com/site/lophocphank57vnua/ 100 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nơng Nghiệp Việt Nam Chương PHÂN TÍCH HỒI QUI 5.1 Hồi quy tuyến tính đơn Có cặp số liệu sau Bảng: Kết theo dõi số dảnh (x) suất lúa (y) N suất (kg/ha) Mơ hình tuyến tính là: Số dảnh (dảnh/m2) 4862 160 5244 175 5128 192 5052 195 5298 238 5410 240 5234 252 5608 282 Y = ax + b Vào số liệu IRRISTAT dạng sau: Các bước phân tích: Chọn Analysis Regression Làm tương tự phân tích phương sai có ảnh 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……… https://sites.google.com/site/lophocphank57vnua/ 101 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Chọn x Y ñưa sang hộp Analysis Variates, sau ñó chọn Regression Model ñể có ảnh sau Chọn New ñể có ảnh 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……… https://sites.google.com/site/lophocphank57vnua/ 102 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam ðưa Y vào hộp Response Variate; x vào hộp Independent Variate chọn Add to the Model hộp Selection để có ảnh sau: Kích chuột vào OK để chạy mơ hình, có kết 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……… https://sites.google.com/site/lophocphank57vnua/ 103 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Từ kết cho thấy hồi quy lập có độ tin cậy cao với phương trình tuyến tính bậc xác định sau: Y = 4.555X + 4241.12 5.2 Hồi quy bôi tuyến tính Hồi quy bội tuyến tính biểu diễn dạng phương trình tổng quát sau: Y = a1X1 + a2X2 + … + anXn + B Trong đó: X1, X2, … Xn biến ñộc lập a1, a2, … hệ số biến ñộc lập B số Y biến phụ thuộc Ví dụ ta có tập hợp số liệu theo dõi suất phân tích hàm lượng N, P, K hạt từ 45 mảnh ruộng sau: PLOT Y14 G 1678 4265 2431 2431 4461 3110 4469 N% 0.9849 1.1714 1.0756 1.0435 1.2101 1.2084 1.1643 P% 0.0901 0.0926 0.0886 0.0807 0.0852 0.0845 0.0972 K% 0.3987 0.3814 0.4134 0.4027 0.3851 0.3746 0.3794 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……… https://sites.google.com/site/lophocphank57vnua/ 104 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 4194 3379 7132 4359 3646 6917 6692 7028 7196 4220 7250 4915 4185 7463 7367 7860 7624 2830 3705 3280 2906 4041 3295 4147 4102 3509 5087 4353 3915 5122 4660 5150 5730 3744 5363 4628 4212 5063 1.2369 1.0247 1.4624 0.9954 1.0404 1.5309 1.5001 1.6112 1.6009 0.9392 1.6738 0.9525 0.8468 1.602 1.6678 1.5521 1.5707 0.8732 1.0743 1.0124 0.9288 1.0279 1.1412 1.1594 1.0733 0.8955 1.2949 1.0541 0.9689 1.4547 1.2749 1.3373 1.2832 0.9977 1.7496 1.0091 0.8938 1.6833 0.085 0.2901 0.2467 0.2872 0.2487 0.2868 0.2776 0.2852 0.3104 0.2558 0.1985 0.2712 0.2474 0.3213 0.1793 0.3058 0.2969 0.2582 0.2396 0.2665 0.2705 0.2492 0.243 0.2867 0.2624 0.2395 0.2548 0.2604 0.2469 0.3239 0.2454 0.3057 0.3059 0.2103 0.1617 0.2208 0.2008 0.2609 0.4373 0.3108 0.2545 0.3095 0.3137 0.3071 0.2506 0.3091 0.3168 0.3039 0.2474 0.3116 0.3051 0.3144 0.2295 0.3296 0.3078 0.3204 0.239 0.3016 0.3078 0.2925 0.2475 0.2924 0.2736 0.3815 0.2978 0.3597 0.3806 0.3195 0.2556 0.3158 0.3104 0.2943 0.2353 0.2993 0.2919 0.2983 Vào số liệu bảng IRRSTAT dạng 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……… https://sites.google.com/site/lophocphank57vnua/ 105 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Tạo thêm biến ñể ghi giá trị dư thừa giá trị phù hợp sau tính hồi qui để có ảnh sau: Tạo xong, save file lại ñể chạy hồi qui 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……… https://sites.google.com/site/lophocphank57vnua/ 106 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Chọn Analysis Regression, viết lệnh cho Comand File xong Open ñể mở File số liệu ảnh dưới: Từ ảnh trên, ñưa biến GY14,N%, P%, K% vào hộp Analysis Variates Xong chọn Regression Model, chọn New chuyển biến từ Analysis Variates vào hộp Response Variate Independent Variates, chọn Add to the Model hộp Selection để có ảnh 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……… https://sites.google.com/site/lophocphank57vnua/ 107 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Chọn Option, Correclation, Compute and Test Residuals ảnh xong kích chuột vào OK để chạy mơ hình Kết chạy mơ hình: CORRELATION MATRIX AND STANDARD ERRORS FILE HQUYBOI 14/ 9/ 11:27 :PAGE SECTION 1 CORR AND RESID SES 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……… https://sites.google.com/site/lophocphank57vnua/ 108 Lớp Học Phần VNUA - Khoa Nông Học - Học Viện Nông Nghiệp Việt Nam Y14G N% P% K% 0.1584525E+04 0.8156641E+00 0.2676436E+00 0.4406525E+00 0.1723501E+00 0.7490746E-01 -0.4226677E+00 -0.3386678E+00 -0.5344221E+00 0.5054769E-01 Y14G N% P% K% REGRESSION OF Y14G (2 ) ON VARIATES FILE HQUYBOI 14/ 9/ 11:27 :PAGE TERM NO CONSTANT N% P% K% SOURCE COEFFICIENT -2284.148 4515.239 6554.451 38.24514 SS REGRESSION RESIDUAL TOTAL 0.8374954E+08 0.2672206E+08 0.1104716E+09 STD ERROR F-VALUE 1484.60 483.325 1922.42 2982.48 DF 41 44 2.367 87.274 11.625 0.000 MS 0.2791651E+08 651757.7 PRBF TOLERANCE 0.128 0.000 0.002 0.987 0.6571E-02 0.8852 0.7143 0.6517 F PRBF 42.833 0.0000 R-SQUARED - 75.8%, (ADJUSTED FOR D.F.- 74.0%) RECORD RECORD 42 HAS INFLUENTIAL X-VAL OBS= 45 HAS RESIDUAL= 2.6 SES OBS= 5363 5063 RES= RES= BOX PLOT OF STUDENTIZED RESIDUALS FROM LPLT= -2.570 NO.UPLT I -0 0.516 Kết cho ta phương trình tuyến tính với ba biến độc lập sau: Y = -2284,148 + 4515,239 N% + 6554,451 P% + 38,245 K% Với ñộ tin cậy phương trình hồi quy đạt mức cao với xác suất nhận Ho nhỏ (0,000) Chứng tỏ phương trình có ý nghĩa cho việc suy đốn từ yếu tố N, P, K suất Y 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……… https://sites.google.com/site/lophocphank57vnua/ 109 ... PAGE Phan tich ANOVA ket qua thi nghiem thiet ke kieu CRD VARIATE V003 NSUAT Tinh bang ta/ha LN SOURCE OF VARIATION DF SUMS OF MEAN F RATIO PROB ER SQUARES SQUARES LN =============================================================================... sơ đồ thí nghiệm phần mềm IRRISTAT 4.0 với cách thi? ??t kế khác ñể từ sơ ñồ thi? ??t kế máy tính người thi? ??t kế có thi? ??t kế hồn tồn khách quan, đảm bảo tính ngun tắc xác thi? ??t kế thí nghiệm, đồng... trình IRRISTAT để bạn đọc có số ý niệm sơ lược khái quát IRRISTAT Sách không nhằm giới thi? ??u tổng quát IRRISTAT mà tập trung vào phần: (1) Giới thi? ??u cách dùng IRRISTAT ver 4.0 (2) Dùng IRRISTAT