PCA SOLE + alq.pmq.V1519 nou.eurm2.V1519 usat.eurm2.V1519
Variables
merged2 <- merged %>%
remove_rownames %>%
column_to_rownames(var="Nom_Barri") %>%
select("n.tot","pc.esp","pc.ext","pc.ue27-esp","pc.20.34","2019-2014","var20192014",
"hotel2019",
"alq.mq","alq.num",#"alq.pmq",
"tot_ann","pc.ent","pc.priv","pc.shared","pc.hotel","pm_ent","pm_priv","pm_sha","pm_hot",
"RFD.2017",
"tot.comp","perc.nou.comp","perc.prot.comp","perc.usat.comp",
"tot.m2","nou.m2","prot.m2","usat.m2",
"tot.eur","nou.eur","usat.eur",
#"tot.eurm2","nou.eurm2","usat.eurm2",
"rest1614",
"alq.pmq.V1519", "nou.eurm2.V1519", "usat.eurm2.V1519",
"pmedio.V1519","pm_ent.V1519","pm_priv.V1519")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
## 22607.6438356 80.7810477 13.4049860 7.1143474
## pc.20.34 2019-2014 var20192014 hotel2019
## 18.9882789 658.3424658 5.2145625 276.0547945
## alq.mq alq.num tot_ann pc.ent
## 69.9929674 167.3698630 279.8356164 33.4375693
## pc.priv pc.shared pc.hotel pm_ent
## 60.9146924 1.7345421 1.1734702 160.8299728
## pm_priv pm_sha pm_hot RFD.2017
## 44.5267676 46.9704709 172.4968876 93.6712329
## tot.comp perc.nou.comp perc.prot.comp perc.usat.comp
## 185.1917808 7.6875336 2.2687579 91.9241064
## tot.m2 nou.m2 prot.m2 usat.m2
## 77.6799087 59.2428082 15.4760274 77.3446347
## tot.eur nou.eur usat.eur rest1614
## 309.5945205 260.8150685 302.6876712 1.1917808
## alq.pmq.V1519 nou.eurm2.V1519 usat.eurm2.V1519 pmedio.V1519
## 35.3976621 13.8070105 448.2268587 -87.3018000
## pm_ent.V1519 pm_priv.V1519
## 0.8279582 491.4042876pca_result$scale## n.tot pc.esp pc.ext pc.ue27-esp
## 14753.886904 9.009606 6.720753 7.320092
## pc.20.34 2019-2014 var20192014 hotel2019
## 3.873115 1490.085971 10.477120 490.230782
## alq.mq alq.num tot_ann pc.ent
## 14.165945 159.168980 443.819033 18.180666
## pc.priv pc.shared pc.hotel pm_ent
## 20.110450 6.114537 1.943317 201.699239
## pm_priv pm_sha pm_hot RFD.2017
## 16.003633 143.964167 403.002786 42.709775
## tot.comp perc.nou.comp perc.prot.comp perc.usat.comp
## 148.403509 10.382977 17.543273 10.563166
## tot.m2 nou.m2 prot.m2 usat.m2
## 18.930482 47.275532 32.108977 18.206886
## tot.eur nou.eur usat.eur rest1614
## 183.338837 440.142831 171.131367 26.158522
## alq.pmq.V1519 nou.eurm2.V1519 usat.eurm2.V1519 pmedio.V1519
## 15.633887 67.962390 3428.920989 28.320141
## pm_ent.V1519 pm_priv.V1519
## 57.451326 1916.511004#pesos asignados a cada componente
pca_result$rotation[,1:5]## PC1 PC2 PC3 PC4
## n.tot -0.2051495626 -0.125978537 0.029816416 -0.0892154232
## pc.esp 0.1104881449 0.306460717 -0.035313476 -0.0096166625
## pc.ext 0.0027815703 -0.312580572 0.051189615 0.0113166232
## pc.ue27-esp -0.1926948624 -0.207820723 0.010888160 0.0342473536
## pc.20.34 -0.1277782227 -0.309738115 -0.054978568 0.1002890221
## 2019-2014 -0.0293338728 0.154155191 -0.034895556 -0.0905019621
## var20192014 0.0471933008 0.171689013 -0.089417170 0.0914727421
## hotel2019 -0.2519723627 -0.163852725 -0.092009903 -0.0107172204
## alq.mq -0.2215140365 0.215085480 0.043523484 0.0343311639
## alq.num -0.2598796661 -0.134172770 -0.044685807 -0.0520924273
## tot_ann -0.2423682585 -0.248540431 -0.081586267 -0.0004533682
## pc.ent -0.2146327672 0.030318516 -0.197701487 0.0010862455
## pc.priv 0.1334119584 -0.028397307 0.175335845 -0.0549894251
## pc.shared 0.0931516258 -0.062303432 -0.071668860 -0.1093704062
## pc.hotel -0.2047397741 -0.102993998 -0.102490428 -0.0783620212
## pm_ent -0.0623895232 0.102705068 -0.331849405 0.3125930157
## pm_priv -0.2065038233 -0.008134245 -0.006627291 -0.1369364885
## pm_sha -0.1243445575 0.009742234 -0.100681778 -0.1250742622
## pm_hot -0.1778265891 -0.068076159 -0.146555507 -0.0634410023
## RFD.2017 -0.2544910829 0.224575941 0.052752250 -0.0381902471
## tot.comp -0.1533774479 -0.153707794 0.110997072 0.0204836405
## perc.nou.comp -0.0297092453 0.007829454 0.271321353 0.3699641784
## perc.prot.comp 0.0441826223 0.033950351 -0.341064039 0.3514029089
## perc.usat.comp 0.0242941300 -0.010246232 -0.225063728 -0.4199955994
## tot.m2 -0.2212498765 0.266952014 0.093299608 -0.0334073040
## nou.m2 -0.1773635551 0.004862997 0.265869598 0.1937137950
## prot.m2 -0.0840809946 -0.087984636 0.072275907 0.1567023976
## usat.m2 -0.2055563620 0.288382242 0.037808993 -0.0075760769
## tot.eur -0.2546826734 0.225795177 0.068787967 -0.0344627311
## nou.eur -0.1737899205 0.072359262 0.222551001 0.1573357596
## usat.eur -0.2537756281 0.233466786 0.043895184 -0.0522135821
## rest1614 0.0350809686 0.158815809 -0.070417957 -0.0440618874
## alq.pmq.V1519 -0.0930931865 -0.031863053 0.173600711 -0.1327706918
## nou.eurm2.V1519 -0.0564059525 -0.116022333 0.226348623 0.3364032415
## usat.eurm2.V1519 0.0007377175 0.042865596 -0.115207285 -0.0315944321
## pmedio.V1519 -0.1591793844 -0.144605215 -0.180658048 -0.0602476234
## pm_ent.V1519 -0.1490852196 0.051182289 -0.284250771 0.1008629889
## pm_priv.V1519 0.0259294101 0.048257568 -0.354480195 0.3530861440
## PC5
## n.tot 0.3238910668
## pc.esp 0.2397182282
## pc.ext -0.1696558068
## pc.ue27-esp -0.2932247087
## pc.20.34 -0.2329006112
## 2019-2014 0.3496079548
## var20192014 0.0614715245
## hotel2019 0.0842221897
## alq.mq -0.1102428574
## alq.num 0.2074380198
## tot_ann -0.0071229271
## pc.ent -0.1175703158
## pc.priv 0.1739112285
## pc.shared 0.0358639195
## pc.hotel 0.1987867592
## pm_ent 0.0867802244
## pm_priv 0.0418153925
## pm_sha -0.0061570432
## pm_hot 0.0606703253
## RFD.2017 -0.0770962174
## tot.comp 0.2661823650
## perc.nou.comp -0.0002000059
## perc.prot.comp 0.0422771727
## perc.usat.comp -0.0137326701
## tot.m2 -0.1147041079
## nou.m2 0.1199575977
## prot.m2 0.2305163799
## usat.m2 -0.1143445759
## tot.eur -0.1348865250
## nou.eur -0.0245895461
## usat.eur -0.1343300867
## rest1614 0.1973596462
## alq.pmq.V1519 0.2468700534
## nou.eurm2.V1519 0.0670834639
## usat.eurm2.V1519 -0.2317432097
## pmedio.V1519 0.0222115669
## pm_ent.V1519 0.1054465744
## pm_priv.V1519 0.0577348574#results
pca_result$x[,1:5]## PC1 PC2
## el Raval -3.339869715 -7.400224313
## el Barri Gòtic -4.782148841 -6.991777719
## la Barceloneta 0.186460822 -4.288199867
## Sant Pere, Santa Caterina i la Ribera -3.825988288 -4.927347617
## el Fort Pienc -3.846182653 -0.135045234
## la Sagrada FamÃlia -2.020426133 -2.045753148
## la Dreta de l'Eixample -9.756166690 -1.919615168
## l'Antiga Esquerra de l'Eixample -5.741525150 -1.410040432
## la Nova Esquerra de l'Eixample -2.676678401 -1.894451947
## Sant Antoni -4.006412940 -1.943567071
## el Poble Sec -2.138378915 -3.365605782
## la Marina del Prat Vermell 3.334918269 1.727819625
## la Marina de Port 1.860869702 -0.092799737
## la Font de la Guatlla 0.675046867 -0.418498314
## Hostafrancs -0.495846023 -0.129327292
## la Bordeta 2.154898405 0.044392725
## Sants - Badal 1.291870685 -0.257949367
## Sants -1.058219682 -0.877487181
## les Corts -2.980966959 -0.031895518
## la Maternitat i Sant Ramon 0.006605866 1.327876185
## Pedralbes -5.598167600 6.763507087
## Vallvidrera, el Tibidabo i les Planes -0.638737466 3.947209162
## Sarrià -2.663994782 4.258939035
## les Tres Torres -3.956368415 6.204168034
## Sant Gervasi - la Bonanova -3.760193226 4.710795837
## Sant Gervasi - Galvany -5.455427109 2.686776599
## el Putxet i el Farró -2.917096152 2.169607254
## Vallcarca i els Penitents 0.041522769 1.373031769
## el Coll 1.935065159 0.126013399
## la Salut 0.639261670 1.015815051
## la Vila de Grà cia -6.071601601 -2.501072406
## el Camp d'en Grassot i Grà cia Nova -1.951126864 -0.044309126
## el Baix Guinardó 0.745911983 0.002656247
## Can Baró 1.790448897 0.124239453
## el Guinardó 0.590833244 -0.256907561
## la Font d'en Fargues 1.667239518 2.438586478
## el Carmel 0.534448863 -1.195890550
## la Teixonera 2.264146107 -0.154211737
## Sant GenÃs dels Agudells 1.474447005 0.628361412
## Montbau 2.677194060 -0.087601277
## la Vall d'Hebron 2.465785149 0.782698944
## la Clota 2.891399280 1.571424227
## Horta 1.482155040 -0.025863686
## Vilapicina i la Torre Llobeta 2.654028480 0.200064459
## Porta 2.198200135 -0.205649319
## el Turó de la Peira 1.913221539 -1.426081185
## Can Peguera 5.867921278 -2.368635640
## la Guineueta 2.260611289 0.937888235
## Canyelles 3.467588233 1.197419046
## les Roquetes 2.974867479 -1.507653819
## Verdun 2.846750116 -1.248080458
## la Prosperitat 3.333986435 -0.857644487
## la Trinitat Nova 3.833573464 -1.641884676
## Torre Baró 3.414160057 0.481919329
## Ciutat Meridiana 3.758636842 -2.425311583
## Vallbona 3.759702300 0.939422668
## la Trinitat Vella 2.829649616 -1.417193793
## Baró de Viver 0.378985378 2.003979532
## el Bon Pastor 2.473716964 -0.080474357
## Sant Andreu 0.106585600 0.337422247
## la Sagrera 1.592324288 0.013771495
## el Congrés i els Indians 1.410734441 0.469814381
## Navas 0.739269247 0.666343460
## el Camp de l'Arpa del Clot -1.587138160 -0.996729109
## el Clot -0.303526954 -0.083651699
## el Parc i la Llacuna del Poblenou 0.736627325 0.150101356
## la Vila OlÃmpica del Poblenou -1.503460455 3.521087609
## el Poblenou -2.015995485 0.599529269
## Diagonal Mar i el Front MarÃtim del Poblenou -2.008076651 3.679901815
## el Besòs i el Maresme -0.151772667 -1.943069590
## Provençals del Poblenou 0.296752695 1.536941084
## Sant Martà de Provençals 1.495008012 0.378119758
## la Verneda i la Pau 2.198063402 -0.420142506
## PC3 PC4
## el Raval 0.01353081 0.25194005
## el Barri Gòtic 1.84221142 2.32130260
## la Barceloneta -0.58683734 0.15063739
## Sant Pere, Santa Caterina i la Ribera 0.14439605 0.58456889
## el Fort Pienc 2.48423480 1.99250916
## la Sagrada FamÃlia -2.37976571 -2.00618319
## la Dreta de l'Eixample -2.01404804 0.12031924
## l'Antiga Esquerra de l'Eixample -1.21023753 -0.31435283
## la Nova Esquerra de l'Eixample 0.02660752 -0.05775211
## Sant Antoni -2.76999932 -0.44213580
## el Poble Sec -0.58467953 -0.62350326
## la Marina del Prat Vermell -10.25494730 8.86627153
## la Marina de Port 0.74703938 -1.55101914
## la Font de la Guatlla -0.61786221 -0.57370032
## Hostafrancs -0.30179108 0.76740415
## la Bordeta -0.65108428 -1.65854837
## Sants - Badal -0.63778398 -1.36987869
## Sants -0.53261005 -0.96220634
## les Corts 0.82320344 0.08107355
## la Maternitat i Sant Ramon -1.19181034 -1.82782114
## Pedralbes 3.87431928 0.72875830
## Vallvidrera, el Tibidabo i les Planes -0.53553276 -1.05818632
## Sarrià 1.52571595 -0.78825269
## les Tres Torres -1.23025127 -2.49145432
## Sant Gervasi - la Bonanova -0.45091125 1.91360118
## Sant Gervasi - Galvany 1.68291181 0.52623694
## el Putxet i el Farró -0.27563455 1.47371667
## Vallcarca i els Penitents 1.26198720 0.79383579
## el Coll 1.34973733 1.38323438
## la Salut 0.75285391 -0.58364071
## la Vila de Grà cia -2.23704522 -1.12079758
## el Camp d'en Grassot i Grà cia Nova -0.80854230 -1.76468282
## el Baix Guinardó -0.93105585 -2.00095513
## Can Baró -0.62002733 -1.02733839
## el Guinardó 0.49925008 -0.51271773
## la Font d'en Fargues -0.99969326 0.87662278
## el Carmel 1.30590641 2.42138618
## la Teixonera 2.55734986 2.08249213
## Sant GenÃs dels Agudells 1.38469835 1.21997967
## Montbau 2.46000065 2.39573558
## la Vall d'Hebron 1.73725075 2.12696386
## la Clota 2.24360343 4.55734996
## Horta 2.34389360 1.07353779
## Vilapicina i la Torre Llobeta -0.34612986 -1.64041893
## Porta 0.24221942 -0.81808354
## el Turó de la Peira 0.55587343 -1.16555728
## Can Peguera -2.32943682 -1.75832568
## la Guineueta 0.54752747 -0.53456703
## Canyelles -0.50903653 -0.72687419
## les Roquetes -0.10706979 -1.13129337
## Verdun 0.77608026 -0.38931196
## la Prosperitat -0.09062517 -1.43566144
## la Trinitat Nova 0.16678643 -1.00703037
## Torre Baró -1.50700602 -0.12120713
## Ciutat Meridiana 0.45823523 -1.01890434
## Vallbona -0.44848915 -0.38292654
## la Trinitat Vella -0.22536543 -1.44075995
## Baró de Viver -3.45459908 -0.78426962
## el Bon Pastor 0.56804808 -0.03305654
## Sant Andreu 0.67085363 0.05766363
## la Sagrera 0.03519687 -1.48686443
## el Congrés i els Indians -0.15851423 -1.59408430
## Navas 1.54600791 -0.08412682
## el Camp de l'Arpa del Clot 0.69639451 1.15114278
## el Clot 1.17266241 -0.21681899
## el Parc i la Llacuna del Poblenou -0.99527193 -1.14035041
## la Vila OlÃmpica del Poblenou -1.93631759 -1.57839381
## el Poblenou -0.13473138 0.85018107
## Diagonal Mar i el Front MarÃtim del Poblenou -2.13502649 -2.16317605
## el Besòs i el Maresme 2.52738083 2.75368458
## Provençals del Poblenou 1.95965420 0.62535060
## Sant Martà de Provençals 1.87983061 1.18290673
## la Verneda i la Pau 1.33631666 0.05678239
## PC5
## el Raval -0.9840557750
## el Barri Gòtic -7.0097941326
## la Barceloneta -3.4521817603
## Sant Pere, Santa Caterina i la Ribera -1.9322029693
## el Fort Pienc 0.6753701099
## la Sagrada FamÃlia 0.3805792231
## la Dreta de l'Eixample 2.8930557364
## l'Antiga Esquerra de l'Eixample 0.5859033313
## la Nova Esquerra de l'Eixample 1.7820531224
## Sant Antoni 1.7722108674
## el Poble Sec -0.4487344424
## la Marina del Prat Vermell 0.8164007740
## la Marina de Port 1.4015166819
## la Font de la Guatlla -0.7659555211
## Hostafrancs -0.5134902755
## la Bordeta 0.4402107399
## Sants - Badal 0.3143058201
## Sants 1.7418918363
## les Corts 2.7708467273
## la Maternitat i Sant Ramon -0.0126297805
## Pedralbes -2.6122986915
## Vallvidrera, el Tibidabo i les Planes -2.5848633410
## Sarrià -0.5475192774
## les Tres Torres -2.2285610459
## Sant Gervasi - la Bonanova 0.3644251449
## Sant Gervasi - Galvany 0.2497080571
## el Putxet i el Farró 0.0702529843
## Vallcarca i els Penitents -0.6140337071
## el Coll -0.2078989061
## la Salut -0.4428799710
## la Vila de Grà cia 0.8238485463
## el Camp d'en Grassot i Grà cia Nova 0.8452990728
## el Baix Guinardó -0.1241032348
## Can Baró -0.7664789315
## el Guinardó 1.0433808144
## la Font d'en Fargues 0.0213758989
## el Carmel 1.6091199388
## la Teixonera 0.2741543084
## Sant GenÃs dels Agudells -0.3304093147
## Montbau 0.5384495228
## la Vall d'Hebron 0.0004907381
## la Clota -1.7395122408
## Horta 1.5078492388
## Vilapicina i la Torre Llobeta 0.3136392754
## Porta 1.4861992828
## el Turó de la Peira 1.2746966055
## Can Peguera 0.3516617663
## la Guineueta 0.1879712379
## Canyelles -0.7856683774
## les Roquetes -0.4020516835
## Verdun -0.0094299205
## la Prosperitat 0.5403632899
## la Trinitat Nova -0.2823559934
## Torre Baró -0.7337970239
## Ciutat Meridiana -0.5609108412
## Vallbona -1.9504194326
## la Trinitat Vella -1.2195255477
## Baró de Viver -5.0632574518
## el Bon Pastor 1.0924618523
## Sant Andreu 3.2333957737
## la Sagrera 0.6971283653
## el Congrés i els Indians 0.0061138564
## Navas 0.1492374592
## el Camp de l'Arpa del Clot 1.5371581199
## el Clot 1.3776611546
## el Parc i la Llacuna del Poblenou 0.0411740079
## la Vila OlÃmpica del Poblenou -1.5651374921
## el Poblenou 1.4661085081
## Diagonal Mar i el Front MarÃtim del Poblenou 0.3715648903
## el Besòs i el Maresme 1.1686047288
## Provençals del Poblenou 0.4961278875
## Sant Martà de Provençals 0.3239457658
## la Verneda i la Pau 0.8522440197Var Exp
VE <- pca_result$sdev^2
PVE <- VE / sum(VE)
round(PVE, 2)## [1] 0.23 0.15 0.09 0.08 0.07 0.06 0.04 0.04 0.03 0.03 0.02 0.02 0.02 0.02 0.02
## [16] 0.01 0.01 0.01 0.01 0.01 0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
## [31] 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
)