收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

非线性混合效应模型及其在林业上应用

符利勇  
【摘要】:非线性混合效应模型(Nonlinear Mixed Effects Models,简称NLMEMs)是依据回归函数依赖于固定效应参数和随机效应参数的非线性关系而建立的。它是分析纵向数据、多水平数据及重复调查数据的近代统计学方法,既能反映总体的平均变化趋势,又能刻画个体差异的统计分析工具。近几年,NLMEMs受到越来越多的学者关注,并且被引入到多门学科,如医学、工学、农学及林学等。经过30多年的发展历程,学者已提出了单水平和嵌套多水平NLMEMs的多种参数估计方法,典型的计算软件有SAS和S-Plus。但在实际应用中发现2个问题,一是主流软件SAS和S-Plus计算经常不收敛,尤其是模型中待估参数较多时最为明显,二是现有的混合模型没有包括所有类型的随机效应的组合(例如交互作用,这在林业中是常用的),因而限制了模型的应用。本文的目的在于(1)提出一种收敛性好的计算方法来分析单水平和嵌套多水平NLMEMs;(2)提出一个包括所有随机效应类型的NLMEMs的统一表达式,并且给出一种参数计算方法;(3)完成上述2个算法的程序代码并在ForStat上实现;(4)应用我们的程序解决一个林学上的实际问题,该问题用已有的程序是无法解决的。 本研究实现了上述4个目的。具体内容如下:1)根据一阶条件期望线性化—数学期望极大方法(First-order conditional expectation linearization–expectation maximation,简称FOCE-EM)的理论推导出计算单水平和嵌套多水平NLMEMs的计算公式、设计了计算流程。2)本研究提出了一种正态NLMEMs的标准表达式,它包括了正态NLMEMs的所有随机效应类型,给出该模型相应的一种参数估计方法,即线性逼近—逐步2次规划算法。3)发现并通过实例说明SAS中mixed模块不能保证随机效应参数方差为非负定矩阵,所以在本文提出采用线性逼近—逐步2次规划算法,同时给出十种方差类型满足正定或半正定的条件,因而算法可以保证不出现类似SAS中的错误。4)模型的标准表达式可以处理固定效应和随机效应参数分级(即数量化问题),并且指出了S-Plus中nlme函数分级差法计算有缺陷,而线性逼近—逐步2次规划算法克服了此问题。5)首次利用带有交互作用的两因素(林分密度和地位级指数)NLMEMs分析了落叶松树高—直径模型。在此基础上进一步分析了随机效应与海拔的关系。 通过本项研究,可以得出以下主要结论:1)本研究提出了一种正态NLMEMs的标准表达式,它包含多种类型的非线性混合效应模型(随机效应参数服从正态分布),例如单水平NLMEMs、逐级嵌套多水平NLMEMs、只含主效应的多因素NLMEMs、包括主效应和交互效应的NLMEMs以及某几种类型组合的一般性NLMEMs等。模型中固定效应参数和随机效应参数可以考虑分级(即数量化)。同时还把正态NLMEMs的标准表达式推广到参数方差与某些因素(称为组变量)有关的NLMEM(s考虑组变量的NLMEMs)。因此模型比传统的NLMEMs表达式更为一般化,具有更广的用途。2)在计算单水平或逐级嵌套的多水平NLMEMs时,FOCE-EM算法与SAS及S-Plus提供的Lindstrom andBates(LB)算法计算精度非常接近,数值实例表明,至少四位小数相同。但FOCE-EM算法从理论上保证了线性步计算收敛,从而使得该算法计算收敛性明显好于LB算法。3)本研究给出线性逼近—逐步2次规划算法计算一般类型的NLMEMs。为保证方差非负定采用逐步2次规划是必要的。该算法同样也能计算考虑组变量的NLMEMs。4)从计算速度上讲,FOCE-EM算法比线性逼近—逐步2次规划算法要快,因此建议在计算单水平和逐级嵌套多水平NLMEMs时,使用FOCE-EM算法,而计算其它类型的NLMEMs时,使用线性逼近—逐步2次规划算法。5)通过两因素NLMEMs对落叶松树高—直径模型研究得出,考虑林分密度与地位级的交互作用能明显提高模型的预测精度,而且把海拔高度作为组变量时,还可以进一步提高模型预测精度。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 ;三天制成落叶松树胶,今后可节省外汇二百万元[J];林业实用技术;1958年09期
2 徐丽,张庆国;Holling圆盘方程参数估计方法的研究[J];安徽农业大学学报;1993年02期
3 ;土面增温剂在落叶松育苗上的试验[J];吉林林业科技;1977年01期
4 杨立群;从进境落叶松树苗上检出落叶松枯梢病[J];植物检疫;2003年05期
5 田志诚,于瑞安,刘新田,张晓放;兴安落叶松树冠发生及生长特点的研究[J];林业科技通讯;1997年02期
6 王增辉;窦森;贾红震;;农药残留函数y=ae~(bt)参数估计的一种新方法[J];农业环境保护;1997年01期
7 许俊利;何学凯;;林木生物量模型研究概述[J];河北林果研究;2009年02期
8 史景钊;杨星钊;陈新昌;;3参数威布尔分布参数估计方法的比较研究[J];河南农业大学学报;2009年04期
9 王怀明;高立刚;;对落叶松人工林森林经营技术的探讨[J];林业勘查设计;2011年02期
10 延晓冬,赵士洞;崔-Lawson和Logistic方程参数的优化估计方法[J];应用生态学报;1991年03期
11 张连翔;曹振喜;刘宗相;;Weibull分布的参数估计及在油松人工林直径分布规律研究中的应用[J];辽宁林业科技;1992年02期
12 黄俊洋,唐建忠,李洪瑞;塑料薄膜围扎落叶松树干防治落叶松毛虫技术[J];防护林科技;2001年01期
13 王福林;王吉权;;非对称养分效应模型参数估计的一种优化回归组合法[J];农业工程学报;2006年02期
14 万飚;高仕春;付湘;纪昌明;;水文频率分析的线性矩法实用化研究[J];中国农村水利水电;2009年06期
15 ;落叶松人工幼林生长与气溫湿度降水土壤含水量关系的初步观測[J];东北林业大学学报;1975年02期
16 唐谦,顾万春;几种落叶松树种叶绿体DNA的分化及其系统演化关系的含义[J];林业科学;1995年04期
17 张元跃;;一种获得非线性生长函数参数初始值的新方法:杂种遗传算法[J];生物数学学报;2010年04期
18 高翔;王若平;夏长高;朱茂桃;;随机变量多重Weibull统计模型及参数最优估计[J];农业机械学报;2006年11期
19 朱敏;冯仲科;;树木生长方程参数估计方法的比较研究[J];北京林业大学学报;2007年S2期
20 蔡细平;沈月琴;;持续性分析用于条件价值评估法数据统计分析——以浙江省临安市CVM调查数据为例[J];东北林业大学学报;2009年01期
中国重要会议论文全文数据库 前10条
1 郭仲英;陶宏;史彦江;;天山地区落叶松种和种源选择试验研究[A];中国林木遗传育种进展[C];1991年
2 文新辉;牛明洁;;基于神经网络理论的一种新的时间序列模型的参数估计方法[A];全国青年管理科学与系统科学论文集(第2卷)[C];1993年
3 潘永泉;;基于神经网络的图像序列参数估计方法[A];1994中国控制与决策学术年会论文集[C];1994年
4 李政;罗飞路;;基于匹配追踪的超声回波参数估计方法[A];2009中国仪器仪表与测控技术大会论文集[C];2009年
5 刘桂英;李传庆;邵联合;;一种递推最小二乘估计简化算法[A];第十届全国电工数学学术年会论文集[C];2005年
6 杨世永;李宏伟;何水明;;乘性和加性噪声中二维谐波的参数估计[A];第十二届全国信号处理学术年会(CCSP-2005)论文集[C];2005年
7 赵明旺;吕勇哉;;随机连续线性系统参数估计的Laguerre多项式方法[A];1993中国控制与决策学术年会论文集[C];1993年
8 葛新科;胡保生;;参数估计及自校正控制的Systolic算法[A];1995中国控制与决策学术年会论文集[C];1995年
9 刘浩;吴季;;雷达目标特征提取的数学模型与RELAX参数估计方法[A];第十三届全国遥感技术学术交流会论文摘要集[C];2001年
10 郭垒;王向东;王智峰;;基于参数估计器的一类时滞系统自适应模糊控制[A];2005中国控制与决策学术年会论文集(下)[C];2005年
中国博士学位论文全文数据库 前10条
1 符利勇;非线性混合效应模型及其在林业上应用[D];中国林业科学研究院;2012年
2 丁建明;车辆动力学性能参数估计方法研究[D];西南交通大学;2012年
3 王兰美;极化阵列的参数估计和滤波方法研究[D];西安电子科技大学;2005年
4 凌丹;威布尔分布模型及其在机械可靠性中的应用研究[D];电子科技大学;2011年
5 周颖颖;住房抵押贷款强度定价模型研究[D];大连理工大学;2009年
6 孙永梅;α稳定分布参数估计与谱分析理论及应用研究[D];大连理工大学;2006年
7 李粤;引文网络的可调优先粘贴模型及其应用[D];清华大学;2007年
8 孙晓东;混沌噪声背景下谐波参数估计方法研究[D];吉林大学;2009年
9 臧嘉捷;Meta分析中时间效应量参数估计方法的建立及应用[D];第二军医大学;2012年
10 于晓辉;噪声背景下chirp信号参数估计理论与方法研究[D];吉林大学;2007年
中国硕士学位论文全文数据库 前10条
1 朱孝武;扩谱信号参数估计与同步算法研究[D];电子科技大学;2004年
2 米亚静;贝叶斯数据挖掘算法在反垃圾邮件中的研究[D];辽宁工程技术大学;2007年
3 吕晓蕊;Alpha稳定分布的模型仿真及参数估计[D];华中科技大学;2008年
4 方明;有相似产品信息情形下的可靠性评定方法[D];延边大学;2007年
5 朱玲;地、海杂波特征建模与参数估计[D];电子科技大学;2009年
6 石屹然;有色噪声背景下极化阵列信号参数估计方法研究[D];兰州理工大学;2010年
7 赵雪骞;电力系统抗差状态估计与参数估计方法研究[D];重庆大学;2012年
8 项凤雏;基于多时段量测信息的独立线路参数估计方法研究[D];重庆大学;2010年
9 詹旻;电网支路参数估计方法及其应用研究[D];哈尔滨工业大学;2011年
10 曹雪;正弦非线性切普键控技术及参数估计方法[D];哈尔滨工程大学;2013年
中国重要报纸全文数据库 前10条
1 本报记者 纳兰海贝 通讯员 孙胜利;马凤玉的五年维权路[N];西部时报;2006年
2 本报通讯员  李瑞雪 本报记者  段丽茜 徐华;红枣局长的“甜蜜事业”[N];河北日报;2006年
3 杜建立;尖山:承包造林利场利民[N];黑河日报;2008年
4 吴庚秀;岫岩百名残疾人栽种万亩自强林[N];鞍山日报 ;2008年
5 贾宏谱 孙阁;陈德廷:为了绿色阜平 我愿付出全部[N];中国绿色时报;2006年
6 辛语 刘玉荣;扑火人员分赴北部原始森林火场[N];内蒙古日报(汉);2006年
7 吴庚秀;岫岩2万农户水上林下赚大钱[N];鞍山日报 ;2008年
8 韩彦 建刚;15年前的承包合同有瑕疵咋办?[N];丹东日报;2005年
9 记者 范洪庄;林权改革给我市农村家庭经营带来巨变[N];抚顺日报;2007年
10 李江天;林改切中抚顺农民生产经营“兴奋点”[N];辽宁日报;2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978