基于非线性混合模型的东北兴安落叶松树高和直径生长模拟
Height and Diameter Growth Modeling of Dahurian Larch Based on Nonlinear Mixed Model in Northeastern China
-
摘要: 以黑龙江省带岭林业局大青川林场80株人工兴安落叶松解析木数据为例,采用Richards生长模型作为基础模型,利用S-PLUS软件中的NLME过程,分别拟合非线性树高和直径生长模型。采用AIC、BIC、对数似然值和似然比检验等模型评价统计指标对不同模型的精度进行比较分析。结果表明:当对树高-年龄关系进行拟合时,b1、b3 同时作为混合参数时模型拟合最好;当对直径-年龄关系进行拟合时,b1 、b3 同时作为混合参数时模型拟合最好。把相关性结构包括一阶自回归结构AR(1)、一阶移动平均结构MA(1)及一阶自回归与移动平均结构 加入到树高和直径最优混合模型中,一阶自回归结构AR(1)显著提高了树高混合模型的拟合精度,一阶移动平均结构MA(1)显著提高了直径混合模型的拟合精度。模型检验结果表明:混合模型通过校正随机参数值能提高模型的预测精度。因此,混合模型在应用上不但能反映树高和直径的平均预测趋势,还能用方差协方差结构和误差相关性结构校正随机参数来反映个体之间的差异。Abstract: With the stem analysis data based on 80 sample trees from dahurian larch (Larix gmelinii Rupr.) plantations located in Daqingchuan Forestry Center of Dailing Forest Bureau in Heilongjiang Province, the Richards growth model was selected to model the height/age and diameter/age relationships of dahurian larch using NLME procedure of S-PLUS. Evaluation statistics, such as AIC, BIC, Log Likelihood and likelihood ratio test were used for precision analyzing and comparing of different models. The results indicated that Richards model with parameters b1 and b3 as mixed effects showed the best performance for both height-age relationships and diameter-age relationships. Correlation structures including first-order autoregressive correlation structure AR(1), moving average correlation structure MA(1) and autoregressive-moving average correlation structure were incorporated into the optimal height and diameter mixed models. AR(1) significantly improved the precision of mixed height model and MA(1) significantly improved the precision of mixed diameter model. Model validation showed that the mixed model with calibration of random parameters could improve the precision of prediction. Therefore, the application of the mixed model showed not only the mean trends of height and diameter prediction, but also the individual difference by calibrating random parameters using variance-covariance structure and correlation structure.
-
[1] 邓红兵, 郝占庆, 王庆礼, 等. 红松单木高生长模型的研究[J]. 生态学杂志, 1999, 18 (3): 19-22 [2] Huang S M, Titus S J. An individual tree height increment model for mixed white spruce-aspen stands in Alberta, Canada [J]. Forest Ecology and Management, 1999, 123:41-53 [3] Wang G G, Huang S M. Height growth pattern of white spruce in natural subregions in Alberta, Canada [J]. Forest Ecology and Management, 1999, 134:271-279 [4] 黄家荣,万兆溟. 马尾松人工林与距离有关的单木模型研究[J]. 山地农业生物学报, 2000, 19(1):10-15 [5] 刘兆刚,李凤日,于金成. 落叶松人工林单木模型的研究[J]. 植物研究,2003, 23(2):237-244 [6] 孟宪宇, 张 弘. 闽北杉木人工林单木模型[J]. 北京林业大学学报, 1996, 18(2):1-8 [7] 吕 勇. 杉木人工林生长率模型的研究[J]. 林业科学,2002, 38(1):146-149 [8] Fang Z, Bailey R L. Nonlinear mixed effects modeling for slash pine dominant height growth following intensive silvicultural treatments [J]. Forest Science, 2001, 47:287-300 [9] Calegario N, Daniels R F, Maestri R, et al. Modeling dominant height growth based on nonlinear mixed-effects model: a clonal Eucalyptus plantation case study [J]. Forest Ecology and Management, 2005, 204:11-20 [10] Nothdurft A, Kublin E, Lappi J. A non-linear hierarchical mixed model to describe tree height growth [J]. European Journal Forest Research, 2006, 125:281-289 [11] 李永慈,唐守正. 用Mixed 和Nlmixed 过程建立混合生长模型 [J]. 林业科学研究,2004, 17(3):279-283 [12] 李春明,张会儒. 利用非线性混合模型模拟杉木林优势木平均高 [J]. 林业科学,2010, 46(3):89-95 [13] Pinheiro J C, Bates D M. Mixed-effects Models in S and S-PLUS [M]. Springer, New York, 2000 [14] Calama R, Montero G. Interregional nonlinear height-diameter model with random coefficients for stone pine in Spain[J]. Canadian Journal Forest Research, 2004, 34:150-163 [15] Jordan L, Daniels R F, Clark A, et al. Multilevel nonlinear mixed effects models for the modeling of earlywood and latewood microfibril angle [J]. Forest Science, 2005, 51(4):357-371 [16] Vonesh E F, Chinchilli V M. Linear and Nonlinear Models for the Analysis of Repeated Measurements [M]. Marcel Dekker, New York, 1997
计量
- 文章访问数: 3573
- HTML全文浏览量: 155
- PDF下载量: 1908
- 被引次数: 0