r语言主成分分析_从0开始入门R语言—主成分分析PCA

683c41dbd9afa909faf54eb92bcca822.png
#读取文件信息,此数据在文末
orgData<-read.csv("data.csv") 
head(orgData)
#方法一:作主层分首先要标准化
orgData<-orgData[,2:5]
apply(orgData,2,sd)
#方法二:采用相关系数矩阵
pr1<-princomp(orgData,cor=TRUE)
summary(pr1) #查看主层分分析结果

# Importance of components:
#                           Comp.1    Comp.2    Comp.3     Comp.4
# Standard deviation     1.5812709 0.9923751 0.7137466 0.07307394
# Proportion of Variance 0.6251044 0.2462021 0.1273585 0.00133495
# Cumulative Proportion  0.6251044 0.8713065 0.9986650 1.00000000

summary(pr1,loadings=TRUE) #查看因子载荷矩阵
# Importance of components:
#                           Comp.1    Comp.2    Comp.3     Comp.4
# Standard deviation     1.5812709 0.9923751 0.7137466 0.07307394
# Proportion of Variance 0.6251044 0.2462021 0.1273585 0.00133495
# Cumulative Proportion  0.6251044 0.8713065 0.9986650 1.00000000
# 
# Loadings:
#          Comp.1 Comp.2 Comp.3 Comp.4
# cnt_call -0.111  0.990              
# cnt_msg  -0.510 -0.127  0.810 -0.262
# cnt_wei  -0.579        -0.559 -0.593
# cnt_web  -0.627        -0.157  0.762

# 因子分析,上述主成分解释性差,因此,通过因子旋转增强主成分的解释能力
fc2<-principal(orgData,nfactor=3,rotate="varimax",score=TRUE)
fc2
# Principal Components Analysis
# Call: principal(r = orgData, nfactors = 3, rotate = "varimax", scores = TRUE)
# Standardized loadings (pattern matrix) based upon correlation matrix
#          RC1   RC2  RC3   h2     u2 com
# CNT_TBM 0.04  0.99 0.10 0.99 0.0092 1.0
# CNT_ATM 0.78 -0.05 0.19 0.65 0.3481 1.1
# CNT_POS 0.84  0.10 0.04 0.71 0.2871 1.0
# CNT_CSC 0.16  0.11 0.98 0.99 0.0113 1.1
# 
#                        RC1  RC2  RC3
# SS loadings           1.34 1.00 1.00
# Proportion Var        0.34 0.25 0.25
# Cumulative Var        0.34 0.59 0.84
# Proportion Explained  0.40 0.30 0.30
# Cumulative Proportion 0.40 0.70 1.00
# 
# Mean item complexity =  1.1
# Test of the hypothesis that 3 components are sufficient.
# 
# The root mean square of the residuals (RMSR) is  0.14 
#  with the empirical chi square  22611.62  with prob <  NA 
# 
# Fit based upon off diagonal values = 0.6

#变量聚类
#orgData<-orgData[,2:5]
names(orgData)
orgData<-orgData[,2:5]
#install.packages("ClustOfVar")
library(ClustOfVar)
tree<-hclustvar(orgData)
plot(tree)

#choice of the number of clusters
stability(tree,B=40)
part<-cutreevar(tree,3,matsim=T)
#print(part)
summary(part)
part$var
part$sim
head(part$scores)
test<-cbind(orgData[,3:4],part$scores)
cor(test)
#选择的最优变量:与自己组cluster越相关越好,与别的组cluster越不想关越好
cnt_wei_1_R2<-(1-0.98*0.98)/(1-0.51*0.51)#所以选择cnt_wei_1_R2
cnt_web_1_R2<-(1-0.98*0.98)/(1-0.739*0.739)

以下是数据~~~~~~~~~~

