多元线性回归计算程序——SHARP PC—1500计算机在林业上的应用

一、数学方法简述 1.多元线性回归的数学模型 x,,x2.……x*,为自变量(i二1,2……p),y为因变量。x;与y之间线性回归方程为: Y为因变量原始数据标准化矩阵。 又知xXT乘积结果为自变量相关矩阵,这样就得到r,j,xy乘积求得r三,,故建立起正规方程组(1)式。 用全主元高斯消元法解正规方程组(l)式,求到了b:y=b。+b lx,+bZx:+……+b*x;为:回归常数,b,、bZ……b;为回归系数。3.求多元回归系数b0根据最小二乘法,选择这样的b。,b,,bZ…… /\b;,使艺(y一y)’=万〔y一(b。+b,xl+…+b;x;)〕2达到极小。为此将上式分别对b。,b:……,b:求偏导数,令期等于0,…     

林下微生境对草坪草叶绿素含量的影响

对林下微生境和全光照下的5个草坪草品种叶绿素含量及叶绿素a/b值进行了对比试验研究.结果表明:林下微生境使得草坪草叶绿素含量及叶绿素a/b值均发生规律性的变化,叶绿素含量林下的明显高于全光照下的,而叶绿素a/b值则相反.根据这两个指标进行分析,普通、欧宝和优异对林下微生境的适应性较强,派尼和高羊茅猎狗5适应性中等,午夜较差.     

基于线性混合模型的落叶松枝条长度和角度模型

以黑龙江省五营林业局丽林林场30株人工落叶松2190个枝条长度和角度数据为例,利用逐步回归技术建立落叶松枝条长度和角度模型:BL=b1+b2DINC+b3DINC2+b4DBH·DINC2,BA=b1+b2DINC+b3DINC2+b4DBH·DINC。利用S-PLUS软件中的LME模块,考虑树木效应拟合线性枝条长度和角度模型。采用AIC、BIC、对数似然值和似然比检验等模型评价统计指标对不同模型的拟合效果进行比较分析。结果表明:当拟合枝条长度和角度模型时,b1,b2,b3同时作为混合参数时模型拟合最好。为了描述混合模型构建过程中产生的异方差现象,把幂函数和指数函数加入到枝条长度和角度混合模型中。指数函数显著提高了枝条长度混合模型的拟合效果,幂函数显著提高了角度混合模型的拟合效果,并且消除了异方差现象。模型检验结果表明:混合模型通过校正随机参数值能提高模型的预测精度。因此,混合模型在应用上不但能反映总体枝条长度和角度预测,而且能通过方差协方差结构校正随机参数来反映树木之间的差异。     

基于混合模型的落叶松微纤丝角模型

以黑龙江省七台河市林业局金沙林场9株人工落叶松2 790个样品数据为例,选择6个常用方程进行非线性回归分析,把拟合精度最高的修正Logistic模型作为微纤丝角基础模型y=b1/[1+ exp(b2x)]+b3,然后,利用S-PLUS软件中的NLME过程,拟合非线性微纤丝角模型.采用AIC、BIC、对数似然值和似然比检验等模型评价统计指标对不同模型的精度进行比较分析.结果表明:当对微纤丝角-年龄关系进行拟合时,b1,b2,b3同时作为混合参数时模型拟合效果最好.把相关性结构包括复合对称结构(CS)、一阶自回归结构AR(1)、一阶移动平均结构MA(1)及一阶自回归与移动平均结构[ARMA(1,1)]加入到微纤丝角最优混合模型中,一阶自回归与移动平均模型[ ARMA(1,1)]显著提高了微纤丝角混合模型的拟合精度.模型检验结果表明:混合模型通过校正随机参数值能提高模型的预测精度.因此,混合模型在应用上不仅能反映总体微纤丝角预测,而且能通过方差协方差结构和误差相关性结构校正随机参数来反映个体微纤丝角差异.     

基于线性混合模型的落叶松枝条基径模型

以黑龙江省五营林业局丽林林场30株人工落叶松2 190个枝条基径数据为例,利用逐步回归技术建立了落叶松枝条基径模型:BD=b1+b2DINC+b3DINC2+b4DBH.DINC2。然后,利用S-PLUS软件中的LME过程,拟合线性枝基径模型。采用AIC、BIC、对数似然值和似然比检验等模型评价统计指标对不同模型的拟合效果进行比较分析。结果表明:当拟合枝条基径模型时,b1、b2、b3同时作为混合参数时模型拟合最好。为了矫正混合模型构建过程中产生的异方差现象,把幂函数和指数函数加入到枝条基径混合模型中。指数函数显著提高了枝条基径混合模型的拟合效果,并且消除了异方差现象。模型模拟表明:对于大小相同树木,枝条基径随着着枝深度(DINC)的增加而增大,对于大小不同的树木,枝条基径随着胸径(DBH)的增加而增大。林木的胸径变量很好地反映了不同大小树木的枝条基径的变化。在不知道详细林分信息的前提下,可以利用树木变量合理地预测兴安落叶松人工林的枝条基径的变化规律。     

