PCA 1bprot
Variables
merged2 <- merged %>%
remove_rownames %>%
column_to_rownames(var="Nom_Barri") %>%
select("n.tot","pc.esp","pc.ue27-esp","pc.20.34","2019-2014","n.esp.V1419",#"var20192014","pc.ext",
"hotel2019","rest1614",
"RFD.2017",
"tot_ann","pmedio","pmedio.V1519","pm_ent.V1519","pm_priv.V1519","pm_sha.V1519",
#"pc.ent","pc.priv","pc.shared","pc.hotel","pm_ent","pm_priv","pm_sha","pm_hot",
#"alq.mq",
"alq.num","alq.pm","alq.pm.V1519","alq.num.V1519",
"tot.comp","tot.eur","perc.nou.comp","perc.prot.comp","tot.comp.V1519",
"nou.eur.V1519","usat.eur.V1519",
#"perc.prot.comp",
#"tot.m2","nou.m2","prot.m2","usat.m2",
#"tot.eur","nou.eur","usat.eur",
#"nou.eurm2","usat.eurm2",
)Standardizar los datos
df <- scale(merged2)Calcular eigenvalues and eigenvectors
pca_result <- prcomp(merged2, scale = TRUE)
pca_result$center## n.tot pc.esp pc.ue27-esp pc.20.34 2019-2014
## 22607.643836 80.781048 7.114347 18.988279 658.342466
## n.esp.V1419 hotel2019 rest1614 RFD.2017 tot_ann
## -1.234187 276.054795 1.191781 93.671233 279.835616
## pmedio pmedio.V1519 pm_ent.V1519 pm_priv.V1519 pm_sha.V1519
## 98.551822 35.607051 59.533227 24.523801 113.220750
## alq.num alq.pm alq.pm.V1519 alq.num.V1519 tot.comp
## 167.369863 858.019940 34.605289 50.827284 185.191781
## tot.eur perc.nou.comp perc.prot.comp tot.comp.V1519 nou.eur.V1519
## 309.594521 7.687534 2.268758 55.060726 21.814261
## usat.eur.V1519
## 49.557724pca_result$scale## n.tot pc.esp pc.ue27-esp pc.20.34 2019-2014
## 14753.886904 9.009606 7.320092 3.873115 1490.085971
## n.esp.V1419 hotel2019 rest1614 RFD.2017 tot_ann
## 5.196927 490.230782 26.158522 42.709775 443.819033
## pmedio pmedio.V1519 pm_ent.V1519 pm_priv.V1519 pm_sha.V1519
## 94.233739 74.042253 111.076238 47.613605 594.339954
## alq.num alq.pm alq.pm.V1519 alq.num.V1519 tot.comp
## 159.168980 357.725402 15.669486 48.467147 148.403509
## tot.eur perc.nou.comp perc.prot.comp tot.comp.V1519 nou.eur.V1519
## 183.338837 10.382977 17.543273 83.625492 96.249820
## usat.eur.V1519
## 35.928599#pesos asignados a cada componente
pca_result$rotation[,1:5]## PC1 PC2 PC3 PC4 PC5
## n.tot -0.287220281 0.002943826 -0.099299664 0.34839672 -0.14168355
## pc.esp 0.224832762 0.304182810 -0.120436897 0.16630550 -0.14087554
## pc.ue27-esp -0.279224758 -0.177719023 0.049471579 -0.33015783 0.14311607
## pc.20.34 -0.251961994 -0.254633603 0.226626769 -0.19027998 0.04289683
## 2019-2014 0.016976469 0.212240412 -0.196988885 0.31221363 -0.10061582
## n.esp.V1419 0.183514599 0.101782872 -0.064601780 -0.13236273 -0.44177168
## hotel2019 -0.317006647 -0.020573233 0.013715354 0.02636191 -0.17376357
## rest1614 0.108735783 0.199376689 -0.077532799 0.13917449 0.10907129
## RFD.2017 -0.182270935 0.267079718 -0.318143377 -0.21021650 0.07235050
## tot_ann -0.343775615 -0.116239885 0.068258335 -0.05001297 -0.09159847
## pmedio -0.153359132 0.353288370 0.289489312 -0.05520862 -0.02837579
## pmedio.V1519 -0.178969201 0.307384497 0.229499933 0.11668365 -0.03106963
## pm_ent.V1519 -0.100998072 0.337635360 0.246527106 0.07457729 -0.03603154
## pm_priv.V1519 -0.113755742 0.143643893 0.328394549 0.07602683 0.18197344
## pm_sha.V1519 -0.144937172 0.009081225 0.013210707 0.08769194 -0.17437311
## alq.num -0.335347053 0.026211443 -0.057735860 0.16600464 -0.16620658
## alq.pm -0.230202541 0.171478392 -0.325534620 -0.06453845 0.24544063
## alq.pm.V1519 -0.188761155 -0.090441353 -0.174834808 0.23409718 0.25161392
## alq.num.V1519 -0.018471070 0.256365258 0.010937530 -0.28035983 -0.01729249
## tot.comp -0.232848968 -0.080035782 -0.048855619 0.31500264 -0.18376302
## tot.eur -0.178997702 0.241323631 -0.313425118 -0.20692872 0.09458074
## perc.nou.comp 0.002718912 -0.048014065 -0.104727340 -0.17452269 -0.50432349
## perc.prot.comp 0.018362966 0.211285332 0.433080019 -0.07338559 -0.05203209
## tot.comp.V1519 0.162040913 -0.164933519 0.100976644 0.26690769 0.06239378
## nou.eur.V1519 -0.105459003 -0.148883516 0.002805124 -0.08305581 -0.38245442
## usat.eur.V1519 -0.052103038 -0.107950023 0.110284327 0.25913149 0.09362090#results
pca_result$x[,1:5]## PC1 PC2
## el Raval -5.90561995 -4.07934267
## el Barri Gòtic -6.23724439 -5.35899385
## la Barceloneta -2.11460234 -3.98120373
## Sant Pere, Santa Caterina i la Ribera -4.81598723 -3.36674170
## el Fort Pienc -2.48714435 -0.42566354
## la Sagrada FamÃlia -3.01163315 -0.16581107
## la Dreta de l'Eixample -7.27835491 1.67726495
## l'Antiga Esquerra de l'Eixample -4.79805593 0.45749927
## la Nova Esquerra de l'Eixample -3.24223250 0.03562211
## Sant Antoni -3.21362123 1.59633524
## el Poble Sec -3.65888076 -2.32850245
## la Marina del Prat Vermell 1.02614915 6.97921952
## la Marina de Port 0.98530565 -0.84809507
## la Font de la Guatlla 0.79050214 -0.22057842
## Hostafrancs -0.23982586 0.06857753
## la Bordeta 1.34874388 0.18623323
## Sants - Badal 0.51578136 -0.07323730
## Sants -1.14622925 0.12132066
## les Corts -1.44336982 0.76486086
## la Maternitat i Sant Ramon -0.44521678 2.56960757
## Pedralbes -1.27974379 2.86448874
## Vallvidrera, el Tibidabo i les Planes 1.23951929 2.13095949
## Sarrià -0.83457461 1.96011385
## les Tres Torres -0.22646151 3.52854626
## Sant Gervasi - la Bonanova -2.28965936 5.84769687
## Sant Gervasi - Galvany -2.61253691 1.35407475
## el Putxet i el Farró -1.97407036 2.71623283
## Vallcarca i els Penitents 0.68567387 0.26406676
## el Coll 1.79456975 -0.71645872
## la Salut 0.86855688 -0.37597686
## la Vila de Grà cia -6.01179173 0.03424970
## el Camp d'en Grassot i Grà cia Nova -1.14232632 0.55503701
## el Baix Guinardó -0.29376329 0.05906344
## Can Baró 1.42605239 -0.34709279
## el Guinardó 0.01577852 -0.30392776
## la Font d'en Fargues 1.74547949 1.86872640
## el Carmel -0.12119942 -0.07509188
## la Teixonera 1.68680961 -1.41558129
## Sant GenÃs dels Agudells 2.40840685 -1.23144860
## Montbau 2.42262338 -1.08199791
## la Vall d'Hebron 2.72918727 0.54823941
## la Clota 4.78613615 -1.43972750
## Horta 1.14123563 -0.82326788
## Vilapicina i la Torre Llobeta 1.44386120 -0.11946506
## Porta 1.12163444 0.26581908
## el Turó de la Peira 1.63882563 -1.07839135
## Can Peguera 4.16641342 -1.80490922
## la Guineueta 2.62748227 -0.15151812
## Canyelles 2.93199821 1.18142578
## les Roquetes 1.78638655 -2.42584354
## Verdun 1.44297324 -1.17242651
## la Prosperitat 1.74112980 -0.96792553
## la Trinitat Nova 2.32484867 -2.21431794
## Torre Baró 3.19126177 0.23705004
## Ciutat Meridiana 1.86431152 -2.00671682
## Vallbona 4.37443130 -1.07858400
## la Trinitat Vella 1.52946387 -2.05106101
## Baró de Viver 3.09435136 -0.65126751
## el Bon Pastor 2.54803448 -1.02718693
## Sant Andreu -0.57713964 0.69010690
## la Sagrera 0.68759751 -0.46532119
## el Congrés i els Indians 0.99253684 0.25346258
## Navas 0.71084638 -0.98719010
## el Camp de l'Arpa del Clot -1.74929193 0.02945872
## el Clot -0.63640685 -0.14501094
## el Parc i la Llacuna del Poblenou 0.51460028 0.88746672
## la Vila OlÃmpica del Poblenou 0.05598277 3.14164798
## el Poblenou -1.33353499 1.31727153
## Diagonal Mar i el Front MarÃtim del Poblenou 0.18583443 3.03471334
## el Besòs i el Maresme -0.58287879 -0.99869703
## Provençals del Poblenou 0.79185953 0.12868411
## Sant Martà de Provençals 1.20839208 -0.34845780
## la Verneda i la Pau 1.11182915 -1.00211165
## PC3 PC4
## el Raval 1.2369419143 -0.47876710
## el Barri Gòtic 2.7334233513 -6.05415611
## la Barceloneta 1.0765289168 -1.60745976
## Sant Pere, Santa Caterina i la Ribera 0.9753829311 -2.17025709
## el Fort Pienc -0.5814350627 -0.79571515
## la Sagrada FamÃlia 0.3355486263 0.99336699
## la Dreta de l'Eixample -0.5548669866 0.91504290
## l'Antiga Esquerra de l'Eixample -0.1675745428 -0.01868104
## la Nova Esquerra de l'Eixample -0.5709912752 0.67542979
## Sant Antoni 1.3685472194 0.63345314
## el Poble Sec 1.0630025287 1.69496738
## la Marina del Prat Vermell 10.5959405683 -1.32022718
## la Marina de Port -0.4150145163 1.88025343
## la Font de la Guatlla 0.6597945225 -0.64452455
## Hostafrancs 0.7159235569 0.25557682
## la Bordeta 0.2467487708 0.59264656
## Sants - Badal -0.0406490991 1.19110010
## Sants -0.5667599168 1.60236785
## les Corts -1.2256398285 1.05952295
## la Maternitat i Sant Ramon 0.3814445787 -0.39020044
## Pedralbes -4.3803361033 -2.57168733
## Vallvidrera, el Tibidabo i les Planes -0.8846830462 -3.11753141
## Sarrià -3.2276660093 -0.53549835
## les Tres Torres -2.3783668404 -1.52655858
## Sant Gervasi - la Bonanova 0.0471417549 -0.09441515
## Sant Gervasi - Galvany -2.8169008942 -0.47529205
## el Putxet i el Farró -0.0008837824 0.31294479
## Vallcarca i els Penitents -0.4776194535 -0.81285558
## el Coll -0.6654379948 -1.12073302
## la Salut -1.0436099522 -0.43082630
## la Vila de Grà cia 0.5904053415 1.13934699
## el Camp d'en Grassot i Grà cia Nova -0.0990753147 0.71675799
## el Baix Guinardó -0.2198270649 0.64942722
## Can Baró 0.1290144703 -0.33873084
## el Guinardó -0.3921962930 1.28348703
## la Font d'en Fargues 0.4864387054 -0.29030551
## el Carmel 1.4112508357 1.97155671
## la Teixonera -0.5230122968 -0.12535500
## Sant GenÃs dels Agudells 0.5960511102 1.12526403
## Montbau -0.3766522977 -0.76612745
## la Vall d'Hebron -0.8049349460 -1.93220624
## la Clota -1.3320754709 -3.25940268
## Horta -0.8623268108 0.56140814
## Vilapicina i la Torre Llobeta 0.1562062211 1.19665345
## Porta 0.0519756784 1.55022719
## el Turó de la Peira 0.6049761149 1.12008093
## Can Peguera 1.5958948321 0.64788858
## la Guineueta -0.2182425486 -0.24552157
## Canyelles 0.7290622743 -1.39689962
## les Roquetes 0.4122873878 1.44487326
## Verdun 0.3055769537 0.23872583
## la Prosperitat -0.2031716912 0.74986445
## la Trinitat Nova 0.1186800139 0.28994306
## Torre Baró 1.9023409125 -0.84889911
## Ciutat Meridiana 2.2183959416 1.45939973
## Vallbona -0.3733044591 -2.22845572
## la Trinitat Vella -0.1918081162 0.40278821
## Baró de Viver 0.2856529751 -0.97689596
## el Bon Pastor -0.0074718762 1.40106972
## Sant Andreu -0.7340653905 3.58001271
## la Sagrera -0.3411435730 1.17756535
## el Congrés i els Indians 0.9090384083 0.76015710
## Navas -0.7577857795 0.29607368
## el Camp de l'Arpa del Clot 0.1272173534 1.23643172
## el Clot -0.1500397603 1.32302808
## el Parc i la Llacuna del Poblenou -0.1793493575 -0.88053556
## la Vila OlÃmpica del Poblenou -0.6110950289 -1.58471484
## el Poblenou -0.6405152631 0.68082673
## Diagonal Mar i el Front MarÃtim del Poblenou -2.6107342152 -0.95242922
## el Besòs i el Maresme -0.4719549381 0.03184803
## Provençals del Poblenou -1.6442638389 -0.07624608
## Sant Martà de Provençals -0.4031338451 0.08456265
## la Verneda i la Pau 0.0797807101 1.14217032
## PC5
## el Raval -1.156772733
## el Barri Gòtic 1.192036180
## la Barceloneta 2.457863656
## Sant Pere, Santa Caterina i la Ribera 0.882413702
## el Fort Pienc -1.522964066
## la Sagrada FamÃlia 0.886895527
## la Dreta de l'Eixample -1.878745949
## l'Antiga Esquerra de l'Eixample -1.092665868
## la Nova Esquerra de l'Eixample -0.615414753
## Sant Antoni 0.068055581
## el Poble Sec 0.713987251
## la Marina del Prat Vermell -0.677460272
## la Marina de Port 0.990566132
## la Font de la Guatlla 0.998361369
## Hostafrancs 0.737360515
## la Bordeta 0.667052645
## Sants - Badal 1.409408804
## Sants 0.116271315
## les Corts -0.658508396
## la Maternitat i Sant Ramon 1.299696569
## Pedralbes 2.001701860
## Vallvidrera, el Tibidabo i les Planes 1.027490786
## Sarrià 1.288850656
## les Tres Torres 1.283630965
## Sant Gervasi - la Bonanova -0.596043175
## Sant Gervasi - Galvany -1.409169953
## el Putxet i el Farró -0.855618506
## Vallcarca i els Penitents -0.353800530
## el Coll -0.359142726
## la Salut 0.742176719
## la Vila de Grà cia -2.336339768
## el Camp d'en Grassot i Grà cia Nova 0.860193810
## el Baix Guinardó 0.198762003
## Can Baró 0.948993397
## el Guinardó -0.221371248
## la Font d'en Fargues 0.007789659
## el Carmel -0.999048494
## la Teixonera -1.025388503
## Sant GenÃs dels Agudells -0.212973675
## Montbau -0.810912403
## la Vall d'Hebron -2.126440026
## la Clota -6.364259230
## Horta -1.096209868
## Vilapicina i la Torre Llobeta 1.085084724
## Porta 0.096870481
## el Turó de la Peira 0.976060424
## Can Peguera 0.084840804
## la Guineueta -0.227983218
## Canyelles 0.226729443
## les Roquetes 0.773321614
## Verdun 0.831218571
## la Prosperitat 0.282952958
## la Trinitat Nova 1.061890317
## Torre Baró -0.992333339
## Ciutat Meridiana 1.166402516
## Vallbona -1.155217966
## la Trinitat Vella 1.176937668
## Baró de Viver -0.228141593
## el Bon Pastor 0.550312886
## Sant Andreu -1.710826847
## la Sagrera 0.323022307
## el Congrés i els Indians 1.300902694
## Navas -0.412524030
## el Camp de l'Arpa del Clot -0.974545562
## el Clot 0.297736201
## el Parc i la Llacuna del Poblenou 0.858569867
## la Vila OlÃmpica del Poblenou 2.140351344
## el Poblenou -0.959397386
## Diagonal Mar i el Front MarÃtim del Poblenou 1.876925846
## el Besòs i el Maresme -1.862569355
## Provençals del Poblenou 0.083008336
## Sant Martà de Provençals -0.639174104
## la Verneda i la Pau -0.440734558Var Exp
VE <- pca_result$sdev^2
PVE <- VE / sum(VE)
round(PVE, 2)## [1] 0.25 0.15 0.11 0.08 0.07 0.05 0.04 0.04 0.03 0.03 0.03 0.02 0.02 0.02 0.01
## [16] 0.01 0.01 0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00vis eigen
fviz_eig(pca_result)
Plot
biplot(pca_result, scale = 0)
graph pca
fviz_pca_ind(pca_result,
col.ind = "cos2",
gradient.cols = c("#00AFBB", "#E7B800", "#FC4E07"),
repel = TRUE
)
graph var
fviz_pca_var(pca_result,
col.var = "contrib",
gradient.cols = c("#00AFBB", "#E7B800", "#FC4E07"),
repel = TRUE
)
graph biplot
fviz_pca_biplot(pca_result, repel = TRUE,
col.var = "#2E9FDF", # Variables color
col.ind = "#696969" # Individuals color
)