Kết quả nghiên cứu
4.5.1. Thăm dò mối quan hệ giữa khả năng Keo lai bị gãy ngang thân với D 1.3, Hvn, Hdc, Dtvà biến phân cành
Biến phân cành được ký hiệu là PC và nhận các giá trị 0 hoặc 1. Trong đó, 0- Cây ít phân cành hoặc phân cành nhánh nhỏ chỉ có một ngọn chính; 1- Cây phân cành lớn (ảnh 1) bao gồm những cây mà trong quá trình phát triển đến một chiều cao nhất định rồi phân ra làm 2 hoặc nhiều hơn 2 cành lớn bằng nhau, mỗi cành phát triển thành một ngọn cây.
Số liệu thu thập trên 2 ô tiêu chuẩn đem gộp chung lại, lọc bỏ những cây bị gãy ngang thân mà không còn nguyên hiện trạng, chỉ để lại những cây còn đủ chỉ tiêu về D1.3, Hvn, Hdc và Dt. Tiến hành xử lý số liệu, thăm dò yếu tố có liên quan đến sự gãy gập của Keo lai theo mô hình hồi qui logistic (trong đó, pi là khả năng Keo lại bị gãy ngang thân), kết quả cho ở bảng 4.9.
Bảng 4.9.Kết quả thăm dò yếu tố có liên quan đến khả năng gãy ngang thân Keo lai bằng mô hình hồi qui logistic (3.16)
Biến độc lập Tham số Sai tiêu chuẩn Sig.
D13 -1,021 0,399 0,010 Hvn 0,009 0,437 0,983 Hdc -0,327 0,356 0,359 Dt 5,554 1,791 0,002 Phân cành 3,232 1,409 0,022 Hằng số tự do -7,566 4,609 0,101
Kết quả tính toán ở bảng 4.9 cho thấy, các tham số bo, b2 (của biến Hvn), b3 (biến Hdc) không tồn tại do giá trị của Sig.>0,05, có nghĩa các nhân tố chiều cao vút ngọn, chiều cao dưới cành ảnh hưởng không rõ đến sự gãy gập của Keo lai. Các tham số b1 (của biến D1.3), b4 (biến Dt), b5(biến Phân cành) là tồn tại trong tổng thể do Sig.<0,05, điều này chứng tỏ các nhân tố đường kính ngang ngực, đường kính tán, phân cành nhánh lớn có ảnh hưởng rõ rệt đến khả năng gãy của Keo lai.
Sau khi loại bỏ các biến không có liên quan, tiến hành xử lý và lập hồi qui logistic với các biến còn lại cho kết quả ghi trong bảng 4.10.
Bảng 4.10. Tổng hợp các tham số trong phân tích hồi qui Logistic dạng Logit(pi) = bo + b1.D1.3+ b2.Dt+ b3.PC
Các biến Tham số Sai tiêu chuẩn Sig. Exp(B)
D13 -1,148 0,340 0,001 0,317
Dt 5,774 1,812 0,001 321,890
PC 3,907 1,342 0,004 49,742
Hằng số tự do -9,766 3,832 0,011 0,000
Phương trình hồi qui cụ thể có dạng như sau:
Logit (pi) = -9,766 - 1,148.D1.3 + 5,774.Dt+ 3,907.PC (4.8) Với R2 = 0,734.
Từ phương trình 4.8 cho thấy, hệ số của D1.3là một số âm, có nghĩa cây có đường kính lớn thì khả năng bị gãy sẽ giảm đi. Hệ số của biến Dtvà phân cành là số dương, có nghĩa cây có đường kính tán và phân cành lớn thì khả năng bị gãy cũng lớn hơn. Căn cứ vào tính chất của hồi qui logistic, chúng ta có thể lượng hoá được khả năng bị gãy của từng cây. Cụ thể, từ các giá trị tương ứng với từng biến ở cột cuối cùng của bảng 4.10 cho biết, nếu 2 cây có cùng DTvà D1.3thì cây phân cành lớn có nguy cơ bị gãy lớn hơn gấp exp[3,907(1 - 0)]49,74 lần so với cây
phân cành nhánh nhỏ (hay không phân nhánh); nếu 2 cây có D1.3 và tình hình phân nhánh như nhau thì cây có đường kính tán 4m có khả năng bị gãy tăng lên là exp[5,775(4-3)]321,89 lần so với cây có đường kính tán 3m; Còn một cây có đường kính ngang ngực tăng lên 1cm (D1.3 = 15cm) có nguy cơ bị gãy giảm đi exp[-1,148(15-14)]0,317 lần (31,7%) so với cây khác (D1.3= 14cm) nếu chúng có cùng Dtvà giống nhau về tình hình phân cành.
Kết quả thăm dò mối quan hệ giữa khả năng Keo lai bị gãy ngang thân với một số nhân tố sinh trưởng bằng mô hình hồi qui logistic còn cho ta bảng phân loại những cây sống, cây gãy ngoài thực tế và theo lý thuyết như bảng 4.11.
Bảng 4.11.Kết quả phân loại
Phân loại Cây bình thường Cây gãy Tỷ lệ phân loại đúng (%)
Cây bình thường 176 1 99,4
Cây gãy 3 10 76,9
Tỷ lệ phân loại đúng tính chung cho cả 2 mẫu (%) 97,9
Từ bảng 4.11 cho thấy, trong số 177 cây không bị gãy ngoài thực tế thì có 1 cây nhảy sang bị gãy theo lý thuyết, kết quả phân loại đúng 99,4%; trong số 13 cây bị gãy ngoài thực tế thì có 3 cây nhảy sang không bị gãy theo lý thuyết, kết quả phân loại dúng 76,9%. Tính chung, phân loại đúng đạt 97,9%.
4.5.2. Thăm dò mối quan hệ giữa khả năng Keo lai bị gãy ngang thân vớiD1.3, Hvn, Dt