PCA 1prot
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.pmq","alq.pmq.V1519","alq.num.V1519",
"tot.comp","tot.eurm2","perc.nou.comp","perc.prot.comp","tot.comp.V1519",
"nou.eurm2.V1519","usat.eurm2.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
## 22607.6438356 80.7810477 7.1143474 18.9882789
## 2019-2014 n.esp.V1419 hotel2019 rest1614
## 658.3424658 -1.2341867 276.0547945 1.1917808
## RFD.2017 tot_ann pmedio pmedio.V1519
## 93.6712329 279.8356164 98.5518223 -87.3018000
## pm_ent.V1519 pm_priv.V1519 pm_sha.V1519 alq.num
## 0.8279582 491.4042876 41.4094048 167.3698630
## alq.pmq alq.pmq.V1519 alq.num.V1519 tot.comp
## 12.5311589 35.3976621 50.8272836 185.1917808
## tot.eurm2 perc.nou.comp perc.prot.comp tot.comp.V1519
## 3722.8534247 7.6875336 2.2687579 55.0607260
## nou.eurm2.V1519 usat.eurm2.V1519
## 13.8070105 448.2268587pca_result$scale## n.tot pc.esp pc.ue27-esp pc.20.34
## 14753.886904 9.009606 7.320092 3.873115
## 2019-2014 n.esp.V1419 hotel2019 rest1614
## 1490.085971 5.196927 490.230782 26.158522
## RFD.2017 tot_ann pmedio pmedio.V1519
## 42.709775 443.819033 94.233739 28.320141
## pm_ent.V1519 pm_priv.V1519 pm_sha.V1519 alq.num
## 57.451326 1916.511004 119.838533 159.168980
## alq.pmq alq.pmq.V1519 alq.num.V1519 tot.comp
## 4.260785 15.633887 48.467147 148.403509
## tot.eurm2 perc.nou.comp perc.prot.comp tot.comp.V1519
## 1266.180075 10.382977 17.543273 83.625492
## nou.eurm2.V1519 usat.eurm2.V1519
## 67.962390 3428.920989#pesos asignados a cada componente
pca_result$rotation[,1:5]## PC1 PC2 PC3 PC4 PC5
## n.tot -0.280062380 0.008117929 -0.160473278 0.32698560 -0.03256051
## pc.esp 0.235556864 -0.204715304 -0.254227351 0.16710037 0.15878587
## pc.ue27-esp -0.288805582 0.091154636 0.179514851 -0.33948927 -0.05034289
## pc.20.34 -0.254248642 0.117065324 0.331411296 -0.14338015 -0.09298993
## 2019-2014 0.019585930 -0.148396038 -0.316223359 0.27356934 0.03254738
## n.esp.V1419 0.181789452 -0.081502924 -0.040558108 -0.03566676 0.43005983
## hotel2019 -0.322626502 -0.048295843 0.065456820 0.16355594 0.13484305
## rest1614 0.107262592 -0.153585359 -0.184327926 0.12109288 -0.01325183
## RFD.2017 -0.171484814 -0.190702916 -0.302642783 -0.25926538 0.17665077
## tot_ann -0.353209198 0.018639012 0.142927681 0.05567696 0.04065077
## pmedio -0.110855954 -0.473800392 0.101270132 0.02634232 -0.07227003
## pmedio.V1519 -0.209731118 -0.017755636 0.146814909 0.24809920 0.19819412
## pm_ent.V1519 -0.135790237 -0.342487593 0.024992530 0.12859507 -0.08492712
## pm_priv.V1519 0.042899443 -0.420599552 0.273707108 0.06174934 -0.09911554
## pm_sha.V1519 -0.141929400 -0.054957840 0.004460411 -0.07194796 0.16336815
## alq.num -0.329283914 -0.053864522 -0.079060140 0.20189147 0.04289234
## alq.pmq -0.228398126 0.016170513 -0.339248184 -0.18140551 -0.21457929
## alq.pmq.V1519 -0.122986265 0.089769050 -0.297134133 0.10811163 -0.28013659
## alq.num.V1519 -0.007078160 -0.285496363 -0.063418295 -0.25713965 0.07568486
## tot.comp -0.229057878 0.085282732 -0.062732955 0.32017142 0.02225100
## tot.eurm2 -0.230799444 -0.103334105 -0.253002901 -0.29113026 0.14963436
## perc.nou.comp -0.002779745 0.074918539 0.002539556 -0.04810266 0.49827702
## perc.prot.comp 0.059390010 -0.391490831 0.297926824 0.05927161 -0.08792449
## tot.comp.V1519 0.155527212 0.195832951 0.102661394 0.28924209 -0.12822983
## nou.eurm2.V1519 -0.087292694 0.119969467 0.136313432 0.15675518 0.40126491
## usat.eurm2.V1519 0.035926859 0.019655803 0.093696362 0.03041667 0.22056150#results
pca_result$x[,1:5]## PC1 PC2
## el Raval -6.46821553 2.250889004
## el Barri Gòtic -6.53637853 3.221482488
## la Barceloneta -2.48367557 2.426571931
## Sant Pere, Santa Caterina i la Ribera -5.34077887 1.663974502
## el Fort Pienc -2.18803263 0.006734010
## la Sagrada FamÃlia -3.19893206 -0.536957883
## la Dreta de l'Eixample -8.21089275 -2.148630842
## l'Antiga Esquerra de l'Eixample -4.82320699 -1.038768991
## la Nova Esquerra de l'Eixample -3.40746963 -0.029197514
## Sant Antoni -3.25957220 -1.842107690
## el Poble Sec -3.65878151 1.032608473
## la Marina del Prat Vermell 3.35459105 -11.428822507
## la Marina de Port 1.19628113 1.043947060
## la Font de la Guatlla 0.74455924 0.055643142
## Hostafrancs -0.13399290 -0.269963384
## la Bordeta 1.30087556 0.310431074
## Sants - Badal 0.60725536 0.045190232
## Sants -1.34741543 -0.305257670
## les Corts -1.59943924 -0.328146080
## la Maternitat i Sant Ramon 0.08540222 -1.620437848
## Pedralbes -0.41354772 -0.961640918
## Vallvidrera, el Tibidabo i les Planes 0.44320770 -1.777742339
## Sarrià -0.44134508 -0.959325675
## les Tres Torres 0.14462410 -2.095454786
## Sant Gervasi - la Bonanova -1.46513122 -3.814839891
## Sant Gervasi - Galvany -2.59048846 -0.686445388
## el Putxet i el Farró -1.46549523 -2.003915310
## Vallcarca i els Penitents 0.66578106 0.020967388
## el Coll 1.69474726 0.728972910
## la Salut 0.68353362 0.419933299
## la Vila de Grà cia -5.48998208 -0.807943440
## el Camp d'en Grassot i Grà cia Nova -1.16896771 -0.255698096
## el Baix Guinardó 0.01249128 0.084403090
## Can Baró 1.32729576 0.513157510
## el Guinardó 0.03471537 0.506674599
## la Font d'en Fargues 1.94272014 -1.205200352
## el Carmel 0.26022155 0.869643573
## la Teixonera 1.72325507 1.695273515
## Sant GenÃs dels Agudells 2.65251511 1.296083879
## Montbau 2.24958379 1.015980442
## la Vall d'Hebron 2.23146275 0.078466908
## la Clota 4.25696706 0.756892907
## Horta 1.06703932 0.972019044
## Vilapicina i la Torre Llobeta 1.66466322 0.512349688
## Porta 1.13959338 0.368005224
## el Turó de la Peira 1.75519162 1.355734141
## Can Peguera 3.52958715 1.516112388
## la Guineueta 2.32355321 0.128085115
## Canyelles 3.13080168 -0.762857281
## les Roquetes 1.61411421 1.736321816
## Verdun 1.48785222 1.128836208
## la Prosperitat 1.62130220 0.975710191
## la Trinitat Nova 1.89084332 1.306053931
## Torre Baró 3.16861153 -0.280060431
## Ciutat Meridiana 2.23867147 2.198456584
## Vallbona 3.98462867 0.433207864
## la Trinitat Vella 1.51703360 1.246212368
## Baró de Viver 2.21261016 0.624236916
## el Bon Pastor 2.61303859 1.085798864
## Sant Andreu -0.32004088 0.001300455
## la Sagrera 0.57249931 0.682722212
## el Congrés i els Indians 1.42102484 0.124403224
## Navas 0.68739682 1.202966824
## el Camp de l'Arpa del Clot -1.61911029 0.244725242
## el Clot -0.88103124 0.497068898
## el Parc i la Llacuna del Poblenou 0.28886202 -0.650692932
## la Vila OlÃmpica del Poblenou 0.09754504 -2.589573311
## el Poblenou -1.68439839 -1.213716070
## Diagonal Mar i el Front MarÃtim del Poblenou 0.02762610 -2.299980115
## el Besòs i el Maresme -0.52835350 0.997522038
## Provençals del Poblenou 0.77148531 0.377356521
## Sant Martà de Provençals 1.07067018 0.940059934
## la Verneda i la Pau 1.21634429 1.214189122
## PC3 PC4
## el Raval 2.95065447 0.797414247
## el Barri Gòtic 5.89858088 -5.044109079
## la Barceloneta 2.54450043 -2.792647457
## Sant Pere, Santa Caterina i la Ribera 2.57478885 -1.979276546
## el Fort Pienc -0.16807207 -0.088184861
## la Sagrada FamÃlia -0.23880212 0.508602865
## la Dreta de l'Eixample 0.14745817 3.368641244
## l'Antiga Esquerra de l'Eixample 0.08902364 0.857092583
## la Nova Esquerra de l'Eixample -0.40042329 1.139531236
## Sant Antoni 0.41937891 1.133559225
## el Poble Sec 0.98525521 0.536964670
## la Marina del Prat Vermell 7.75939461 0.933247105
## la Marina de Port -0.88043167 1.313418886
## la Font de la Guatlla 0.65609621 -1.099419797
## Hostafrancs 0.23800908 -0.362453155
## la Bordeta -0.78801638 0.081722298
## Sants - Badal -0.83094141 0.043988326
## Sants -1.02315214 1.179575001
## les Corts -1.70669997 1.158480339
## la Maternitat i Sant Ramon -1.37226091 -1.221780171
## Pedralbes -2.71791548 -2.824672622
## Vallvidrera, el Tibidabo i les Planes -0.52107026 -3.493748514
## Sarrià -2.76444382 -1.330374343
## les Tres Torres -2.69029977 -1.940787326
## Sant Gervasi - la Bonanova -2.18903825 -0.274087237
## Sant Gervasi - Galvany -2.14675161 -0.008187678
## el Putxet i el Farró -0.87380891 0.365946646
## Vallcarca i els Penitents -0.53932097 -0.961430527
## el Coll 0.12690416 -0.844172724
## la Salut -0.67685338 -1.083728492
## la Vila de Grà cia 0.56702637 0.990016302
## el Camp d'en Grassot i Grà cia Nova -1.01625270 0.174080245
## el Baix Guinardó -0.76725095 -0.374027069
## Can Baró -0.34697832 -0.572721002
## el Guinardó -0.51929693 0.967793061
## la Font d'en Fargues -0.18666669 -0.310451239
## el Carmel 0.27007497 1.968421972
## la Teixonera -0.05780385 0.310715479
## Sant GenÃs dels Agudells 0.83580657 1.174436721
## Montbau 0.30356814 -0.417728644
## la Vall d'Hebron -0.78862575 -1.207700482
## la Clota 1.45152255 -1.292069487
## Horta -0.69177761 0.658325476
## Vilapicina i la Torre Llobeta -0.40981121 0.348995726
## Porta -1.11870374 1.115561152
## el Turó de la Peira 0.09386486 0.923850648
## Can Peguera 3.22873191 2.689073020
## la Guineueta -0.10044259 -0.021024874
## Canyelles 0.52128106 -1.204558689
## les Roquetes 0.65935451 0.615722422
## Verdun 0.32579187 -0.056688466
## la Prosperitat -0.17224054 0.606983047
## la Trinitat Nova 0.64614569 0.093262081
## Torre Baró 2.07261675 0.053785820
## Ciutat Meridiana 1.62702134 1.167690019
## Vallbona 1.76014752 -0.680631442
## la Trinitat Vella 0.84020880 -0.340260650
## Baró de Viver 2.55395782 0.575040239
## el Bon Pastor 0.25232012 1.162172609
## Sant Andreu -1.68450146 3.248440534
## la Sagrera -1.04111509 0.634764919
## el Congrés i els Indians -0.45722094 -0.266330014
## Navas -0.36514331 0.497738181
## el Camp de l'Arpa del Clot -0.22476469 1.449803083
## el Clot -0.81405126 0.517030579
## el Parc i la Llacuna del Poblenou -0.97049417 -1.201182006
## la Vila OlÃmpica del Poblenou -2.02905313 -2.622277860
## el Poblenou -1.34638693 0.484312037
## Diagonal Mar i el Front MarÃtim del Poblenou -3.36316765 -1.572839959
## el Besòs i el Maresme 0.10002718 0.624336995
## Provençals del Poblenou -1.64793558 -0.483994857
## Sant Martà de Provençals -0.08645695 0.369347625
## la Verneda i la Pau 0.23493181 1.133662610
## PC5
## el Raval 0.042981200
## el Barri Gòtic -0.180342701
## la Barceloneta -1.601590571
## Sant Pere, Santa Caterina i la Ribera -0.624217571
## el Fort Pienc 0.316754340
## la Sagrada FamÃlia -1.214926700
## la Dreta de l'Eixample 2.869834777
## l'Antiga Esquerra de l'Eixample 0.446943690
## la Nova Esquerra de l'Eixample 0.367711853
## Sant Antoni -0.554419558
## el Poble Sec -1.026179371
## la Marina del Prat Vermell -1.469146492
## la Marina de Port -1.299071441
## la Font de la Guatlla -0.518424290
## Hostafrancs -0.541478217
## la Bordeta -1.394976605
## Sants - Badal -1.604146852
## Sants 0.332732243
## les Corts 0.753565412
## la Maternitat i Sant Ramon -0.953074965
## Pedralbes 0.731461995
## Vallvidrera, el Tibidabo i les Planes 1.661974120
## Sarrià 0.177160765
## les Tres Torres -0.186114496
## Sant Gervasi - la Bonanova 0.487804383
## Sant Gervasi - Galvany 1.764739573
## el Putxet i el Farró 0.791233951
## Vallcarca i els Penitents 0.969432475
## el Coll 0.548811795
## la Salut -0.176400471
## la Vila de Grà cia 0.540933745
## el Camp d'en Grassot i Grà cia Nova -0.984108086
## el Baix Guinardó -1.020164899
## Can Baró -1.065816903
## el Guinardó 0.096956967
## la Font d'en Fargues -0.393013986
## el Carmel 0.252684081
## la Teixonera 1.001452298
## Sant GenÃs dels Agudells 0.234744562
## Montbau 1.179706847
## la Vall d'Hebron 2.522424415
## la Clota 6.908442231
## Horta 0.780957363
## Vilapicina i la Torre Llobeta -1.450575393
## Porta -0.960009755
## el Turó de la Peira -1.540646629
## Can Peguera 0.006958327
## la Guineueta -0.080484073
## Canyelles 0.300856828
## les Roquetes -1.187881990
## Verdun -1.135279854
## la Prosperitat -0.929807290
## la Trinitat Nova -1.774172627
## Torre Baró 0.706196477
## Ciutat Meridiana -1.698402351
## Vallbona 0.935155264
## la Trinitat Vella -1.406439909
## Baró de Viver 3.418088457
## el Bon Pastor -0.623768043
## Sant Andreu 0.263872612
## la Sagrera -0.740764712
## el Congrés i els Indians -0.997746018
## Navas 0.393385409
## el Camp de l'Arpa del Clot 0.054084355
## el Clot -0.035938441
## el Parc i la Llacuna del Poblenou -1.197899313
## la Vila OlÃmpica del Poblenou -1.106959845
## el Poblenou 0.816142547
## Diagonal Mar i el Front MarÃtim del Poblenou -1.269965014
## el Besòs i el Maresme 1.103784958
## Provençals del Poblenou 0.278235023
## Sant Martà de Provençals 0.949347547
## la Verneda i la Pau -0.063197451Var Exp
VE <- pca_result$sdev^2
PVE <- VE / sum(VE)
round(PVE, 2)## [1] 0.25 0.13 0.12 0.08 0.07 0.07 0.05 0.04 0.03 0.03 0.03 0.02 0.02 0.01 0.01
## [16] 0.01 0.01 0.00 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
)