R语言主成分分析的案例.docxVIP

R 语言主成分分析的案例 R 语言也介绍到案例篇了，也有不少同学反馈说还是不是特别明白一些基础的东西，希望能够有一些比较浅显的可以操作的入门。其实这些之前 SPSS 实战案例都不少，老实说一旦用上了开源工具就好像上瘾了，对于以前的SAS、clementine 之类的可视化工具没有一点感觉了。本质上还是觉得要装这个、装那个的比较麻烦，现在用R 或者 python 直接简单安装下，导入自己需要用到的包，活学活用一些命令函数就可以了。以后平台上集成R、python 的开发是趋势，包括现在 BAT 公司内部已经实现了。 今天就贴个盐泉水化学分析资料的主成分分析和因子分析通过R 语言数据挖掘的小李 子： 有条件的同学最好自己安装下R，操作一遍。 今有 20 个盐泉，盐泉的水化学特征系数值见下表.试对盐泉的水化学分析资料作主成分分析和因子分析.（数据可以自己模拟一份） 其中x1:矿化度(g/L) 导入数据保存在对象 saltwell 导入数据保存在对象 saltwell 中 saltwell-read.table(c:/saltwell.txt,header=T) saltwell x2:Br?103/Cl x3:K?103/Σ 盐 x4:K?103/Cl x5:Na/K x6:Mg?102/Cl x7:εNa/εCl. 1. 1.数据准备 2. 2.数据分析 标准误、方差贡献率和累积贡献率 arrests.pr- prcomp(saltwell, scale = TRUE) summary(arrests.pr，loadings=TRUE) 每个变量的标准误和变换矩阵 prcomp(saltwell, scale = TRUE) 查看对象 arests.pr 中的内容 str(arrests.pr) 利用主成分的标准误计算出主成分的累积方差比例 cumsum(arrests.pr$sdev^2)/7 [1]0.6067060 0.7850968 0.9165341 0.9790524 0.9954128 0.9999024 1.0000000 各个化学成分占主成分的得分 arrests.pr$x f.C I PC2 m FO` 阳25 ”, I-I 心?”“°心－L l 5 3 钝？5 ?`8 l 297比·O1 =(I且2 l S华妇斗-O.09 55938j [2 ) ? i .(iS9丑 0巧 心．9 J 0分 89 I. I6 6840沁 0 . 仇S7牛9Ol 52 O.OSI 叩 1 2 [3,l 4 川983世5d ! ·,O l ] ,l $617 ..2 3 67611 四 0 －0.68730“ 4 口 19 765, l ［斗J - I..38 5的 03中．I．忒9砒03 3., 4 9?，如心5 -I)正6?15 185 工l_ i4炽64? [5,I4 ，名S中 立 5斗书 ，3邸 究 i, I.6?i91?6-Ol, 4 0 . 5? i12 758 O.i 31 90 邸 l [6.J ,-tl1, i;l6j ,Sil35 -G.!)J. 茂 59 2 -3 . 030 6 幻 c--0 l ..0.OJ,4l24 l!n 02 0Cl斗gs3 仅14 ．勹5 迁．D64七 2.，士拓4r ? I _2 J.:J.314c+oo 1_1心化忧“； ° 0 .，g? ? C!Si 固 I 心，心吕63 ,}65 ] 心 ，71 !i印;22 3.81 l69 沁气Ol 庄，785 1 ! O49 心．＂叩 733 2 [9J l {11_2 200 B9-4 - I .0 39122_5, .7S774S 91 0.“0159ll 0 . ］2 089601 [ IOJ l 化2914？9 打 ·0j.9 i, 51 936 －624148c4 1 心斗 3:S7593 4 . E” :S2加？ [ 11, I 仇2 伲 0 15的 ．庄？l 3立嘈6 I _0叩 511c+00 O_73 的 5 n n ...10_2” i43351 [ 12,I 0-.7 75J152 l - OJ. 78氐2 6?O -3 _08 9？? 8巳 O) 03 n 6247?？－0. l 万 05749 [ l 3 , I 仆1．l16 68? r - LOO？ 斗 9 －I2 4 8生5？点 0 O另6 09位 56 -0”.， “”“ [ I勺，I l.7706 li61 O．知 957.2 -6i o649-I8c4 1－仆＿128 93651 心，位14 钻 65 [ 15. J l.,