Ảnh hưởng của các nhân tố sinh thái đến quần thể trà hoa vàng (camellia inusitata orel, curry luu) tại vườn quốc gia bidoup núi bà, tỉnh lâm đồng

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƢỜNG ĐẠI HỌC LÂM NGHIỆP ĐINH VĂN TÝ ẢNH HƢỞNG CỦA CÁC NHÂN TỐ SINH THÁI ĐẾN QUẦN THỂ TRÀ HOA VÀNG (Camellia inusitata Orel, Curry & Luu), TẠI VƢỜN QUỐC GIA BIDOUP NÚI BÀ, TỈNH LÂM DỒNG LUẬN VĂN THẠC SỸ LÂM NGHIỆP Đồng Nai, tháng 12 năm 2018 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƢỜNG ĐẠI HỌC LÂM NGHIỆP ĐINH VĂN TÝ ẢNH HƢỞNG CỦA CÁC NHÂN TỐ SINH THÁI ĐẾN QUẦN THỂ TRÀ HOA VÀNG (Camellia inusitata Orel, Curry & Luu) TẠI VƢỜN QUỐC GIA BIDOUP NÚI BÀ, TỈNH LÂM DỒNG CHUYÊN NGÀNH LÂM HOC MÃ NGÀNH: 60.62.60 LUẬN VĂN THẠC SỸ LÂM NGHIỆP NGƢỜI HƢỚNG DẪN KHOA HỌC TS Phạm Văn Hƣờng Đồng Nai, tháng 12 năm 2018 LỜI CAM ĐOAN Tơi cam đoan, cơng trình nghiên cứu riêng Các số liệu, kết nêu luận văn trung thực chƣa đƣợc cơng bố cơng trình nghiên cứu khác Đinh Văn Tý i LỜI NÓI ĐẦU Luận văn đƣợc hồn thành theo chƣơng trình đào tạo cao học khóa K24B Phân hiệu trƣờng Đại học Lâm nghiệp Việt Nam Để hồn thành khóa học luận văn này, xin gửi lời cảm ơn chân thành đến Ban giám hiệu trƣờng Đại học Lâm nghiệp, khoa Đào tạo sau đại học, thầy cô giáo giảng dạy thời gian qua Tôi xin đặc biệt cảm ơn TS Phạm Văn Hƣờng ngƣời trực tiếp hƣớng dẫn khoa học, truyền đạt kiến thức, kinh nghiệm q báu, tận tình giúp đỡ tơi hồn thành luận văn Nhân dịp tơi xin tỏ lịng biết ơn tới Ban lãnh đạo Vƣờn quốc gia Bidoup – Núi Bà, nhà khoa học, toàn thể đồng nghiệp bạn bè có đóng góp động viên giúp đỡ tơi hồn thành luận văn Mặc dù làm việc với tất nỗ lực, nhƣng hạn chế trình độ thời gian, nên luận văn khơng tránh khỏi thiếu sót định Tôi mong nhận đƣợc ý kiến đóng góp xây dựng quý báu nhà khoa học bạn bè đồng nghiệp Tôi xin chân thành cảm ơn./ Đồng Nai, ngày tháng 12 năm 2018 Tác giả Đinh Văn Tý ii TÓM TẮT Đề tài luận văn: "Ảnh hƣởng nhân tố sinh thái đến quần thể Trà hoa vàng (camellia inusitata Orel, Curry& Luu), Vƣờn Quốc Gia Bidoup – Núi Bà, tỉnh Lâm Đồng" đƣợc tiến hành nghiên cứu đối tƣợng ảnh hƣởng yếu tố sinh thái chủ đạo đến quần thể Trà hoa vàng phân bố tự nhiên đai độ cao khác nhau, thuộc kiểu phụ rừng lùn núi, Vƣờn Quốc Gia Bidoup – Núi bà Bằng phƣơng pháp luận dựa lý thuyết sinh thái, kết điều tra OTC đai cao khác nhau, nơi Trà hoa vàng phân bố để mô tả trạng thái rừng nơi Trà hoa vàng phân bố Đồng thời tiến hành thiết lập 150 OTC thứ cấp (100m ), để tiến hành xem xét ảnh hƣởng yếu tố sinh thái đến đặc điểm tầng tái sinh Trà hoa vàng Kết hợp với phƣơng pháp xử lý phân tích số liệu khác Kết nghiên cứu , Vùng phân bố Trà hoa vàng đặc hữu Lâm Đồng có diện tích cƣ trú 89 ha, dƣới tán rừng rộng thƣờng xanh khu vực rừng lùn đỉnh núi, độ cao từ 1500m đến 1900 m Vật hậu học Trà hoa vàng, kỳ chồi non Trà hoa vàng diễn liên tục từ cuối tháng đến đầu tháng 9; kỳ nọn chuyển màu thành già trung bình xuất từ gữa tháng 9, kết thúc khoảng cuối tháng hàng năm Kỳ nụ hoa xuất vào đầu tháng kết thúc vào cuối tháng 1, Kỳ nụ, hoa kết diễn liên tục từ tháng 12 đến tháng 5; Thời gian rụng bao gồm thời gian già chín tách hạt để phát tán, thời gian kéo dài trung bình 39 ngày Đặc điểm tổ thành lồi đai độ cao nơi Trà hoa vàng phân bố, đai độ cao 1500m có tổng cộng 71 loài xuất hiện, đai độ cao 1700 m , số lƣợng loài xuất 56 loài Ở độ cao 1900m có 42 lồi xuất Phân bố khoảng cách phù hợp để mô phân bố N-D tầng gỗ lâm phần đai độ cao có Trà hoa vàng phân bố Phần lớn số tập trung cỡ đƣờng kính 11-12 cm Phân bố thực nghiệm N/H cho thấy, hầu hết rừng tập trung cỡ chiều cao 7-8 m Trà hoa vàng phân bố đai cao 1700 m cao với mật độ trung bình 947 cây/ha, thấp đai độ cao 1900m Trà hoa vàng cấp tuổi khác có mật độ phân bố khác Nhìn chung Trà hoa vàng trƣởng thành (cấp tuổi 5, với Hvn> 200cm) chiếm tỷ lệ mật độ cao so với cấp tuổi khác Trà hoa vàng cấp tuổi khác mật độ khác nhau, xu hƣớng chung cấp tuổi tăng mật độ iii Trà hoa vàng có xung hƣớng giảm dần, nhiên quy luật phân bố không rõ nét, khác đai độ cao Phân bố Trà hoa vàng theo nguồn gốc đai cao có nguồn gốc từ chồi phân bố cao có nguồn gốc từ hạt, phân bố Trà hoa vàng theo phẩm chất tốt, trung bình phân bố tăng dần từ đai cao 1500 m đến 1700 m giảm dần xuống 1900 m Trà hoa vàng xuất nhiều hƣớng Đơng hƣớng Đơng bắc, hƣớng lại lần lƣợt theo thứ tự từ cao đến thấp hƣớng Nam > Đông Nam > Đơng Bắc Trong đó, mật độ phân bố nhiều cấp tuổi hai hƣớng chủ yếu hƣớng Đơng Nam hƣớng Nam, hai hƣớng cịn lại phân bố thấp Cấp tuổi Trà hoa vàng o o xuất nhiều nơi có độ dốc từ 19 đến 20 nơi xuất nơi có độ dốc 0 từ 11 đến 12 , xuất Trà hoa vàng độ dốc lại lần lƣợt theo thứ tự từ cao đến thấp 20 > 17-18 > 13-14 > 15-16 Tối ƣu sinh thái Trà hoa vàng tuổi dao động độ pH từ 5.87 đến 6.33, giới hạn tối ƣu sinh thái cấp tuổi độ pH từ 5.07±6.67 đến 5.67±6.99, biên độ sinh thái tối thích độ pH từ 4.26±7.48 đến 5.02±7.64 Trà hoa vàng xuất khu vực có độ pH từ 4.26 đến 7.94 Tối ƣu sinh thái cấp tuổi độ ẩm 85.12%, cấp tuổi 2, 3, 4, tƣơng ứng 75.18%, 76.64%, 72.15%, 71.23% Tối ƣu sinh thái tàn che cấp tuổi độ tàn che 80.21% cao tất cấp tuổi 2, 3, 4, tƣơng ứng 65.11%, 63.65%, 62.01% 61.3%, cấp tuổi tăng lên tối ƣu sinh thái tàn che giảm dần, đồng thời biên độ sinh thái tàn che tăng lên, Trà hoa vàng thích nghi phân bố nơi có độ tàn che thấp, lồi chiu bóng Độ tàn che thảm cỏ tăng mật độ Trà hoa vàng giảm Trà hoa vàng tái sinh chiu ảnh hƣởng chiều cao thảm cỏ, chiều cao thảm cỏ tăng lên mật độ tái sinh giảm dần Mật độ phân bố Trà hoa vàng phù hợp với đƣờng cong hồi quy, có dạng phân bố giảm, độ đầy thảm cỏ tăng mật độ Trà hoa vàng giảm Khi độ dầy thảm khô tăng mật độ Trà hoa vàng giảm, điều có nghĩa nơi có độ dầy thảm mục cao thi tái sinh khó để rễ tiếp xúc với đất khả tái sinh thấp Khối lƣợng thảm khơ tăng mật độ Trà hoa vàng giảm iv MỤC LỤC LỜI CAM ĐOAN i LỜI NÓI ĐẦU ii TÓM TẮT iii MỤC LỤC v DANH LỤC TỪ VIẾT TẮT viii DANH MỤC BẢNG BIỂU ix DANH LỤC HÌNH ẢNH, ĐỒ THỊ x ĐẶT VẤN ĐỀ Chƣơng TỔNG QUAN NGHIÊN CỨU 1.1 Các nghiên cứu đặc điểm lâm học, cấu trúc rừng 1.1.1 Nghiên cứu cấu trúc rừng 1.1.2 Nghiên cứu tái sinh rừng 1.2 Các nghiên cứu chi Trà Cmellia 1.2.1 Đặc điểm phân loại 1.2.2 Đặc điểm hình thái 1.3 Các nghiên cứu Trà hoa vàng 1.3.1 Các nghiên cứu giới 1.3.2 Các nghiên cứu Việt Nam 1.3.3 Các nghiên cứu VQG Bi doup – Núi bà 10 1.4 Các nghiên cứu ảnh hƣởng yếu tố sinh thái đến Trà hoa vàng 11 1.4.1 Ảnh hƣởng yếu tố phi sinh vật 11 1.4.2 Ảnh hƣởng yếu tố sinh vật 11 1.5 Thảo luận chung 11 Chƣơng ĐỐI TƢỢNG, MỤC TIÊU, NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 13 2.1 Đối tƣợng, phạm vi nghiên cứu 13 2.1.1 Đối tƣợng nghiên cứu 13 2.1.2 Phạm vi nghiên cứu 13 2.2 Mục tiêunghiên cứu 13 2.2.1 Mục tiêu tổng quát 13 v 2.2.2 Mục tiêu cụ thể 13 2.2.3 Ý nghĩa đề tài 14 2.3 Nội dung – phƣơng pháp nghiên cứu 14 2.3.1 Nội dung nghiên cứu 14 2.3.2 Phƣơng pháp nghiên cứu 15 2.3.2.1 Quan điểm phƣơng pháp luận 15 2.3.2.2 Phƣơng pháp thu thập số liệu 16 Chƣơng ĐẶC ĐIỂM ĐIỀU KIỆN TỰ NHIÊN KHU VỰC NGHIÊN CƢU 24 3.1 Khái quát điều kiện tự nhiên 24 3.2 Vị trí địa lý, địa hình 24 3.2.1 Vị trí địa lý 24 2.2 Địa hình 25 3.3 Chế độ khí hậu - thủy văn 25 Chƣơng KẾT QUẢ VÀ THẢO LUẬN 28 4.1 Đặc điểm phân loại, phân bố vật hậu học Trà hoa vàng 28 4.1.1 Đặc điểm phân loại 28 4.1.2 Đặc điểm phân bố 28 4.1.3 Đặc điểm vật hậu học 31 4.2 Đặc điểm lâm học trạng thái rừng nơi Trà hoa vàng phân bố 32 4.2.1 Đặc điểm tổ thành loài 32 4.2.1.1 Đai độ cao 1500m 32 4.2.1.2 Đai độ cao 1700m 33 4.2.1.3 Đai độ cao 1900m 33 4.2.2 Đặc điểm cấu trúc phân bố 35 4.2.2.1 Phân bố N-D 35 4.2.2.2 Qui luật phân bố số theo chiều cao (N/Hvn) 36 4.3 Đặc điểm quần thể Trà hoa vàng 37 4.3.1 Đặc điểm phân bố theo cấp tuổi 37 4.3.2 Đặc điểm phân bố mật độ Trà hoa vàng theo nguồn gốc 40 4.3.3 Đặc điểm phân bố theo phẩm chất 41 4.4 Ảnh hƣởng nhân tố sinh thái đến phân bố Trà hoa vàng 43 vi 4.4.1.Ảnh hƣởng yếu tố địa hình – địa mạo 43 4.4.1.1 Ảnh hƣởng hƣớng phơi 43 4.4.1.2 Ảnh hƣởng độ dốc 45 4.4.2 Ảnh hƣởng yếu tố thổ nhƣỡng 48 4.4.2.1 Ảnh hƣởng độ pH độ ẩm tầng đất mặt 48 4.4.2.2 Độ ẩm tầng đất mặt 49 4.4.2.3 Ảnh hƣởng độ tàn che tán rừng 51 4.4.3 Ảnh hƣởng thảm cỏ 53 4.4.3.1 Ảnh hƣởng độ tàn che thảm cỏ 53 4.4.3.2 Ảnh hƣởng chiều cao thảm cỏ 53 4.4.3.3 Ảnh hƣởng độ đầy thảm cỏ 54 4.4.4 Ảnh hƣởng thảm khô – thảm mục 55 4.4.4.1 Ảnh hƣởng độ dày thảm khô – thảm mục 55 4.4.4.2 Ảnh hƣởng khối lƣợng thảm khô – thảm mục 56 4.5 Giải pháp bảo tồn phát triển bền vững quần thể Trà hoa vàng 56 KẾT LUẬN, TỒN TẠI, KIẾN NGHỊ 58 TÀI LIỆU THAM KHẢO CHÍNH 62 MỘT SỐ HÌNH ẢNH a PHỤ LỤC c vii DANH LỤC TỪ VIẾT TẮT [1] Thứ tự tài liệu tham khảo Lk Loài khác Cn Cồng nhám Lx Luống xƣơng D1,3 Chiều cao ngang ngực (cm) N Mật độ (cây/ha) Dd Dung đen ODB Ô dạng Dg Dẻ gai ÔTC Ô tiêu chuẩn Dl Dung láng Sd Sồi poilane Dr Dẻ rừng Ss Sơn trâm speng Đq Đỗ quyên St Sơn trâm hẹp HDC Chiều cao dƣới cành (m) Su Sụ Hvn Chiều cao vút (m) Tt Thơng tre Kt Kha thụ nhím VQG Vƣờn quốc gia Trv Trà hoa vàng viii 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 1500 1500 1500 1500 1500 1500 1500 1500 1500 1500 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 5.0 6.0 6.0 5.5 5.5 5.5 5.0 6.0 4.0 4.0 6.0 4.0 4.0 4.5 5.0 6.0 6.0 5.5 5.5 5.5 5.0 6.0 4.0 4.0 6.0 4.0 4.0 4.5 5.0 5.5 5.5 5.0 5.5 4.5 4.5 4.0 4.0 4.0 6.0 6.0 6.0 4.0 5.5 4.5 6.0 77.00 79.00 81.00 85.00 75.00 79.00 1.00 80.00 60.00 60.00 76.00 75.00 85.00 65.00 70.00 79.00 82.00 74.00 82.00 79.00 69.00 75.00 65.00 60.00 80.00 77.00 69.00 78.00 70.00 75.00 75.00 80.00 75.00 65.00 65.00 60.00 68.00 72.00 81.00 85.00 86.00 67.00 83.00 65.00 80.00 85 87 86 91 83 84 79 90 69 66 81 77 87 84 78 87 83 85 83 84 79 90 69 66 89 70 69 74 78 81 75 79 86 75 74 68 69 70 89 91 91 69 78 74 89 1 0 1 0 0 1 1 1 0 1 0 0 0 0 1 bb 0 1 0 0 1 0 1 0 0 0 1 1 0 1 0 0 1 0 0 1 0 1 0 0 0 1 1 0 1 0 0 1 0 0 1 1 0 1 0 0 0 1 1 0 1 0 0 1 0 0 1 0 1 0 0 0 1 1 0 1 0 1 1 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1700 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 5.5 5.0 6.0 4.5 4.0 4.5 4.5 6.0 5.0 6.0 5.5 5.0 6.0 4.5 4.0 4.5 4.5 6.0 5.0 6.0 5.5 5.0 6.0 4.5 4.0 4.5 4.5 6.0 5.0 6.0 5.5 5.0 6.0 4.5 4.0 4.5 4.5 6.0 5.0 5.5 5.5 6.0 5.0 5.0 6.0 75.00 70.00 85.00 65.00 70.00 65.00 75.00 88.00 70.00 82.00 75.00 70.00 80.00 65.00 60.00 65.00 65.00 75.00 70.00 81.00 77.00 76.00 84.00 65.00 55.00 65.00 65.00 83.00 70.00 80.00 75.00 70.00 80.00 65.00 66.00 54.00 65.00 80.00 56.00 75.00 79.00 72.00 70.00 70.00 80.00 84 79 88 74 69 74 74 87 79 89 84 79 88 74 69 74 74 87 79 75 79 79 85 74 64 74 74 87 79 89 84 79 78 74 69 74 78 87 77 82 82 73 77 79 85 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 cc 0 0 1 1 0 0 1 0 1 0 0 0 0 0 0 1 0 0 1 1 0 0 1 0 1 0 0 0 0 1 0 1 0 0 1 1 0 0 1 0 1 0 0 0 0 1 0 1 0 0 1 1 0 0 1 0 1 0 0 0 0 1 0 1 1 1 1 0 1 1 1 0 1 1 1 0 0 0 1 1 1 0 1 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 1900 5.5 5.5 6.0 5.0 5.0 6.0 5.5 5.5 6.0 5.0 5.0 6.0 5.5 5.5 6.0 5.0 5.0 6.0 3.5 4.0 5.2 75.00 75.00 65.00 70.00 70.00 80.00 75.00 75.00 65.00 65.00 50.00 79.00 19.00 75.00 82.00 67.00 72.00 77.00 65.00 69.00 72.00 88 83 82 77 79 90 88 85 85 77 79 84 88 72 88 77 79 74 76 81 85 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 1 0 1 1 1 1 0 0 1 0 1 1 1 0 0 Phụ lục Phƣơng trình hồi quy độ pH cấp tuổi Trà hoa vàng Logistic Regression - Tr1 Dependent variable: Tr1 Factors: pH pH^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -47.0159 30.0902 pH 14.7095 9.77157 2.44486E6 pH^2 -1.16193 0.790961 0.312882 Analysis of Deviance Source Deviance Df P-Value Model 3.69336 0.1578 Residual 187.261 147 0.0139 Total (corr.) 190.954 149 Percentage of deviance explained by model = 1.93416 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor ChiDf P-Value Squared pH 2.62722 0.1050 pH^2 2.47715 0.1155 Residual Analysis dd 1 0 1 1 1 0 0 1 0 1 1 1 0 0 n MSE MAE MAPE ME Estimation 150 0.0491995 0.444558 Validation 0.000243425 MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between Tr1 and independent variable(s) The equation of the fitted model is Tr1 = exp(-47.0159 + 14.7095*pH - 1.16193*pH^2)/(1+exp(-47.0159 + 14.7095*pH 1.16193*pH^2)) Logistic Regression - Tr2 Dependent variable: Tr2 Factors: pH pH^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -47.5629 28.7592 pH 15.2448 9.38896 4.17558E6 pH^2 -1.23202 0.763775 0.291704 Analysis of Deviance Source Deviance Df P-Value Model 3.11228 0.2109 Residual 189.198 147 0.0108 Total (corr.) 192.311 149 Percentage of deviance explained by model = 1.61836 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor ChiDf P-Value Squared pH 3.00535 0.0830 pH^2 2.94946 0.0859 Residual Analysis Estimation Validation n 150 MSE 0.0500803 MAE 0.446779 MAPE ME 0.000285065 MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between Tr2 and independent variable(s) The equation of the fitted model is Tr2 = exp(-47.5629 + 15.2448*pH - 1.23202*pH^2)/(1+exp(-47.5629 + 15.2448*pH 1.23202*pH^2)) ee Logistic Regression - Tr3 Dependent variable: Tr3 Factors: pH pH^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -29.0773 25.4046 pH 9.17907 8.31242 9692.15 pH^2 -0.731416 0.677459 0.481227 Analysis of Deviance Source Deviance Df P-Value Model 1.61941 0.4450 Residual 199.445 147 0.0026 Total (corr.) 201.065 149 Percentage of deviance explained by model = 0.805416 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor ChiDf P-Value Squared pH 1.29167 0.2557 pH^2 1.23076 0.2673 Residual Analysis Estimation Validation n 150 MSE 0.0570382 MAE 0.473977 MAPE ME 0.000444527 MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between Tr3 and independent variable(s) The equation of the fitted model is Tr3 = exp(-29.0773 + 9.17907*pH - 0.731416*pH^2)/(1+exp(-29.0773 + 9.17907*pH 0.731416*pH^2)) Logistic Regression - Tr4 Dependent variable: Tr4 Factors: pH pH^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -28.9892 24.5606 pH 9.21269 8.04393 10023.5 pH^2 -0.734063 0.656114 0.479955 Analysis of Deviance Source Deviance Df P-Value Model 1.73835 0.4193 Residual 204.041 147 0.0013 Total (corr.) 205.779 149 ff Percentage of deviance explained by model = 0.844765 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor ChiDf P-Value Squared pH 1.37608 0.2408 pH^2 1.30971 0.2524 Residual Analysis Estimation Validation n 150 MSE 0.0605403 MAE 0.487842 MAPE ME -0.000932106 MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between Tr4 and independent variable(s) The equation of the fitted model is Tr4 = exp(-28.9892 + 9.21269*pH - 0.734063*pH^2)/(1+exp(-28.9892 + 9.21269*pH 0.734063*pH^2)) Logistic Regression - Tr5 Dependent variable: Tr5 Factors: pH pH^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -27.4159 24.5095 pH 8.71751 8.0296 6108.95 pH^2 -0.696 0.655126 0.498576 Analysis of Deviance Source Deviance Df P-Value Model 1.50789 0.4705 Residual 203.762 147 0.0014 Total (corr.) 205.27 149 Percentage of deviance explained by model = 0.734591 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor ChiDf P-Value Squared pH 1.23287 0.2668 pH^2 1.17784 0.2778 Residual Analysis Estimation Validation n 150 MSE 0.0603256 MAE 0.486899 MAPE ME 0.000732105 gg MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between Tr5 and independent variable(s) The equation of the fitted model is Tr5 = exp(-27.4159 + 8.71751*pH - 0.696*pH^2)/(1+exp(-27.4159 + 8.71751*pH - 0.696*pH^2)) Phụ lục 01: Phƣơng trình hồi quy độ tàn che cấp tuổi Trà hoa vàng Logistic Regression - TrC1 Dependent variable: TrC1 Factors: TC TC^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -19.8887 15.9301 TC 0.486633 0.436593 1.62683 TC^2 -0.00303347 0.00297281 0.996971 Analysis of Deviance Source Deviance Df P-Value Model 7.31336 0.0258 Residual 183.641 147 0.0217 Total (corr.) 190.954 149 Percentage of deviance explained by model = 3.8299 Adjusted percentage = 0.68779 Likelihood Ratio Tests Factor Chi-Squared Df P-Value TC 1.96245 0.1612 TC^2 1.31419 0.2516 Residual Analysis Estimation n MSE MAE MAPE ME MPE Vali dati on 150 0.0481856 0.445134 -0.000371137 The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C1 and independent variable(s) The equation of the fitted model is TrC1 = exp(-19.8887 + 0.486633*TC - 0.00303347*TC^2)/(1+exp(-19.8887 + 0.486633*TC - 0.00303347*TC^2)) Logistic Regression - TrC2 hh Dependent variable: TrC2 Factors: TC TC^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -4.97964 5.15113 TC 0.137607 0.147235 1.14752 TC^2 -0.00105669 0.00106053 0.998944 Analysis of Deviance Source Deviance Df P-Value Model 2.59054 0.2738 Residual 189.72 147 0.0101 Total (corr.) 192.311 149 Percentage of deviance explained by model = 1.34706 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value TC 2.58302 0.1080 TC^2 2.34802 0.1254 Residual Analysis Estimation n MSE MAE MAPE ME MPE Valid ation 150 0.0510647 0.451531 -0.0000401259 The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C2 and independent variable(s) The equation of the fitted model is C2 = exp(-4.97964 + 0.137607*TC - 0.00105669*TC^2)/(1+exp(-4.97964 + 0.137607*TC 0.00105669*TC^2)) Logistic Regression - TrC3 Dependent variable: TrC3 Factors: TC TC^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -4.58419 4.39543 TC 0.137533 0.126435 1.14744 ii TC^2 -0.00108045 0.000920441 0.99892 Analysis of Deviance Source Deviance Model 3.25185 Residual 198.652 Total (corr.) 201.904 Df P-Value 0.1967 147 0.0029 149 Percentage of deviance explained by model = 1.6106 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value TC 3.2518 0.0713 TC^2 3.02336 0.0821 Residual Analysis Estimation n 150 MSE 0.0578368 MAE 0.480077 MAPE ME 0.0000929459 MPE Validation The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C3 and independent variable(s) The equation of the fitted model is TrC3 = exp(-4.58419 + 0.137533*TC - 0.00108045*TC^2)/(1+exp(-4.58419 + 0.137533*TC - 0.00108045*TC^2)) Logistic Regression - TrC4 Dependent variable: TrC4 Factors: TC TC^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -5.35403 5.17951 TC 0.174579 0.148049 1.19074 TC^2 -0.00140762 0.00106433 0.998593 Analysis of Deviance Source Deviance Model 4.97885 Residual 200.8 Total (corr.) 205.779 Df P-Value 0.0830 147 0.0021 149 Percentage of deviance explained by model = 2.41951 Adjusted percentage = 0.0 jj Likelihood Ratio Tests Factor Chi-Squared Df P-Value TC 4.88294 0.0271 TC^2 4.89675 0.0269 Residual Analysis Estimation n 150 MSE 0.0596861 MAE 0.488024 MAPE ME 0.0000778536 MPE Validation The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C4 and independent variable(s) The equation of the fitted model is C4 = exp(-5.35403 + 0.174579*TC - 0.00140762*TC^2)/(1+exp(-5.35403 + 0.174579*TC 0.00140762*TC^2)) Logistic Regression - TrC5 Dependent variable: TrC5 Factors: TC TC^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -6.03236 5.98428 TC 0.199413 0.170349 1.22069 TC^2 -0.00162304 0.00121506 0.998378 Analysis of Deviance Source Deviance Model 5.96189 Residual 199.308 Total (corr.) 205.27 Df P-Value 0.0507 147 0.0027 149 Percentage of deviance explained by model = 2.90442 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value TC 5.6707 0.0172 TC^2 5.94947 0.0147 Residual Analysis Estimation n 150 MSE 0.058639 MAE 0.484234 MAPE Validation kk ME MPE 0.000193967 The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C5 and independent variable(s) The equation of the fitted model is TrC5 = exp(-6.03236 + 0.199413*TC - 0.00162304*TC^2)/(1+exp(-6.03236 + 0.199413*TC - 0.00162304*TC^2)) Phụ Lục 2: Phƣơng trình hồi quy độ ẩm cấp tuổi Logistic Regression - TrC1 Dependent variable: TrC1 Factors: W W^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -36.3724 25.3801 W 0.845602 0.631048 2.32938 W^2 -0.00496697 0.00390384 0.995045 Analysis of Deviance Source Deviance Model 5.09309 Residual 185.861 Total (corr.) 190.954 Df P-Value 0.0784 147 0.0166 149 Percentage of deviance explained by model = 2.66718 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value W 1.94852 0.1627 W^2 1.74076 0.1870 Residual Analysis Estimation n 150 MSE 0.0487801 MAE 0.443822 MAPE ME-0.000304593 MPE Validation The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C1 and independent variable(s) The equation of the fitted model is C1 = exp(-36.3724 + 0.845602*W - 0.00496697*W^2)/(1+exp(-36.3724 + 0.845602*W 0.00496697*W^2)) ll where eta = -36.3724 + 0.845602*W - 0.00496697*W^2 Logistic Regression - TrC2 Dependent variable: TrC2 Factors: W W^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -12.0632 21.7835 W 0.307579 0.548839 1.36013 W^2 -0.00204575 0.00343746 0.997956 Analysis of Deviance Source Deviance Model 0.960388 Residual 191.35 Total (corr.) 192.311 Df P-Value 0.6187 147 0.0082 149 Percentage of deviance explained by model = 0.499394 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value W 0.319091 0.5722 W^2 0.360176 0.5484 Residual Analysis Estimation n 150 MSE 0.0510245 MAE 0.44803 MAPE ME 0.000232354 MPE Validation The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C2 and independent variable(s) The equation of the fitted model is C2 = exp(-12.0632 + 0.307579*W - 0.00204575*W^2)/(1+exp(-12.0632 + 0.307579*W 0.00204575*W^2)) Logistic Regression - TrC3 Dependent variable: TrC3 Factors: W W^2 mm Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -19.8214 21.2903 W 0.5124 0.536275 1.66929 W^2 -0.00334299 0.00335801 0.996663 Analysis of Deviance Source Deviance Model 1.8226 Residual 200.081 Total (corr.) 201.904 Df P-Value 0.4020 147 0.0024 149 Percentage of deviance explained by model = 0.902709 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value W 0.9343 0.3337 W^2 1.0154 0.3136 Residual Analysis Estimation Validation n 150 MSE 0.057669 MAE 0.476257 MAPE ME 0.00030799 MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C3 and independent variable(s) The equation of the fitted model is TrC3 = exp(-19.8214 + 0.5124*W - 0.00334299*W^2)/(1+exp(-19.8214 + 0.5124*W - 0.00334299*W^2)) Logistic Regression - C4 Dependent variable: C4 Factors: W W^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -9.23967 20.6598 W 0.25541 0.520346 1.29099 W^2 -0.0017699 0.00325749 0.998232 Analysis of Deviance Source Deviance Df Model 1.66671 P-Value 0.4346 nn Residual 204.112 Total (corr.) 205.779 147 0.0013 149 Percentage of deviance explained by model = 0.809952 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value W 0.242096 0.6227 W^2 0.296933 0.5858 Residual Analysis Estimation n 150 MSE 0.0606831 MAE 0.488035 MAPE ME 0.000152827 MPE Validation The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C4 and independent variable(s) The equation of the fitted model is C4 = exp(-9.23967 + 0.25541*W - 0.0017699*W^2)/(1+exp(-9.23967 + 0.25541*W 0.0017699*W^2)) Logistic Regression - TrC5 Dependent variable: TrC5 Factors: W W^2 Estimated Regression Model (Maximum Likelihood) Standard Estimated Parameter Estimate Error Odds Ratio CONSTANT -9.09344 20.729 W 0.254721 0.522241 1.2901 W^2 -0.00178802 0.00327035 0.998214 Analysis of Deviance Source Deviance Df P-Value Model 2.03883 0.3608 Residual 203.231 147 0.0015 Total (corr.) 205.27 149 Percentage of deviance explained by model = 0.993244 Adjusted percentage = 0.0 Likelihood Ratio Tests Factor Chi-Squared Df P-Value W 0.239058 0.6249 W^2 0.300722 0.5834 Residual Analysis Estimation Validation oo n 150 MSE 0.0599975 MAE 0.485405 MAPE ME 0.000127905 MPE The StatAdvisor The output shows the results of fitting a logistic regression model to describe the relationship between C5 and independent variable(s) The equation of the fitted model is TrC5 = exp(-9.09344 + 0.254721*W - 0.00178802*W^2)/(1+exp(-9.09344 + 0.254721*W - 0.00178802*W^2)) pp