统计建模与R软件（薛毅）第八章答案.docVIP

下载本文档

88
0
约9.53千字
约 9页
2018-01-09 发布于重庆
举报
版权申诉

统计建模与R软件（薛毅）第八章答案.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

统计建模与R软件（薛毅）第八章答案

第八章 8.1 ####按数据框形式输入数据 classX1-data.frame( x1=c(1.9,-6.9,5.2,5.0,7.3,6.8,0.9,-12.5,1.5,3.8), x2=c(3.2,10.4,2.0,2.5,0.0,12.7,-15.4,-2.5,1.3,6.8) ) classX2-data.frame( x1=c(0.2,-0.1,0.4,2.7,2.1,-4.6,-1.7,-2.6,2.6,-2.8), x2=c(0.2,7.5,14.6,8.3,0.8,4.3,10.9,13.1,12.8,10.0) ) testx- c(8.1,2.0) ####调用已编好的程序source(discriminiant.distance.R) ####方差相等的情形 discriminiant.distance(classX1, classX2,testx, var.equal=TRUE) 1 blong 1 ####方差不等的情形 discriminiant.distance(classX1, classX2,testx) 1 blong 1 ####Bayes判别 source(discriminiant.bayes.R) discriminiant.bayes(class1,class2,rate=1,testx,var.equal=TRUE) 1 blong 1 discriminiant.bayes(classx1,classx2,rate=1,testx) 1 blong 1 source(discriminiant.fisher.R) discriminiant.fisher(class1, class2,testx) 1 blong 1 由综述上述结果可得结论：明天下雨 8.2 ####输入数据 -data.frame( x1=c(8.11,9.36,9.85,2.55,6.01,9.64,4.11,8.90,7.71,7.51,8.06,6.80, 8.68, 5.67, 8.10, 3.71, 5.37, 9.89,5.22, 4.71, 4.71, 3.36, 8.27), x2=c(261.01,185.39,249.58,137.13,231.34,231.38,260.25,259.91,273.84,303.59,231.03,308.90, 258.69, 355.54, 476.69, 316.32, 274.57, 409.42,330.34, 331.47,352.50, 347.31, 189.56), x3=c(13.23,9.02,15.61,9.21,14.27,13.03,14.72,14.16,16.01,19.14,14.41, 15.11,14.02,15.13,7.38,17.12,16.75,19.47,18.19,21.26,20.79,17.90,12.74), x4=c(7.36,5.99,6.11,4.35,8.79,8.53,10.02,9.79,8.79,8.53,6.15,8.49, 7.16,9.43,11.32,8.17,9.67,10.49,9.61,13.72,11.00,11.19,6.94) ) G - factor(rep(c(1,2,3),c(11,7,5))) p-c(11/23,7/23,5/23) ####p是先验概率：11/23,7/23,5/23 source(distinguish.distance.R) ####打开多类距离判别算法 distinguish.distance(heart,G,var.equal=F) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 blong 2 1 1 1 1 1 1 1 1 2 1 2 1 2 2 2 1 2 2 3 3 3 3 ####由结果可知，当协方差不一致的时候，第1,10,13,17,19号样本错判，回代的判别正确率为23-5=18,18/23=78.26% distinguish.distance(heart,G,var.equal=T) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 blong 1 1 1 1 1 1 3 1 1 2 1 2 1 2