应用多元统计分析作业.docxVIP

PAGE  PAGE 14 多元统计分析 实 验 报 告 实验课程名称 多元统计分析 实验项目名称 多元统计理论的计算机实现 年 级 2013 专 业 应用统计学 学生姓名 侯杰 成 绩 理 学 院 实验时间： 2015 年 05 月07 日 学生所在学院：理学院 专业：应用统计学 班级：9131137001 姓 名侯杰 学 号913113700124实验组 实验时间 16.05.07指导教师李建军实验项目名称多元统计理论的计算机实现 实验目的及要求： 目的：熟悉R（或SPSS）软件，掌握多元统计分析中多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容。 要求：程序要有注释，尽量体现多元统计分析多元正态分布均值向量和协差阵的检验，判别方法，聚类分析，主成分分析，因子分析，相应分析内容内容的基本原理。  实验硬件及软件平台： 计算机、 R、 网络 实验内容（包括实验具体内容、算法分析、源代码等等）： 指导教师意见： 签名： 年 月 日 代码及运行结果分析 均值检验 问题重述：某医生观察了16名正常人的24小时动态心电图，分析出早晨3小时各小时的低频心电频谱值（LF）、高频心电频谱值（HF），数据见压缩包，试分析这两个指标的各次重复测定均值向量是否有显著差异。 代码如下： Tsq.test-function(data,alpha=0.05){ data-as.matrix(read.table(ch37.csv,header=TRUE,sep=,)) #读取数据 xdat-data[,2:4]; xbar-apply(xdat,2,mean); #计算LF指标的均值 ydat-data[,5:7]; ybar-apply(ydat,2,mean); #计算HF指标数据 xcov-cov(xdat); #计算LF样本协差阵 ycov-cov(ydat); #计算HF样本协差阵 sinv-solve(xcov+ycov);#求逆矩阵 Tsq-(16+16-2)*t(sqrt(16*16/(16+16)*(xbar-ybar)))%*%sinv%*%sqrt(16*16/(16+16)*(xbar-ybar)); #计算T统计量 Fstat-((16+16-2)-3+1)/((16+16-2)*3)*Tsq; #计算F统计量 pvalue-as.numeric(1-pf(Fstat,3,16+16-3-1)); cat(p值=,pvalue,\n); if(pvalue0.05) #结果输出 cat(均值向量不存在差异) else cat(均值向量存在差异); } 运行结果及分析： 通过运行程序，我们可以得到如下结果： Tsq.test() p值= 1.632028e-14 均值向量存在差异 即LF与HF这两个指标的各次重复测定均值向量存在显著差异。 判别分析 问题重述：银行的贷款部门需要判别每个客户的信用好坏（是否未履行还贷责任），以决定是否给予贷款。可以根据贷款申请人的年龄（）、受教育程度（）、现在所从事工作的年数（）、未变更住址的年数（）、收入（）、负债收入比例（）、信用卡债务（）、其它债务（）等来判断其信用情况。数据见压缩包。⑴根据样本资料分别用距离判别法、Bayes判别法和Fisher判别法建立判别函数和判别规则。⑵某客户的如上情况资料为（53，1，9，18，50，11.20，2.02，3.58），对其进行信用好坏的判别。 代码如下： #距离判别法 discrim.dist-function(x)