3种红豆杉光合特性及叶绿素荧光差异比较

为红豆杉资源的合理开发和人工栽培提供理论基础,以4年生的东北红豆杉Taxus cuspidata、南方红豆杉T.chinensis ziar.mairei和曼地亚红豆杉T.media为试验材料,比较了这3种红豆杉叶色素含量、光合及叶绿素荧光特性的差异.结果表明:3种红豆杉中叶绿素a、叶绿素b、总叶绿素和类胡萝卜素含量,...     

盆景环境对对节白蜡叶片生理及显微结构的影响

以对节白蜡盆景为试材,地栽对节白蜡为对照,测定了叶绿素a、叶绿素b、类胡萝卜素、可溶糖和可溶蛋白为生理指标,同时分析了其显微结构.结果表明:盆景叶片叶绿素a、b与类胡萝卜素含量无明显变化,但可溶性总糖与可溶性蛋白含量显著提高;同时叶片更为厚实,栅栏组织排列紧密,层数较多;盆景环境(修剪与盆栽)下的对节白蜡叶片的生理与结...     

‘84K’杨组氨酸激酶基因PaHK3b的克隆及功能分析

[目的]本研究克隆了银腺杨‘84K’(Populus alba×P. glandulosa ‘84K’)组氨酸激酶基因PaHK3b启动子及编码区,并对其表达进行检测及功能鉴定,为深入研究PaHK3b基因在杨树生长发育的调控作用提供线索,为杨树分子育种及品种改良奠定基础。[方法]根据毛果杨(P. trichocarpa Torr.Gray)基因组信息,设计引物克隆‘84K’杨组氨酸激酶基因PaHK3b启动子及CDS序列,并对其保守结构域和启动子顺式作用元件进行分析;同时,对‘84K’杨进行植物激素处理(10μmol·L~(-1)ABA、10μmol·L~(-1)6-BA、10μmol·L~(-1) IBA、10μmol·L~(-1)GA3及10μmol·L~(-1)水杨酸(SA))及非生物胁迫处理(42℃高温、0℃低温、200 mmol·L~(-1) NaCl和5%PEG6000),利用实时定量PCR(qRT-PCR)方法检测PaHK3b基因表达情况与表达响应差异,并采用原核表达方法初步确定PaHK3b基因的生物学功能。[结果]PaHK3b基因编码框区长度为3 060 bp,编码1 019个氨基酸,PaHK3b蛋白具有CHASE、HisKA和REC等典型的细胞分裂素受体结构域。PaHK3b基因启动子序列中不仅含有大量TATA框和CAAT框常见核心元件,还包含低温响应元件LTR、防御与胁迫响应元件TC-rich repeats、赤霉素响应元件GARE-motif、水杨酸响应元件TCA-element等顺式作用元件,这些元件与杨树的激素响应和逆境胁迫响应密切相关。qRT-PCR分析表明:PaHK3b基因在叶片中表达最高,根部中等,茎中最少;另外,与正常条件下相比,在高温、低温、NaCl及PEG处理时,PaHK3b基因表达量与对照明显增高,分别为对照的2.67、2.61、2.28、1.87倍;用IBA诱导处理时,基因表达量与对照相比差异不大,而在6-BA、ABA、GA3及SA处理时,基因表达量与对照相比均呈下调表达;在添加5%PEG6000的LB液体培养基中,转入PaHK3b基因原核表达载体的大肠杆菌菌株生长速度显著高于对照,在添加50～150 mmol·L~(-1)NaCl的LB固体培养基上,转入PaHK3b基因原核表达载体的大肠杆菌菌株单克隆生长均好于对照。[结论]‘84K’杨PaHK3b基因启动子含有逆境和激素响应元件,表明PaHK3b基因与杨树植物激素信号及非生物胁迫信号响应密切相关。经非生物胁迫处理、激素处理及原核表达证实,杨树PaHK3b基因参与杨树植物激素信号响应,并在其抗逆境胁迫过程中发挥重要调控作用。     

油茶种子数量相关基因EPFL2的鉴别

