4.3.1. M i quan h gi a nh ng b ph n sinh kh i khô c a cây tràm
K t qu phân tích ma tr n t ng quan gi a nh ng b ph n sinh kh i khô (Ph l c 1b) trên m t t c a cây tràm c ghi l i b ng 4.15. T ó cho th y:
+ Gi a nh ng b ph n sinh kh i khô (TSKk, SKTk, SKCk và SKLk) và Dcv có m i quan h r t ch t ch v i nhau (r = 0,860 – 0,922; P < 0,001). T ng t , gi a nh ng b ph n sinh kh i khô (TSKk, SKTk, SKCk và SKLk) và H c ng t n t i m i quan h r t ch t ch v i nhau (r = 0,854 – 0,883; P < 0,001).
+ Gi a nh ng b ph n sinh kh i khô (TSKk, SKTk, SKCk và SKLk) c a cây tràm cajuputi c ng có m i quan h r t ch t ch v i nhau (r = 0,868 – 0,998; P <
0,001).
T nh ng m i quan h ch t ch này cho phép xác nh nhanh nh ng thành ph n sinh kh i khô (TSKk, SKTk, SKCk và SKLk) khó o c tr c ti p ngoài tr i thông qua m i quan h gi a chúng v i Dcv(cm) và H(m).
0 10 20 30 40 50 60 70 80
0 10 20 30 40 50 60 70 80
TSK(t)thc nghim
TSK(t) oán
0 10 20 30 40 50 60
0 10 20 30 40 50 60
SKT(t)thc nghim
SKT(t) oán
0 2 4 6 8 10 12
0 2 4 6 8 10 12
SKC(t)thc nghim
SKC(t) oán
0 1 2 3 4 5 6
0 1 2 3 4 5 6
TSK(t)thc nghim
SKL(t) oán
Hình 4.13. th mô t sinh kh i t i th c nghi m và d oán
B ng 4.15. T ng quan gi a nh ng b ph n sinh kh i khô c a cây tràm Ch tiêu Th ng kê Dcv(cm) H(m) SKT(k) SKC(k) SKL(k)
H(m) r 0,964 1
P 0,0000 ,
N 34 34
SKT(k) r 0,913 0,876 1
P 0,0000 0,0000 ,
N 34 34 34
SKC(k) r 0,918 0,859 0,944 1
P 0,0000 0,0000 0,0000 ,
N 34 34 34 34
SKL(k) r 0,860 0,854 0,903 0,868 1
P 0,0000 0,0000 0,0000 0,0000 ,
N 34 34 34 34 34
TSK(t) r 0,922 0,883 0,998 0,959 0,915
P 0,0000 0,0000 0,0000 0,0000 0,0000
N 34 34 34 34 34
4.3.2. Xây d ng mô hình sinh kh i khô c a cây tràm theo c p Dcv 4.3.2.1. Xây d ng mô hình t ng sinh kh i khô theo c p Dcv
Nh ng tính toán cho th y (Ph l c 10), n u s d ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative) mô t m i quan h gi a TSKt c a cây tràm v i c p Dcv, thì k t qu c ghi l i b ng 4.16-4.18 và hình 4.16.
B ng 4.16.Mô hình mô t t ng sinh kh i khô c a cây tràm theo c p Dcv
TT Hàm Mô hình
(1) (2) (3)
1 Gompertz TSK(k) = 61,24018*exp(-11,57389*exp(-0,30892*Dcv)) 2 Drakin-Vuevski TSK(k) = 63,56878*(1-exp(-0,27840*Dcv))^8,87102 3 Schumacher TSK(k) = 175,79899*exp(-22,75976/Dcv^1,15959) 4 Korf TSK(k) = 175,79899*exp(-22,75976*Dcv^-1,15959) 5 S-Curve TSK(k) = exp(5,50725 – 19,09841/Dcv)
6 Multiplicative TSK(k) = 0,11388*Dcv^2,51746
B ng 4.17.Nh ng th ng kê sai l ch c a 6 hàm dùng mô t t ng sinh kh i khô theo c p Dcv c a cây tràm
Th t Hàm r2 ±Se MAE MAPE SSR
(1) (2) (3) (4) (5) (6) (7)
1 Gompertz 96,44 2,40 1,61 35,95 178,99
2 Drakin-Vuevski 96,31 2,44 1,67 39,59 185,37
3 Schumacher 96,07 2,52 1,71 40,56 197,35
4 Korf 96,07 2,52 1,71 40,56 197,35
5 S-Curve 96,06 2,48 1,71 39,82 197,84
6 Multiplicative 95,44 2,67 1,69 23,98 228,94
T s li u c a b ng 4.4 và hình 4.1 cho th y, hàm Gompertz có h s R2max (96,44%), còn hàm Multiplicative có R2min (95,44%). Hàm Gompertz nh n Semin (2,40); còn Semax là hàm Multiplicative (2,67). Giá tr MAEmin là hàm Gompertz (1,61), còn MAEmax là ba hàm Schumacher, Korf và S-Curve (1,71). Giá tr MAPEmin là hàm Multiplicative (23,98%), còn MAPEmax là hàm Schumacher và Korf (40,56%). Giá tr SSRmin thu c v hàm Gompertz (178,99), còn SSRmax là hàm Multiplicative (228,94).
Phân tích sai l ch c a 6 hàm trên ây c ng cho th y, khi d oán TSKk theo c p Dcv c a cây tràm, thì c 6 hàm này u có sai khác l n hai c p Dcv u (1 và 2cm) và hai c p Dcv cu i (11 và 12cm)(B ng 4.5; Hình 4.17).
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40 TSKk (kg/cây)
Dcv(cm) Drakin-Vuevski
Hình 4.16. th t ng sinh kh i khô c a cây tràm nh ng c p Dcv khác nhau c mô t b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative).
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40 TSKk (kg/cây)
Dcv(cm) Schumacher
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40 TSKk (kg/cây)
Dcv(cm) Gompertz
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40 TSKk (kg/cây)
Dcv(cm) Korf
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40
Multiplicative TSKk (kg/cây)
Dcv(cm)
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40
S-Curve TSKk (kg/cây)
Dcv(cm)
B ng 4.18. T ng sinh kh i khô c a cây tràm cajuputi theo c p Dcv c mô t b ng 6 hàm khác nhau
D(cm) Gompertz Schumacher Korf Drakin S-Curve Multiplicative
(1) (2) (3) (4) (5) (6) (7)
1 0,0 0,0 0,0 0,0 0,1 0,0
2 0,1 0,0 0,0 0,0 0,7 0,0
3 0,6 0,3 0,3 0,4 1,8 0,4
4 2,1 1,8 1,8 1,9 3,7 2,1
5 5,2 5,2 5,2 5,0 6,5 5,4
6 10,0 10,2 10,2 10,0 10,4 10,2
7 16,2 16,2 16,2 16,3 15,3 16,1
8 23,0 22,8 22,8 23,1 21,4 22,6
9 29,9 29,6 29,6 29,9 28,8 29,5
10 36,2 36,4 36,4 36,1 37,5 36,5
11 41,6 42,9 42,9 41,6 47,7 43,4
12 46,1 49,1 49,1 46,2 59,3 50,2
T nh ng phân tích th ng kê sai l ch c a 6 hàm cho th y, n u s d ng tiêu chu n R2max và SSRmin ch n hàm mô t TSKk theo c p Dcv c a cây tràm, thì hàm Gompertz là hàm phù h p nh t. Theo ó, mô hình mô t TSKk theo c p Dcv c a cây tràmcajuputi b ng hàm Gompertz có d ng (Hình 4.18):
0.0 10.0 20.0 30.0 40.0 50.0 60.0
1 2 3 4 5 6 7 8 9 10 11 12
Gompertz Schumacher Korf Drakin S-Curve Multiplicative
TSKk (kg/cây)
Dcv(cm)
Hình 4.17. th t ng sinh kh i khô c a cây tràm nh ng c p Dcv
khác nhau c mô t b ng 6 hàm (Schumacher, Korf, Drakin- Vuevsi, Gompertz, S-Curve và Multiplicative).
TSK(k) = 61,24018*exp(-11,57389*exp(-0,30892*Dcv)) (4.9) R2 = 96,44%; Se = ±2,40; MAE = 1,61; MAPE = 35,95%.
4.3.2.2. Xây d ng mô hình sinh kh i thân khô theo c p Dcv
K t qu tính toán cho th y (Ph l c 11), n u s d ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative) mô t m i quan h gi a SKTk c a cây tràm v i c p Dcv, thì k t qu c ghi l i b ng 4.19-4.20 và hình 4.19.
B ng 4.19.Mô hình mô t sinh kh i thân khô c a cây tràm theo c p Dcv
TT Hàm Mô hình
(1) (2) (3)
1 Gompertz SKT(k) = 44,32462*exp(-15,41488*exp(-0,35859*Dcv)) 2 Drakin-Vuevski SKT(k) = 44,83933*(1-exp(-0,33996*Dcv))^12,96299 3 Schumacher SKT(k) = 85,52350*exp(-41,89655/Dcv^1,58881) 4 Korf SKT(k) = 85,52350*exp(-41,89655*Dcv^-1,58881) 5 S-Curve SKT(k) = exp(5,39043 – 20,07379/Dcv)
6 Multiplicative SKT(k) = 0,072465*Dcv^2,61975
1 2 3 4 5 6 7 8 9 10
0 10 20 30 40
TSKk (kg/cây)
p Dcv(cm)
Hình 4.18. th mô t t ng sinh kh i khô c a cây tràm theo c p Dcv b ng hàm Gompertz.
B ng 4.20.Nh ng th ng kê sai l ch c a 6 hàm dùng mô t sinh kh i thân khô theo c p Dcv c a cây tràm
Th t Hàm r2 ±Se MAE MAPE SSR
(1) (2) (3) (4) (5) (6) (7)
1 Gompertz 96,18 2,01 1,31 39,93 125,87
2 Drakin-Vuevski 96,07 2,04 1,33 41,20 128,87
3 Schumacher 95,82 2,11 1,37 42,25 137,56
4 Korf 95,82 2,10 1,37 42,25 137,56
5 S-Curve 95,70 2,10 1,39 41,06 141,42
6 Multiplicative 94,82 2,12 1,49 25,51 170,44
Phân tích s li u c a b ng 4.20 và hình 4.19 cho th y, hàm Gompertz có h s R2max (96,18%), còn hàm Multiplicative có R2min (94,82%). Giá tr Semin thu c v hàm Gompertz (2,01); còn Semax là hàm Multiplicative (2,12). Hàm Gompertz nh n giá tr MAEmin (1,31), còn hàm Multiplicative nh n MAEmax (1,49). Giá tr MAPEmin là hàm Gompertz (125,87%), còn MAPEmax là hàm S-Curve (38,97%). Hàm Gompertz có SSRmin (339,43), còn SSRmax thu c v hàm Multiplicative (170,44).
Hình 4.19. th SKTk c a cây tràm nh ng c p Dcv khác nhau c mô b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative).
1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30 SKTk (kg/cây)
Dcv(cm) Schumacher
1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30 SKTk (kg/cây)
Dcv(cm) Gompertz
1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30 SKTk (kg/cây)
Dcv(cm) Drakin-Vuevski
1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30
S-Curve SKTk (kg/cây)
Dcv(cm) 1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30
Multiplicative SKTk (kg/cây)
Dcv(cm)
1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30 SKTk (kg/cây)
Dcv(cm) Korf
Nh ng tính toán c ng cho th y, khi d oán SKTk theo c p Dcv c a cây tràm, thì 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative) u có sai sai l n hai c p Dcv u (1 và 2cm) và hai c p Dcv cu i (11 và 12cm)(B ng 4.21; Hình 4.20).
T nh ng th ng kê sai l ch c a 6 hàm cho th y, n u s d ng tiêu chu n R2max và SSRmin ch n hàm mô t SKTk theo c p Dcv, thì hàm Gompertz là hàm phù h p nh t. Vì th , hàm Gompertz c ch n mô t SKTk t ng ng v i nh ng c p Dcv. Theo ó, hàm Gompertz có d ng (Hình 4.21):
SKT(k)= 44,32462*exp(-15,41488*exp(-0,35859*Dcv)) (4.10) R2 = 96,18%; Se = ±2,01; MAE = 1,31; MAPE = 39,93%.
B ng 4.21. Sinh kh i thân khô c a cây tràmcajuputi theo c p Dcv c mô t b ng 6 hàm khác nhau
D(cm) Gompertz Schumacher Korf Drakin S-Curve Multiplicative
(1) (2) (3) (4) (5) (6) (7)
1 0,0 0,0 0,0 0,0 0,0 0,1
2 0,0 0,0 0,0 0,0 0,0 0,4
3 0,2 0,1 0,1 0,1 0,3 1,3
4 1,1 0,8 0,8 1,0 1,5 2,7
5 3,4 3,3 3,3 3,3 4,0 4,9
6 7,4 7,5 7,5 7,4 7,7 7,9
7 12,7 12,8 12,8 12,7 12,5 11,9
8 18,5 18,3 18,3 18,5 17,8 16,8
9 24,0 23,9 23,9 24,1 23,6 22,9
10 28,9 29,0 29,0 28,9 29,5 30,2
11 32,9 33,8 33,8 32,8 35,4 38,8
12 36,0 38,1 38,1 35,9 41,2 48,7
4.3.2.3. Xây d ng mô hình sinh kh i cành khô theo c p Dcv
K t qu làm phù h p 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative) v i SKCk nh ng c p Dcv khác nhau c a cây tràm c ghi l i b ng 4.22-4.23; hình 4.22 và ph l c 12.
0.0 10.0 20.0 30.0 40.0 50.0 60.0
1 2 3 4 5 6 7 8 9 10 11 12
Gompertz Schumacher Korf Drakin S-Curve Mult iplicat ive
SKTk (kg/cây)
Dcv(cm)
Hình 4.20. th SKTk c a cây tràm nh ng c p Dcv khác nhau c mô t b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative).
1 2 3 4 5 6 7 8 9 10
0 5 10 15 20 25 30
SKTk (kg/cây)
Dcv(cm)
Hình 4.21. th SKTk c a cây tràm nh ng c p Dcv khác nhau c mô t b ng hàm Gompertz.
B ng 4.22.Nh ng mô hình mô t sinh kh i cành khô c a cây tràm theo c p Dcv
TT Hàm Mô hình
(1) (2) (3)
1 Gompertz SKC(k) = 11,26295*exp(-7,71210*exp(-0,23137*Dcv)) 2 Drakin-Vuevski SKC(k) = 15,53585*(1-exp(-0,14986*Dcv))^4,28091 3 Schumacher SKC(k) = 733,92382*exp(-12,76688/Dcv^0,41324) 4 Korf SKC(k) = 733,92382*exp(-12,76688*Dcv^-0,41324) 5 S-Curve SKC(k) = exp(3,35060 - 17,01081/Dcv)
6 Multiplicative SKC(k) = 0,02606*Dcv^2,31525
B ng 4.23.Nh ng th ng kê sai l ch c a 6 hàm dùng mô t sinh kh i cành khô theo c p Dcv c a cây tràm
Th t Hàm R2 ±Se MAE MAPE SSR
(1) (2) (3) (4) (5) (6) (7)
1 Gompertz 92,55 0,51 0,33 42,97 7,98
2 Drakin-Vuevski 92,43 0,51 0,34 36,84 8,10
3 Schumacher 92,34 0,51 0,34 37,91 8,19
4 Korf 92,34 0,51 0,34 37,91 8,19
5 S-Curve 92,21 0,35 45,60 8,34
6 Multiplicative 92,18 0,51 0,33 44,55 8,37
Phân tích s li u b ng 4.23 cho th y, hàm Gompertz có h s R2max (92,55%), còn hàm Multiplicative có R2min (92,18%). Giá tr Se c a 6 hàm là không khác nhau (0,51). Giá tr MAEmin thu c v hàm Gompertz và hàm Multiplicative (0,33), còn MAEmax là 3 hàm S-Curve (0,35). Hàm Drakin-Vuevski nh n giá tr MAPEmin (36,84%), còn MAPEmax (45,60%) thu c v hàm S-Curve. Hàm Gompertz có SSRmin (7,98), còn SSRmax thu c v hàm Multiplicative (8,37).
Nh ng tính toán c ng cho th y, c 6 hàm (Schumacher, Korf, Drakin- Vuevsi, Gompertz, S-Curve và Multiplicative) u có sai s l n hai c p Dcv u (1 và 2cm) và hai c p Dcv cu i (11 và 12cm)(B ng 4.24; Hình 4.23).
Hình 4.22. th bi u di n SKCk nh ng c p Dcv khác nhau c a cây tràm c mô t b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S- Curve và Multiplicative).
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 SKCk (kg/cây)
Dcv(cm) Schumacher
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 SKCk (kg/cây)
Dcv(cm) Gompertz
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 SKCk (kg/cây)
Dcv(cm) Drakin-Vuevski
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6
S-Curve SKCk (kg/cây)
Dcv(cm) 1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6
Multiplicative SKCk (kg/cây)
Dcv(cm)
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 SKCk (kg/cây)
Dcv(cm) Korf
B ng 4.24. Sinh kh i cành khô c a cây tràmcajuputi theo c p Dcv c mô t b ng 6 hàm khác nhau
D(cm) Gompertz Schumacher Korf Drakin S-Curve Multiplicative
(1) (2) (3) (4) (5) (6) (7)
1 0,0 0,0 0,0 0,0 0,0 0,0
2 0,1 0,1 0,1 0,0 0,0 0,1
3 0,2 0,2 0,2 0,2 0,1 0,3
4 0,5 0,5 0,5 0,5 0,4 0,6
5 1,0 1,0 1,0 1,0 0,9 1,1
6 1,6 1,7 1,7 1,7 1,7 1,7
7 2,4 2,4 2,4 2,5 2,5 2,4
8 3,4 3,3 3,3 3,3 3,4 3,2
9 4,3 4,3 4,3 4,3 4,3 4,2
10 5,3 5,3 5,3 5,3 5,2 5,4
11 6,1 6,4 6,4 6,2 6,1 6,7
12 7,0 7,6 7,6 7,2 6,9 8,2
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0
1 2 3 4 5 6 7 8 9 10 11 12
Gompertz Schumacher Korf Drakin S-Curve Multiplicat ive
SKCk (kg/cây)
Dcv(cm)
Hình 4.23. th bi u di n SKC(k) theo c p Dcv c a cây tràm c mô t b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative).
T nh ng phân tích th ng kê sai l ch c a 6 hàm cho th y, n u s d ng tiêu chu n R2max và SSRmin ch n hàm mô t SKCk theo c p Dcv, thì hàm Gompertz là hàm phù h p nh t. Theo ó, SKCk theo c p Dcv c a cây tràmcajuputi c bi u th b ng hàm Gompertz có d ng (Hình 4.24):
SKC(k) = 11,26295*exp(-7,71210*exp(-0,23137*Dcv)) (4.11) R2 = 92,55%; Se = ±0,51; MAE = 0,33; MAPE = 42,97%.
4.3.2.4. Xây d ng mô hình sinh kh i lá khô theo c p Dcv
K t qu làm phù h p 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative) v i SKL(k) nh ng c p Dcv khác nhau c a cây tràm
c ghi l i b ng 4.25-4.26, hình 4.25 và ph l c 13.
B ng 4.25.Nh ng mô hình mô t sinh kh i lá khô c a cây tràm theo c p Dcv
TT Hàm Mô hình
(1) (2) (3)
1 Gompertz SKL(k) = 3,72228*exp(-5,20861*exp(-0,20409*Dcv)) 2 Drakin-Vuevski SKL(k) = 9,93823*(1-exp(-0,06220*Dcv))^2,14752 3 Schumacher SKL(k) = 72463,02*exp(-14,99518/Dcv^0,15294) 4 Korf SKL(k) = 72463,02*exp(-14,99518*Dcv^-0,15294) 5 S-Curve SKL(k) = exp(1,76202 – 11,550/Dcv)
6 Multiplicative SKL(k) = 0,03704*Dcv^1,71569
1 2 3 4 5 6 7 8 9 10
0 2 4 6 8 10 12
SKCt (kg/cây)
Dcv(cm)
Hình 4.24. th mô t SKCt theo c p Dcv c a cây tràm b ng hàm Gompertz.
B ng 4.26.Nh ng th ng kê sai l ch c a 6 hàm dùng mô t sinh kh i lá khô theo c p Dcv c a cây tràm
Th t Hàm R2 ±Se MAE MAPE SSR
(1) (2) (3) (4) (5) (6) (7)
1 Gompertz 76,17 0,35 0,21 72,38 3,85
2 Drakin-Vuevski 76,13 0,35 0,21 41,96 3,85
3 Schumacher 76,11 0,35 0,21 40,39 3,86
4 Korf 76,11 0,35 0,21 40,39 3,86
5 S-Curve 75,32 0,35 0,23 44,98 3,99
6 Multiplicative 76,07 0,35 0,21 57,05 3,86
Phân tích s li u b ng 4.26 cho th y, hàm Gompertz có h s R2max (76,17%), còn hàm Multiplicative có R2min (76,07%). Giá tr Se c a 6 hàm không có
Hình 4.25. th bi u di n SKL(k) c a cây tràm nh ng c p Dcv khác nhau c mô t b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative).
1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3 SKLk (kg/cây)
Dcv(cm) Schumacher
1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3 SKLk (kg/cây)
Dcv(cm) Gompertz
1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3 SKLk (kg/cây)
Dcv(cm) Drakin-Vuevski
1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3
S-Curve SKLk (kg/cây)
Dcv(cm) 1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3
Multiplicative SKLk (kg/cây)
Dcv(cm)
1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3 SKLk (kg/cây)
Dcv(cm) Korf
s khác nhau (0,35). T ng t , giá tr MAEmin c a 6 hàm c ng không khác nhau rõ r t (0,21-0,23). Hàm Gompertz nh n giá tr MAPEmax (72,38%), còn ah2m Schumacher và Korf có MAPEmin (40,39%). Hai hàm Gompertz và Drakin-Vuevski có SSRmin (3,85), còn SSRmax thu c v hàm S-Curve (3,99).
Nh ng tính toán c ng cho th y, c 6 hàm (Schumacher, Korf, Drakin- Vuevsi, Gompertz, S-Curve và Multiplicative) u có sai s l n hai c p Dcv u (1 và 2cm) và hai c p Dcv cu i (11 và 12cm)(B ng 4.27; Hình 4.26).
B ng 4.27. Sinh kh i lá khô c a cây tràmcajuputi theo c p Dcv c mô t b ng 6 hàm khác nhau
D(cm) Gompertz Schumacher Korf Drakin S-Curve Multiplicative
(1) (2) (3) (4) (5) (6) (7)
1 0,1 0,0 0,0 0,0 0,0 0,0
2 0,1 0,1 0,1 0,1 0,0 0,1
3 0,2 0,2 0,2 0,2 0,1 0,2
4 0,4 0,4 0,4 0,4 0,3 0,4
5 0,6 0,6 0,6 0,6 0,6 0,6
6 0,8 0,8 0,8 0,8 0,8 0,8
7 1,1 1,1 1,1 1,1 1,1 1,0
8 1,3 1,3 1,3 1,3 1,4 1,3
9 1,6 1,6 1,6 1,6 1,6 1,6
10 1,9 1,9 1,9 1,9 1,8 1,9
11 2,1 2,2 2,2 2,2 2,0 2,3
12 2,4 2,6 2,6 2,5 2,2 2,6
T nh ng phân tích th ng kê sai l ch c a 6 hàm cho th y, n u s d ng tiêu chu n R2max và SSRmin ch n hàm mô t SKLt theo c p Dcv, thì hàm Schumacher là hàm phù h p. Theo ó, mô hình mô t SKLt theo c p Dcv c a cây tràmcajuputi b ng hàm Schumacher có d ng (Hình 4.27):
SKL(k) = 3,72228*exp(-5,20861*exp(-0,20409*Dcv)) (4.12) R2 = 76,17%; Se = ±0,35; MAE = 0,21; MAPE = 72,38%.
0.0 0.5 1.0 1.5 2.0 2.5 3.0
1 2 3 4 5 6 7 8 9 10 11 12
Gompert z Schumacher Korf Drakin S-Curve Multiplicative
SKL(k) (kg/cây)
Dcv(cm)
Hình 4.26. th bi u di n SKL(k) theo c p Dcv c a cây tràm c mô t b ng 6 hàm (Schumacher, Korf, Drakin-Vuevsi, Gompertz, S-Curve và Multiplicative).
1 2 3 4 5 6 7 8 9 10
0 0.5 1 1.5 2 2.5 3
SKLk (kg/cây)
Dcv(cm)
Hình 4.27. th mô t SKL(k) theo c p Dcv c a cây tràm b ng hàm Gompertz.
4.3.3. Xây d ng mô hình sinh kh i khô c a cây tràm theo c p Dcv và H
K t qu phân tích h i quy t ng quan t ng b c (Ph l c 14-17) cho th y, khi s d ng hai nhân t Dcv(cm) và H(m) d oán TSK(k), SKT(k), SKC(k) và SKL(k), thì nh ng mô hình phù h p có d ng (Hình 4.28):
- i v i t ng sinh kh i khô c a cây tràmcajuputi
TSK(k) = 0,02494*D^1,22020*H^1,93104 (4.13) R2 = 97,8%; Se = ±1,88; MAE = 1,21; MAPE = 17,6%; P < 0,001.
- i v i sinh kh i thân khô c a cây tràmcajuputi
SKT(k) = 0,01218*D^1,13112*H^2,23340 (4.14) R2 = 97,9%; Se = ±1,49; MAE = 1,0; MAPE = 21,3%; P < 0,001.
- i v i sinh kh i cành khô c a cây tràmcajuputi
SKC(k) = 0,03154*D^2,49879*H^-0,26289 (4.15) R2 = 92,2%; Se = ±0,52; MAE = 0,34; MAPE = 51,2%; P < 0,001.
- i v i sinh kh i lá khô c a cây tràmcajuputi
SKL(k) = 0,00237*D^-0,31717*H^3,19088 (4.16) R2 = 83,8%; Se = ±0,29; MAE = 0,19; MAPE = 42,2%; P < 0,001.
0 5 10 15 20 25 30 35 40
0 5 10 15 20 25 30 35 40
TSK(tkthc nghim
TSK(t) oán
0 5 10 15 20 25 30 35 40
0 5 10 15 20 25 30 35 40
SKT(k)thc nghim
SKT(t) oán
0 1 2 3 4 5 6
0 1 2 3 4 5 6
SKC(k)thc nghim
SKC(t) oán
0 0.5 1 1.5 2 2.5 3
0 0.5 1 1.5 2 2.5 3
TSK(t)thc nghim
SKL(k) oán
Hình 4.28. th mô t sinh kh i khô th c nghi m và d oán