Ứng dụng tin học trong công nghệ sinh học (2)

3/6/2011 Bố trí hồn tồn ngẫu nhiên Completely Randomized Design (CRD) CRD kiểu bố trí mà nghiệm thức bố trí cách hồn tồn ngẫu nhiên Do vậy, đơn vị thí nghiệm có hội nhận nghiệm thức  Trong CRD, nghiệm thức đ ợc bố trí lập lại tất số liệu đ ợc ghi nhận điều kiện nói hồn tồn giống Ví dụ nh tất nghiệm thức đ ợc bố trí lập lại m t khu đất, điều kiện môi tr ờng nuôi cấy, tuổi cây,  Bất kỳ khác đơn vị thí nghiệm nhận m t nghiệm thức đ ợc xem nh sai số thí nghiệm  Do vậy, CRD thích hợp cho thí nghiệm có đơn vị thí nghiệm đồng nhất, nh thí nghiệm phịng Lab, nơi mà ảnh h ởng yếu tố môi tr ờng đ ợc kiểm sốt dễ dàng Vì vậy, mà số lần lập lại nghiệm thức khác  Thí nghiệm ngồi đồng, có nhiều biến đ ng l n lơ thí nghiệm, CRD đ ợc sử dụng Thí dụ: Bố trí thí nghiệm gồm nghiệm thức (A,B,C,D,E) v i lần lặp lại đ ợc xếp theo kiểu hoàn toàn ngẫu nhiên Các b c thực B c 1: Xác định số đơn vị thí nghiệm (n) Có t = 5, r = n = Số nghiệm thức (t) x số lần lặp lại (r) = t x r = x = 20 B c 2: Phân lơ thí nghiệm (từ đến n) Trong thí nghiệm này, số lơ thí nghiệm đ ợc tính từ 1- 20 lô B c 3: Phân nghiệm thức vào lơ thí nghiệm ph ơng pháp sau: Dựa theo bảng phân bố ngẫu nhiên Lấy ngẫu nhiên m t số đến 19 số kế theo Cho nghiệm thức vào theo A- số đầu B- số kế E- số kế C- số kế D- số kế Xếp theo thứ tự 1- 20 => Suy vị trí lơ Dựa theo phiếu rút thăm + Chuẩn bị giấy đủ t nghiệm thức (5) + Làm lần lặp lại + Tr n ngẫu nhiên + Rút thăm theo lơ từ 1- 20 (đã bố trí sẵn) 3/6/2011 10 11 12 13 14 15 16 17 18 19 20 Ví dụ: Khảo sát ảnh h ởng m t số thuốc diệt cỏ lên suất giống lúa IR1 thu đ ợc kết sau: A E D D 10 Bây giờ, bố trí Randomized A 14 CÁC B D E 16 C 19 E B 12 B C 18 A Dol-Mix (1kg) Dol-Mix (2kg) DDT + -BHC Azodrin Dimecron-Boom Dimecron-Knap Đối chứng D B Dol-Mix (1kg) Dol-Mix (2kg) DDT + -BHC Azodrin Dimecron-Boom Dimecron-Knap Đối chứng 2.537 3.366 2.536 2.387 1.997 1.796 1.401 2.069 2.591 2.459 2.453 1.679 1.704 1.516 2.104 2.211 2.827 1.556 1.649 1.904 1.270 Tổng nghiệm thức (Ti) 1.797 2.544 2.385 2.116 1.859 1.320 1.077 Tổng chung (G) Trung bình chung CÁC B B 8.507 10.712 10.207 8.512 7.184 6.724 5.264 Trung bình nghiệm thức 2.127 2.678 2.552 2.128 1.796 1.681 1.316 57 57.110 2.040 2.040 C PHÂN TÍCH c 3: Tính bậc tự (df) + Đ tự tổng c ng (dfT0) = r x t - + Đ tự nghiệm thức (dfT) = t - + Đ tự sai số (dfE) = t (r- 1) Đ tự sai số đ ợc tính: dfE = dfT0 - dfT 2.069 2.591 2.459 2.453 1.679 1.704 1.516 2.104 2.211 2.827 1.556 1.649 1.904 1.270 1.797 2.544 2.385 2.116 1.859 1.320 1.077 Nên sử dụng loại thuốc diệt cỏ để đạt suất lúa l n nhất? (α=0.05) CÁC B C PHÂN TÍCH Năng suất (kg/ha) 2.537 3.366 2.536 2.387 1.997 1.796 1.401 20 B c 1: Tính tổng Ti nghiệm thức, tổng tồn b G giá trị trung bình nghiệm thức Nghiệm thức Năng suất (kg/ha) 15 A 17 Nghiệm thức C E 11 C 13 B B Nguồn biến thiên C PHÂN TÍCH c 2: Đ a bảng phân tích ph ơng sai Độ tự (df) Tổng bình phương (SS) Trung bình bình phương (MS) F bảng F tính 5% 1% Nghiệm thức (T) Sai số (E) Tổng cộng (T0) CÁC B C PHÂN TÍCH B c 4: Tính số hiệu chỉnh (CF) loại tổng bình ph ơng (SS) CF  G2 n SST0   xij2  CF r i 1 SST   t i 1 Ti  CF r SSE  SST0  SST 3/6/2011 CÁC B C PHÂN TÍCH CÁC B B c 5: Tính trung bình bình ph ơng (MS) cho nguồn biến đ ng C PHÂN TÍCH B c 6: Tính trị số Ftính để kiểm tra mức ý nghĩa khác công thức Ftn  MST  SST dfT SSE MSE  dfE CÁC B B Nghiệm thức (T) Trung bình bình phương (MS) 5587174 931196 Sai số (E) 21 1990238 Tổng cộng (T0) 27 7577412 Nguồn biến thiên Độ tự (df) Tổng bình phương (SS) c 7: Tìm giá trị F(α,dfT, dfE) bảng F CÁC B C PHÂN TÍCH B c 8: Điền giá trị tính đ ợc vào bảng phân tích ph ơng sai F bảng F tính 9,83** 5% 1% 2,57 3,81 94773 Chú ý: - Bất F tỏ khơng có ý nghĩa, ng ời nghiên cứu cần xem lại đ l n sai số thí nghiệm khác trung bình cơng thức Nếu giá trị l n nên làm lại thí nghiệm tìm cách làm cẩn thận để giảm sai số thí nghiệm để thấy rõ khác công thức Mặc khác, hai giá trị nhỏ chứng tỏ khác công thức nhỏ nên không cần thí nghiệm thêm - Kiểm tra F cho biết sai khác chung công thức không cho biết khác hay không đôi công thức, muốn biết phải so sánh trung bình đơi m t MST MSE C PHÂN TÍCH B c 9: So sánh giá trị Ftn F tra bảng mức ý nghĩa khác - Nếu Ftn>Ftra(0,01) α = 0,01 có khác l n nghiệm thức thí nghiệm (ký hiệu **) - Nếu Ftra(0,01)≥Ftn>Ftra(0,05) có khác nghiệm thức mức α=0,05 (ký hiệu *) - Nếu Ftn≤Ftra(0,05) khác cơng thức khơng có ý nghĩa (ký hiệu ns) CÁC B C PHÂN TÍCH B c 10: Tính trung bình tồn b hệ số biến đ ng CV% X G n CV %  MSE x100% X 3/6/2011 CV% đ xác so sánh nghiệm thức m t giá trị tin cậy thí nghiệm Nó biểu thị sai số thí nghiệm thơng qua phần trăm trung bình, vậy, giá trị CV cao thí nghiệm có giá trị tin cậy CV% biến đ ng tuỳ theo loại thí nghiệm, loại trồng, đặc tính ghi nhận Giá trị CV% chấp nhận đ ợc v i: - thí nghiệm phịng - thí nghiệm chậu, vại, nhà l ≤ 1% i ≤ 5% - suất lúa 6- 8% - thí nghiệm phân bón 10- 12% - thí nghiệm thuốc trừ sâu - thí nghiệm thuốc trừ cỏ 13- 15% 13- 15% V i kết F ta biết nghiệm thức khác biệt có ý nghĩa hay không, mà so sánh nghiệm thức v i Để xác định đ ợc ta sử dụng ph ơng pháp so sánh cặp Có hai ph ơng pháp đ ợc sử dụng phổ biến là: + Least significant difference (LSD) + Duncan’s multiple range test (DMRT) Thực theo b Least significant difference Test (LSD) LSD ph ơng pháp đơn giản đ ợc sử dụng phổ biến để so sánh cặp Nó đ ợc sử dụng để xác định khác biệt có ý nghĩa hay khơng có ý nghĩa cặp nghiệm thức so sánh Tuy nhiên khơng thích hợp để so sánh cho tất trung bình mà số nghiệm thức nhiều LSD nên sử dụng F test cho thấy có khác biệt có ý nghĩa số nghiệm thức không sáu (6) B c nh sau: c 11: Tính giá trị LSDα * Tr ờng hợp lần lặp lại nhau: LSD  t( ,dfE ) 2MSE r •Tr ờng hợp lần lặp lại khác nhau: 1 1 LSD  t( ,dfE )   ,  MSE  ni ni  So sánh giá trị trung bình cặp nghiệm thức so với LSD để biết mức độ khác biệt có ý nghĩa - Hiển thị kết sau Nghiệm thức Năng suất trung bình (Kg/ha) So sánh cặp trung bình nghiệm thức - Sắp xếp số trung bình nghiệm thức theo thứ tự giảm dần tăng dần - Nhận diện nhóm tất trung bình nghiệm thức khơng khác biệt có ý nghĩa v i T2 T3 T4 T1 T5 T6 T7 Suy ra: Nghiệm thức T2 T3 T4 T1 T5 T6 T7 2.678 2.552 2.128 2.127 1.796 1.681 1.316 a b c d Năng suất trung bình (Kg/ha) 2.678 a 2.552 ab 2.128 bc 2.127 bc 1.796 c 1.681 cd 1.316 d 3/6/2011 Bảng lập lại sau: Nghiệm thức Năng suất trung bình (Kg/ha) Xếp hạng Dol-Mix (1kg) (T1) Dol-Mix (2kg) (T2) DDT + -BHC (T3) Azodrin (T4) Dimecron-Boom (T5) Dimecron-Knap (T6) Đối chứng (T7) 2.127 bc 2.678 a 2.552 ab 2.128 bc 1.796 c 1.681 cd 1.316 d ... số thí nghiệm để thấy rõ khác công thức Mặc khác, hai giá trị nhỏ chứng tỏ khác cơng thức q nhỏ nên khơng cần thí nghiệm thêm - Kiểm tra F cho biết sai khác chung công thức không cho biết khác... + -BHC Azodrin Dimecron-Boom Dimecron-Knap Đối chứng D B Dol-Mix (1kg) Dol-Mix (2kg) DDT + -BHC Azodrin Dimecron-Boom Dimecron-Knap Đối chứng 2.537 3.366 2.536 2.387 1.997 1.796 1.401 2.069... xác so sánh nghiệm thức m t giá trị tin cậy thí nghiệm Nó biểu thị sai số thí nghiệm thơng qua phần trăm trung bình, vậy, giá trị CV cao thí nghiệm có giá trị tin cậy CV% biến đ ng tuỳ theo loại