An Improved Hight Prediction Mixed Effect Model for Quercus variabilis

doi:10.13466/j.cnki.lyzygl.2018.05.010

Abstract

Abstract:

Tree height and DBH are two important tree-building factors in the survey of forest resources.In view of the relatively difficulty for height measurement,in this article we took 94 cork oaks from Xixi Town,Jinhua City as the research object and compared the fitting effects of 10 tree height curve models commonly used in forestry,improved the traditional Gompertz tree height curve model,proposed and constructed a Gompertz mixed effect tree high prediction model with site factors.The results show that:(1) The model fits best when the random parameters b₁ and b₃ are introduced for Gompertz mixed effect model.The model fits best when the random parameters b1 and b4 are introduced for Gompertz mixed effect improvement model.(2) The Gompertz mixed effect improvement model constructed in this paper has a coefficient of determination of 0.779 with 0.553 Gompertz mixed effect model,and with 0.542 for Gompertz model.That is building the model only by the mixed effect method is not obvious for improving the accuracy of the prediction model.(3) Experiments show that the Gompertz mixed effect improvement model constructed in this paper greatly improves the prediction accuracy of the height of Quercus variabilis,and provides a new method for studying the relationship between tree height and diameter at breast height.

Key words: Cork oak, tree height prediction, site factor, Gompertz model, mixed effect

CLC Number:

S758
S792.18

HUANG Feng, XU Aijun, TANG Lihua. An Improved Hight Prediction Mixed Effect Model for Quercus variabilis[J]. FOREST RESOURCES WANAGEMENT, 2018, 0(5): 54-62.

Figures/Tables 11

Tab.1

Tab.2

Tree height curve empirical model

序号	模型形式	作者/文献
(1)	$H = β 0 - β 1 × exp (- β 2 × DBH)$	Banin (2012)/[8]
(2)	$H = β 0 + β 1 / (β 2 + DBH)$	Tang(1991)/[11]
(3)	$H = β 0 × (1 - exp - β 1 × DBH)$	Meyer(1940)/[3]
(4)	H=β₀×exp(-β₁×exp $- β 2 × DBH$ )	Gompertz(1992)/[9]
(5)	$H = β 0 × exp (- β 1 × DB H β 2)$	Lundqvist-Korf(1939)/[12]
(6)	$H = β 0 + β 1 × (1 - exp - β 2 × DBH)$	Mitscherlish(1919)/[13]
(7)	$H = 10 β 0 + β 1 DBH$	Schumacher(1939)/[14]
(8)	$H = β 0 × exp (β 1 / (DBH + β 2))$	Ratkowsky(1990)/[11]
(9)	H=DBH²/(β₀+β₁×DBH+β₂×DBH²)	Curtis(1967)/[11]
(10)	$H = β 0 × exp (- β 1 / DBH)$	寺崎(1915)/[13]

Tab.2

Fig.1

Tab.3

Estimation values and fitting statistics of 10 tree height curve models of Quercus variabilis

