聚类分析与稳健统计方法对CNAS T0402数据处理比较和分析.docVIP

聚类分析与稳健统计方法对CNAS T0402数据处理比较和分析 　　摘要建立了一种基于对能力验证数据进行聚类分析的方法，克服了以往采用稳健统计方法处理的层次性和深度不足的缺陷，可将参加实验室按测试结果及相似性进行分组，便于对检测过程存在问题的查找和分析，也有利于检测实验室对自身测试水平的认识和评价。 　　关键词聚类分析；稳健统计方法；计量；对比 　　中图分类号 O212.1 文献标识码A文章编号 1007-5739（2009）06-0258-06 　　 　　聚类分析（Cluster Analysis）是进行数据统计分析的一类重要工具[1-5]，广泛应用于商业、经济、医药等领域，通过数据分类获得有用的统计信息，是多元统计分析的一个重要分支，根据分类对象的数量指标，定量地确定分类对象之间的相似关系并进行分类。聚类分析中应用广泛的是层次聚类法和K-Means算法。能力验证，作为我国实现量值溯源并确保分析实验室数据可靠、准确一致的重要工具，发挥着日益广泛和深入的作用[6]。在数据处理上一般采用稳健统计方法进行评价[7-11]。 　　稳健统计方法是世界各国计量机构在组织能力验证或国际比对中常用的统计方法，该方法在数据非正态分布或存在离群值时优势明显，Z比分数为其评价的重要指标[12]。其缺点是提供信息较少，缺乏深入和细致的描述，无法描述实验室数据相关性之间的联系。当测试对象为多元素或多组分时，只能单个项目分析，缺乏整体判断和描述，无法适应比对内容与对象不断增加的趋势[13]。 　　本文尝试将聚类分析应用在能力验证数据处理中，结合稳健统计方法，以获得更加细致全面的数据信息。通过把不同实验室测量数据划分为不同的集合和小组，方便对能力验证数据的分析和比对原因的查找，可对不同实验室的共性问题，如测量方法的比较等进行深入的探讨，为化学计量领域理论和实践中深入处理实验数据提供方法基础和探讨依据。通过文献检索，目前国内外尚无类似报道。 　　本文将聚类分析与稳健统计方法结合起来，对CNAS T0402粮食中重金属铅、镉含量测定的数据进行处理。中国合格评定国家认可委员会于2008年组织的全国范围的一次能力验证，由中国计量科学研究院组织实??，其目的是了解我国目前粮食检测机构实际的测量水平和潜在问题。研究表明，该方法可以分层次、多角度对数据进行分析，具有广泛的应用价值。 　　 　　1材料与方法 　　 　　1.1能力验证样品的制备 　　此次能力验证样品选用河南省重金属污染区生长的小麦，经清理除杂，去除糠、麸皮以及子实、胚芽，用鄂式破碎机初碎，将候选物在70℃烘干24h，去除水分。用高铝球磨机研磨36~48h，使样品99%以上通过80目。制备好的样品装入清洁塑料桶中，封口保存。经均匀性检验合格后混匀分装。用Co60辅照灭活，置于干燥阴凉处保存。 　　1.2样品的均匀性和稳定性 　　样品采用等离子体发射光谱法（ICPOES）和等离子体质谱法（ICPMS）进行均匀性和稳定性检验。在均匀性检验中，随机抽取15瓶样品，每瓶在上部和底部取2个样品。每个样品准确称取0.2g，加入5mL浓硝酸，用微波消解炉进行消解。完毕待冷却后，转移到洁净塑料瓶中，加入In或Y内标，上机测试。经F分析，样品中Pb、Cd均匀性良好。经t检验，样品稳定性检验良好。 　　1.3推荐方法 　　此次能力验证推荐方法为GB/T5009.12-2003食品中铅的测定方法和GB/T5009.15-2003食品中镉的测定方法。 　　1.4数据统计分析方法和软件 　　该研究选用SPSS软件进行数据处理。分层聚类分析子模块，聚类法采用组间连接法（Between Groups Linkage），测度方法选择欧氏距离平方（Squared Euclidean Distance），即2项之间的距离是每个变量值之差的平方和[14]。 　　 　　2结果与讨论 　　 　　2.1采用分析方法的统计情况 　　在此次能力验证中，共采用4种测量方法，分别是原子吸收分析方法（AAS）、等离子体发射光谱法（ICPOES）、等离子体质谱法（ICPMS）、原子荧光光谱法（AFS）（见表1）。综合分析，AAS、ICP、ICPMS之间的比例接近7:1:1。这说明原子吸收作为国家标准GB/T5009.12-2003和GB/T5009.15-2003的推荐方法，在实验室分析领域仍占据相当的比重。 　　 　　2.2正态分布检验 　　稳健统计方法适合多种分布情况的数据分析，尤其分析非正态分布数据具有很大的优势。采用SPSS软件进行分析，铅元素和镉元素检验结果分别为Asymy.sig=0.10和Asymy.sig=0.02。进一步采用峰态系数和偏态系数法数据进行检验[15]。将数据汇总后，按从小到大