- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
R语言主成份分析的案例
R语言也介绍到案例篇了,也有不少同学反馈说还是不是特殊明白一些基础的东西,希翼能够有一些比较浅显的可以操作的入门。其实这些之前SPSS实战案例都不少,老实说一旦用上了开源工具就好像上瘾了,对于以前的SAS、clementine之类的可视化工具没有一点感觉了。本质上还是觉得要装这个、装那个的比较麻烦,现在用R或者python直接简单安装下,导入自己需要用到的包,活学活用一些命令函数就可以了。以后平台上集成R、python的开辟是趋势,包括现在BAT公司内部已经实现了。
今天就贴个盐泉水化学分析资料的主成份分析和因子分析通过R语言数据挖掘的小李子:
有条件的同学最好自己安装下R,操作一遍。
今有20个盐泉,盐泉的水化学特征系数值见下表.试对盐泉的水化学分析资料作主成份分析和因子分析.(数据可以自己摹拟一份)
其中x1:矿化度(g/L)
2.数据分析
x2:Br?103/Cl
x3:K?103/Σ盐
x4:K?103/Cl
x5:Na/K
x6:Mg?102/Cl
x7:εNa/εCl.
1.数据准备
导入数据保存在对象saltwell中
saltwell-read.table(c:/saltwell.txt,header=T)
saltwell
1标准误、方差贡献率和积累贡献率
arrests.pr-prcomp(saltwell,scale=TRUE)
summary(arrests.pr,loadings=TRUE)
2每一个变量的标准误和变换矩阵
prcomp(saltwell,scale=TRUE)
3查看对象arests.pr中的内容
str(arrests.pr)
4利用主成份的标准误计算出主成份的积累方差比例
cumsum(arrests.pr$sdev^2)/7
[1]0.60670600.78509680.91653410.97905240.99541280.9999024
1.0000000
5各个化学成份占主成份的得分
arrests.pr$x
6数据分析结果图形表示
screeplot(arrests.pr,main=saltwell)
biplot(arrests.pr)
按第一主成份排序的结果:
data.frame(sort(arrests.pr$x[,1]))
主因子分析
计算数据的相关系数矩阵
saltwell.cor-cor(saltwell)
saltwell.cor
计算特征值和特征向量及因子的贡献率和积累贡献率
saltwell.eigen-eigen(saltwell.cor)
saltwell.eigen
根据主成份分析结果确定公共因子个数.
saltwell.pr-princomp(saltwell,cor=T)
summary(saltwell.pr)
均值
saltwell.pr$center
标准误
saltwell.pr$scale
下面用特征值的平方根乘以相应的特征向量得到因子载荷矩阵.并且只显示前2个因子的结果:
t(sqrt(saltwell.eigen$values)*t(saltwell.eigen$vectors))[,1:2]
用R语言自带的函数factanal()进行分析
saltwell.fa-factanal(saltwell,factors=2)
print(saltwell.fa,cutoff=0.001)
下面用回归方法(regression)计算因子得分并作图,然后对样本进行分类.
saltwell.fa-factanal(saltwell,factors=2,scores=regression)
saltwell.fa$scores
您可能关注的文档
- UV光固化墨水项目工程项目投资估算.docx
- SPC管理作业指导书.docx
- SAP项目PS项目管理模块调研提纲.docx
- SAP QM 物料主数据QM视图维护.docx
- R007-贸易金融授信业务总协议.docx
- ppt课件制作过程.docx
- PCL6143运动控制器的原理及应用.docx
- origin-直接处理XPS数据说明.docx
- NETGEAR城中村宽带网络运营方案.pdf
- MODIS数据介绍、下载及处理.docx
- 新的一年工作展望.docx
- 医生年终个人工作的述职报告(3篇).docx
- 2023年消防设施操作员之消防设备中级技能考前冲刺练习题附答案详解.docx
- 2022-2023年环境影响评价工程师之环评技术导则与标准通关练习题包括详细解答.docx
- 2023年中级注册安全工程师之安全生产管理考前冲刺检测卷和答案.docx
- 2023年中级银行从业资格之中级银行管理考前冲刺测试卷提供答案解析.docx
- 2023年公共营养师之二级营养师通关模拟考试试卷附带答案.docx
- 证券分析师之发布证券研究报告业务考前冲刺模拟题库.docx
- 2022-2023年二级建造师之二建建设工程法规及相关知识综合提升测试卷附答案.docx
- 2023年二级建造师之二建机电工程实务通关模拟考试试卷提供答案解析.docx
文档评论(0)