基于两层次线性混合效应模型的杉木林单木胸径生长量模型.pdfVIP

第 3 期 李春明:基于两层次线性混合效应模型的杉木林单木胸径生民量模型 67 国内一些学者在单木生长模型方面开展了一些 研究(杜纪山， 1999; 黄家荣等， 2000; 如l 兆刚， 2004; 王文斗等，2005) ，所用的建模数据通常来自 于林分内多个样地在一定时间间隔内多次的重复观 测:首先，这些数据之间存在着序列相关性，且相邻 测量时间的数据相关性较大，而离得较远的数据相 关性减小;其次，由于分布区域、竞争环境和立地条 件不同，以及林木个体之间的差异，树木之间存在着 生长上的差异;再次，这些数据之间不满足独立性 要求，丽且存在着异方差性。以往这些问题经常被 忽略，这样就可能会造成很大的估计和预测误差，给 林业生产经营和决策部门提供错误的信息和依据。 混合效应模型方法被广泛地用来解决这类问题，该 方法主要优点在于它不仅考虑固定效应，也考虑随 机效应，使得信息利用更为充分;不仅可以满足常 规要求的重复测量数据，而且还可用于说明分析变 量水平及其变化趋势;对个别观察单位重复测量次 数不完全相同的情况，也同样有效。 国外很多学者研究了混合模型在模拟和预测单 木胸径生长中的应用: Adame 等 (2008 )利用截距效 应混合模型方法研究西班牙西北部比利牛斯烁 ( Quercus pyrenacia) 的胸径生长，考虑了林木大小、 林木竞争、立地条件和生物地理气候等因子; Uzoh 等(2006)利用多层次线性混合模型研究美国西部 美国黄松 (Pinus ponderosa) 同龄林分的单木胸径生 长，与 Adame 等(2008 )不同的是没有考虑生物地理 气候对单木生长的影响; Rafael 等 (2005 )在建立石 松 ( Pinus p阳α) 的胸径生长量模型时，与 Adame (2008 )的方法和考虑的因子相同; Weiskittel 等 (2007 )在分析花旗松( Pseudotsuga 阳nz!es!~ var. menziesii) 、异叶铁杉 ( Tsugα helerophylla )和美国赤 杨 ( Alnω rubra) 3 个树种的胸径生长时，利用混合效 应模型方法，并考虑了气候和经营相互作用的影响; Kieman 等 (2008) 在研究不同采伐方式下糖枫 (Acer saccharum)异龄林的单木胸径生长时，用混合效应 模型来解释重复测量产生的时间序列相关性。这些 学者在利用线性混合效应模型模拟单木胸径生长 时，均认为考虑随机效应后，模型精度会有很大的 提高。 线性混合效应模型包括随机截距效应和随机系 数效应 2 种( LiUell et al忖 2000) ，上述学者利用线性 混合效应模型方法在构建单木胸径生长模型时都采 用了随机截距效应。到目前为止，利用随机系数效 应线性混合模型方法构建单木胸径生长模型的研究 还未见诸于文献。本研究根据所用不同区域和不同 样地数据之间存在着嵌套性的特点，采用 SAS 软件 进行模拟。在模拟过程中利用参数随机效应方法分 别对最套的样地层次、区域层次及同时考虑 2 个层 次进行模拟，并考虑异方差和序列相关性问题，利用 LRT( 似然比检验)进行方差分析，最后利用均方根 误差、平均绝对残差和确定系数对独立的验证数据 进行验证。 1 研究数据和方法 1.1 数据来源 数据来源于江西省的人工杉木( Cunninghamia lanceolata) ( 杉木树种蓄职占样地总蓄视 809毛以上) 同定样地调查。共计 82 个区域 365 块杉木样地，从 中选择了 5 416 株树木(树木起测胸径为 5 cm ，从 1991 年开始每 5 年测量 1 次，共 4 次，4500 株树木 作为模拟数据，其余树木为验证数据)进行研究。 样地面积为 0.08 hm2 D 样地基本调查因子包括单 木胸径、海拔、平均年龄、平均胸径、祥木总株数等。 样地汇总因子统计结果见表 1 。 褒 1 江西省人工杉木 5 年间跚跚单木胸径 生长量统计 (N=16248) Tab.l τbe result of 目皿mary 民atisti岱 w据d to develop th啤 S-year periodic annual diameter increment model zor fir plantation in Jiangxi Provinc，鸣~ (N = 1(248) 测量因子 最大值最小值平均值 标准差 Measure factor Max. Min. Mean SD 海拔 1150 21 394.5 213.8 Eleva