拟南芥(Arabidopsis thaliana)AtEPFL2基因属于EPF/EPFL基因家族中的EPFL1/2/3亚家族,与种子数量有关.对普通油茶(Camellia oleifera)幼叶、未成熟种仁、花蕾、根的转录组数据进行挖掘,得到3条油茶EPFL2基因序列;从油茶基因组3代测序数据中鉴别出含有前述3条序列的基因组DNA序列;最后,得到3条油茶EPFL2基因全长CDS序列,分别命名为CoEPFL2a、CoEPFL2b、CoEPFL2c.用MEGA7构建EPFL1/2/3亚家族蛋白质系统进化树,结果表明:CoEPFL2a、CoEPFL2b、CoEPFL2c属于EPFL2分支,并且CoEPFL2a、CoEPFL2c与AtEPFL2在同一个亚分支.用SignalP-5.0预测蛋白信号肽与酶切位点,结果表明:AtEPFL2、CoEPFL2c含信号肽与酶切位点,CoEPFL2a、CoEPFL2b不含;CoEPFL2c可能与AtEPFL2具有同样的功能,即通过调节胚珠密度影响种子数量.     

甲维盐和阿维菌素对杨树叶片叶绿素含量的影响

在田间,分别用不同浓度的杀虫剂阿维菌素和甲维盐50μL处理3年生107杨树幼苗叶片,24 h后测定幼苗叶片的叶绿素含量。结果表明:阿维菌素和甲维盐均导致杨树幼苗叶片叶绿素a、叶绿素b和叶绿素总量下降,类胡萝卜素含量、叶绿素a/b比值上升。表明2种药剂虽然对杨树幼苗会产生轻微的胁迫作用,但通过改变类胡萝卜素的含量和叶绿素a/b比值,增强对作物的抗逆性,消除杀虫剂所造成的胁迫。     

抗虫基因CryIA（b）植物表达载体的构建

苏云金芽孢杆菌杀虫晶体蛋白基因具有对鳞翅目昆虫的毒杀作用,因此被广泛用于转基因研究,以提高植物的抗虫能力.故以含有CryIA(b)基因的PKUB质粒为模板,利用PCR技术克隆出抗虫基因[CryIA(b)],将其构建到PGEMT-Easy上,获得重组质粒PGEMT-CryIA(b),并测定全部序列,将PGEMT-CryIA(b)用BamHI 和SmaI酶切后插入表达载体PBI121的CaMV35S启动子和NOS终止子之间,构建成准备用于丽江云杉等针叶树种遗传转化CryIA(b)基因的植物表达载体PBI-CryIA(b),采用液氮冻融法将其转化到根癌农杆菌LBA4404和C58中,为下一步的转化工作奠定了基础.     

喀斯特森林植被自然恢复中凋落物现存量 及其碳库特征演化

采用空间代替时间的方法,研究了茂兰退化喀斯特森林自然恢复中凋落物现存量及其碳库特征.结果表明:随群落恢复凋落物现存量呈减少趋势,早期减幅大,中后期减幅小;随群落恢复凋落物现存量演化受地貌因子、群落生活史、土壤微生物的影响较大,受水热条件的影响较小;凋落物现存量生物量与其分解失重率、表层土壤微生物量碳有较强的负相关关系;凋落物分解失重率与表层土壤微生物量碳呈极强的正相关关系;随群落恢复凋落物现存量含碳率变化不显著,凋落物现存量碳密度变化规律与凋落物现存量变化规律一致,单位面积凋落物碳储量表现为碳源效应,且早期碳源效应较强、固碳能力不稳定,中后期碳源效应较弱、固碳能力稳定;凋落物碳密度(Y)与其现存生物量(x)的关系:Y=b.+b1x(b0、b1为常数).     

利用2种非线性混合效应模型(2水平)对杉木林胸径生长量的分析

选择2种2水平非线性混合模型对杉木林胸径生长量进行分析,其中模型1为一般的2水平非线性混合模型,模型2在模型1的基础上进一步考虑固定效应参数随某一特定因子水平变化而变化。本文通过对这2种模型的分析,首先确定构建2水平非线性混合模型的基础模型,然后对模型1衍生出的665种模型及模型2衍生出的2703种模型进行计算和比较:对于模型1,有57种模型计算收敛,当形参b0同时考虑区组和样地效应、而b4和b5只考虑区组效应时,模型拟合效果最好,因此把该模型作为模型1的最佳拟合模型;对于模型2,有24种模型计算收敛,当形参b5同时考虑区组和样地效应、b1只考虑区组效应并且固定效应b0的取值与各区组水平有关时,模型拟合效果最好,因此把该模型作为模型2的最佳拟合模型。最后对传统的非线性回归模型、模型1及模型2进行比较:模型1和模型2的拟合效果都比传统的非线性回归模型好,且模型2的拟合效果最好。     

