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.usat.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.usat.comp tot.comp.V1519 nou.eur.V1519
## 309.594521 7.687534 91.924106 55.060726 21.814261
## usat.eur.V1519
## 49.557724
pca_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.usat.comp tot.comp.V1519 nou.eur.V1519
## 183.338837 10.382977 10.563166 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.286729978 -0.021582928 0.037099456 -0.15160979 0.35633104
## pc.esp -0.224096321 -0.324517636 -0.025709562 -0.03909388 0.19606773
## pc.ue27-esp 0.278681729 0.189581923 -0.063903008 0.04523500 -0.34239330
## pc.20.34 0.251950783 0.295611899 0.028181945 0.18988729 -0.17605079
## 2019-2014 -0.016812157 -0.245405653 0.002103398 -0.21137615 0.31121494
## n.esp.V1419 -0.183649506 -0.100611476 -0.308803300 0.16169830 0.09725100
## hotel2019 0.316905437 0.024468289 -0.025251290 0.04234416 0.07566150
## rest1614 -0.108198845 -0.215264516 0.054770071 -0.07766618 0.08763702
## RFD.2017 0.181946334 -0.312753532 -0.232169293 -0.19125044 -0.19144836
## tot_ann 0.343545577 0.128898575 -0.002226149 0.06345491 -0.02254553
## pmedio 0.156686781 -0.296936240 0.121347840 0.36732419 -0.03309403
## pmedio.V1519 0.181118705 -0.273967764 0.166372265 0.33053608 0.11677566
## pm_ent.V1519 0.103522595 -0.298150522 0.147710212 0.35564038 0.08878587
## pm_priv.V1519 0.116432951 -0.087783179 0.276418759 0.23525151 -0.01396179
## pm_sha.V1519 0.144977640 -0.007982933 0.020751873 0.02564709 0.11349893
## alq.num 0.335117003 -0.036005488 -0.004120047 -0.05269583 0.19677526
## alq.pm 0.229431099 -0.224143855 -0.117598160 -0.31009116 -0.13479742
## alq.pm.V1519 0.187607553 0.060364905 0.017397730 -0.30201086 0.14695906
## alq.num.V1519 0.019814893 -0.248764813 -0.060487626 0.11973594 -0.26093002
## tot.comp 0.232166803 0.072917836 -0.013974427 -0.08441307 0.36858922
## tot.eur 0.178681878 -0.286534472 -0.209765595 -0.21531435 -0.20335988
## perc.nou.comp -0.004244247 0.050256998 -0.526793071 0.21268804 0.16436467
## perc.usat.comp 0.005286657 -0.034061879 0.504567161 -0.26941706 -0.16260407
## tot.comp.V1519 -0.162571476 0.172643678 0.138167684 0.02366186 0.22804270
## nou.eur.V1519 0.104639311 0.161712956 -0.261615134 0.12198748 0.12629338
## usat.eur.V1519 0.052140682 0.123132496 0.132205624 0.01539091 0.22044024
#results
pca_result$x[,1:5]
## PC1 PC2
## el Raval 5.89440408 4.288845098
## el Barri Gòtic 6.21613977 5.934295689
## la Barceloneta 2.10238722 4.122595554
## Sant Pere, Santa Caterina i la Ribera 4.80571465 3.545316628
## el Fort Pienc 2.48150888 0.380396642
## la Sagrada FamÃlia 3.01559898 0.143368839
## la Dreta de l'Eixample 7.27983143 -1.692371488
## l'Antiga Esquerra de l'Eixample 4.79716342 -0.458410072
## la Nova Esquerra de l'Eixample 3.23936438 -0.146186542
## Sant Antoni 3.22921445 -1.417342223
## el Poble Sec 3.65560124 2.437102211
## la Marina del Prat Vermell -0.81362269 -4.003233661
## la Marina de Port -0.98824516 0.691328492
## la Font de la Guatlla -0.78655395 0.264340361
## Hostafrancs 0.23842434 0.042057884
## la Bordeta -1.34424879 -0.229360356
## Sants - Badal -0.51446906 -0.025711760
## Sants 1.14628665 -0.265455922
## les Corts 1.43908814 -0.963803988
## la Maternitat i Sant Ramon 0.46239244 -2.559126920
## Pedralbes 1.26937290 -3.534008480
## Vallvidrera, el Tibidabo i les Planes -1.23012280 -2.288319965
## Sarrià 0.83001985 -2.507799156
## les Tres Torres 0.23302916 -3.931009441
## Sant Gervasi - la Bonanova 2.30852820 -5.778909973
## Sant Gervasi - Galvany 2.59584949 -1.722528286
## el Putxet i el Farró 1.98046347 -2.679889005
## Vallcarca i els Penitents -0.69446695 -0.293182571
## el Coll -1.80859250 0.672523120
## la Salut -0.87406375 0.159386611
## la Vila de Grà cia 6.01170654 0.059507241
## el Camp d'en Grassot i Grà cia Nova 1.14752448 -0.630853033
## el Baix Guinardó 0.29567806 -0.170487639
## Can Baró -1.42411588 0.303229461
## el Guinardó -0.01759316 0.187645928
## la Font d'en Fargues -1.73487867 -1.834189642
## el Carmel 0.11382380 0.354062620
## la Teixonera -1.71149550 1.425431275
## Sant GenÃs dels Agudells -2.42108465 1.321281161
## Montbau -2.44050919 1.132799513
## la Vall d'Hebron -2.74757463 -0.546050849
## la Clota -4.82799474 1.504798759
## Horta -1.15599925 0.720219648
## Vilapicina i la Torre Llobeta -1.44118724 0.039875367
## Porta -1.11860558 -0.335055465
## el Turó de la Peira -1.63889881 1.084058293
## Can Peguera -4.16549163 1.937361863
## la Guineueta -2.62962340 0.064473853
## Canyelles -2.92000906 -1.123540128
## les Roquetes -1.79325827 2.390698997
## Verdun -1.44576540 1.155267637
## la Prosperitat -1.74405604 0.848505782
## la Trinitat Nova -2.33215187 2.132610245
## Torre Baró -3.17956963 0.010258195
## Ciutat Meridiana -1.86075521 2.253986091
## Vallbona -4.37904993 0.960992383
## la Trinitat Vella -1.53763274 1.927536179
## Baró de Viver -3.09370300 0.628347897
## el Bon Pastor -2.55243938 0.949464398
## Sant Andreu 0.57211279 -0.843067607
## la Sagrera -0.68831753 0.336888598
## el Congrés i els Indians -0.98378540 -0.175926399
## Navas -0.71722483 0.831720167
## el Camp de l'Arpa del Clot 1.74593002 -0.007226575
## el Clot 0.63676640 0.075173596
## el Parc i la Llacuna del Poblenou -0.50909372 -0.966248309
## la Vila OlÃmpica del Poblenou -0.04108845 -3.283877781
## el Poblenou 1.32903617 -1.375199913
## Diagonal Mar i el Front MarÃtim del Poblenou -0.18033458 -3.499889674
## el Besòs i el Maresme 0.56021641 1.054446603
## Provençals del Poblenou -0.80744722 -0.355496170
## Sant Martà de Provençals -1.22041152 0.305965934
## la Verneda i la Pau -1.11764609 0.965594181
## PC3 PC4
## el Raval 0.08842016 0.320716992
## el Barri Gòtic -1.84822212 3.076267623
## la Barceloneta 0.75164269 -0.400321973
## Sant Pere, Santa Caterina i la Ribera -0.12635947 0.239993364
## el Fort Pienc -1.38321099 0.004565837
## la Sagrada FamÃlia 1.62107396 -0.567720895
## la Dreta de l'Eixample -0.81952285 0.302833897
## l'Antiga Esquerra de l'Eixample -0.46091333 0.284792376
## la Nova Esquerra de l'Eixample 0.06322903 -0.695639668
## Sant Antoni 1.61065119 1.407860501
## el Poble Sec 1.62925513 -0.358054987
## la Marina del Prat Vermell 3.61174980 8.076641508
## la Marina de Port 1.25121653 -1.607751395
## la Font de la Guatlla 0.99956198 0.293943568
## Hostafrancs 0.23813921 0.868460831
## la Bordeta 1.23574835 -0.264192942
## Sants - Badal 1.40805153 -0.874021440
## Sants 0.74387993 -1.189507325
## les Corts -0.51164028 -1.043786228
## la Maternitat i Sant Ramon 1.41034394 0.381555432
## Pedralbes -2.73781891 -3.189694417
## Vallvidrera, el Tibidabo i les Planes -0.31914282 -0.299626053
## Sarrià -1.02843241 -2.999391334
## les Tres Torres -0.72622544 -1.528704419
## Sant Gervasi - la Bonanova -0.64995153 2.434001717
## Sant Gervasi - Galvany -2.95869535 -1.372614661
## el Putxet i el Farró -0.58635432 1.335474693
## Vallcarca i els Penitents -1.38160690 0.337801076
## el Coll -1.81597421 0.054053193
## la Salut -0.08263057 -1.354956595
## la Vila de Grà cia 0.32687317 0.834680939
## el Camp d'en Grassot i Grà cia Nova 1.17792766 -0.728644520
## el Baix Guinardó 1.00897057 -0.773397066
## Can Baró 0.82484409 -0.394340156
## el Guinardó 0.59346852 -0.946132587
## la Font d'en Fargues 0.62177832 1.043357304
## el Carmel -0.41341077 2.183920129
## la Teixonera -2.56075134 0.510100404
## Sant GenÃs dels Agudells -0.46095787 0.920123656
## Montbau -2.36821671 0.598811100
## la Vall d'Hebron -3.53663869 1.477228797
## la Clota -7.53182668 2.905238196
## Horta -1.43939439 -0.444431702
## Vilapicina i la Torre Llobeta 1.58273263 -0.690447644
## Porta 1.12904833 -0.442819799
## el Turó de la Peira 1.48273133 -0.373804746
## Can Peguera 1.86741846 0.716717297
## la Guineueta 0.13169804 -0.227046613
## Canyelles 0.97958892 0.839115980
## les Roquetes 1.31863077 -0.901131738
## Verdun 0.88272587 -0.455129390
## la Prosperitat 0.87569878 -0.973279503
## la Trinitat Nova 1.00327009 -0.979225872
## Torre Baró 1.09049573 2.205308865
## Ciutat Meridiana 2.46197593 0.881664851
## Vallbona -0.46663312 -0.195299949
## la Trinitat Vella 0.88068951 -1.316039966
## Baró de Viver 0.45891565 0.138115677
## el Bon Pastor 0.86918257 -0.743490589
## Sant Andreu 0.11703255 -0.658593411
## la Sagrera 0.92631257 -1.140493130
## el Congrés i els Indians 1.79228187 0.057143949
## Navas -0.21775049 -1.116933125
## el Camp de l'Arpa del Clot -0.11117413 0.340734704
## el Clot 0.81752820 -0.782483754
## el Parc i la Llacuna del Poblenou 0.57692678 -0.249133724
## la Vila OlÃmpica del Poblenou 0.66793154 -0.247098951
## el Poblenou -1.10500196 0.291553387
## Diagonal Mar i el Front MarÃtim del Poblenou -0.19650793 -2.437919615
## el Besòs i el Maresme -2.82145232 0.784936756
## Provençals del Poblenou -1.72380430 -1.059592945
## Sant Martà de Provençals -1.09204118 0.148474534
## la Verneda i la Pau 0.35262141 -0.273294308
## PC5
## el Raval -0.081486501
## el Barri Gòtic -5.373189620
## la Barceloneta -2.429958889
## Sant Pere, Santa Caterina i la Ribera -2.274372142
## el Fort Pienc 0.055084007
## la Sagrada FamÃlia 0.183408957
## la Dreta de l'Eixample 1.603241314
## l'Antiga Esquerra de l'Eixample 0.402375957
## la Nova Esquerra de l'Eixample 0.670404140
## Sant Antoni 0.290964491
## el Poble Sec 1.108797003
## la Marina del Prat Vermell -0.705539299
## la Marina de Port 1.117701092
## la Font de la Guatlla -1.181363659
## Hostafrancs 0.174784759
## la Bordeta 0.009024956
## Sants - Badal 0.325534474
## Sants 1.234007506
## les Corts 1.238100642
## la Maternitat i Sant Ramon -1.182000202
## Pedralbes -2.769640806
## Vallvidrera, el Tibidabo i les Planes -3.499988478
## Sarrià -1.002506943
## les Tres Torres -1.961786825
## Sant Gervasi - la Bonanova 0.511221385
## Sant Gervasi - Galvany 0.503663793
## el Putxet i el Farró 0.876224175
## Vallcarca i els Penitents -0.290674812
## el Coll -0.498896696
## la Salut -0.795000948
## la Vila de Grà cia 1.583964561
## el Camp d'en Grassot i Grà cia Nova 0.028851358
## el Baix Guinardó 0.119664666
## Can Baró -0.881574802
## el Guinardó 1.052435719
## la Font d'en Fargues -0.399259593
## el Carmel 2.745576220
## la Teixonera 1.034466944
## Sant GenÃs dels Agudells 1.561471646
## Montbau 0.241832194
## la Vall d'Hebron -0.225878749
## la Clota 0.964936510
## Horta 1.246712613
## Vilapicina i la Torre Llobeta 0.360565444
## Porta 1.159736528
## el Turó de la Peira 0.440269046
## Can Peguera 0.144867840
## la Guineueta -0.320278477
## Canyelles -1.787801405
## les Roquetes 0.848066226
## Verdun -0.279159666
## la Prosperitat 0.287713394
## la Trinitat Nova -0.368544274
## Torre Baró -0.698007708
## Ciutat Meridiana 0.667116529
## Vallbona -2.000227365
## la Trinitat Vella -0.292187023
## Baró de Viver -1.108167422
## el Bon Pastor 0.959623875
## Sant Andreu 3.954251710
## la Sagrera 0.692201455
## el Congrés i els Indians -0.081314457
## Navas 0.368580710
## el Camp de l'Arpa del Clot 1.601366244
## el Clot 0.931099820
## el Parc i la Llacuna del Poblenou -1.379270263
## la Vila OlÃmpica del Poblenou -2.459994422
## el Poblenou 1.313403660
## Diagonal Mar i el Front MarÃtim del Poblenou -1.804631687
## el Besòs i el Maresme 1.508870309
## Provençals del Poblenou 0.283659964
## Sant Martà de Provençals 0.575432800
## la Verneda i la Pau 1.151426496
Var Exp
VE <- pca_result$sdev^2
PVE <- VE / sum(VE)
round(PVE, 2)
## [1] 0.25 0.15 0.10 0.09 0.08 0.06 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.00
vis 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
)