ID	cnt_call	cnt_msg	cnt_wei	cnt_web
1964627	46	90	36	31
3107769	53	2	0	2
3686296	28	24	5	8
3961002	9	2	0	4
4174839	145	2	0	1
5068087	186	4	3	1
5574319	34	34	31	16
5733390	61	3	0	2
6148334	87	14	14	7
6190226	103	2	0	2
6794044	22	2	0	1
8806327	58	2	0	1
9337857	56	30	45	22
9545675	106	2	0	2
12830089	54	5	0	2
13088420	28	9	30	12
13830644	110	54	28	22
14321378	89	50	31	22
17432437	13	9	1	2
18297289	28	40	21	14
18633203	3	2	0	0
19006275	90	30	10	9
19501849	57	2	0	1
22213441	75	34	46	23
23746234	15	13	7	7
24574017	144	33	90	37
31203185	84	8	1	2
33523674	79	16	7	6
38470850	144	2	0	2
40524929	12	10	1	2
42571470	14	2	0	2
44444171	60	2	0	1
44733230	10	5	0	3
46056236	59	5	12	4
49342221	163	19	4	5
49853771	56	42	11	13
51873745	45	2	0	1
54493626	40	2	0	1
56918870	92	34	8	10
57335986	41	5	0	2
59319033	19	11	2	3
64837074	7	2	0	0
66008620	111	19	41	18
67695886	80	2	0	0
70393919	48	52	22	18
71533065	110	14	6	6
71601391	50	7	1	1
76205684	89	2	0	1
76856670	50	20	4	5
80175682	33	2	0	0
80532151	45	34	40	21
81616872	14	2	0	0
82578429	68	2	0	3
83350197	30	42	17	14
83644879	328	2	0	1
85307522	9	23	5	9
87241863	58	88	55	35
87469631	135	21	5	7
89616802	190	16	53	21
90120796	33	2	23	8
92749578	44	10	30	12
93567222	36	4	14	5
95011593	79	2	6	2
96766333	13	2	0	0
98409940	124	9	1	4
98988350	14	2	0	0
100050907	57	2	0	1
101124040	63	5	0	2
104011739	43	37	10	9
105023634	39	2	0	0
105080681	13	4	0	1
105103356	23	14	15	10
107977170	110	32	8	8
108669317	7	29	14	10
115103248	65	2	0	0
118563063	101	25	30	17
119043769	144	2	0	1
119597006	75	2	0	0
120129156	46	23	13	10
121518058	70	2	0	1
122042785	133	19	105	39
125419152	4	6	0	4
128254222	10	2	0	0
128326127	8	6	0	3
129176318	120	2	0	2
130055921	12	2	0	2
130078026	81	2	8	4
136274207	7	11	2	5
136322860	174	21	5	8
137042222	286	2	0	0
138730564	116	30	19	10
140385122	24	20	4	7
140463210	25	7	1	4
142437323	38	29	29	16
143920827	45	18	12	8
147221477	16	6	0	2
148303754	40	7	8	4
151423050	49	34	9	9
152309015	17	9	48	18
153564854	152	24	5	8
154501820	259	21	5	6
155357211	38	2	0	1
156868289	69	6	53	19
157261209	165	17	3	4
157775478	7	4	0	0
162645333	47	25	6	7
162886701	56	15	3	5
163529785	85	16	4	5
164763253	79	2	0	4
165804830	10	3	0	0
168754193	20	5	0	2
168764830	61	37	57	26
168840229	137	24	75	31
169886881	24	30	94	39
170126716	176	40	11	13
170871624	30	16	36	14
171753874	39	2	35	13
175386308	7	11	2	6
176105780	16	2	0	2
176376649	60	115	159	77
177129638	9	12	4	5
177710649	105	10	2	3
178008696	36	15	4	5
178468702	6	2	0	1
181051663	28	15	4	4
182141494	13	2	0	0
185172771	14	28	22	12
187902501	52	2	0	1
189628669	98	37	31	18
190351908	119	2	0	1
191185562	98	2	0	1
191244383	74	45	12	13
192635000	11	5	0	2
192646586	15	2	0	1
192687201	10	10	2	2
193272768	24	25	18	10
193419148	62	26	50	23
193450876	18	2	0	1
193772712	107	3	0	3
198425150	15	11	2	4
198477840	81	2	0	0
198659389	140	2	0	1
199251899	147	44	12	12
202657878	17	2	0	0
203042801	47	2	12	5
204368269	40	74	19	20
204765109	47	2	0	2
206703691	62	12	2	3
206889690	49	2	0	1
211761132	50	2	0	0
215449998	4	5	0	1
217335228	121	11	13	9
218022716	70	55	19	17
218216440	34	43	25	15
218251830	30	12	2	3
220242626	63	7	1	1
221362115	7	10	2	3
224253857	72	31	36	19
224500170	41	44	46	24
225095226	12	2	0	0
226461814	51	27	18	12
227798415	58	20	5	7
229325525	28	2	0	2
229499324	20	2	0	1
229677653	118	24	7	6
231038265	110	24	54	23
232278654	14	23	5	8
232830575	46	11	2	2
233994349	11	2	0	0
234322379	93	14	7	4
234605761	38	7	2	3
235629625	47	27	55	23
235698973	5	4	0	2
238134152	33	20	7	7
239169161	77	11	19	8
240288690	31	2	0	1
240762769	36	17	23	12
242498032	113	74	60	34
243525808	47	52	58	29
243531708	135	7	73	25
244879228	148	27	17	13
246056082	210	2	0	1
250384566	166	50	77	36
251764291	83	34	17	14
251810754	66	39	27	18
252853220	35	29	7	8
254353144	36	14	8	4
257225967	150	2	0	2
261742482	67	2	0	0
264225572	5	2	0	1
264455199	30	73	19	20
270648342	57	33	14	11
271215695	58	43	41	22
271218445	56	31	48	22
273711983	46	2	0	2
275008090	143	40	148	58
276726625	167	8	57	20
277283908	73	28	102	41
282983900	152	23	30	15
284897498	133	25	6	10
286319548	85	3	0	1
289503865	171	14	10	8
292721037	13	3	0	1
294943851	7	2	0	2
295884239	123	28	11	11
295940261	70	25	6	7
297994041	49	20	4	5
304076405	729	10	2	2
308167734	7	2	0	1
310001595	52	2	0	1
311092829	54	8	7	5
313475031	81	12	12	7
314932328	12	10	2	3
315970539	15	18	4	6
317030510	25	34	8	9
322758859	103	4	44	16
324687072	101	5	0	3
325877779	17	51	13	14
326634395	66	15	20	11
327575126	16	2	0	0
330265163	21	35	122	47
330265949	70	46	11	12
333203708	68	18	51	20
336836057	144	30	7	8
340569949	34	72	34	27
343108714	48	6	61	23
344319344	35	44	30	19
345444952	100	2	0	1
345720692	213	2	0	2
346020003	68	2	0	2
347613964	32	9	8	4
348882956	32	11	21	12
349213426	11	2	0	1
350548517	97	31	20	14
350665021	20	6	0	1
351078587	16	3	0	2
351312836	68	12	2	3
354698510	92	10	2	3
354716830	157	34	11	9
354792280	55	26	50	20
356459709	10	2	0	0
357111303	20	19	91	34
358162030	251	23	5	7
360666444	32	2	0	0
361871527	120	11	2	2
362874964	105	30	9	10
365146902	15	2	0	0
365860229	28	2	0	2
367173963	84	28	40	20
375688431	78	31	26	16
376505122	35	40	11	10
376692017	102	29	20	12
377766421	73	22	38	19
378472849	129	2	0	1
383242049	9	28	6	9
384492406	179	7	18	6
384891198	4	3	0	0
385805472	29	2	0	0
387948195	118	10	12	8
392984608	103	31	17	11
393764674	11	40	10	13
394287936	87	13	3	5
395892189	130	54	133	54
396279836	119	29	8	8
399877282	94	8	2	2
401379469	122	18	31	13
402449552	80	42	92	38
402717068	11	50	30	20
402780161	46	2	0	1
403294078	65	42	10	13
403805130	77	2	0	0
406201146	19	10	2	6
410070392	49	40	14	11
411510105	11	2	0	1
413373958	16	2	0	1
417232208	335	2	0	2
418274123	92	60	15	17
419810729	115	32	45	21
423858722	119	50	49	26
430353698	49	23	10	8
430790027	155	2	0	2
435087249	100	18	35	14
435910345	15	22	5	7
438426230	22	2	0	0
438910180	39	3	0	2
442577005	98	76	20	22
444222959	12	2	0	0
445028776	14	2	0	2
447043685	37	19	4	5
452867822	42	2	0	1
455085174	65	2	0	2
455924458	13	2	0	0
456059645	124	46	11	13
458629065	84	2	0	1
458642206	18	2	0	0
462779857	86	2	0	2
467902384	76	26	28	15
468498662	44	60	44	26
468708662	39	92	75	43
470437598	19	23	8	7
470781859	21	12	7	3
471636216	104	15	4	6
472408093	48	2	0	0
474198174	39	24	7	9
476597091	299	2	0	1
477657737	18	2	0	1
477873550	39	80	31	24
477922043	100	2	0	0
478906375	71	50	25	19
479197700	77	2	2	2
479401506	56	20	4	5
479552306	38	22	11	8
480257379	17	8	1	0
482466956	13	2	0	1
482811141	23	23	16	11
483118700	63	9	22	9
483738000	48	44	36	20
483982930	21	42	10	11
486255587	59	10	1	4
488916278	43	7	3	3
489523596	10	2	0	0
489651183	51	19	11	8
489920977	13	2	0	2
490710942	90	9	1	3
491389162	133	2	0	1
492595299	38	15	20	9
493652994	376	25	49	20
493735234	158	19	76	31
494809813	33	29	7	9
495722705	13	2	0	2
497007887	66	42	46	24
497462799	98	29	7	8
497767311	45	34	10	11
497776544	35	2	0	1
502423820	21	4	0	2
505129228	34	28	6	7
505681500	74	7	5	3
509749667	47	13	2	4
519925068	160	40	11	13
520541916	105	17	3	4
523082179	43	2	0	2
523390654	61	12	9	7
524655821	486	37	45	24
524713690	11	19	16	10
527976247	52	10	54	20
529020378	22	2	0	2
532125100	25	13	27	13
532411315	28	58	162	66
535123721	59	19	4	5
535721744	20	20	51	21
535922130	135	19	4	7
536556002	45	10	1	5
536664276	96	49	12	15
539593685	69	47	12	15
540291272	51	2	0	2
543674226	156	10	1	5
543746587	38	2	1	3
545009339	101	26	22	12
545234790	152	2	0	2
547802603	13	18	4	4
548056973	61	4	0	0
552062896	11	3	0	1
552369098	40	2	0	1
554317772	2	3	0	1
554324963	13	9	1	3
555379264	10	2	0	1
555552154	49	34	51	24
558036968	181	45	55	27
558306285	12	2	0	1
559638324	13	14	14	7
561003868	10	20	7	7
563161555	213	2	0	2
566719682	67	19	16	10
568531799	99	2	0	2
569695490	43	10	1	2
569877596	43	12	11	7
571867042	34	7	1	2
571976021	48	2	0	2
574205255	35	32	19	14
575994348	48	3	9	5
577313418	10	2	0	0
578412313	11	27	6	8
579031491	55	26	50	23
579973017	14	6	0	4
580981235	82	2	0	0
581773639	24	56	14	13
583123681	141	9	1	1
584807224	69	20	158	58
587609179	106	16	3	5
589507892	121	10	3	3
591362700	52	8	10	6
591577054	36	2	0	0
592430068	49	25	30	14
593201235	46	16	7	6
593264201	58	12	14	8
596244596	39	28	62	25
598138224	143	21	8	7
598883375	58	63	46	28
600250936	23	34	8	9
600913566	93	40	17	16
601784559	88	10	38	14
602281615	81	2	0	2
604350016	10	2	0	1
611798503	15	2	0	2
614200055	42	12	2	4
615270207	34	63	30	22
617553313	68	2	0	0
619192689	91	4	0	0
620105739	67	4	0	0
621715163	77	2	0	3
625142619	72	32	8	9
626475754	13	2	0	1
627404797	20	2	0	1
629042564	63	13	35	15
629547279	4	6	0	0
630012022	82	2	21	9
630255521	12	2	0	2
634066419	31	19	4	7
634439388	43	14	4	4
635277388	64	29	17	10
635688238	19	108	42	34
635881317	47	2	0	0
636890074	18	7	1	2
639051396	77	25	10	8
639964508	62	32	9	8
640337161	55	30	28	15
641265316	79	2	0	0
641677722	7	10	1	3
642636074	39	53	22	18
644021697	3	2	0	2
644796571	100	50	33	21
645586747	16	14	4	4
645794972	165	35	94	38
646561461	36	11	3	5
659900820	113	7	1	2
660194097	25	22	15	10
660604731	10	3	0	0
661877968	30	63	16	17
662977430	127	2	0	0
664579271	173	2	0	1
666676324	402	2	0	0
666756456	78	95	52	36
666761558	209	25	8	9
671053761	88	22	9	8
671845060	9	8	1	3
674932407	62	22	38	16
676426784	39	16	74	29
676467048	110	32	31	18
676509327	15	5	0	1
678947954	64	37	13	10
682126197	11	20	4	6
684888661	83	2	0	0
685022257	129	15	3	4
685376162	69	2	0	1
685798745	43	20	18	10
686969673	16	2	0	0
688647604	11	2	0	1
689084566	15	5	0	2
692402551	56	25	10	8
694380312	48	20	32	14
696932273	43	17	6	7
697539758	58	58	53	28
698082185	36	26	7	8
699582886	83	9	1	3
702031302	10	26	6	6
702435137	70	27	6	8
703124286	19	68	17	19
704600057	78	52	149	59
705396948	83	2	0	3
706020739	54	10	3	3
710384631	27	2	0	0
712451029	13	2	0	0
713738344	58	73	20	20
714552701	21	2	0	0
714775856	12	8	3	5
714797559	88	15	28	11
714955919	86	20	34	16
715277139	89	2	0	0
718259860	9	40	23	15
718955394	86	31	56	27
719473522	21	2	0	1
722896108	55	7	5	2
723864727	184	2	0	0
724228916	61	44	54	27
727393472	29	22	20	10
729115343	69	24	30	17
733122336	23	49	12	13
737289482	11	2	0	1
737355159	78	6	0	3
738494206	35	2	0	0
741931256	23	2	0	2
742103587	4	39	11	11
744625870	106	3	1	2
745257584	28	28	6	8
747018468	31	16	6	5
750405765	78	13	24	10
750876578	122	24	70	29
751081613	125	2	0	0
753744433	68	10	7	6

相关文章

QQ ZHOU:从零开始入门R语言-软件下载及安装​zhuanlan.zhihu.com
QQ ZHOU:从零开始入门R语言—Rstudio下载与安装​zhuanlan.zhihu.com
QQ ZHOU:从0开始入门R语言—DataFrame及聚类算法​zhuanlan.zhihu.com
8e8b46d5b5b050abe5fe97d0baba98af.png
QQ ZHOU:从0开始入门R语言—决策树之C5.0​zhuanlan.zhihu.com
55cd3fd5a8885c15c09dff47ce79c56a.png
QQ ZHOU:R语言实现推荐算法—基于用户协同过滤​zhuanlan.zhihu.com
041c537abc56ddf7690ef35fb42c8211.png
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值