BỘ GIÁO DỤC VÀ ĐÀO TAO TRƯỜNG ĐẠI HỌC NÔNG LÂM TP HỒ CHÍ MINH **************** NGUYỄN HỮU TRUNG XÁC ĐỊNH NHỮNG MƠ HÌNH PHÙ HỢP NHẤT ĐỂ MƠ TẢ Q TRÌNH SINH TRƯỞNG CỦA RỪNG TRỒNG THƠNG BA LÁ (Pinus keysia Royle ex Gordon) Ở KHU VỰC ĐỨC TRỌNG TỈNH LÂM ĐỒNG LUẬN VĂN THẠC SĨ KHOA HỌC NÔNG NGHIỆP TP Hồ Chí Minh, tháng 12/2011 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NÔNG LÂM THÀNH PHỐ HỒ CHÍ MINH ****************** NGUYỄN HỮU TRUNG XÁC ĐỊNH NHỮNG MƠ HÌNH PHÙ HỢP NHẤT ĐỂ MƠ TẢ Q TRÌNH SINH TRƯỞNG CỦA RỪNG TRỒNG THƠNG BA LÁ (Pinus keysia Royle ex Gordon) Ở KHU VỰC ĐỨC TRỌNG TỈNH LÂM ĐỒNG Chuyên ngành: Lâm Học Mã số : 60.62.60 LUẬN VĂN THẠC SĨ KHOA HỌC NÔNG NGHIỆP Hướng dẫn khoa học: PGS.TS NGUYỄN VĂN THÊM Thành phố Hồ Chí Minh Tháng 12/2011 XÁC ĐỊNH NHỮNG MƠ HÌNH PHÙ HỢP NHẤT ĐỂ MƠ TẢ Q TRÌNH SINH TRƯỞNG CỦA RỪNG TRỒNG THÔNG BA LÁ (Pinus keysia Royle ex Gordon) Ở KHU VỰC ĐỨC TRỌNG, TỈNH LÂM ĐỒNG NGUYỄN HỮU TRUNG Hội đồng chấm luận văn: Chủ tịch: TS NGƠ AN Trường Đại học Nơng Lâm TP.HCM Thư ký: TS BÙI VIỆT HẢI Trường Đại học Nông Lâm TP.HCM Phản biện 1: TS GIANG VĂN THẮNG Hội khoa học lâm nghiệp TP.HCM Phản biện 2: TS PHẠM TRỌNG THỊNH Phân viện Điều tra quy hoạch rừng Nam Bộ Ủy viện: TS PHẠM TRỊNH HÙNG Trường Đại học Nông Lâm TP.HCM LÝ LỊCH CÁ NHÂN Tôi tên Nguyễn Hữu Trung sinh ngày 12 tháng 12 năm 1977 huyện Chương Mỹ thành phố Hà Nội, tốt nghiệp Đại học ngành Lâm nghiệp xã hội, hệ quy Đại học lâm nghiệp thành phố Hà Nội Từ năm 2002 đến cán Ban quản lý rừng phòng hộ Đại Ninh Từ tháng năm 2009 theo học cao học ngành Lâm nghiệp trường đại học Nơng Lâm Thành phố Hồ Chí Minh Địa liên lạc: Nguyễn Hữu Trung Ban quản lý rừng phòng hộ Đại Ninh thơn Phú Thịnh, xã Phú Hội, huyện Đức Trọng, tỉnh Lâm Đồng Điện thoại di động: 01682036875 LỜI CAM ĐOAN Tôi cam đoan cơng trình nghiên cứu tơi Các số liệu, kết nêu luận văn trung thực Chưa công bố cơng trình khác Người cam đoan Nguyễn Hữu Trung LỜI CẢM TẠ Luận văn hoàn thành theo chương trình đào tạo thạc sỹ lâm nghiệp, khóa 2009-2011 Trường Đại học Nơng Lâm TP Hồ Chí Minh Trong trình học tập thực luận văn, tác giả nhận quan tâm, giúp đỡ tạo điều kiện thuận lợi Ban Giám hiệu Thầy – Cô Khoa lâm nghiệp - Đại Học Nơng Lâm TP Hồ Chí Minh, Ban lãnh đạo Ban quản lý rừng phòng hộ Đại Ninh huyện Đức Trọng tỉnh Lâm Đồng Nhân dịp tác giả bày tỏ lòng biết ơn sâu sắc trước quan tâm, giúp đỡ qúy báu Kết luận văn tách rời dẫn nhiệt tình giáo viên hướng dẫn khoa học: PGS TS Nguyễn Văn Thêm – Trưởng Bộ môn Lâm sinh – Khoa Lâm nghiệp - Đại Học Nông Lâm TP Hồ Chí Minh Tác giả xin bày tỏ lòng biết ơn sâu sắc kính trọng dẫn chân tình thầy hướng dẫn Trong trình làm luận văn, tác giả nhận giúp đỡ cổ vũ chân tình Bố - Mẹ, vợ con, anh chị em gia đình, bạn đồng nghiệp quan khóa học Tác giả xin chân thành cảm ơn ghi nhớ giúp đỡ cổ vũ vô tư TP Hồ Chí Minh, ngày 30 tháng 12 năm 2011 Nguyễn Hữu Trung TÓM TẮT Đề tài “Xác định mơ hình phù hợp để mơ tả q trình sinh trưởng rừng trồng Thơng ba (Pinus keysia Royle ex Gordon)” thực Đức Trọng tỉnh Lâm Đồng Thời gian nghiên cứu từ tháng đến tháng năm 2011 Mục tiêu đề tài xây dựng mơ hình thống kê phù hợp để mơ tả phân tích quy luật sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khác khu vực Đức Trọng tỉnh Lâm Đồng Kết nghiên cứu rằng: (1) Nếu sử dụng hàm (Korf, Schumacher, Gompertz, Korsun-Strand Drakin-Vuevski) để làm phù hợp với số liệu thực nghiệm trình sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khu vực Đức Trọng tỉnh Lâm Đồng, hàm Korf hàm phù hợp (2) Nếu sử dụng mô hình khác nhau, đặc trưng sinh trưởng đường kính, chiều cao thể tích thân Thơng ba cấp đất báo cáo khác (3) Cùng mơ hình sinh trưởng, sử dụng phương pháp tiêu chuẩn khác để kiểm định mơ hình phù hợp, đặc trưng sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất báo cáo khác (4) Sinh trưởng đường kính thân Thông ba 30 năm đầu thay đổi tùy theo tuổi cấp đất Thời điểm xuất ZDmax cấp đất cấp tuổi Thời điểm xuất Dmax cấp đất I-IV tuổi 16, cấp đất V cấp tuổi 18; trung bình tuổi 16 Tốc độ sinh trưởng đường kính giai đoạn 30 tuổi cấp đất I lớn cấp đất II, III, IV V tương ứng 1,07; 1,14; 1,22 1,32 lần (5) Sinh trưởng chiều cao thân Thông ba 30 năm đầu thay đổi tùy theo tuổi cấp đất Thời điểm xuất ZHmax cấp đất I, II III cấp tuổi 4, cấp đất IV V cấp tuổi 6; trung bình cấp tuổi Thời điểm xuất Hmax cấp đất I, II III cấp tuổi 8, cấp đất IV V cấp tuổi 10; trung bình cấp tuổi Tốc độ sinh trưởng chiều cao giai đoạn 30 tuổi cấp đất I lớn cấp đất II, III, IV V tương ứng 1,05; 1,10; 1,17 1,23 lần (6) Thể tích thân Thông ba 30 năm đầu thay đổi tùy theo tuổi cấp đất Tốc độ sinh trưởng thể tích thân giai đoạn 30 tuổi cấp đất I lớn cấp đất II, III, IV V tương ứng 1,14; 1,28; 1,44 1,61 lần Từ kết nghiên cứu, đề tài đề xuất mơ hình phù hợp để dự đốn sinh trưởng Thông ba cấp đất khu vực Đức Trọng tỉnh Lâm Đồng SUMMARY The thesis “Identification the most appropriate models to illustrate the growing process of Pinus keysia Royle ex Gordon artificial forest at Duc Trong District, Lam Dong province” The research period was from April to September 2011 The research objectives was to build the most appropriate models to exemplify and analyze the growing rules of Pinus keysia diameter, height and treetrunk volume on different site indexes at Duc Trong district, Lam Dong province To describe the Pinus keysia diameter, height, and tree-trunk volume growing process, the author used five statistic models – they are Korf, Schumacher, Gompertz, Korsun-Strand and Drakin-Vuevski model The parameters of five models are identified by two method – they are minimum square method and nonlinear regression method The suitable statistic model is tested by five criteria – they are R2max , Minimum standard error of estimation, minimum mean absolute error, minimum mean absolute percent error and minimum risidual sum of square The research results showed that: (1) If we use five models (Korf, Schumacher, Gompertz, Korsun-Strand and Drakin-Vuevski) to make suitability of experimental data of the Pinus diameter, height, and tree-trunk volume growing process on five Site Indexes in Duc Trong District, Lam Dong Province, the Korf Model is the most suitable (2) If we use different models, the Pinus keysia diameter, height, and tree-trunk volume growing characteristics on five Site Indexes are also different (3) In the same growing model, if we use different methods and criteria to test the appropriate model, the Pinus keysia diameter, height, and tree-trunk volume growing characteristics on five Site Indexes are also different (4) Diameter grows of Pinus keysia in the first 30 years change depending on age and site index The time of ZDmax appearance at all five site indexes is at years old-age The time of Dmax appearance of site index I-IV is at 16 years old-age, while site index V is at 18 years old-age, on average at 16 years old-age The diameter growing speed in 30 years old-age on site index I greater than site index II, III, IV, V is 1,07; 1,14; 1,22 and 1,32 time respectively (5) Height grows of Pinus keysia in the first 30 years change depending on age and site index The time of ZHmax appearance on site index I, II, and III is at years old-age, while index IV and V is at year old-age; on average is at years oldage The time of ZHmax appearance on site index I, II, and III is at years oldage, while index IV and V is at 10 year old-age; on average is at years oldage The height growing speed in 30 years old-age on site index I greater than site index II, III, IV, V is 1,05; 1,10; 1,17 and 1,23 time respectively (6) Tree-trunk volume grows of Pinus keysia in the first 30 years change depending on age and site index The tree-trunk volume growing speed in the time of 30 years old on site index I greater than site index II, III, IV, V is 1,14; 1,28; 1,44 1,61 time respectively 10 Phụ lục 21 Số liệu chiều cao thực nghiệm Thông ba 26 tuổi cấp đất V khu vực Đức Trọng tỉnh Lâm Đồng A(năm) (1) 10 12 14 16 18 20 22 24 26 Chiều cao thân mẫu (m): Cây Cây Cây Bình quân (2) (3) (4) (5) 1,9 1,6 2,0 1,8 3,7 3,1 3,7 3,5 5,1 4,3 5,2 4,9 6,5 5,6 6,5 6,2 7,9 6,6 7,9 7,5 9,0 7,4 9,1 8,5 9,9 8,3 10,0 9,4 10,8 9,3 10,8 10,3 11,7 9,7 11,8 11,1 12,3 10,3 12,6 11,7 12,9 10,7 13,6 12,4 13,8 11,5 14,0 13,1 Tăng trưởng: ZH (m/năm) H(m/năm) (6) (7) 0,46 0,83 0,71 0,66 0,63 0,52 0,44 0,44 0,38 0,34 0,34 0,35 HV (m) 0,46 0,58 0,61 0,62 0,62 0,61 0,59 0,57 0,55 0,53 0,52 0,50 Ph(%) (8) 25,0 23,7 14,5 10,7 8,5 6,1 4,7 4,3 3,5 2,9 2,7 2,7 ZH H 14 0.9 12 0.8 0.7 10 0.6 0.5 0.4 0.3 0.2 0.1 0.0 H 10 12 14 16 18 20 22 24 26 ZH A(năm) AH Đồ thị mơ tả q trình sinh trưởng chiều cao thân Thông ba thực nghiệm cấp đất V khu vực Đức Trọng tỉnh Lâm Đồng 170 Phụ lục 22 Phân tích hồi quy tương quan H-A rừng Thông ba cấp đất V theo mơ hình Korf, Schumacher, Gompertz, Korsun-Strand Drakin-Vuevski 22.1 Mơ hình Korf Source DF Sum of Squares Mean Square Regression 2957.04057 985.68019 Residual 33 22.04943 66816 Uncorrected Total 36 2979.09000 (Corrected Total) 35 460.72306 R squared = - Residual SS / Corrected SS = 95214 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y B C 35.968052443 17.431880339 502625230 71.433479657 6.544703566 1.054284645 4.399745328 8.689661805 572212948 187443095 190857104 953568792 22.2 Mơ hình Schumacher Source DF Sum of Squares Mean Square Regression 2957.04057 985.68019 Residual 33 22.04943 66816 Uncorrected Total 36 2979.09000 (Corrected Total) 35 460.72306 R squared = - Residual SS / Corrected SS = 95214 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 35.968053507 17.431944837 502495072 71.433611943 B 6.544703508 1.054281689 4.399751284 8.689655731 C 572212937 187443070 190857144 953568730 22.3 Mơ hình Gompertz Source DF Sum of Squares Mean Square Regression 2955.72802 985.24267 Residual 33 23.36198 70794 Uncorrected Total 36 2979.09000 (Corrected Total) 35 460.72306 R squared = - Residual SS / Corrected SS = 94929 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y B C 14.380501220 893471274 12.562720244 16.198282196 3.066738788 347311756 2.360127708 3.773349867 125879124 017400867 090476794 161281454 22.4 Mô hình Korsun-Strand 171 Polynomial Regression Analysis Dependent variable: H Standard T Parameter Estimate Error Statistic P-Value CONSTANT 4.87665 2.38906 2.04125 0.0493 A 0.710738 0.358499 1.98254 0.0558 A^2 0.0429288 0.0117096 3.66613 0.0009 Analysis of Variance Source Sum of Squares Df Mean Square F-Ratio P-Value Model 6972.46 3486.23 396.88 0.0000 Residual 289.874 33 8.78406 Total (Corr.) 7262.33 35 R-squared = 96.0085 percent R-squared (adjusted for d.f.) = 95.7666 percent Standard Error of Est = 2.96379 Mean absolute error = 2.30695 Durbin-Watson statistic = 0.318554 The StatAdvisor The output shows the results of fitting a second order polynomial model to describe the relationship between H and A The equation of the fitted model is H = 4.87665 + 0.710738*A + 0.0429288*A^2 22.5 Mơ hình Drakin-Vuevski Source DF Sum of Squares Mean Square Regression 2956.77187 985.59062 Residual 33 22.31813 67631 Uncorrected Total 36 2979.09000 (Corrected Total) 35 460.72306 R squared = - Residual SS / Corrected SS = 95156 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 16.240605209 2.048227705 12.073454610 20.407755808 K 077038482 023693631 028833427 125243536 M 1.564721525 317707113 918341544 2.211101507 172 Phụ lục 23 Số liệu chiều cao thực nghiệm rừng Thông ba 26 tuổi khu vực Đức Trọng tỉnh Lâm Đồng A(năm) 10 12 14 16 18 20 22 24 26 H(m) 3,4 6,0 8,2 10,2 12,0 13,5 14,9 16,0 17,2 18,1 19,0 20,0 ZH (m/năm) H(m/năm) 0,85 1,30 1,11 0,99 0,91 0,76 0,66 0,58 0,59 0,47 0,44 0,47 0,85 1,00 1,03 1,02 1,00 0,97 0,93 0,89 0,86 0,82 0,79 0,77 H(m) PH% 25,0 21,7 13,6 9,7 7,6 5,6 4,5 3,6 3,4 2,6 2,3 2,4 ZH H 25 1.40 1.20 20 1.00 15 0.80 H(m) ZH 0.60 10 DH 0.40 0.20 0.00 10 12 14 16 18 20 22 24 26 A(năm) Đồ thị mơ tả q trình sinh trưởng chiều cao thực nghiệm rừng Thông ba khu vực Đức Trọng tỉnh Lâm Đồng 173 Phụ lục 24 Phân tích hồi quy tương quan H-A rừng Thơng ba theo mơ hình Korf, Schumacher, Gompertz, Korsun-Strand Drakin-Vuevski 24.1 Mơ hình Korf Source DF Sum of Squares Mean Square Regression 13778.90218 4592.96739 Residual 60 840.86782 14.01446 Uncorrected Total 63 14619.77000 (Corrected Total) 62 2729.44857 R squared = - Residual SS / Corrected SS = 69193 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y B C 122.69817717 293.73212835 -464.8535594 710.24991377 5.746743024 757831012 4.230855302 7.262630746 360090958 387540288 -.415105036 1.135286952 24.2 Mơ hình Schumacher Source DF Sum of Squares Mean Square Regression 13778.90218 4592.96739 Residual 60 840.86782 14.01446 Uncorrected Total 63 14619.77000 (Corrected Total) 62 2729.44857 R squared = - Residual SS / Corrected SS = 69193 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 122.70015379 293.84623852 -465.0798371 710.48014472 B 5.746747347 758104666 4.230312235 7.263182460 C 360088445 387539144 -.415105261 1.135282150 24.3 Mơ hình Gompertz Source DF Sum of Squares Mean Square Regression 13763.28486 4587.76162 Residual 60 856.48514 14.27475 Uncorrected Total 63 14619.77000 (Corrected Total) 62 2729.44857 R squared = - Residual SS / Corrected SS = 68621 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y B C 24.872045954 4.211872682 16.447046201 33.297045706 2.586989029 495168296 1.596504965 3.577473093 100711066 033366800 033967528 167454603 174 24.4 Mơ hình Korsun-Strand Polynomial Regression Analysis Dependent variable: H Standard T Parameter Estimate Error Statistic P-Value CONSTANT 0.5705 3.4786 0.164003 0.8703 A 0.913949 0.49924 1.83068 0.0721 A^2 0.0145699 0.0155905 0.934534 0.3538 Analysis of Variance Source Sum of Squares Df Mean Square F-Ratio P-Value Model 6310.81 3155.41 90.93 0.0000 Residual 2081.99 60 34.6998 Total (Corr.) 8392.8 62 R-squared = 75.1932 percent R-squared (adjusted for d.f.) = 74.3663 percent Standard Error of Est = 5.89065 Mean absolute error = 4.43957 Durbin-Watson statistic = 0.0845843 The StatAdvisor The output shows the results of fitting a second order polynomial model to describe the relationship between H and A The equation of the fitted model is H = 0.5705 + 0.913949*A + 0.0145699*A^2 24.5 Mơ hình Drakin-Vuevski Source DF Sum of Squares Mean Square Regression 13775.23711 4591.74570 Residual 60 844.53289 14.07555 Uncorrected Total 63 14619.77000 (Corrected Total) 62 2729.44857 R squared = - Residual SS / Corrected SS = 69058 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 33.133028689 19.814509579 -6.501891666 72.767949045 K 041238899 050621644 -.060019465 142497263 M 1.115261045 482816265 149484722 2.081037368 175 Phụ lục 25 Xây dựng mối quan hệ hình số với tuổi Thơng ba Regression Analysis - Multiplicative model: Y = a*X^b Dependent variable: F Independent variable: A Standard T Parameter Estimate Error Statistic P-Value Intercept -0.0999003 0.102654 -0.973176 0.3450 Slope -0.16691 0.0370222 -4.50839 0.0004 NOTE: intercept = ln(a) Analysis of Variance Source Sum of Squares Df Mean Square F-Ratio P-Value Model 0.183492 0.183492 20.33 0.0004 Residual 0.144443 16 0.00902766 Total (Corr.) 0.327935 17 Correlation Coefficient = -0.748023 R-squared = 55.9539 percent Standard Error of Est = 0.095014 The StatAdvisor The output shows the results of fitting a multiplicative model to describe the relationship between F and A The equation of the fitted model is F = 0.904928*A^-0.16691 or ln(F) = -0.0999003 - 0.16691*ln(A) Since the P-value in the ANOVA table is less than 0.01, there is a Phụ lục 26 Thể tích thân Thơng ba cấp đất A Thể tích thân Thơng ba cấp đất: I II III IV V V(Bình quân) 0,00206 0,00128 0,00070 0,00038 0,00018 0,00117 0,01557 0,00979 0,00573 0,00318 0,00150 0,00649 0,05009 0,03183 0,01927 0,01084 0,00523 0,01891 10 0,10908 0,06986 0,04313 0,02450 0,01203 0,04031 12 0,19156 0,12348 0,07717 0,04417 0,02201 0,07162 14 0,29458 0,19088 0,12028 0,06927 0,03496 0,11311 16 0,41457 0,26982 0,17100 0,09899 0,05052 0,16470 18 0,54812 0,35809 0,22787 0,13250 0,06828 0,22603 20 0,69217 0,45370 0,28958 0,16904 0,08786 0,29663 22 0,84415 0,55494 0,35499 0,20792 0,10891 0,37597 24 1,00190 0,66036 0,42313 0,24859 0,13113 0,46347 26 1,16366 0,76878 0,49322 0,29056 0,15425 0,55857 28 1,32799 0,87922 0,56459 0,33344 0,17804 0,66073 30 1,49375 0,99089 0,63674 0,37690 0,20233 0,76942 176 Phụ lục 27 Xây dựng mơ hình V-A Thông ba cấp đất theo hàm Korf 27.1 Mơ hình V-A Thơng ba cấp đất I Source DF Sum of Squares Mean Square Regression 8.15487 2.71829 Residual 12 3.859108E-07 3.215923E-08 Uncorrected Total 15 8.15487 (Corrected Total) 14 3.72750 R squared = - Residual SS / Corrected SS = 1.00000 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 23.020976941 126628283 22.745077614 23.296876269 B 21.875041307 030915824 21.807681512 21.942401101 C 611322803 000995927 609152865 613492742 V(I) = 23.020976*exp(-21.875041*A^-0.611322) 27.2 Mơ hình V-A Thông ba cấp đất II Source DF Sum of Squares Mean Square Regression 3.55450 1.18483 Residual 12 1.561080E-07 1.300900E-08 Uncorrected Total 15 3.55450 (Corrected Total) 14 1.63709 R squared = - Residual SS / Corrected SS = 1.00000 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 16.369897527 089814998 16.174207456 16.565587598 B 22.035728201 029936320 21.970502563 22.100953838 C 606097949 000964328 603996858 608199040 V(II) = 16.369897*exp(-22.035728*A^-0.606098) 27.3 Mơ hình V-A Thơng ba cấp đất III Source DF Sum of Squares Mean Square Regression 1.46023 48674 Residual 12 6.691340E-08 5.576117E-09 Uncorrected Total 15 1.46023 (Corrected Total) 14 67709 R squared = - Residual SS / Corrected SS = 1.00000 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 9.920963692 053734689 9.803885863 10.038041521 B 22.793732296 032959338 22.721920068 22.865544524 C 622245305 000994529 620078411 624412199 V(III) = 9.920963*exp(-22.793732*A^-0.622245) 27.4 Mơ hình V-A Thơng ba cấp đất IV Source DF Sum of Squares Mean Square Regression 50610 16870 Residual 12 1.928800E-08 1.607334E-09 177 Uncorrected Total 15 50610 (Corrected Total) 14 23668 R squared = - Residual SS / Corrected SS = 1.00000 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 6.353440126 032638341 6.282327290 6.424552962 B 22.982645297 030255324 22.916724609 23.048565985 C 616359884 000912374 614371992 618347776 V(IV) = 6.353440*exp(-22.982645*A^-0.616359) 27.5 Mơ hình V-A Thơng ba cấp đất V Source DF Sum of Squares Mean Square Regression 14230 04743 Residual 12 4.697442E-09 3.914535E-10 Uncorrected Total 15 14230 (Corrected Total) 14 06778 R squared = - Residual SS / Corrected SS = 1.00000 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y B C 4.227358274 022425464 4.178497386 4.276219163 23.244367618 028027262 23.183301461 23.305433775 598152075 000859423 596279553 600024596 V(V) = 4.227358*exp(-23.244367*A^-0.598152) 27.6 Mơ hình V-A bình qn Thơng ba Source DF Sum of Squares Mean Square Regression 1.88288 62763 Residual 12 8.369572E-09 6.974644E-10 Uncorrected Total 15 1.88288 (Corrected Total) 14 93678 R squared = - Residual SS / Corrected SS = 1.00000 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std Error Lower Upper Y 294.41612296 1.662825073 290.79313836 298.03910756 B 20.694475397 002561923 20.688893447 20.700057347 C 366626025 000311559 365947197 367304853 V = 294.416122*exp(-20.694475*A^-0.366626) 178 Phụ lục 28 Kiểm định khác biệt đường kính thân Thơng ba cấp đất 28.1 Mơ hình kiểm định D-A cấp đất Multiple Regression Analysis -Standard T Parameter Estimate Error Statistic P-Value -CONSTANT 2.17361 0.00644068 337.482 0.0000 LnA -0.626604 0.00234228 -267.519 0.0000 Cap dat=2 -0.0303383 0.00910849 -3.33077 0.0015 Cap dat=3 -0.0375067 0.00910849 -4.11777 0.0001 Cap dat=4 -0.0552889 0.00910849 -6.07004 0.0000 Cap dat=5 -0.062366 0.00910849 -6.84702 0.0000 LnA*Cap dat=2 0.0408829 0.00331249 12.3421 0.0000 LnA*Cap dat=3 0.0750383 0.00331249 22.6532 0.0000 LnA*Cap dat=4 0.114845 0.00331249 34.6702 0.0000 LnA*Cap dat=5 0.153527 0.00331249 46.3479 0.0000 -Analysis of Variance Source Sum of Squares Df Mean Square F-Ratio P-Value Model 8.47592 0.941769 35363.12 0.0000 Residual 0.00159788 60 0.0000266314 Total (Corr.) 8.47752 69 R-Squared = 99.9812 percent R-Squared (adjusted for d.f.) = 99.9783 percent Standard Error of Est = 0.00516056 Mean absolute error = 0.00390762 Durbin-Watson statistic = 0.58247 Residual Analysis Estimation Validation n 70 MSE 0.0000266314 MAE 0.00390762 MAPE 1.25838 ME 6.3937E-17 MPE -0.538006 28.2 Kiểm định dộ dốc Further ANOVA for Variables in the Order Fitted Source Sum of Squares Df Mean Square F-Ratio P-Value LnA 7.33515 7.33515 275432.41 0.0000 Intercepts 1.07026 0.267566 10047.01 0.0000 Slopes 0.0705105 0.0176276 661.91 0.0000 Model 8.47592 179 Phụ lục 29 Kiểm định khác biệt chiều cao thân Thông ba cấp đất 29.1 Mơ hình kiểm định H-A cấp đất Multiple Regression Analysis -Standard T Parameter Estimate Error Statistic P-Value -CONSTANT 1.66247 0.0170282 97.6306 0.0000 LnA -0.427957 0.00619265 -69.1072 0.0000 Cap dat=2 0.0298142 0.0240815 1.23805 0.2205 Cap dat=3 0.0717975 0.0240815 2.98143 0.0041 Cap dat=4 0.112215 0.0240815 4.65978 0.0000 Cap dat=5 0.170856 0.0240815 7.0949 0.0000 LnA*Cap dat=2 0.0212534 0.00875773 2.42682 0.0183 LnA*Cap dat=3 0.0394415 0.00875773 4.50362 0.0000 LnA*Cap dat=4 0.0616752 0.00875773 7.04237 0.0000 LnA*Cap dat=5 0.0810569 0.00875773 9.25547 0.0000 -Analysis of Variance Source Sum of Squares Df Mean Square F-Ratio P-Value Model 4.9755 0.552833 2969.78 0.0000 Residual 0.0111692 60 0.000186153 Total (Corr.) 4.98667 69 R-Squared = 99.776 percent R-Squared (adjusted for d.f.) = 99.7424 percent Standard Error of Est = 0.0136438 Mean absolute error = 0.0108404 Durbin-Watson statistic = 0.406812 Residual Analysis Estimation Validation n 70 MSE 0.000186153 MAE 0.0108404 MAPE 1.72398 ME -4.44089E-17 MPE 0.179091 29.2 Kiểm định dộ dốc Further ANOVA for Variables in the Order Fitted Source Sum of Squares Df Mean Square F-Ratio P-Value LnA 3.64013 3.64013 19554.53 0.0000 Intercepts 1.31545 0.328861 1766.62 0.0000 Slopes 0.0199258 0.00498146 26.76 0.0000 Model 4.9755 180 Phụ lục 30 Kiểm định khác biệt thể tích thân Thơng ba cấp đất 30.1 Mơ hình kiểm định V-A cấp đất Multiple Regression Analysis -Standard T Parameter Estimate Error Statistic P-Value -CONSTANT 3.44479 0.0353877 97.3443 0.0000 LnA -0.931633 0.0128695 -72.391 0.0000 Cap dat=2 -0.138059 0.0500458 -2.75865 0.0077 Cap dat=3 -0.216311 0.0500458 -4.32227 0.0001 Cap dat=4 -0.28285 0.0500458 -5.65182 0.0000 Cap dat=5 -0.314272 0.0500458 -6.27969 0.0000 LnA*Cap dat=2 0.115078 0.0182002 6.32291 0.0000 LnA*Cap dat=3 0.201281 0.0182002 11.0593 0.0000 LnA*Cap dat=4 0.282849 0.0182002 15.541 0.0000 LnA*Cap dat=5 0.353101 0.0182002 19.401 0.0000 -Analysis of Variance Source Sum of Squares Df Mean Square F-Ratio P-Value Model 17.1774 1.9086 2374.00 0.0000 Residual 0.0482377 60 0.000803962 Total (Corr.) 17.2257 69 R-Squared = 99.72 percent R-Squared (adjusted for d.f.) = 99.678 percent Standard Error of Est = 0.0283542 Mean absolute error = 0.0166332 Durbin-Watson statistic = 0.351318 Residual Analysis Estimation Validation n 70 MSE 0.000803962 MAE 0.0166332 MAPE 2.42745 ME -1.34733E-15 MPE -0.886613 30.2 Kiểm định dộ dốc Further ANOVA for Variables in the Order Fitted Source Sum of Squares Df Mean Square F-Ratio P-Value LnA 13.3328 13.3328 16583.89 0.0000 Intercepts 3.47053 0.867634 1079.20 0.0000 Slopes 0.374075 0.0935188 116.32 0.0000 Model 17.1774 181 Phụ lục 31: Tọa độ ô tiêu chuẩn Stt Tọa độ ô tiêu chuẩn X Y 0573726 1305698 0572868 1304874 0561120 1300247 0573310 1306328 0561982 1301382 182 Cấp đất MỘT SỐ HÌNH ẢNH ĐO ĐẾM VÒNG NĂM 183 184 ... cam đoan cơng trình nghiên cứu tơi Các số liệu, kết nêu luận văn trung thực Chưa công bố cơng trình khác Người cam đoan Nguyễn Hữu Trung LỜI CẢM TẠ Luận văn hồn thành theo chương trình đào tạo thạc... I, II III cấp tuổi 4, cấp đất IV V cấp tuổi 6; trung bình cấp tuổi Thời điểm xuất Hmax cấp đất I, II III cấp tuổi 8, cấp đất IV V cấp tuổi 10; trung bình cấp tuổi Tốc độ sinh trưởng chiều cao... Tổng bình phương sai lệch (Sum of Squared Residuals) Sai số tuyệt đối trung bình (Mean Absolute Error) Sai số tuyệt đối trung bình theo phần trăm (Mean Absolute Percent Error) Trị lớn nhỏ 15