参数、统计量/指标	模型编号
参数、统计量/指标	(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$β 0$	18.656	22.376	18.642	18.350	0.005	0.073	1.340	22.253	0.001	21.859
$β 1$	18.584	-4.582	3.465	1.543	-8.280	18.583	-0.077	-0.196	0.007	0.178
$β 2$	3.449	0.146	—	4.325	0.045	3.449	—	0.021	0.047	—
RMSE	2.906	2.945	2.906	2.887	3.113	2.906	2.912	2.942	2.938	2.912
SSE	264.4	268.0	264.4	262.7	283.2	264.4	267.8	267.8	267.4	267.8
R²	0.539	0.533	0.539	0.542	0.506	0.539	0.533	0.533	0.534	0.533

Tab.3

Fig.2

Fig.3

Tab.4

Tab.5

Gompertz mixed effect model simulation

模拟次数	随机参数	参数个数	AIC	BIC	R²
基本模型	无	3	0.184	-3.375	0.542
模拟1	$b 1$	4	0.159	-3.500	0.553
模拟2	$b 2$	4	0.159	-3.500	0.553
模拟3	$b 1, b 3$	5	0.159	-3.500	0.553
模拟4	$b 2$ , $b 3$	5	-0.257	-3.816	0.468
模拟5	$b 1, b 2, b 3$	6	-0.256	-3.816	0.468

Tab.5

Tab.6

Improved Gompertz mixed-effect model simulation

模拟次数	随机参数	参数个数	AIC	BIC	R²
基本模型	无	4	22.570	-0.569	0.596
模拟1	$b 3$	5	26.430	3.291	0.777
模拟2	$b 4$	5	26.881	3.742	0.709
模拟3	$b 1, b 3$	6	26.931	3.793	0.703
模拟4	$b 1, b 4$	6	26.471	3.333	0.779
模拟5	$b 2, b 3$	6	26.465	3.327	0.726
模拟6	$b 3, b 4$	6	22.589	-0.550	0.598
模拟7	$b 1, b 2, b 3$	7	26.921	3.782	0.702
模拟8	$b 1, b 3, b 4$	7	21.994	-1.144	0.512
模拟9	$b 2, b 3, b 4$	7	26.921	3.782	0.702
模拟10	$b 1, b 2, b 3, b 4$	8	26.454	3.315	0.778

Tab.6

Tab.7

Model parameters,variance estimates and fitting statistics

参数/ 统计量	参数	Gompertz 模型	Gompertz混合效应模型	本文改进的Gompertz 混合效应模型
固定参数	$β 1$	18.35	-33.039	3.853
	$β 2$	1.543	30.984	13.263
	$β 3$	4.325	-168.259	4.643
	$β 4$	—	—	0.053
拟合统计量	Bias	-0.0009	0.0001	0.4542
	RMSE	2.887	2.911	1.118
	R²	0.542	0.553	0.779

Tab.7

Fig.4

References 28

[1]	李令, 郑道爽. 栓皮栎的特征特性及栽培技术[J].现代农业科技, 2010(5):189-189.
[2]	王冬至, 张冬燕, 张志东, 等. 基于非线性混合模型的针阔混交林树高与胸径关系[J]. 林业科学, 2016,52(01):30-36.
[3]	韦雪花, 王佳, 冯仲科. 北京市13个常见树种胸径估测研究[J]. 北京林业大学学报, 2013,35(5):56-63.
[4]	Ram P.Sharma J B. Modeling height-diameter relationships for Norway spruce,Scots pine,and downy birch using Norwegian national forest inventory data[J]. Forest Science and Technology, 2015,11(1):44-53. doi: 10.1080/21580103.2014.957354
[5]	Duan Guangshuang, Gao Zhigang, Wang Qiuyan, et al. Comparison of Different Height-Diameter Modelling Techniques for Prediction of Site Productivity in Natural Uneven-Aged Pure Stands[J]. Forests, 2018,9(2):63. doi: 10.3390/f9020063
[6]	董云飞. 福建杉木人工林单木生长模型的研究[D]. 北京:北京林业大学, 2015.
[7]	Saunders M R, Wagner R G. Height-diameter models with random coefficients and site variables for tree species of Central Maine[J]. Annals of Forest Science, 2008,65(2):203-203. doi: 10.1051/forest:2007086
[8]	Kearsley E, Moonen P C, Hufkens K, et al. Model performance of tree height-diameter relationships in the central Congo Basin[J]. Annals of Forest Science, 2017,74(1):7. doi: 10.1007/s13595-016-0611-0
[9]	Marcos A S S, Luciana F A, Simone A V, et al. Height-diameter relationships of tropical Atlantic moist forest trees in southeastern Brazil[J]. Scientia Agricola, 2012,69(1):26-37. doi: 10.1590/S0103-90162012000100005
[10]	朱光玉, 罗小浪. 湖南栎类天然混交林优势木树高曲线哑变量模型研究[J].林业资源管理, 2017(4):22-29.
[11]	Fang Zixing, Bailey R L. Height-diameter models for tropical forests on Hainan Island in southern China[J]. Forest Ecology & Management, 1998,110:315-327.
[12]	Krisnawati H, Ades P K. Generalized Height-Diameter Models for Acacia mangium Willd.Plantations In South Sumatra[J]. Praxis, 2010,98(21):1219-24. pmid: 19844880
[13]	李婉婷. 基于混合效应的兴安落叶松树高曲线模型[D]. 哈尔滨:东北林业大学, 2014.
[14]	李政龙, 孟京辉. 基于Cross-Validation的福建杉木人工林最优树高曲线模型筛选研究[J]. 西北林学院学报, 2016,31(1):92-96.
[15]	樊伟, 许崇华, 崔珺, 等. 基于混合效应的大别山地区杉木树高-胸径模型比较[J]. 应用生态学报, 2017,28(9):2831-2839.
[16]	祖笑锋, 李秋实, 倪成才, 等. 非线性混合效应生长模型的拟合、随机效应预测和应变量预测间对应关系[J]. 林业科学, 2016,52(10):72-79. doi: 10.11707/j.1001-7488.20161009
[17]	王明初, 孙玉军. 基于混合效应模型及EBLUP预测杉木树高生长过程[J]. 浙江农林大学学报, 2017,34(5):782-790.
[18]	高慧淋, 董利虎, 李凤日. 基于混合效应的人工落叶松树冠轮廓模型[J]. 林业科学, 2017,53(3):84-93.
[19]	Pinheiro J C, Bates D M. Mixed-Effects Models in S and S-PLUS[M] //Mixed-effects models in S and S-PLUS.Springer, 2000: 1135-1136.
[20]	Fang Zheng. Nonlinear mixed effects modeling for slash pine dominant height growth following intensive silvicultural treatments[J]. Forest Science, 2001,47(47):287-300.
[21]	Peng Changhui, Zhang Lianjun, Zhou Xiaolu, et al. Developing and evaluating tree height-diamter models at three geographic scales for black spruce in Ontario.[J]. Northern Journal of Applied Forestry, 2004,21(2):83-92. doi: 10.1093/njaf/21.2.83
[22]	Özçelik R, Yavuz H, Karatepe Y, et al. Development of ecoregion-based height-diameter models for 3 economically important tree species of southern Turkey[J]. Turkish Journal of Agriculture & Forestry, 2014,38(3):399-412.
[23]	许崇华, 崔珺, 黄兴召, 等. 基于线性混合效应模型的杉木树高-胸径模型[J]. 西北农林科技大学学报:自然科学版, 2017,45(6):53-60.
[24]	Yuan Lianzhen, Yang Bin, Shi Fuqiang, et al. Effect of Different Slope Directions on Aluns nepalensis and Betula alnoides Afforestation[J]. Journal of Sichuan Forestry Science & Technology, 2015.
[25]	Temesgen H, Zhang Chunyu, Zhao Xiuhai. Modelling tree height-diameter relationships in multi-species and multi-layered forests:A large observational study from Northeast China[J]. Forest Ecology & Management, 2014,316(316):78-89.
[26]	Lin Jinhua. Type,Aspect and Slope Position of Seedlings of Cunninghamia lanceolatain Early Stage[J]. Protection Forest Science & Technology, 2015.
[27]	RupåiYs P. Height-diameter models with stochastic differential equations and mixed-effects parameters[J]. Journal of Forest Research, 2015,20(1):9-17. doi: 10.1007/s10310-014-0454-1
[28]	Ryan T P. Modern Regression Methods[M] //Modern regression methods.Wiley, 2009: 1315-1315.

变量	最小值	最大值	标准差	均值	偏度	峰度
坡度/(°)	20.73	42.73	3.21	36.05	-1.50	5.76
坡向	1	3	1.01	1.98	0.04	-2.04
坡位	1	3	0.70	1.51	1.02	-0.25
海拔/m	169	521	41.70	475.89	-4.49	30.99
树高(H)/m	8	19	2.48	14.50	-0.28	-0.63
胸径(DBH)/cm	25	100	16.40	47.17	1.19	1.37

检验指标	(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
ME	0.0001	-0.0003	0.0002	0.0057	-0.3003	-0.00008	-0.0390	1.4600	0.0070	-0.0036
MAE	1.1996	1.2168	1.1995	1.1907	1.2804	1.1996	1.2108	1.9769	1.3922	1.2152
M%E	-0.0162	-0.0164	-0.0162	-0.0161	-0.0389	-0.0162	-0.0191	-0.0156	-0.0245	-0.0166
MA%E	0.0910	0.0920	0.0910	0.0904	0.0985	0.0910	0.0920	0.0920	0.0927	0.0920