实验判别分析、聚类分析.pptVIP

一、判别分析 判别分析的基本思想是根据一批分类明确的样本在若干指标上的观察值，建立一个关于指标的判别函数和判别准则，然后根据这个判别函数和判别准则对新的样本进行分类，并且根据回代判别的准确率评估它的实用性。 Fisher 判别法是以距离为判别准则来分类，仅适用于两类判别 Bayes 判别法是以概率为判别准则来分类，适用于多类判别。SAS软件采用Bayes 判别法 1. 判别分析过程 PROC DISCRIM; CLASS G; VAR X1-X5； PRIORS PROP； RUN; DISCRIM过程后选项： METHOD=NORMOL参数法(默认) 、NPAR非参数法 POOL=YES (默认) 、NO、TEST先检验，方差齐则合并 ANOVA、MANOVA 进行一元、多元方差分析检验 DISTANCE 输出类内及类间的平均距离 LISTERR 输出回代错判的个体及后验概率 CROSSLISTERR 输出回代结果中原分类与现分类的列联表 OUT= 运算结果输出到指定数据集 OUTSTAT= 中间统计量及二次方程系数输出到指定数据集 SIMPLE 可输出各自变量的简单统计描述 NOPRINT 禁止屏幕输出任何结果 PRIORS 语句指定先验概率，有3个选项： EQUAL 指定各类的先验概率相等，为默认值 PROR 各类的先验概率等于各类样本频率 ‘1’=P1 ‘2’=P2 ‘3’=P3… 指定各类的先验概率为P1、P2….Pi， (P1+P2+…+Pi=1) 组外回代 DATA TEST; INPUT X1 X2@@; CARDS; PROC DISCRIM DATA=A1 TESTDATA=TEST TESTLIST POSTERR; CLASS G; VAR X1 X2; RUN; 练习1： 通过已确诊的100名四类病人的血清生化指标测定值建立判别函数式 PA: serum protein PA AG: serum protein ?1-AG AT: serum protein ?1-AT HP: serum protein HP GROUP: 健康=1； 肝癌(APF+)=2； 肝癌(APF-)=3； 肝硬化=4 程序： d:\mv_sas\unit4\data\discrim1.sas 数据： d:\mv_sas\unit4\data\disc1.xls 练习2： 组外回代 调查50-59岁的女性冠心病人15名和正常人16名，试用判别分析法建立判别冠心病与正常人的判别函数式，并作组外病例回代。 X1 舒张压 X2 血浆胆固醇 G 病人=1；正常人=2 程序： d:\mv_sas\unit4\data\discrim2.sas 2. 逐步判别分析过程 仅适用于类内为多元正态分布，并具有相同方差协方差矩阵的数据（POOL=TEST时Pα）。 PROC STEPDISC; CLASS G; VAR X1-X5； RUN; STEPDISC过程后可加选项： METHOD=STEPWISE(默认) FORWARD、BACKWARD SLE=0.15 SLS=0.15 练习3： 调查心肌梗塞和健康人各30名，分别测得其6项血脂指标，试建立判别函数式 TC 总胆因醇 TG 甘油三酯 HDLC 高密度脂蛋白胆固醇 LDLC 低高密度脂蛋白胆固醇 apoA 载脂蛋白A apoB 载脂蛋白B GROUP 心肌梗塞=1；正常人=2 程序： d:\mv_sas\unit4\data\stepdisc.sas 数据： d:\mv_sas\unit4\data\stepdisc.xls 练习4： 通过三组人群的5项心电图指标，建立判别函数 X1、X2、X3、X4、X5为5项心电图指标 GROUP: 健康=1；硬化症=2；冠心病=3 设各类先验概率相等 程序： d:\mv_sas\unit4\data\discrim3.sas 数据： d:\mv_sas\unit4\data\disc3.x