学位论文 > 优秀研究生学位论文题录展示
混合效应模型在森林生长模拟研究中的应用
作 者: 李春明
导 师: 唐守正
学 校: 中国林业科学研究院
专 业: 森林经理学
关键词: 优势木平均高 林分断面积 林分蓄积 单木直径生长量 线性混合效应模型 非线性混合效应模型 联立方程组系统
分类号: F224
类 型: 博士论文
年 份: 2010年
下 载: 130次
引 用: 0次
阅 读: 论文下载
内容摘要
混合效应模型包括线性混合效应模型和非线性混合效应模型两种,目前这一方法在医学、农业、经济、林业及其他领域有了广泛的应用。森林在生长发育过程中,受地理或物理环境、立地条件、气候变化、树种特性及林分本身的影响,在地域上或样地上存在着生长上的差异。以往在模拟森林生长及收获预估时,建立的模型基本上都没有考虑样地之间或地域之间存在的差异性。森林生长及收获模型采用的重复测量数据是林业研究中十分常见的数据类型。针对这种数据的分析方法,经过多年的发展,已经有了很大改善,尤其是通过拟合混合效应模型,相对于用一般的多变量分析而言,模型建立更为灵活,结果更为可信。但是,一般以往文献中应用混合效应模型多是针对单一反应变量重复测量的情况,数据之间的相关性也仅存在于同一个体的多次测量结果之间。本研究是在了解混合效应模型原理方法的基础上,运用分布在江西省不同地区的杉木数据,研究混合效应模型在重复测量资料分析中的应用,建立与完善其模型,并进行实例分析。结果认为该方法不仅可以更加有效、深入的挖掘该类数据中蕴藏的信息,得到固定效应和随机效应的估计值,同时可以得到反应变量本身之间的相关系数和变量多次重复测量间的相关系数,对数据间的相关性进行更加细致的分解和研究。另外,也可拟合不同的混合效应模型,通过比较拟合统计量的大小和协方差参数的多少,从中选择针对该数据最切合实际解释的统计分析模型。本文首先总结了混合效应模型在森林生长模型中的应用,然后考虑区域效应、样地效应及两层效应,对混合效应模型在模拟江西省杉木林分的优势木平均高、林分单位断面积、林分单位蓄积及单木直径生长量中的应用进行研究。本文主要的研究内容、结论和创新如下:(1)优势木平均高:以江西省杉木样地数据为例,利用4个常用的Richards和Schumacher模型形式,分别考虑样地效应、区域效应以及同时考虑样地和区域效应,在此基础上考虑模型的异方差及多次观测的时间序列相关性。研究结果表明,无论是考虑样地效应或区域效应,还是同时考虑样地效应和区域效应,模型的拟合效果要好于传统的最小二乘模拟方法。当考虑样地效应时,在这4个方程中,Fang (2001)提出的修改后的Richards模型模拟的效果最好。在考虑模型的误差方差协方差矩阵时,指数函数形式的异方差结构及AR(1)的时间序列相关结构矩阵形式模拟效果最好。当考虑区域效应时,Fang提出的修改后的Richards模型模拟的效果最好。在考虑模型的误差方差协方差矩阵时,指数函数形式的异方差结构及AR(1)的时间序列相关结构矩阵形式模拟效果最好。同时考虑样地效应和区域效应时,结果为与单独考虑一个效应相比,提高了模拟的精度,并且差异显著,在考虑误差方差协方差结构时,指数函数的异方差方程和AR(1)的时间序列相关结构矩阵形式模拟效果最好。利用新的验证数据进行优势木平均高的预测也证实了上述结论。(2)林分断面积:选择常用的Richards和Schumacher方程来模拟林分的断面积,结果发现自变量包括林分公顷株数、林龄及优势木平均高的Schumacher型断面积模型模拟精度最高。因此在此模型基础上构建非线性混合效应模型。当考虑林分地貌特征时能够提高模型的估计精度,因此在构建断面积模型时地貌特征作为定性因子加入到模型中。在构建混合模型时首先考虑所有参数都分别作为混合参数,然后多个参数作为混合参数进行模拟。研究结果表明,在考虑样地效应时,G=exp(b0+b01*dm+b1/t)Nb2+b4(?)Hb3+u(?)+b31*dm+(b5+u(?))(?)+εi模拟效果最好。当考虑区域效应时,G=exp(b0+b01*dm+b1/t)Nb2+u(?)+b4/tHb3+u2i+b31*dm+b5/t+εi的模拟效果最好。同时考虑样地效应和区域效应时,G=exp(b0+b01*dm+b1/t)Nb2+b4/tHb3+u(?)+v(?)ij+b31*dm+(b5+v2ij)/t+εi的模拟效果最好。无论是考虑样地效应、区域效应还是两者同时考虑,当考虑方程的异方差结构时,幂函数形式模拟效果最好;当考虑时间序列相关性时,AR(1)结构的模拟效果最好。最后进行方差分析表明,同时考虑异方差和时间序列相关性时,比单纯考虑异方差或时间序列相关性模拟效果都好,且差异显著。最后,对3种效应的模拟效果和最小二乘方法的模拟精度进行验证并选择新的样地进行预测。(3)林分蓄积:分两种情况进行模拟,一种是单纯混合模型进行模拟,一种是基于混合效应模型的优势木平均高、林分断面积和林分蓄积的联立方程组模型进行模拟。结果为,优势木平均高是联立方程组中最基本的组成部分,其准确预测是联立方程组的关键之一。优势木平均高和林分断面积是林分蓄积估计的主要误差来源。在模拟中考虑三个因变量的自相关性,在优势木平均高和林分断面积模型中考虑参数的随机效应后则蓄积的随机效应可忽略。最后对未来林分的优势木平均高、林分断面积和林分蓄积进行预测。(4)林分单木直径生长量:采用线性混合效应模型方法进行研究,在构建模型时,考虑单木大小、立地条件、林分因子及竞争对单木直径生长量的影响。在构建混合模型时分两种情况,一种是截距的样地效应、区域效应及两者的结合;一种是随机参数的样地效应、区域效应及两者的结合。在构建线性混合模型时两种情况都考虑了异方差和时间序列相关性的误差方差协方差结构。最后结果表明:①在考虑截距效应或随机参数效应时,无论考虑样地效应、区域效应或两种效应的结合,模拟效果都要好于传统最小二乘方法。②在考虑样地的截距效应或考虑区域的截距效应时,选择指数函数作为异方差形式及选择ARMA(1,1)的时间序列相关结构时,模拟效果最好。同时考虑样地和区域的截距效应时,选择幂函数作为异方差形式及选择ARMA(1,1)的时间序列相关结构时,模拟效果最好。③在考虑随机效应时:分别考虑样地效应或区域效应的影响时,则包括截距在内的所有参数同时作为混合参数时的模拟精度最高。在同时考虑样地效应和区域效应影响时,的模型模拟效果最好。最后利用验证数据分不同情况进行验证。
|
全文目录
摘要 5-8 ABSTRACT 8-13 第一章 绪论 13-31 1.1 引言 13-29 1.1.1 研究的目的和意义 13-14 1.1.2 混合效应模型简介 14-20 1.1.3 混合模型在森林生长模型中的应用 20-29 1.2 研究的主要目的和内容 29-30 1.2.1 主要目的 29 1.2.2 主要内容 29-30 1.3 研究的技术路线 30-31 第二章 研究材料和方法 31-38 2.1 研究区域概况 31 2.2 杉木简介 31-32 2.3 数据来源 32-38 2.3.1 样地调查因子 32-33 2.3.2 样地主要因子的计算 33-36 2.3.3 样地数据分布 36-38 第三章 基于混合模型的杉木林优势木平均高研究 38-68 3.1 方法 39-42 3.1.1 基础模型 39-40 3.1.2 非线性混合效应模型 40-42 3.1.3 模型精度比较 42 3.2 结果与分析 42-66 3.2.1 样地效应模拟结果 42-50 3.2.2 区域效应模拟结果 50-58 3.2.3 两层次模拟结果 58-63 3.2.4 随机参数估计 63-64 3.2.5 模型预测 64-66 3.3 小结 66-68 第四章 基于混合模型的杉木林分断面积模拟研究 68-89 4.1 断面积模型的选择 69 4.2 传统最小二乘模拟结果 69-71 4.3 非线性混合效应模型的模拟 71-88 4.3.1 样地效应 71-77 4.3.2 区域效应 77-82 4.3.3 两层效应 82-86 4.3.4 模拟预测 86-88 4.4 小结 88-89 第五章 基于联立方程混合效应模型的蓄积研究 89-103 5.1 单纯蓄积模型的构建 89-93 5.1.1 传统最小二乘模拟结果 90-91 5.1.2 混合效应模型模拟结果 91-92 5.1.3 样地内方差协方差矩阵 92-93 5.1.4 模型验证 93 5.2 基于非线性混合效应的蓄积模型的联立系统 93-98 5.2.1 联立模型方程组的确定 93-94 5.2.2 混合参数的确定 94-96 5.2.3 考虑时间序列相关性 96-98 5.3 模型预测 98-101 5.4 小结 101-103 第六章 线性混合模型在模拟单木胸径生长量中的应用 103-128 6.1 研究方法 104-106 6.1.1 影响单木生长的因子 104-106 6.1.2 模型自变量的选择 106 6.1.3 模型验证 106 6.2 结果与分析 106-126 6.2.1 传统回归方法结果 106-107 6.2.2 截距效应 107-116 6.2.3 随机参数效应 116-125 6.2.4 模型验证 125-126 6.3 小结 126-128 第七章 结论与讨论 128-132 7.1 结论 128-130 7.1.1 优势木平均高 129 7.1.2 断面积 129-130 7.1.3 蓄积联立方程组 130 7.1.4 单木直径生长量 130 7.2 讨论 130-132 参考文献 132-139 在读期间的学术研究 139-141 致谢 141
|
相似论文
- 批准兽药数据库的建立与恩诺沙星的回顾性药动数据分析,S851.66
- 群体PK中的非线性混合效应模型及SDE模型研究,R911
- 中国健康人群及造血干细胞移植患者他克莫司的群体药动学研究,R96
- 异位妊娠患者甲氨蝶呤的群体药动学及治疗结果预测因素分析,R969.1
- 左旋多巴群体药动学模型探讨及多巴丝肼合用多巴胺受体激动剂治疗帕金森病临床疗效观察,R969.1
- 环孢素的临床群体药物动力学研究,R96
- 丙戊酸在癫痫患儿中的群体药代动力学研究,R742.1
- 卫生Ⅷ项目效果评价,R173
- 左旋多巴与普拉克索的临床药学研究及帕金森统一评分量表的改进,R96
- 广义半参数模型的参数估计与影响分析,O212.1
- 中国肝移植患者他克莫司(FK506)的群体药物动力学研究,R96
- 重复测量线性混合模型在医学研究中的应用,R195
- GLMM与考虑重复效应的ANOVA在重复测量数据中应用,R195
- 面向故障过程的多设设可靠性分析与维修决策,TH186
- 番鸭体重生长非线性混合效应模型评价的研究,S834
- 利奈唑胺在中国人群的群体药代/药效动力学研究,R96
- 冠心病心绞痛血瘀证的药效学解构与动态演变的多维度诠释,R259
- 群体药代/药效动力学数据管理和分析的方法学研究,R96
- 非正态及非线性重复测量资料分析模型及其医学应用,R181.3
- K公司计划及预测改进对于合理库存配置的研究,F224
- 我国上市公司董事长更换与业绩的实证研究,F224
中图分类: > 经济 > 经济计划与管理 > 经济计算、经济数学方法 > 经济数学方法
© 2012 www.xueweilunwen.com
|