4实验指导四.doc

数学实验指导四实验名称：多元统计分析软件SAS的应用实验项目： 判别分析专业班级： 姓名：学号：指导教师：成绩：日期：一、实验目的熟悉sas软件中的discrim过程，学会利用该过程对未知类别的样品分类二、实验内容及要求盐泉含钾性判别某地区经勘探证明，A盆地是一个钾盐矿区，B盆地是一个钠盐（不含钾）矿区，其他盐盆地是否含钾盐有待作出判断。今从A和B两盆地各抽取5个盐泉样品；从其他盆地抽得3个盐泉样品。试对后3个待判盐泉进行含钾性判别。（用Fisher 判别法）盐泉类别序号X1X2X3X4类别号第一类：含钾盐泉（A盆地）1234513.8522.3128.8215.2928.792.794.674.633.544.907.8012.3116.187.5016.1249.6047.8062.1543.2058.1011111第二类：含钠盐泉（B盆地）6789102.183.8511.403.6612.101.060.800.002.420.001.224.063.502.145.6820.6047.100.0015.100.0022222待判盐泉1237.903.1912.402.403.205.104.301.434.4333.209.9024.60?要求：写出判别函数、训练样本误判率、待判样本判别结果。三、实验结果Sas程序中过程discrim基本句法及注释: Proc discrim (判别分析过程名) data=训练样本数据集 method=normal(使用正态总体方法) pooled=yes(合并方差) crossvalidate(输出刀切法判别结果) list(列出回代判别结果) crosslisterr(误判结果) testdata=待判样本数据集 testlist (列出待判样本判别结果) testout=(待判数据集判别结果) ;Class 分类变量名;Var 判别分析使用的指标变量名;Run;提示：更详细说明可查看帮助文件，学习更多功能，需要一定的统计理论基础。书上108页实例程序:data a;input g$ n$ country$ x1-x3;cards;11美国 7699537412日本 79.599535913瑞士 7899537214阿根廷72.195.9524215阿联酋73.877.7537026保加利亚 71.293425027古巴 75.394.9341228巴拉圭7091.2339029格鲁吉亚 72.8992300210南非 62.980.63799;proc print;data b;input n$ country$ x1-x3;cards;11中国68.379.3195012罗马尼亚69.996.9284013希腊77.693.8523314哥伦比亚69.390.35158;proc discrim data=a testdata=b method=normal pool=yes crosslisterr list testout=q;