生物量建模过程中剔除扰动数据方法探讨.pdfVIP

2012年 3月 华东森林经理 Mar．2012 第26卷 第 1期 EastchinaForestManagement Vo1．26NO．1 生物量建模过程中剔除扰动数据方法探讨 陈文灿 张伟东 申屠惠良 郑春茂 (国家林业局华东林业调查规划设计院 浙江杭州 310019) 摘要：采集生物量建模数据过程中，因季节、天气变化、立地条件，以及外业调查人员 对数据采集过程的把握，诸多因素均会对采集到的数据产生影响，导致部分数据异常。 异常数据称为扰动数据，若参与建模，会引起建模结果的异常，应该剔除。从残差分析 的角度，对如何从众多数据中过滤出扰动数据提出了解决方案。 关键词：扰动数据；残差分析；三倍标准差 中图分类号：TP391 文献标识码 ：B 文章编号：1004—7743(2012)01-0061-04 1数据准备 为样木树高； 为样木材积。 表 1 生物量建模样木特征数据 生物量建模数据的准备分为外业数据采集 和 内业数据处理两个阶段。外业数据采集阶段， 按树根、树干、树枝、树叶四个部位，采集样 木的鲜重数据 ，以及样木各部位样 品的鲜重数 据。内业数据处理阶段，主要是测定样木各部 位样品干重与含水率，根据样木鲜重、样品鲜 2．1建立回归方程 重与样 品干重，计算样木整体干重，即生物量 。 使用普通最小二乘法 (OLS)，对生物量建 此处，样木指调查主体；样品指样木的一 模样木数据进行多元非线性回归计算，获得模 部分，用来测定生物量。生物量的具体计算过 程如下： 型的参数a、b，由此建立回归方程，根据此 4 回归方程，计算模型估计值。 ． M=∑ x(z／wi) i=l 2．2残差分析 其中 为样木总体生物量； 为样木各 部位鲜重；w 为样木各部位样 品鲜重；Z为 生物量建模数据中，样木实测生物量称为 样木各部位样 品干重。 观察值，根据回归模型计算的生物量称为估计 经过外业调查和 内业数据处理，获得某个 值，观察值和估计值之差，称为残差。为了消 树种的生物量总体建模数据，实际样木数量 51 除负数对分析过程 的影响，将残差的平方作为 统计量。 株，各样本平均值、最大值、最小值的特征数 据见表 1。 表达式为：E =(M 一M )，其中：E 为单株样木的残差平方； 为生物量的观察 2扰动数据分析 值 ； 为生物量的估计值 。 选择生物量模型形式：M =a(D 1 ， 残差平方的平均数： ： E ／力 J_ ‘ 其中a、b为模型参常数；D 为样木胸径； i=l