R语言判别分析实验报告.pdf

R R 语语言言判判别别分分析析实实验验报报告告集团标准化小组：VVOPPT-JOPP28-JPPTL98-LOPPNNR R语语言言判判别别分分析析实实验验报报告告班级：应数 1201学号：姓名：麦琼辉时间：2014 年 11 月 28 号1 1 实验目的及要求实验目的及要求1)了解判别分析的目的和意义；2)熟悉 R 语言中有关判别分析的算法基础。2 2 实验设备（环境）及要求实验设备（环境）及要求个人计算机一台，装有 R 语言以及 RStudio 并且带有 MASS 包。3 3 实验内容实验内容企业财务状况的判别分析4 4 实验主要步骤实验主要步骤1)数据管理：数据管理：实验对 21 个破产的企业收集它们在前两年的财务数据，对 25个财务良好的企业也收集同一时期的数据。数据涉及四个变量：CF_TD（现金/总债务）；NI_TA（净收入/总资产）；CA_CL（流动资产/流动债务）；CA_NS（流动资产/净销售额），一个分组变量：企业现状（1：非破产企业，2：破产企业）。2)调入数据：调入数据：对数据复制，然后在 RStudio 编辑器中执行如下命令。case5=read.table(clipboard,head=T)head(case5)3)FisherFisher 判别效果（等方差，线性判别判别效果（等方差，线性判别 ldalda）：）：采用 Bayes 方式，即先验概率为样本例数，相关的 RStudio 程序命令如下所示。library(MASS)ld=lda(G.,data=case5);ld#线性判别ZId=predict(ld)addmargins(table(case5$G,ZId$class)4)FisherFisher 判别效果（异方差，非线性判别判别效果（异方差，非线性判别-二次判别二次判别qdaqda）：）：再次采用Bayes 方式，相关的 RStudio 程序命令如下所示。library(MASS)qd=qda(G.,data=case5);qd#二次判别Zqd=predict(qd)addmargins(table(case5$G,Zqd$class)5 5 实验结果实验结果表 1 线性判别 lda 效果新分类原分类12 合计124125231821合计271946符合率91.30%由表 1 和表次 判 别-非 线比 lda（一次判表 2 二次判别 qda 效果新分类原分类12 合计124125221921合计262046符合率93.50%2 可知，qda（二性判别）的效果别）要好。6 6 实验小结实验小结通过本次实析 的 目 的 和 意验了解了判别分义，并熟悉 R 语言中有关判别分析的算法基础。