由于格鲁布斯Grubbs检验法使用了所有数据的平均值和标准偏差课件

地球化学数据处理地球化学数据处理 地球化学数据处理1第一节第一节 地球化学数据的误差分析地球化学数据的误差分析1 系统误差 特点：影响准确度，不影响精密度影响准确度，不影响精密度 （1 1）对分析结果的影响比较恒定，可以测定对分析结果的影响比较恒定，可以测定 和校正；和校正；（2 2）在同一条件下，重复测定，重复出现；在同一条件下，重复测定，重复出现；（3 3）影响准确度，不影响精密度；）影响准确度，不影响精密度；（4 4）可以消除。）可以消除。一、误差的分类一、误差的分类第一节 地球化学数据的误差分析1 系统误差一、误差的分类2 产生的原因：产生的原因：（1）方法误差选择的方法不够完善例例：重重量量分分析析中中沉沉淀淀的的溶溶解解损损失失，滴滴定定分分析析中中指指示示剂剂选选择不当择不当 （2）试剂误差所用试剂有杂质例例：去离子水不合格；试剂纯度不够去离子水不合格；试剂纯度不够 （3）仪器误差仪器本身的缺陷 例：天平两臂不等，砝码未校正；滴定管容量瓶未校正例：天平两臂不等，砝码未校正；滴定管容量瓶未校正 （4）主观误差操作人员主观因素造成例：对指示剂颜色辨别偏深或偏浅；滴定管读数不准例：对指示剂颜色辨别偏深或偏浅；滴定管读数不准 产生的原因：32 2 随机误差随机误差(偶然误差偶然误差)特点特点特点特点:（1 1）不恒定不恒定,无法校正；无法校正；（2 2）服从正态分布服从正态分布规律规律：大小相近的正误差和负误差出现的几率相等大小相近的正误差和负误差出现的几率相等;小误差出现的频率较高，而大误差出现的频率较低，小误差出现的频率较高，而大误差出现的频率较低，很大误差出现的几率近于零。很大误差出现的几率近于零。产生的原因产生的原因产生的原因产生的原因:（1 1）偶然因素偶然因素(室温，气压的微小变化室温，气压的微小变化)；（2 2）个人辩别能力个人辩别能力(滴定管读数滴定管读数)注意注意注意注意：过失误差属于不应有的过失。过失误差属于不应有的过失。过失误差属于不应有的过失。过失误差属于不应有的过失。二、误差的减免二、误差的减免（一）（一）（一）（一）系统误差的减免系统误差的减免系统误差的减免系统误差的减免2 随机误差(偶然误差)4 1.1.1.1.方法误差方法误差方法误差方法误差采用标准方法作对照试验 2.2.2.2.仪器误差仪器误差仪器误差仪器误差校准仪器 3.3.3.3.试剂误差试剂误差试剂误差试剂误差作空白试验（二）（二）（二）（二）随机误差的减免随机误差的减免随机误差的减免随机误差的减免 增加平行 测定的次数，取其平均值,可以减少随 机误差。正态分布曲线正态分布曲线 1.方法误差采用标准方法作对照试验 正态分布曲线5第一节第一节 地球化学数据的误差分析地球化学数据的误差分析三、误差的来源三、误差的来源 1 1取样产生误差取样产生误差 取样产生的误差是影响地球化学数据质量的最重要的也是最不取样产生的误差是影响地球化学数据质量的最重要的也是最不易被发现的误差。产生这种误差的途径可能有：易被发现的误差。产生这种误差的途径可能有：1 1）错误地将同一空间产出的不同单元作为同一单元。）错误地将同一空间产出的不同单元作为同一单元。2 2）没有正确确定不同空间的岩石是同一系列的同一单元还是不同）没有正确确定不同空间的岩石是同一系列的同一单元还是不同系列的不同单元。系列的不同单元。3 3）不能准确区分蚀变与原岩的界线。）不能准确区分蚀变与原岩的界线。4 4）判别风化岩石和新鲜岩石的界线不准确。）判别风化岩石和新鲜岩石的界线不准确。5 5）没有查明污染层的广度与深度。）没有查明污染层的广度与深度。6 6）土壤剖面分层不准确。）土壤剖面分层不准确。7 7）采集水样的时间或季节不当。）采集水样的时间或季节不当。8 8）对矿化与非矿石岩石的区分不当。）对矿化与非矿石岩石的区分不当。9 9）未能消除因断层等构造对原岩中元素活化迁移的影响。）未能消除因断层等构造对原岩中元素活化迁移的影响。1010）取样工具、器具及取样方法不当。）取样工具、器具及取样方法不当。第一节 地球化学数据的误差分析三、误差的来源6三、误差的来源三、误差的来源n2制备样品产生误差制备样品产生误差 1 1）没有将整个样品全部粉碎过筛，而仅粉碎了其易）没有将整个样品全部粉碎过筛，而仅粉碎了其易碎部分，难碎部分弃掉而造成粉碎的样品不代表整个碎部分，难碎部分弃掉而造成粉碎的样品不代表整个样品。样品。2 2）样品分选不彻底。）样品分选不彻底。3 3）筛网的原料不当，如不少样品不能用铜网或铁网）筛网的原料不当，如不少样品不能用铜网或铁网筛。筛。4 4）碎样工具和筛具清理不净，造成样品污染。）碎样工具和筛具清理不净，造成样品污染。5 5）样品缩分方法不当，未按缩分样品程序进行缩分）样品缩分方法不当，未按缩分样品程序进行缩分而随意分取一部分样品。而随意分取一部分样品。三、误差的来源2制备样品产生误差7n3分析测试产生误差分析测试产生误差 1 1）分析方法本身的原因，即分析的精度和灵敏度；）分析方法本身的原因，即分析的精度和灵敏度；2 2）分析测试人员因操作的原因而产生误差；）分析测试人员因操作的原因而产生误差；3 3）所用器具清洗不净产生误差；）所用器具清洗不净产生误差；4 4）仪器设备的精密度和分析结果的再现性达不到有）仪器设备的精密度和分析结果的再现性达不到有 关要求而产生误差；关要求而产生误差；5 5）化学试剂达不到分析要求或试剂选取不当产生误）化学试剂达不到分析要求或试剂选取不当产生误差；差；6 6）标准达不到要求产生误差；）标准达不到要求产生误差；7 7）分析测试环境产生误差。）分析测试环境产生误差。三、误差的来源三、误差的来源3分析测试产生误差三、误差的来源8n四、分析质量的控制四、分析质量的控制 1 1 内检内检 2 2 外检外检 3 3 双样品分析双样品分析四、分析质量的控制9第二节第二节 地球化学数据的预处理地球化学数据的预处理n一、基本概念一、基本概念 1总体、个体、样本 总体总体是指我们要研究的对象的全体；个体个体（或叫样品）是指总体中的一个单位；样本样本（或叫子样）是指总体中取出的一部分个体。样本所包含的个体数目。叫做样本的大小（或 叫样本容量样本容量）。第二节 地球化学数据的预处理一、基本概念10一、基本概念一、基本概念n2均值（1）算术平均值（2）加权平均值（3）几何平均值一、基本概念2均值11一、基本概念一、基本概念n3极值与方差 （1）数据的极小值与极大值 极小值 A min x1，x2，xn 极大值 B maxx1，x2，xn （2）方差、标准差n4变异系数Cv100%很大起伏。一、基本概念3极值与方差CvF表，说明两组数据的精密度存在显著性差异 若F计F表，说明两组数据的精密度无显著性差异，再用t检验法检验两组数据的准确度有无显著性差异。1 1 1 1）检验法检验法检验法检验法（1）计算两个样本的方差计算两个样本的方差S S 2 2（2）计算计算值值：3 显著性检验显著性检验 （3）查表（表），比较：1）检验法3 显著性检验24由于格鲁布斯Grubbs检验法使用了所有数据的平均值和标准偏差课件252 2 2 2）t t 检验法检验法检验法检验法 分析方法准确度的检验分析方法准确度的检验系统误差的判断系统误差的判断(1)(1)平均值与标准值平均值与标准值()的比较的比较 a.a.计算计算t 值值 b.由要求的置信度和测定次数,查表得到:t t表 c.比较t t计与t t表,若t t计 t t表,表示有显著性差异,存在系统误差,被检验方法需要改进。若t t计 t t表表,表示有显著性差异 t t计计 t t表表,表示无显著性差异(2)两组数据的平均值比较（同一试样）计算值：查表27数据的检验解决两类问题:1.1.可疑数据的取舍可疑数据的取舍过失误差的判断过失误差的判断可疑值检验可疑值检验：用数理统计方法检验测定数据是否存在应剔除的值 方法方法：Q检验法和格鲁布斯检验法 结论：结论：确定某个数据是否可用2.2.2.2.分析方法的准确性分析方法的准确性分析方法的准确性分析方法的准确性系统误差的判断(对照试验是检查分析过程中有无系统误差的最有效方法)显著性检验显著性检验显著性检验显著性检验：用数理统计方法检验被处理的数据是否存在统计上的显著性差异方法方法方法方法：t t 检验法和F F 检验法结论：结论：结论：结论：确定某种方法是否可用数据的检验解决两类问题:可疑值检验：用数理统计方法检验测定数284 异常值的检验异常值的检验n异常值是指分析数据中特别高或特别低的含量值。n特异值产生的原因：n（1）异常值：与特殊地质背景、土壤类型或用地类型相关，往往有一定的分布规律；n（2）分析数据错误：常呈单点出现，特高或特低。4 异常值的检验异常值是指分析数据中特别高或特别低的含量值。294 异常值的检验异常值的检验n可以使用Excel的数据分析功能绘制直方图、累计频率图等图件，由直方图可以较清晰的分辨出母体的多重性，以及离群值的分布状况。直方图的绘制与简要读图直方图的绘制与简要读图4 异常值的检验可以使用Excel的数据分析功能绘制直方图、305 质量分析质量分析n项目收到成果资料后，应将重复样品的分析成果筛选出来，并按样品类型进行排列，即第一次采样第一次分析、第一次采样第二次分析，第二次采样第一次分析、第二次采样第二次分析，分别计算分析的偏差。计算公式为：1）重复样品分析成果处理）重复样品分析成果处理5 质量分析项目收到成果资料后，应将重复样品的分析成果筛选出315 质量分析质量分析n设因素A具有n个不同的水平，如有n个不同的采样点A1，A2An。因素B具有m个不同的水平，如每个采样点重复采样m次或每个样重复分析m次，或有m个单位各分析一次等，即B1，B2Bm。n两因素方差分析的数学模型为：nXij=+i+j+ij （i=1,2,n），（j=1,2,m）n式中Xij第i个采样点第j次分析结果（或第j次采样）；某元素含量的总平均值（即数学期望值）；i第i个采样点的真值与实测平均值之差；j第j次分析（或第j次采样）间的系统误差；ij第i个采样点第j次分析（或第j次重复采样）的偶然误差。2 2）两因素方差分析）两因素方差分析5 质量分析设因素A具有n个不同的水平，如有n个不同的采样点32（1）首先计算各离差平方和）首先计算各离差平方和nSS总 称为总离差平方和，其自由度f总=nm1nSSA称为A因素离差平方和（如采样点间的离差平方和，它反映的是元素的自然变化，即与重复采样和样品多次分析无关，其自由度fA=n1nSSB称为B因素离差平方和（如多次重复采样的离差平方和或多次重复分析的离差平方和或多个分析单位间离差平方和等，它与采样的具体地点无关）。其自由度 fB=m1nSSE称为剩余离差平方和（常常反映偶然误差平方和或分析及偶然误差离差平方和，即除A和B两个因素之外的因素）。其自由度 fE=f总fAfB=（n1）（m1）计算步骤计算步骤（1）首先计算各离差平方和SS总 称为总离差平方和，其自由33（2）计算样本方差）计算样本方差（2）计算样本方差 34（3）作）作F检验检验（1）作统计量F1，F2（3）作F检验（1）作统计量F1，F235由于格鲁布斯Grubbs检验法使用了所有数据的平均值和标准偏差课件36由于格鲁布斯Grubbs检验法使用了所有数据的平均值和标准偏差课件375 质量分析质量分析3 3）三因素方差分析）三因素方差分析5 质量分析3）三因素方差分析38计计算算步步骤骤计算步骤39由于格鲁布斯Grubbs检验法使用了所有数据的平均值和标准偏差课件40由于格鲁布斯Grubbs检验法使用了所有数据的平均值和标准偏差课件416 变量的均匀化变换变量的均匀化变换n（1 1）均值计量变换）均值计量变换 均值计量变换即用绝对量的平均值作为均匀化均值计量变换即用绝对量的平均值作为均匀化因子，去除该变量每个样品的原始观测数据。因子，去除该变量每个样品的原始观测数据。n（2 2）极差变换）极差变换（又称正规化或规格化变换）（又称正规化或规格化变换）极差变换就是将原始数据减去该变量的最小值极差变换就是将原始数据减去该变量的最小值然后，以极差作为均匀化因子去除然后，以极差作为均匀化因子去除 。n（3 3）标准化变换）标准化变换 标准化变换是将每一变量的原始观测值减去平标准化变换是将每一变量的原始观测值减去平均值，并以该变量的标准差作为均匀化因子去均值，并以该变量的标准差作为均匀化因子去除。除。6 变量的均匀化变换（1）均值计量变换42第三节第三节 背景值和异常下限值的确定背景值和异常下限值的确定n一、图解法一、图解法1剖面图解法第三节 背景值和异常下限值的确定一、图解法43第三节第三节 背景值和异常下限值的确定背景值和异常下限值的确定一一 图解法图解法 2 直方图法Mo的横坐标值，即为所求的背景值c。（或背景值的对数值）。由频率（或频数）的极大值的0.6倍处，作一平行横坐标的直线，与曲线一侧相交，其横坐标长度即为均方差Sx。由Mo向右量取23倍的Sx，该处所指示的含量（或其对数值）即为异常下限值（或其对数值）。第三节 背景值和异常下限值的确定一 图解法Mo的横坐标值，44第三节第三节 背景值和异常下限值的确定背景值和异常下限值的确定一一 图解法图解法3概率格纸图解法概率格纸图解法也是建立在元素在地质体中呈正态分布（或对数正态分布）的基础上的。直线的斜率即为均方差 累计频率50处与频率曲线交点的横坐标（中位数Me）即为背景值C0（或背景值的对数值lgC0）。累计频率84.1或15.9与频率曲线的交点的横坐标与中位数的差值即为Sx（或其对数值），以累计频率97.7处在横坐标上的对应值为异常下限（或异常下限的对数值）。第三节 背景值和异常下限值的确定一 图解法概率格纸图解法也45第三节第三节 背景值和异常下限值的确定背景值和异常下限值的确定n二、计算法二、计算法 这种方法确定背景值和异常下限值也是建立在地质体中元素呈正态分布（或对数正态分布）的基础上。它是通过计算来确定背景值和异常下限值，以均值或众数值来代替背景值，以高出背景值的23倍为异常下限值。第三节 背景值和异常下限值的确定二、计算法46有效数字及其运算规则有效数字及其运算规则一、有效数字一、有效数字一、有效数字一、有效数字:指实际上能测量到的数字指实际上能测量到的数字。有效数字有效数字有效数字有效数字 =各位确定数字各位确定数字各位确定数字各位确定数字 +最后一位可疑数字最后一位可疑数字最后一位可疑数字最后一位可疑数字。1 1 1 1实验过程中常遇到两类数字实验过程中常遇到两类数字实验过程中常遇到两类数字实验过程中常遇到两类数字：（1 1）表示数目表示数目(非测量值非测量值):):如测定次数如测定次数；倍数倍数；系数系数；分分数数（2 2）测量值或计算值测量值或计算值。数据的位数与测定的准确度有关数据的位数与测定的准确度有关。记录的数字不仅表示数量的大小记录的数字不仅表示数量的大小，还要正确地反映还要正确地反映测量的精确程度测量的精确程度。结果 绝对误差 相对误差 有效数字位数 0.32400 0.00001 0.002%5 0.3240 0.0001 0.002%4 0.324 0.001 0.2%3有效数字及其运算规则一、有效数字:指实际上能测量到的数字。47 2 2数字零在数据中具有双重作用：数字零在数据中具有双重作用：（1 1）若作为普通数定使用，是有效数字）若作为普通数定使用，是有效数字 如如 0.3180 4 0.3180 4位有效数字位有效数字 3.180 3.180 10 10-1-1 （2 2）若只起定位作用，不是有效数字。）若只起定位作用，不是有效数字。如如 0.0318 3 0.0318 3位有效数字位有效数字 3.18 3.18 10 10-2-2 3 3改变单位不改变有效数字的位数：改变单位不改变有效数字的位数：如如 19.02mL 19.02mL为为19.0219.02 10 10-3-3 L L 2数字零在数据中具有双重作用：48有效数字的运算规则1.加减运算：几个数据相加或相减时，它们的和或差的有效数字的保留，应依小数点后位数最少的数据为根据，即取决于绝对误差最大的那个数据。例例例例：0.0122 绝对误差绝对误差绝对误差绝对误差：0.0001 25.64 0.01 1.051 0.001 25.7032 0.0121+25.64+1.057=25.70有效数字的运算规则492.乘除运算：几个数据的乘除运算中，所得结果的有效数字的位数取决于有效数字位数最少的那个数，即相对误差最大的那个数。例例例例：（：（0.0325 0.0325 5.103 5.103）/139.8=0.00119/139.8=0.00119 相对误差相对误差相对误差相对误差：0.0325 0.0001/0.0325 0.0325 0.0001/0.0325 100%=0.3%100%=0.3%5.103 0.001/5.103 5.103 0.001/5.103 100%=0.02%100%=0.02%9.8 0.1/139.8 9.8 0.1/139.8 100%=0.07%100%=0.07%几个数据的乘除运算中，所得结果的有效数字的位数取决于有503整化原则：（在取舍有效数字位数时，应注意以下几点）（1）在计算中，经常会遇到一些分数、整数、倍数等，这些数可视为足够有效。（2）若某一数据第一位有效数字等于或大于8，则有效数字的位数可多算一位。如：9.98，按4位算。（3）在计算结果中，可根据四舍五入原则（最好采用“四舍六入五留双”原则）进行整化。（4）有的浓度计算中，一般保留二位或三位有效数字。pH值的小数部分才为有效数字，一般保留一位或二位有效数字。例如,H+=5.210-3 molL-1，则pH=2.28（5）表示误差时，取一位有效数字已足够，最多取二位。3整化原则：（在取舍有效数字位数时，应注意以下几点）51