资源描述：
《应用统计方法电子课件工科研究生应用数理统计作业佘源琦.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、CHINAUNIVERSITYOFPETROLEUM成绩：2009—2010学年第二学期《应用统计方法》期末课程论文题目逐步回归分析和动态聚类分析在白石泉Ni硫化物矿床研究中的应用专业：矿物学、岩石学、矿床学姓名：佘源琦学号：S090010030班级：一班完成日期：2010年5月4日摘要新疆东天山的白石泉Ni硫化物矿床一直受到众多构造学家、岩石学家和矿床学家的关注,但是数学地质在此地的研究尚属空白。文中在前人的基础查基础上,进行了系统岩石地球化学分析,应用多元线性回归和逐步回归分析相结合的方法,借助SAS数学统计软件系统的研究Ni变量与造岩元素变量的关系,建
2、立Ni变量与其他变量的回归方程,验证与解释白石泉Ni硫化物矿床的成因、岩浆来源等问题,对指导找矿勘探有重要作用,对应用数学地质方法有指引作用。关键词:白石泉;Ni硫化物矿床;;逐步回归分析;动态聚类；SAS新疆东天山是中国重要的Cu-Au-Ni-Fe成矿区带【1.8】,20世纪80年代以来,特别是“十五”期间,国家和相关部门在东天山地区实施一系列重点找矿与基础地质研究,使该区的找矿实践取得突破性进展。已发现的矿床有黄山东、香山、黄山、葫芦、图拉尔根、马蹄和土墩基性-超基性岩体与Cu-Ni矿床。新疆东天山的东段哈密地区因发育众多的晚古生代基性-超基性岩及与其伴
3、生的Cu-Ni矿而受到众多构造学家、岩石学家和矿床学家的关注[4,6-7]。但是数学地质在此地的研究尚属空白,本文在前人研究的基础上,进行了系统岩石地球化学分析,应用回归分析和动态聚类两种研究Ni变量与其他变量的关系,建立Ni变量与其他变量的回归方程,旨在预测铜镍矿化的强度与变化趋势,揭示矿化与蚀变的关系,以指导找矿与勘探工作。1矿区地质特征概况白石泉铜镍矿床位于哈密市东南170km处,位于中天山地块北缘,阿齐克库都克—沙泉子断裂带南侧(图1)。其中阿齐克库都克—沙泉子断裂是中天山地块与觉罗塔格构造带的分界断裂,伴生有硅化等蚀变的逆冲断裂破碎带。矿区位于阿齐
4、克库都克—沙泉子深大断裂南侧1~3km处。北部为元古宇卡瓦布拉克群(Ptk)的一套绢云母石英片岩、绿泥石石英片岩、石英岩、大理岩;南部为元古宇星星峡组中亚组(Ptx)的一套(绢、黑)云母斜长石英片岩、斜长角闪片岩、麻粒岩,以F2断层与北侧卡瓦布拉克群接触,其次在低洼地带及现代冲沟中有第四系冲洪积物覆盖。图1新疆东天山大地构造分区与矿床分布图矿区的侵入岩为华力西期产物,以中性(偏基性)岩为主,在其边部及中部有基性-超基性岩分布,并构成白石泉杂岩体;脉岩从基性到酸性均有分布,但一般规模较小(图2)图2白石泉铜镍矿区地质图2 逐步回归方法的基本思想在多元线性回归分
5、析中,我们知道回归方程中所包含的自变量越多,那么回归平方和U就越大,则剩余平方和Q就越小,一般情况下剩余标准差(S剩余)也随之减少,回归方程效果越好,而精度也越高。在“最优”回归方程中总希望包括尽可能多的自变量,特别是对因变量y有显著影响的自变量不能遗漏。但回归方程所包括的自变量太多,也带来不利的一面,首先,若要求自变量多,则在预测时必须测定许多量,并且计算也不方便;其次,如果在回归方程中包括有对因变量y不起作用或作用极小的自变量,那么剩余平方和Q也不会由于自变量的增加而减少,相反由于Q的自由度的减少,反而使剩余标准差增大,这就影响回归方程的精度;第三,由于
6、存在着对因变量y影响不显著的自变量,以致影响回归方程的稳定性,使预测效果下降。因而,在“最优”回归方程中,又希望不包括对因变量y影响不显著的自变量[9]。综上所述,所谓最优回归方程,就是在回归方程中包括所有对因变量y有显著影响的自变量,而不包括对因变量y影响不显著的自变量的回归方程。最优回归方程的建立,是采取将自变量逐个引入的方法。引入自变量的条件是:该自变量的偏回归平方和经检验是所有自变量中最显著的。同时,每引入一个新变量后,要求对已引入的自变量逐个进行检验,将偏回归平方和变得不显著的自变量及时剔除。由于每步都作检验,因而保证了最后所建立的回归方程中所有自
7、变量都是显著的。上述这种建立最优回归方程的理论和方法,称为逐步回归分析[10-12]。3 白石泉铜镍硫化物矿床的回归分析采用组合变量多元线性回归分析和逐步回归分析,综合研究岩体中Ni变量与其他变量关系(表1)[13-14],然后利用SAS数学统计软件计算,得到相关系数矩阵及检验结果、结合地质实际结果,检验和分析白石泉铜镍硫化物矿床的岩浆来源、矿化因素等。表1岩石化学分析结果SiO2TiO2Al2O3Fe2O3FeOMnOMgOCaONa2OK2OP2O5NiCu40.350.437.115.378.060.1828.674.090.480.050.07658
8、.1095.2040.420.596.386.019
