PCA
Variables
merged2 <- merged %>%
remove_rownames %>%
column_to_rownames(var="Nom_Barri") %>%
select("n.tot","pc.esp","pc.ext","pc.ue27-esp","pc.20.34","var20192014",
"hotel2019",
"alq.mq","alq.num","alq.pmq",
"tot_ann","pm_priv",
"RFD.2017",
"tot.comp",
"tot.m2","nou.m2","usat.m2",
"tot.eur","nou.eur","usat.eur",
"tot.eurm2","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.ext pc.ue27-esp pc.20.34 var20192014
## 22607.643836 80.781048 13.404986 7.114347 18.988279 5.214562
## hotel2019 alq.mq alq.num alq.pmq tot_ann pm_priv
## 276.054795 69.992967 167.369863 12.531159 279.835616 44.526768
## RFD.2017 tot.comp tot.m2 nou.m2 usat.m2 tot.eur
## 93.671233 185.191781 77.679909 59.242808 77.344635 299.213356
## nou.eur usat.eur tot.eurm2 nou.eurm2 usat.eurm2
## 256.620662 289.042123 3592.717123 2409.018721 3335.016172pca_result$scale## n.tot pc.esp pc.ext pc.ue27-esp pc.20.34 var20192014
## 14753.886904 9.009606 6.720753 7.320092 3.873115 10.477120
## hotel2019 alq.mq alq.num alq.pmq tot_ann pm_priv
## 490.230782 14.165945 159.168980 4.260785 443.819033 16.003633
## RFD.2017 tot.comp tot.m2 nou.m2 usat.m2 tot.eur
## 42.709775 148.403509 18.930482 47.275532 18.206886 187.686169
## nou.eur usat.eur tot.eurm2 nou.eurm2 usat.eurm2
## 538.160472 175.932336 1372.568549 3007.385279 1505.311684#pesos asignados a cada componente
pca_result$rotation[,1:5]## PC1 PC2 PC3 PC4 PC5
## n.tot -0.17838738 -0.169021341 0.484169158 0.084205999 -0.031903240
## pc.esp 0.08233687 0.362767661 0.328285661 0.105742874 0.005352715
## pc.ext 0.02663318 -0.350202283 -0.285095180 -0.012472212 -0.010525100
## pc.ue27-esp -0.17197728 -0.264862348 -0.304640158 -0.182013449 0.024671473
## pc.20.34 -0.08719817 -0.370724288 -0.240287630 -0.102202130 -0.065737727
## var20192014 0.06671246 0.176340666 -0.077276983 -0.029330665 -0.530222775
## hotel2019 -0.18723776 -0.213358021 0.209002898 -0.093136111 -0.386849077
## alq.mq -0.22910657 0.176560882 -0.097855171 -0.098625874 -0.199811017
## alq.num -0.21161636 -0.183782072 0.362578153 -0.008583847 -0.193079137
## alq.pmq -0.23497673 -0.013902941 0.095793125 -0.160491692 0.471616566
## tot_ann -0.17887933 -0.308012681 0.072343067 -0.135396364 -0.284719017
## pm_priv -0.21634936 -0.042382924 0.066269454 -0.276404247 0.221067207
## RFD.2017 -0.26686666 0.192546635 -0.062548428 -0.139621732 -0.105023289
## tot.comp -0.13682506 -0.190442616 0.338597331 0.214280471 0.054566833
## tot.m2 -0.23823558 0.244825879 -0.089336314 -0.070422563 -0.134830464
## nou.m2 -0.19542927 -0.023893791 -0.059424805 0.438263309 0.033519501
## usat.m2 -0.22069617 0.271153699 -0.080149642 -0.101744084 -0.179044558
## tot.eur -0.28728223 0.170592147 -0.092035995 -0.033583853 0.012489313
## nou.eur -0.18911195 0.020882204 -0.222492402 0.494979207 -0.051756363
## usat.eur -0.28502491 0.169877547 -0.033910066 -0.122241141 0.065621841
## tot.eurm2 -0.29317610 0.052574443 -0.008499432 -0.104323640 0.162757980
## nou.eurm2 -0.20814255 -0.051151091 -0.139152997 0.488534978 -0.030566276
## usat.eurm2 -0.30051144 0.004282554 -0.037892608 0.121459150 0.183384384#results
pca_result$x[,1:5]## PC1 PC2
## el Raval -2.39259310 -7.561672500
## el Barri Gòtic -4.97223394 -6.902052471
## la Barceloneta -0.64400131 -4.413215022
## Sant Pere, Santa Caterina i la Ribera -3.16587204 -5.080370682
## el Fort Pienc -5.34489941 -0.664233088
## la Sagrada FamÃlia -1.03994823 -2.341868997
## la Dreta de l'Eixample -7.21823043 -2.506783610
## l'Antiga Esquerra de l'Eixample -4.54155376 -1.727786740
## la Nova Esquerra de l'Eixample -2.49859533 -1.863289099
## Sant Antoni -1.60764603 -1.668119816
## el Poble Sec -1.46929077 -3.543187403
## la Marina del Prat Vermell 4.45526628 0.687051288
## la Marina de Port 1.53108120 0.041761048
## la Font de la Guatlla 0.71343344 -0.024385023
## Hostafrancs -0.53510405 -0.772028072
## la Bordeta 2.00861149 0.190859599
## Sants - Badal 1.06723475 -0.540259364
## Sants -0.82596213 -1.310020367
## les Corts -3.07744281 0.218693340
## la Maternitat i Sant Ramon -0.46171170 1.210940047
## Pedralbes -9.50782413 5.384602829
## Vallvidrera, el Tibidabo i les Planes -2.17107773 3.840419133
## Sarrià -3.97207109 3.550209644
## les Tres Torres -4.51362143 5.397679930
## Sant Gervasi - la Bonanova -4.62785235 3.403763788
## Sant Gervasi - Galvany -6.05139661 2.209670237
## el Putxet i el Farró -3.06771775 1.434703187
## Vallcarca i els Penitents -0.87697698 1.353616640
## el Coll 1.75907274 0.542962643
## la Salut -0.59317185 1.161268867
## la Vila de Grà cia -4.07639428 -2.335485579
## el Camp d'en Grassot i Grà cia Nova -1.63114027 -0.009170747
## el Baix Guinardó 0.63560527 0.067612959
## Can Baró 2.10325745 0.423133141
## el Guinardó -0.07207559 -0.308426480
## la Font d'en Fargues 1.28312720 2.524930747
## el Carmel 1.21541878 -0.875397235
## la Teixonera 1.92522314 0.348352199
## Sant GenÃs dels Agudells 1.50576926 0.907579623
## Montbau 2.27536388 0.785211908
## la Vall d'Hebron 1.73329120 1.327096110
## la Clota 4.15775901 1.806233177
## Horta 0.75382533 0.488554912
## Vilapicina i la Torre Llobeta 2.05161014 0.372908340
## Porta 2.21396964 -0.263139117
## el Turó de la Peira 2.48786674 -0.747452221
## Can Peguera 6.92051551 -0.504780475
## la Guineueta 1.33219444 1.326715141
## Canyelles 3.26189238 1.856365710
## les Roquetes 2.94287145 -1.181933303
## Verdun 2.34199740 -0.839042244
## la Prosperitat 3.01911021 -0.252934738
## la Trinitat Nova 4.03368217 -1.076499492
## Torre Baró 5.44700726 0.364842845
## Ciutat Meridiana 3.94192381 -2.000711835
## Vallbona 4.43005315 1.537562118
## la Trinitat Vella 2.76586120 -1.414912418
## Baró de Viver 5.14415072 1.098254991
## el Bon Pastor 2.55400888 0.278298998
## Sant Andreu -0.27936379 -0.126029631
## la Sagrera 1.27974918 0.280890108
## el Congrés i els Indians 1.01168622 0.487101780
## Navas 0.02842782 0.901402620
## el Camp de l'Arpa del Clot -1.34159161 -1.058459923
## el Clot -1.15212941 0.066672567
## el Parc i la Llacuna del Poblenou 0.72589484 -0.205039093
## la Vila OlÃmpica del Poblenou -2.07476675 2.723739877
## el Poblenou -1.70134639 0.269785514
## Diagonal Mar i el Front MarÃtim del Poblenou -1.53232158 3.201768106
## el Besòs i el Maresme 0.05035407 -2.000097132
## Provençals del Poblenou -0.08297846 1.261384619
## Sant Martà de Provençals 0.40184609 0.678073385
## la Verneda i la Pau 1.61088939 0.106110208
## PC3 PC4
## el Raval -0.81445810 -0.25798285
## el Barri Gòtic -4.48712474 -1.69008706
## la Barceloneta -2.43328609 -0.83309954
## Sant Pere, Santa Caterina i la Ribera -2.36568515 -0.81421896
## el Fort Pienc -2.07982799 5.58748928
## la Sagrada FamÃlia 2.05336549 -1.94549855
## la Dreta de l'Eixample 2.30908737 -1.14916994
## l'Antiga Esquerra de l'Eixample 1.09822391 -0.21367912
## la Nova Esquerra de l'Eixample 2.29617932 -0.09127318
## Sant Antoni 1.21533887 -1.37385002
## el Poble Sec 0.58470301 -0.58355855
## la Marina del Prat Vermell -0.79617463 -0.06196336
## la Marina de Port 0.75142384 -0.51034610
## la Font de la Guatlla -0.87488332 -0.61508994
## Hostafrancs -1.17112908 -0.06156162
## la Bordeta 0.43605950 -0.79204994
## Sants - Badal 0.41875959 -0.29444780
## Sants 1.21112054 0.32275949
## les Corts 2.05802511 1.50012507
## la Maternitat i Sant Ramon 0.85445450 -1.77606227
## Pedralbes -4.15275112 3.49018401
## Vallvidrera, el Tibidabo i les Planes -1.16083793 -3.14648611
## Sarrià 0.08401267 -1.38647459
## les Tres Torres -0.59616040 -3.20001374
## Sant Gervasi - la Bonanova -0.33709164 0.17986796
## Sant Gervasi - Galvany 1.32096458 0.27407393
## el Putxet i el Farró 0.50933577 0.66838631
## Vallcarca i els Penitents -0.28796230 0.07030917
## el Coll -0.65691906 0.35237852
## la Salut -0.55912001 0.65169521
## la Vila de Grà cia 2.18430941 -0.06766014
## el Camp d'en Grassot i Grà cia Nova 1.48687991 -0.81031908
## el Baix Guinardó 1.21457669 -0.83225380
## Can Baró -0.28853998 -0.89016060
## el Guinardó 1.20631276 0.91461725
## la Font d'en Fargues 0.20949113 0.01814993
## el Carmel 1.30275674 1.23637643
## la Teixonera -0.34236654 0.60992154
## Sant GenÃs dels Agudells -0.80621093 0.25057977
## Montbau -0.50242887 0.65772518
## la Vall d'Hebron -0.06572905 0.71435884
## la Clota -1.78496535 1.05606892
## Horta 0.79999371 1.10901522
## Vilapicina i la Torre Llobeta 0.72081753 -0.57408209
## Porta 0.28202714 0.25689192
## el Turó de la Peira -0.37373518 0.12167755
## Can Peguera 0.30425842 0.74891289
## la Guineueta 0.23432940 1.12610746
## Canyelles 0.51841612 -0.23754678
## les Roquetes -0.29332269 -0.40518207
## Verdun -0.63888337 1.26078654
## la Prosperitat 0.86251260 -0.09275148
## la Trinitat Nova -0.66763066 -0.28382975
## Torre Baró -1.24162738 -0.29021315
## Ciutat Meridiana -1.33587906 -0.33938677
## Vallbona -1.01853293 0.40782224
## la Trinitat Vella -1.25077435 -0.42107747
## Baró de Viver -1.00458490 0.30577568
## el Bon Pastor -0.05724731 0.32681038
## Sant Andreu 3.47436077 2.13691476
## la Sagrera 1.30006782 -0.51153132
## el Congrés i els Indians 0.31657059 -0.80661789
## Navas 0.31484049 -0.10754766
## el Camp de l'Arpa del Clot 1.45626185 1.49521966
## el Clot 1.01894195 0.97396775
## el Parc i la Llacuna del Poblenou -0.83772773 -0.93554047
## la Vila OlÃmpica del Poblenou -1.12009189 -2.78121279
## el Poblenou 0.29912189 0.78267718
## Diagonal Mar i el Front MarÃtim del Poblenou -1.01661907 -2.53667939
## el Besòs i el Maresme -0.61054729 1.79738340
## Provençals del Poblenou -0.31958309 -0.21072619
## Sant Martà de Provençals 0.46179448 1.25151950
## la Verneda i la Pau 1.18074371 1.27468317
## PC5
## el Raval -1.11818860
## el Barri Gòtic 0.76420779
## la Barceloneta 2.20690481
## Sant Pere, Santa Caterina i la Ribera -0.07168142
## el Fort Pienc -0.24666067
## la Sagrada FamÃlia -0.80461206
## la Dreta de l'Eixample -3.42031840
## l'Antiga Esquerra de l'Eixample -1.83913335
## la Nova Esquerra de l'Eixample -0.91663594
## Sant Antoni -0.91143966
## el Poble Sec -0.42066628
## la Marina del Prat Vermell -1.78835855
## la Marina de Port 0.39462142
## la Font de la Guatlla 0.54576000
## Hostafrancs 0.60418423
## la Bordeta 0.53144510
## Sants - Badal 0.70957342
## Sants 0.08770557
## les Corts 0.57557978
## la Maternitat i Sant Ramon 0.76377941
## Pedralbes -0.09854298
## Vallvidrera, el Tibidabo i les Planes 0.15483725
## Sarrià 0.31560890
## les Tres Torres -0.10298797
## Sant Gervasi - la Bonanova -0.20825116
## Sant Gervasi - Galvany -0.96326184
## el Putxet i el Farró -0.04298803
## Vallcarca i els Penitents 0.58019598
## el Coll 0.50151186
## la Salut 0.78301304
## la Vila de Grà cia -0.90960959
## el Camp d'en Grassot i Grà cia Nova 0.83404419
## el Baix Guinardó 0.79604500
## Can Baró 0.70251733
## el Guinardó 0.27953816
## la Font d'en Fargues 0.48532502
## el Carmel 1.02551508
## la Teixonera 0.57836662
## Sant GenÃs dels Agudells 0.04502248
## Montbau 0.47674622
## la Vall d'Hebron 1.07893800
## la Clota -4.52244104
## Horta 0.54140839
## Vilapicina i la Torre Llobeta 0.57150089
## Porta 0.03566262
## el Turó de la Peira 0.51142657
## Can Peguera -0.38875355
## la Guineueta 0.79751480
## Canyelles 0.51667870
## les Roquetes 0.72602327
## Verdun 0.74098291
## la Prosperitat 0.55541306
## la Trinitat Nova 0.81752543
## Torre Baró -3.18670975
## Ciutat Meridiana 0.25333430
## Vallbona -2.35428139
## la Trinitat Vella 0.09292959
## Baró de Viver -2.67372591
## el Bon Pastor 0.09551122
## Sant Andreu -0.04173136
## la Sagrera 0.57436522
## el Congrés i els Indians 0.73858845
## Navas 0.49290762
## el Camp de l'Arpa del Clot 0.22621535
## el Clot 0.91223398
## el Parc i la Llacuna del Poblenou -0.09351505
## la Vila OlÃmpica del Poblenou 0.72151828
## el Poblenou 0.17757658
## Diagonal Mar i el Front MarÃtim del Poblenou -0.71800476
## el Besòs i el Maresme 0.30288416
## Provençals del Poblenou 0.38839391
## Sant Martà de Provençals 0.73371800
## la Verneda i la Pau 0.49719933Var Exp
VE <- pca_result$sdev^2
PVE <- VE / sum(VE)
round(PVE, 2)## [1] 0.42 0.22 0.08 0.08 0.05 0.03 0.03 0.02 0.02 0.01 0.01 0.01 0.01 0.00 0.00
## [16] 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
)