Kendall系列非参数检验方法在区域水质年际变化趋势研究中应用-于雷－环境科学年会 .pptVIP

Kendall系列非参数检验方法在区域水环境质量年际变化趋势研究中的应用 环境保护部环境规划院 于雷 2008.12.27 为什么要用统计检验方法？ 传统方法：2001~2002、2002~2003、2003~2004和2004~2005污染物浓度变化趋势分别为升高、降低、升高和降低，降低和升高类数目相同，升or降？ 非参数统计检验：仅考虑数据的相对排列，受随机变量的基本概率分布和未检出值的影响较小 Kendall系列：Mann-Kendall（以下简称MK）、Seasonal-Kendall（以下简称SK）和Regional-Kendall（以下简称RK） MK检验原理和过程 统计检验值计算：假设随机变量X的数据样本与时间独立，没有趋势存在，定义 式中：xj和xk分别为第j年和k年的数值，n为数值的个数，并有 MK检验原理和过程 当统计检验值近似服从正态分布，对进行标准化处理和连续性修正，得到统计学意义上的趋势检验评价值 式中：Var(S)为分布的方差，计算式为 分异性检验：S为正，表明存在上升趋势；S为负，表明存在下降趋势；S为零，则没有趋势。 显著性检验：给定显著性水平 ，从标准正态分布表获得 ，应用双尾趋势检验，若| z |≤ ，无趋势，反之则有趋势存在。 RK /SK检验原理和过程 SK/RK检验首先对每个季节/点位进行MK检验，然后组合这些MK检验结果。 每个季节/点位的值相互之间比较，不与其它季节/点位的值相比较。 RK检验分异性和显著性检验过程与MK检验相同。 Kendall检验显著性评价和检验工具 高度显著 显著 适度显著 不显著 USGS（美国地质调查局）开发的计算机检验程序， 单断面水环境质量年际变化趋势分析 S=－6， =0.084 适度显著降低 MK非正态检验模式 多断面水环境质量年际变化趋势分析 S=－48 ， =0.0003 高度显著降低趋势 采用RK正态检验模式，检验出显著趋势 如果没有显著性趋势，怎么办？ 多断面分组水环境质量年际变化趋势分析 S=－4 ， =0.8162 不显著降低趋势 RK正态检验 RK检验适用正态检验最少样本量分析 时间跨度=5年，n=？，满足正态分布 步骤： 1、求取n=1时，不同Si对应组合的个数 Si=0，22； Si=2，20； Si=4，15； Si=6，9； Si=8，4； Si=10，1。负值组合个数与正值组合个数相对应。 2、求取n=2时，不同不同Si对应组合的个数 Si可能为0，2，4，6，8，10，12，14，16，18，10以及对应负值。 对于Si=0，可能为0+0，2-2，4-4，6-6，8-8，10-10。 对于0+0，组合可能个数为22*22 对于2-2，组合可能个数为20*20*2 求出所有组合累积，即可得Si=0时全部组合数。 同理，得出Si=其他结果的组合数。 3、求取每种组合个数占所有组合个数的比例，即为对应的真实概率，与对应正态概率对比。 4、累积真实概率，求取Si绝对值大于等于某值的累积概率，与正态累积概率相对比。 根据两误差大小判断是否近似认为Si分布符合正态分布。 RK检验适用正态检验最少样本量分析 RK分组检验 MK检验结果，总体呈“两头降低，中间升高”的特征，即上游段（江宁河口-九乡河口）和下游段（高港码头-小李港）断面呈不同程度降低趋势，中游段（小河口上游-青龙山段）断面呈不同程度升高趋势。 进一步的分组RK检验结果表明，上游段呈不显著降低趋势（S=－6，=0.3865），中游段呈不显著升高趋势（S=12， =0.1779），下游段呈不显著降低趋势（S=－10， =0.1190）。 * * 1.505 1.695 1.648 1.722 1.700 年均值 （mg/L） 2005 2004 2003 2002 2001 年份 1.505 1.695 1.648 1.722 1.700 年均值（mg/L