基于L*a*b*颜色空间对木材分类的研究

L*a*b*均匀颜色空间具有等距性和色差高分辨力的特点,非常适合色差较小情况下的颜色测量和比较.木材材色分布范围较窄,利用L*a*b*颜色空间中的颜色特征表示木材表面颜色,有利于木材材色之间的比较和划分.基于L*a*b*颜色空间,提取了东北常见五种树种木材图像的颜色特征进行分类研究,通过仿真试验得到了满意的分类结果.     

红楠播种育苗试验研究

红楠为中国亚热带常绿阔叶林区,处于Ⅱb层即第二林层b亚层的常绿阔叶乔木树种,在侧方遮阴条件下生长良好,被列为生态公益林树种及园林工程景观绿化树种.对红楠随采随播和芽苗移栽进行试验研究,结果表明:从起苗到移栽的间隔时间限于3h之内;从播种开始到移栽时的芽苗生长期30 d为最佳时间;采用50孔穴盘规格54 cm X27 cm x9 cm(深度),生产红楠小苗,较为合适.     

灵石山米槠林优势种不同叶龄叶绿体色素沿海拔梯度的变化

对灵石山米槠林优势种不同叶龄叶绿体色素的差异以及优势种叶绿体色素沿海拔梯度的变化规律进行分析。结果表明:除少数优势种在某海拔梯度以外,不同叶龄叶绿素a、叶绿素b和叶绿素(a+b)存在显著或极显著差异,一般在中高海拔梯度差异最大,且2年生叶片含量高于1年生叶片;而不同叶龄间叶绿素a/b和类胡萝卜素的差异性表现复杂,表现了植物对环境的适应。不同优势种叶绿素[叶绿素a、叶绿素b和叶绿素(a+b)]在海拔梯度间存在显著差异,并且沿海拔梯度的变化呈现一定的规律性,米槠林优势种中有7种优势种不同叶龄叶绿素[叶绿素a、叶绿素b、叶绿素(a+b)]含量均随海拔的升高先升高,一般在A5(442m)或A6(531m)达到最大值,然后随海拔的升高逐渐降低,而叶绿素a/b随海拔升高变化趋势比较多样化。类胡萝卜素作为功能色素由于其功能比较复杂,对海拔的响应规律也不尽相同。海拔作为重要的环境因素对于植物叶绿体色素产生显著影响,不同的优势种沿海拔的变化规律存在差异,说明其对环境的适应对策不同。     

竹子捕光叶绿素a/b结合蛋白基因全长的克隆和序列分析

采用RT-PCR技术与RACE技术,从绿竹中克隆到捕光叶绿素a/b结合蛋白基因全长cDNA,命名为cab-BO1(GenBank登记号:EF061137).该基因长1 102 bp,从第64～861 bp含有1个开放阅读框,编码265个氨基酸.cab-BO1编码的氨基酸中第64～232位包括典型的捕光叶绿素a/b结合蛋白功能域,第9～11位为蛋白激酶C的磷酸化位点,第27～52位为泛素结构功能域信号,第55～58位为酪氨酸激酶(CK2)磷酸化位点.序列相似性分析结果表明:cab-BO1 DNA序列和编码的氨基酸序列与玉米、小果野蕉、小麦、水稻等的cab基因及其氨基酸序列的相似性分别在85%和89%以上,显示了cab家族基因在进化过程中的保守性.     

具有缓和曲线圆曲线的二种测设方法

本文介绍具有缓和曲线圆曲线的2种测设方法:a.分布测设缓和曲线和圆曲线;b.在直缓点或缓直点1次设站连续测设近端缓和曲线和圆曲线。同时给出2种测设方法的切线支距法和偏角法的公式。     

肥料配置对国槐幼苗生长的影响

为探明不同肥料配置对国槐幼苗生长的影响,以不施肥料为对照,设置了6种不同氮磷钾配比,以研究国槐幼苗苗高、地径、冠幅、叶绿素含量和光合参数.结果表明:不同配料配置处理下国槐的苗高、地径和冠幅均显著高于对照,施肥处理的叶绿素a、b含量,净光合速率、气孔导度、胞间二氧化碳浓度和蒸腾速率均显著高于对照,氮磷钾比例为2:2:1时...     

中国产的兰属植物

中国人通常所谓的国兰,指的就是兰属的地生种类,如传统的春兰、蕙兰(九华)、寒兰、墨兰(报岁兰)和建兰(四季兰),和近年比较火爆的春剑和莲瓣兰。全世界约有67种,其中49种见于我国。分属于3个亚属:兰亚属、大花亚属和建兰亚属。1.纹瓣兰Cym b.aloifolium(L.)Sw.产广东、广西、贵