矿统实习讲义实习2回归分析法.doc

班级 学号 姓名 成绩 实习二、用多元线性回归分析法进行矿床统计预测 目的 通过实习，学会使用多元线性回归分析进行矿床统计预测，加深对该方法原理的理解。要求 （1）根据所提供资料，自己动手完成预测计算的各环节，用控制单元数据建立回归模型，对所提供的各未知单元，评价它们含有小型及以上矿床的可能性；按时提交实习报告。（2）复习课程“多元线性回归分析”有关内容。资料 研究区是湖北省某地区一个铁矿成矿带。为在该区进行矿床统计预测，已将研究区划分为500m×500m基本单元408个，并提取了地质变量。为应用回归分析法，选取一种矿床值（储量乘以某个系数得到的矿床价值）为因变量和多个地质因素、找矿标志为自变量()。本次实习为简便只使用其中两个自变量：为单元磁异常值，为单元中闪长岩体出露面积比。表2-1（控制单元数据及回归估值表）最左4列列出了实习所用数据。 表2-1控制单元数据及回归估值表 单元 序 单元 序 6 0.778 0.958 0.01 58 2.553 0.001 8.50 12 2.204 1.456 24.20 59 0.829 2.950 11.0 14 0.477 0.001 0.01 60 0.477 1.491 11.0 20 4.007 1.861 36.25 61 0.477 0.612 0.01 21 0.477 2.301 0.75 63 0.954 0.001 0.75 27 5.385 2.045 6.25 64 0.954 0.001 1.50 28 5.636 3.057 35.21 65 4.625 0.001 5.75 33 0.477 1.513 7.00 66 0.477 0.001 1.50 35 5.763 2.000 22.70 67 0.477 0.001 0.01 36 3.992 1.468 7.25 68 0.954 0.001 0.75 37 0.477 0.001 0.01 76 3.437 0.001 0.01 43 0.477 0.783 0.57 77 0.477 0.001 0.01 45 5.467 1.301 75.00 83 0.778 0.001 0.01 46 5.035 2.672 24.2 87 0.477 0.001 1.00 47 3.464 3.250 75.2 89 2.813 0.001 0.15 56 0.477 1.096 7.00 100 0.477 0.001 0.01 注：y值小于1为矿点，1-3.5为小型矿床，大于3.5为大中型矿床。 65.829 30.829 1477.273 755.116 方法步骤 第一步：分析研究区内铁矿特征及控矿地质条件和找矿标志，划分基本单元，提取地质变量，地质变量赋值并做适当变换（使因变量尽量服从正态分布、使因变量与自变量之间有线性关系）。这些工作已经完成（不必重新做），见表2-1的最左4列。 第二步：建立回归方程。二元回归方程为 （1） 上式中各系数用最小二乘法确定。最小二乘法就是使偏差平方和 （2） 达到极小而求出。上式中为观测样品总数。为此令 （3） 将（1）和（2）代入（3），可得到关于的线性方程组，称正规方程组。解正规方程组可求出各系数，从而得到回归方程。虽然可以采用矩阵形式，但为利于深入细节，本次实习采用分步骤计算的方式。先计算，等值，填满表2-1的最下面一行。然后计算以下各值，结果填入空白处： _________，_________，_________， 66.328-30.8292/32=36.627_______, ___________________________________________, ___________________________________________, _________________________________, _________________________________, _________________________________, ___________________ ___________________ ___________________________