基于非线性混合模型的落叶松木材管胞长度模拟

摘要

以黑龙江省七台河市林业局金沙林场9株人工落叶松6825对早、晚材管胞长度样品数据为例,选择6个常用方程进行非线性回归分析,把拟合精度最高的Richards模型作为早、晚材管胞长度基础模型y=β1[1-exp(-β2x)]β3+ε。基于Richards模型,利用非线性混合模型技术构建落叶松早、晚材管胞长度混合效应模型yij=(β1+b1i){1-exp[-(β2+b2i)t]}β3+b3i+εij。结果表明:当对早材管胞长度进行拟合时,b1i、b2i、b3i同时作为随机参数时早材管胞长度模型拟合最好;当对晚材管胞长度进行拟合时,b1i、b2i、b3i同时作为随机参数时晚材管胞长度模型拟合最好;一阶自回归模型AR(1)能够较好地表达树木内误差相关性;同时考虑随机效应和时间序列相关性结构能够提高落叶松早、晚材管胞长度混合模型的预测精度。