setwd("~/Library/CloudStorage/GoogleDrive-icarounam@gmail.com/Mi unidad/Agrosavia/Ingreso/Evas")
dat<-read.table("eva2.csv", header=T, sep = ",")
library(FactoMineR)
## Warning: package 'FactoMineR' was built under R version 4.1.2
library (factoextra)
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 4.1.2
## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa
dat.CA<-dat[c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12",
"13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
"25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36",
"37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48",
"49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60",
"61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72",
"73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84",
"85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96",
"97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107",
"108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118",
"119", "120", "121", "122", "123", "124", "125", "126", "127") ,
c("Bovinos", "Bufalos", "Caprinos", "Ovinos", "Aguacate", "Arveja", "Avena",
"Cacao", "Cebolla.De.Bulbo", "Cebolla.De.Rama", "Durazno", "Frijol",
"Haba", "Habichuela", "Lechuga", "Limón", "Maíz", "Mora", "Papa", "Piña",
"Plátano", "Remolacha", "Repollo", "Tabaco", "Tomate", "Trigo", "Yuca",
"Zanahoria", "Aptitud")]
dat.CA1<-dat[c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12",
"13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
"25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36",
"37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48",
"49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60",
"61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72",
"73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84",
"85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96",
"97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107",
"108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118",
"119", "120", "121", "122", "123", "124", "125", "126", "127") ,
c("Aguacate", "Arveja", "Avena","Cacao", "Cebolla.De.Bulbo", "Durazno", "Frijol",
"Haba", "Habichuela", "Lechuga", "Limón", "Maíz", "Mora", "Papa", "Piña",
"Plátano", "Remolacha", "Repollo", "Tabaco", "Tomate", "Trigo", "Yuca",
"Zanahoria", "Aptitud")]
dat.CA2<-dat[c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12",
"13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24",
"25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36",
"37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48",
"49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60",
"61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72",
"73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84",
"85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96",
"97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107",
"108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118",
"119", "120", "121", "122", "123", "124", "125", "126", "127") ,
c("Bovinos", "Bufalos", "Caprinos", "Ovinos")]
## Ganado y Plantas
row.names(dat.CA)<-dat$Municipio
res<-CA(dat.CA, ncp=3, row.sup=NULL, col.sup=NULL, graph = T)
get_eigenvalue(res)
## eigenvalue variance.percent cumulative.variance.percent
## Dim.1 1.855562e-01 15.506004044 15.50600
## Dim.2 1.686353e-01 14.092012090 29.59802
## Dim.3 1.523572e-01 12.731735508 42.32975
## Dim.4 1.383623e-01 11.562245075 53.89200
## Dim.5 1.080036e-01 9.025323149 62.91732
## Dim.6 9.946163e-02 8.311513428 71.22883
## Dim.7 9.110032e-02 7.612800558 78.84163
## Dim.8 5.020210e-02 4.195139278 83.03677
## Dim.9 3.335163e-02 2.787029678 85.82380
## Dim.10 3.230339e-02 2.699433562 88.52324
## Dim.11 2.244435e-02 1.875562926 90.39880
## Dim.12 2.095002e-02 1.750688569 92.14949
## Dim.13 1.962538e-02 1.639995295 93.78948
## Dim.14 1.572847e-02 1.314350300 95.10383
## Dim.15 1.211723e-02 1.012576205 96.11641
## Dim.16 1.110707e-02 0.928162522 97.04457
## Dim.17 7.030524e-03 0.587505856 97.63208
## Dim.18 5.779320e-03 0.482948957 98.11503
## Dim.19 4.420220e-03 0.369375776 98.48440
## Dim.20 4.153166e-03 0.347059374 98.83146
## Dim.21 3.227675e-03 0.269720701 99.10118
## Dim.22 3.003822e-03 0.251014478 99.35220
## Dim.23 2.438931e-03 0.203809289 99.55601
## Dim.24 2.343440e-03 0.195829644 99.75184
## Dim.25 2.009364e-03 0.167912563 99.91975
## Dim.26 6.235448e-04 0.052106531 99.97186
## Dim.27 3.164912e-04 0.026447597 99.99830
## Dim.28 2.030812e-05 0.001697049 100.00000
fviz_eig(res, addlabels = TRUE, ylim = c(0, 20))
get_ca_row(res)
## Correspondence Analysis - Results for rows
## ===================================================
## Name Description
## 1 "$coord" "Coordinates for the rows"
## 2 "$cos2" "Cos2 for the rows"
## 3 "$contrib" "contributions of the rows"
## 4 "$inertia" "Inertia of the rows"
get_ca_col(res)
## Correspondence Analysis - Results for columns
## ===================================================
## Name Description
## 1 "$coord" "Coordinates for the columns"
## 2 "$cos2" "Cos2 for the columns"
## 3 "$contrib" "contributions of the columns"
## 4 "$inertia" "Inertia of the columns"
fviz_ca_row(res, repel=T)
## Warning: ggrepel: 119 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
fviz_ca_col(res, repel =T)
## Warning: ggrepel: 17 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
fviz_ca_biplot(res, repel = T)
## Warning: ggrepel: 122 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
## Warning: ggrepel: 17 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
res.hcpc<-HCPC(res ,nb.clust=-1,consol=TRUE,min=3,max=10,graph=TRUE)
res.hcpc$data.clust[,ncol(res.hcpc$data.clust),drop=F]
## clust
## PUERTO WILCHES 2
## BARRANCABERMEJA 2
## SILOS ns 2
## SABANA DE TORRES 2
## TONA 1
## MUTISCUA ns 2
## CHITAGA ns 2
## CIMITARRA 2
## BOLIVAR 2
## PUERTO PARRA 2
## PAMPLONA ns 2
## SIMACOTA 2
## VELEZ 2
## GAMBITA 2
## ALBANIA 2
## OIBA 2
## LANDAZURI 2
## GUAVATA 2
## LA ESPERANZA ns 2
## HERRAN ns 2
## LABATECA ns 2
## EL PEÑON 2
## GRAMALOTE ns 2
## TOLEDO ns 2
## GUADALUPE 2
## LOURDES ns 2
## GUAPOTA 2
## CHARTA 2
## CHARALA 2
## PUENTE NACIONAL 2
## SANTIAGO ns 2
## SALAZAR ns 2
## SANTA HELENA DEL OPON 2
## SUCRE 2
## CACHIRA ns 2
## CHINACOTA ns 2
## BOCHALEMA ns 2
## CUCUTA ns 2
## PALMAS DEL SOCORRO 2
## DURANIA ns 2
## CHIPATA 2
## EL GUACAMAYO 2
## VALLE DE SAN JOSE 2
## LA BELLEZA 2
## RAGONVALIA ns 2
## CUCUTILLA ns 2
## PALMAR 2
## AGUADA 2
## EL ZULIA ns 2
## JESUS MARIA 2
## EL TARRA ns 2
## CONTRATACION 2
## FLORIAN 2
## PAMPLONITA ns 2
## SARDINATA ns 2
## HACARI ns 2
## COROMORO 2
## ARBOLEDAS ns 2
## SAN ANDRES 2
## CHIMA 2
## RIONEGRO 2
## GUACA 2
## CABRERA 2
## VETAS 2
## VILLA CARO ns 2
## SAN CAYETANO ns 2
## ENCINO 2
## BARBOSA 2
## BETULIA 2
## TIBU ns 2
## CACOTA ns 2
## VILLA DEL ROSARIO ns 2
## CALIFORNIA 2
## SANTA BARBARA 2
## ONZAGA 2
## MATANZA 2
## SAN VICENTE DE CHUCURI 3
## CONVENCION ns 2
## CONFINES 2
## OCAMONTE 2
## SURATA 2
## PUERTO SANTANDER ns 2
## MALAGA 2
## HATO 2
## GALAN 2
## OCAÑA ns 2
## BUCARASICA ns 3
## EL CARMEN ns 2
## CERRITO 2
## CARCASI 2
## EL CARMEN DE CHUCURI 3
## SOCORRO 2
## CONCEPCION 2
## PIEDECUESTA 2
## ABREGO ns 2
## ENCISO 2
## SUAITA 2
## EL PLAYON 3
## MOGOTES 2
## LA PAZ 3
## SAN MIGUEL 2
## PINCHOTE 2
## LA PLAYA ns 2
## ZAPATOCA 2
## TEORAMA ns 3
## SAN CALIXTO ns 2
## GÜEPSA 2
## SAN JOSE DE MIRANDA 2
## FLORIDABLANCA 2
## PARAMO 2
## SAN JOAQUIN 2
## MACARAVITA 2
## SAN GIL 2
## BARICHARA 2
## CURITI 2
## LOS PATIOS ns 2
## MOLAGAVITA 2
## LOS SANTOS 2
## BUCARAMANGA 2
## ARATOCA 2
## VILLANUEVA 2
## CAPITANEJO 2
## JORDAN 2
## CEPITA 2
## LEBRIJA 4
## SAN BENITO 4
## GIRON 4
res.hcpc$desc.var
## $`1`
## Intern % glob % Intern freq Glob freq p.value
## Cebolla.De.Rama 20.848970061 1.448864e-01 2500 3632 0.000000e+00
## Papa 2.668668168 3.011020e-01 320 7548 8.681717e-186
## Aguacate 2.518555583 3.090005e-01 302 7746 2.105149e-165
## Ovinos 5.170544575 2.371996e+00 620 59461 1.430671e-68
## Arveja 0.116754232 5.588816e-02 14 1401 1.781457e-02
## Avena 0.000000000 7.978324e-05 0 2 1.908786e-02
## Tabaco 0.000000000 4.188620e-02 0 1050 1.300353e-02
## Aptitud 0.000000000 4.635407e-02 0 1162 7.598542e-03
## Durazno 0.000000000 5.026344e-02 0 1260 4.748508e-03
## Caprinos 1.067467267 1.439529e+00 128 36086 4.053296e-04
## Frijol 0.216829289 4.877150e-01 26 12226 2.795398e-06
## Cebolla.De.Bulbo 0.008339588 1.530243e-01 1 3836 3.945653e-07
## Yuca 0.116754232 5.585625e-01 14 14002 6.920350e-15
## Limón 0.000000000 3.063677e-01 0 7680 1.921613e-16
## Maíz 0.116754232 6.793144e-01 14 17029 4.744649e-20
## Piña 0.000000000 3.854727e-01 0 9663 1.380122e-20
## Plátano 0.133433408 9.301530e-01 16 23317 1.299261e-29
## Bufalos 0.016679176 2.812599e+00 2 70506 1.435472e-144
## Cacao 0.000000000 2.800232e+00 0 70196 1.083456e-148
## Bovinos 66.666666667 8.573970e+01 7994 2149316 0.000000e+00
## v.test
## Cebolla.De.Rama Inf
## Papa 29.069673
## Aguacate 27.409638
## Ovinos 17.500080
## Arveja 2.369450
## Avena -2.343811
## Tabaco -2.483673
## Aptitud -2.669407
## Durazno -2.823617
## Caprinos -3.536589
## Frijol -4.685308
## Cebolla.De.Bulbo -5.071561
## Yuca -7.785928
## Limón -8.226875
## Maíz -9.169671
## Piña -9.301858
## Plátano -11.300869
## Bufalos -25.602373
## Cacao -25.969915
## Bovinos -Inf
##
## $`2`
## Intern % glob % Intern freq Glob freq p.value
## Bufalos 3.060523000 2.812598732 69183 70506 0.000000e+00
## Bovinos 87.627140238 85.739702377 1980808 2149316 0.000000e+00
## Frijol 0.529220136 0.487714976 11963 12226 1.236137e-253
## Caprinos 1.513030768 1.439529087 34202 36086 3.092757e-227
## Maíz 0.723558016 0.679314439 16356 17029 5.049146e-185
## Cebolla.De.Bulbo 0.167131461 0.153024264 3778 3836 2.762948e-99
## Papa 0.319752833 0.301101966 7228 7548 7.244543e-74
## Mora 0.173413269 0.160922805 3920 4034 2.329258e-67
## Durazno 0.055651503 0.050263444 1258 1260 4.724916e-53
## Arveja 0.061358215 0.055888163 1387 1401 1.013035e-43
## Zanahoria 0.016146899 0.014560442 365 365 8.053090e-17
## Lechuga 0.004910427 0.004427970 111 111 2.067278e-05
## Repollo 0.004821951 0.004348187 109 109 2.542331e-05
## Remolacha 0.002344618 0.002114256 53 53 8.327842e-03
## Trigo 0.002167666 0.001954689 49 49 1.259489e-02
## Yuca 0.554878221 0.558562497 12543 14002 1.927063e-02
## Aptitud 0.043309079 0.046354065 979 1162 3.680688e-10
## Plátano 0.899537093 0.930152960 20334 23317 7.812187e-49
## Tabaco 0.033797892 0.041886204 764 1050 2.497337e-57
## Piña 0.065295404 0.385472748 1476 9663 0.000000e+00
## Limón 0.117319385 0.306367660 2652 7680 0.000000e+00
## Cebolla.De.Rama 0.050077505 0.144886373 1132 3632 0.000000e+00
## Cacao 1.209513310 2.800232329 27341 70196 0.000000e+00
## Aguacate 0.197169117 0.309000507 4457 7746 0.000000e+00
## v.test
## Bufalos Inf
## Bovinos Inf
## Frijol 34.017305
## Caprinos 32.182591
## Maíz 29.009117
## Cebolla.De.Bulbo 21.149940
## Papa 18.181411
## Mora 17.340445
## Durazno 15.331292
## Arveja 13.866349
## Zanahoria 8.330457
## Lechuga 4.257498
## Repollo 4.211009
## Remolacha 2.638481
## Trigo 2.495023
## Yuca -2.340254
## Aptitud -6.267001
## Plátano -14.686977
## Tabaco -15.958225
## Piña -Inf
## Limón -Inf
## Cebolla.De.Rama -Inf
## Cacao -Inf
## Aguacate -Inf
##
## $`3`
## Intern % glob % Intern freq Glob freq p.value
## Cacao 23.495789916 2.800232329 41354 70196 0.000000e+00
## Aguacate 1.583468745 0.309000507 2787 7746 0.000000e+00
## Plátano 1.625512767 0.930152960 2861 23317 1.403741e-181
## Remolacha 0.000000000 0.002114256 0 53 4.220553e-02
## Repollo 0.000000000 0.004348187 0 109 7.158466e-04
## Lechuga 0.000000000 0.004427970 0 111 6.188499e-04
## Zanahoria 0.000000000 0.014560442 0 365 5.759140e-12
## Tomate 0.085792530 0.160364322 151 4020 1.378955e-18
## Habichuela 0.001136325 0.035782785 2 897 2.039822e-25
## Yuca 0.372146404 0.558562497 655 14002 1.130968e-30
## Tabaco 0.000000000 0.041886204 0 1050 1.250630e-33
## Durazno 0.001136325 0.050263444 2 1260 1.317135e-36
## Mora 0.056816245 0.160922805 100 4034 5.802065e-38
## Arveja 0.000000000 0.055888163 0 1401 9.870069e-45
## Limón 0.134086338 0.306367660 236 7680 1.199862e-51
## Cebolla.De.Bulbo 0.032385260 0.153024264 57 3836 3.962270e-58
## Cebolla.De.Rama 0.000000000 0.144886373 0 3632 2.430707e-115
## Frijol 0.130677363 0.487714976 230 12226 3.595999e-151
## Piña 0.048293808 0.385472748 85 9663 1.374130e-191
## Papa 0.000000000 0.301101966 0 7548 1.953545e-239
## Maíz 0.126700226 0.679314439 223 17029 1.145042e-275
## Caprinos 0.560776337 1.439529087 987 36086 3.582929e-288
## Bufalos 0.713612036 2.812598732 1256 70506 0.000000e+00
## Bovinos 68.635160165 85.739702377 120802 2149316 0.000000e+00
## v.test
## Cacao Inf
## Aguacate Inf
## Plátano 28.734775
## Remolacha -2.031488
## Repollo -3.383437
## Lechuga -3.423216
## Zanahoria -6.885484
## Tomate -8.799118
## Habichuela -10.418577
## Yuca -11.513278
## Tabaco -12.086114
## Durazno -12.637168
## Mora -12.880418
## Arveja -14.032437
## Limón -15.119756
## Cebolla.De.Bulbo -16.072733
## Cebolla.De.Rama -22.827765
## Frijol -26.188469
## Piña -29.525032
## Papa -33.043113
## Maíz -35.476190
## Caprinos -36.278108
## Bufalos -Inf
## Bovinos -Inf
##
## $`4`
## Intern % glob % Intern freq Glob freq p.value
## Piña 13.89732242 0.38547275 8102 9663 0.000000e+00
## Limón 8.21969502 0.30636766 4792 7680 0.000000e+00
## Tabaco 0.49057445 0.04188620 286 1050 3.792539e-210
## Yuca 1.35508328 0.55856250 790 14002 9.956611e-109
## Tomate 0.43911559 0.16036432 256 4020 7.946393e-45
## Habichuela 0.15437658 0.03578279 90 897 5.009290e-30
## Aptitud 0.17496012 0.04635407 102 1162 4.369354e-29
## Maíz 0.74786875 0.67931444 436 17029 4.656569e-02
## Caprinos 1.31906208 1.43952909 769 36086 1.307855e-02
## Cacao 2.57465823 2.80023233 1501 70196 7.584420e-04
## Zanahoria 0.00000000 0.01456044 0 365 3.721027e-04
## Durazno 0.00000000 0.05026344 0 1260 2.637977e-13
## Arveja 0.00000000 0.05588816 0 1401 9.540189e-15
## Ovinos 1.86109539 2.37199576 1085 59461 2.023284e-17
## Cebolla.De.Rama 0.00000000 0.14488637 0 3632 1.434900e-37
## Cebolla.De.Bulbo 0.00000000 0.15302426 0 3836 1.171954e-39
## Mora 0.00000000 0.16092280 0 4034 1.102159e-41
## Papa 0.00000000 0.30110197 0 7548 1.114339e-77
## Frijol 0.01200707 0.48771498 7 12226 4.165271e-112
## Plátano 0.18182130 0.93015296 106 23317 1.326071e-118
## Bufalos 0.11149419 2.81259873 65 70506 0.000000e+00
## Bovinos 68.11780648 85.73970238 39712 2149316 0.000000e+00
## v.test
## Piña Inf
## Limón Inf
## Tabaco 30.937067
## Yuca 22.152101
## Tomate 14.047800
## Habichuela 11.384258
## Aptitud 11.193870
## Maíz 1.990229
## Caprinos -2.481623
## Cacao -3.367532
## Zanahoria -3.559118
## Durazno -7.311693
## Arveja -7.745239
## Ovinos -8.492449
## Cebolla.De.Rama -12.810348
## Cebolla.De.Bulbo -13.178172
## Mora -13.525749
## Papa -18.656685
## Frijol -22.499844
## Plátano -23.153963
## Bufalos -Inf
## Bovinos -Inf
##
## attr(,"class")
## [1] "descfreq" "list"
res.hcpc$desc.axes
##
## Link between the cluster variable and the quantitative variables
## ================================================================
## Eta2 P-value
## Dim.1 0.6777691 4.159063e-30
## Dim.3 0.6582307 1.531399e-28
## Dim.2 0.6200443 1.003300e-25
##
## Description of each cluster by quantitative variables
## =====================================================
## $`1`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.2 1165.28786 4.6216089 1.137261e-01 0 0.3569451
## Dim.3 275.38344 0.9234894 7.282353e-05 0 0.3094012
## Dim.1 -58.09645 -0.1748147 5.396080e-02 0 0.3633474
## p.value
## Dim.2 0
## Dim.3 0
## Dim.1 0
##
## $`2`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.2 174.0250 0.12427572 1.137261e-01 0.2268727 0.3569451
## Dim.1 -686.0734 0.01162426 5.396080e-02 0.2032272 0.3633474
## Dim.3 -1008.3151 -0.05291071 7.282353e-05 0.1725872 0.3094012
## p.value
## Dim.2 0
## Dim.1 0
## Dim.3 0
##
## $`3`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.3 1245.97203 1.12359200 7.282353e-05 0.32738349 0.3094012
## Dim.1 -15.65215 0.03738607 5.396080e-02 0.04279501 0.3633474
## Dim.2 -418.17329 -0.32129268 1.137261e-01 0.08948808 0.3569451
## p.value
## Dim.3 0.000000e+00
## Dim.1 3.211806e-55
## Dim.2 0.000000e+00
##
## $`4`
## v.test Mean in category Overall mean sd in category Overall sd p.value
## Dim.1 1302.3945 2.1910238 5.396080e-02 0.45333714 0.3633474 0
## Dim.3 -153.4909 -0.2143926 7.282353e-05 0.08275925 0.3094012 0
## Dim.2 -180.5003 -0.1772330 1.137261e-01 0.05118814 0.3569451 0
res.hcpc$desc.ind
## $para
## Cluster: 1
## TONA
## 0
## ------------------------------------------------------------
## Cluster: 2
## GALAN CONFINES SANTA BARBARA PIEDECUESTA SURATA
## 0.02704024 0.05627685 0.05811471 0.06400537 0.07736660
## ------------------------------------------------------------
## Cluster: 3
## EL CARMEN DE CHUCURI SAN VICENTE DE CHUCURI TEORAMA ns
## 0.1013280 0.1194722 0.1379134
## BUCARASICA ns EL PLAYON
## 0.2914150 0.4319223
## ------------------------------------------------------------
## Cluster: 4
## SAN BENITO LEBRIJA GIRON
## 0.5631282 0.5944524 1.1527144
##
## $dist
## Cluster: 1
## TONA
## 4.591288
## ------------------------------------------------------------
## Cluster: 2
## SILOS ns CAPITANEJO PUERTO WILCHES CEPITA JORDAN
## 2.560564 2.191972 1.937743 1.918890 1.913227
## ------------------------------------------------------------
## Cluster: 3
## LA PAZ EL CARMEN DE CHUCURI SAN VICENTE DE CHUCURI
## 2.0709652 1.3528934 1.1649636
## TEORAMA ns BUCARASICA ns
## 1.1528859 0.9642495
## ------------------------------------------------------------
## Cluster: 4
## GIRON SAN BENITO LEBRIJA
## 3.786275 2.078730 2.054302
summary(res, nb.dec = 3, nbelements=10, nbind = 10, ncp = 3, file="")
##
## Call:
## CA(X = dat.CA, ncp = 3, row.sup = NULL, col.sup = NULL, graph = T)
##
## The chi square of independence between the two variables is equal to 2999810 (p-value = 0 ).
##
## Eigenvalues
## Dim.1 Dim.2 Dim.3 Dim.4 Dim.5 Dim.6 Dim.7
## Variance 0.186 0.169 0.152 0.138 0.108 0.099 0.091
## % of var. 15.506 14.092 12.732 11.562 9.025 8.312 7.613
## Cumulative % of var. 15.506 29.598 42.330 53.892 62.917 71.229 78.842
## Dim.8 Dim.9 Dim.10 Dim.11 Dim.12 Dim.13 Dim.14
## Variance 0.050 0.033 0.032 0.022 0.021 0.020 0.016
## % of var. 4.195 2.787 2.699 1.876 1.751 1.640 1.314
## Cumulative % of var. 83.037 85.824 88.523 90.399 92.149 93.789 95.104
## Dim.15 Dim.16 Dim.17 Dim.18 Dim.19 Dim.20 Dim.21
## Variance 0.012 0.011 0.007 0.006 0.004 0.004 0.003
## % of var. 1.013 0.928 0.588 0.483 0.369 0.347 0.270
## Cumulative % of var. 96.116 97.045 97.632 98.115 98.484 98.831 99.101
## Dim.22 Dim.23 Dim.24 Dim.25 Dim.26 Dim.27 Dim.28
## Variance 0.003 0.002 0.002 0.002 0.001 0.000 0.000
## % of var. 0.251 0.204 0.196 0.168 0.052 0.026 0.002
## Cumulative % of var. 99.352 99.556 99.752 99.920 99.972 99.998 100.000
##
## Rows (the 10 first)
## Iner*1000 Dim.1 ctr cos2 Dim.2 ctr
## AGUADA | 0.376 | -0.059 0.002 0.012 | 0.054 0.002
## ALBANIA | 0.740 | -0.098 0.015 0.036 | 0.129 0.028
## ARATOCA | 6.359 | 0.600 0.341 0.100 | 0.492 0.253
## BARBOSA | 1.153 | -0.019 0.000 0.000 | 0.032 0.001
## BARICHARA | 4.361 | 0.257 0.121 0.051 | 0.380 0.290
## BARRANCABERMEJA | 17.977 | -0.216 0.792 0.082 | -0.213 0.842
## BETULIA | 0.371 | -0.018 0.001 0.005 | -0.006 0.000
## BOLIVAR | 3.017 | -0.127 0.259 0.159 | -0.065 0.074
## BUCARAMANGA | 1.860 | 0.447 0.157 0.156 | -0.048 0.002
## CABRERA | 0.776 | -0.028 0.001 0.003 | 0.057 0.006
## cos2 Dim.3 ctr cos2
## AGUADA 0.010 | -0.076 0.005 0.020 |
## ALBANIA 0.064 | -0.088 0.014 0.030 |
## ARATOCA 0.067 | -0.371 0.159 0.038 |
## BARBOSA 0.001 | -0.074 0.004 0.005 |
## BARICHARA 0.112 | -0.187 0.078 0.027 |
## BARRANCABERMEJA 0.079 | -0.423 3.678 0.312 |
## BETULIA 0.001 | 0.117 0.052 0.215 |
## BOLIVAR 0.041 | -0.083 0.133 0.067 |
## BUCARAMANGA 0.002 | 0.490 0.230 0.188 |
## CABRERA 0.013 | -0.107 0.024 0.046 |
##
## Columns (the 10 first)
## Iner*1000 Dim.1 ctr cos2 Dim.2 ctr
## Bovinos | 14.495 | -0.048 1.084 0.139 | -0.020 0.195
## Bufalos | 120.619 | -0.548 4.552 0.070 | -0.540 4.871
## Caprinos | 108.468 | 0.607 2.855 0.049 | 0.785 5.262
## Ovinos | 41.930 | 0.171 0.372 0.016 | 0.352 1.746
## Aguacate | 10.084 | 0.215 0.077 0.014 | 0.195 0.070
## Arveja | 9.630 | 0.049 0.001 0.000 | 0.876 0.254
## Avena | 2.923 | -0.018 0.000 0.000 | 0.496 0.000
## Cacao | 132.696 | 0.159 0.382 0.005 | -0.566 5.315
## Cebolla.De.Bulbo | 61.373 | 0.212 0.037 0.001 | 1.104 1.106
## Cebolla.De.Rama | 152.891 | -0.287 0.064 0.001 | 8.659 64.425
## cos2 Dim.3 ctr cos2
## Bovinos 0.023 | -0.032 0.585 0.062 |
## Bufalos 0.068 | -0.944 16.435 0.208 |
## Caprinos 0.082 | -0.480 2.177 0.031 |
## Ovinos 0.070 | -0.063 0.062 0.002 |
## Aguacate 0.012 | 1.346 3.676 0.555 |
## Arveja 0.045 | 0.008 0.000 0.000 |
## Avena 0.000 | -0.301 0.000 0.000 |
## Cacao 0.068 | 1.960 70.595 0.811 |
## Cebolla.De.Bulbo 0.030 | 0.386 0.150 0.004 |
## Cebolla.De.Rama 0.711 | 1.734 2.860 0.029 |
res$eig
## eigenvalue percentage of variance cumulative percentage of variance
## dim 1 1.855562e-01 15.506004044 15.50600
## dim 2 1.686353e-01 14.092012090 29.59802
## dim 3 1.523572e-01 12.731735508 42.32975
## dim 4 1.383623e-01 11.562245075 53.89200
## dim 5 1.080036e-01 9.025323149 62.91732
## dim 6 9.946163e-02 8.311513428 71.22883
## dim 7 9.110032e-02 7.612800558 78.84163
## dim 8 5.020210e-02 4.195139278 83.03677
## dim 9 3.335163e-02 2.787029678 85.82380
## dim 10 3.230339e-02 2.699433562 88.52324
## dim 11 2.244435e-02 1.875562926 90.39880
## dim 12 2.095002e-02 1.750688569 92.14949
## dim 13 1.962538e-02 1.639995295 93.78948
## dim 14 1.572847e-02 1.314350300 95.10383
## dim 15 1.211723e-02 1.012576205 96.11641
## dim 16 1.110707e-02 0.928162522 97.04457
## dim 17 7.030524e-03 0.587505856 97.63208
## dim 18 5.779320e-03 0.482948957 98.11503
## dim 19 4.420220e-03 0.369375776 98.48440
## dim 20 4.153166e-03 0.347059374 98.83146
## dim 21 3.227675e-03 0.269720701 99.10118
## dim 22 3.003822e-03 0.251014478 99.35220
## dim 23 2.438931e-03 0.203809289 99.55601
## dim 24 2.343440e-03 0.195829644 99.75184
## dim 25 2.009364e-03 0.167912563 99.91975
## dim 26 6.235448e-04 0.052106531 99.97186
## dim 27 3.164912e-04 0.026447597 99.99830
## dim 28 2.030812e-05 0.001697049 100.00000
res$var
## NULL
res$ind
## NULL
res$ind.sup
## NULL
dimdesc(res, axes=1:3)
## $`Dim 1`
## $`Dim 1`$row
## coord
## PUERTO WILCHES -0.435039064
## BARRANCABERMEJA -0.216480413
## SILOS ns -0.207962525
## SABANA DE TORRES -0.197424093
## TONA -0.174814747
## MUTISCUA ns -0.146029094
## CHITAGA ns -0.135221091
## CIMITARRA -0.129135932
## BOLIVAR -0.127394194
## PUERTO PARRA -0.118474778
## PAMPLONA ns -0.107525137
## SIMACOTA -0.103089986
## VELEZ -0.101498350
## GAMBITA -0.100947259
## ALBANIA -0.097957153
## OIBA -0.096480988
## LANDAZURI -0.093247532
## GUAVATA -0.091842225
## LA ESPERANZA ns -0.091464584
## HERRAN ns -0.091381168
## LABATECA ns -0.089318897
## EL PEÑON -0.087935226
## GRAMALOTE ns -0.087072260
## TOLEDO ns -0.086920992
## GUADALUPE -0.086915986
## LOURDES ns -0.086175410
## GUAPOTA -0.085527032
## CHARTA -0.084288407
## CHARALA -0.084123095
## PUENTE NACIONAL -0.083949961
## SANTIAGO ns -0.081906416
## SALAZAR ns -0.079936445
## SANTA HELENA DEL OPON -0.078372973
## SUCRE -0.077352508
## CACHIRA ns -0.076875313
## CHINACOTA ns -0.075741080
## BOCHALEMA ns -0.075599923
## CUCUTA ns -0.074954852
## PALMAS DEL SOCORRO -0.074804034
## DURANIA ns -0.068077542
## CHIPATA -0.067844679
## EL GUACAMAYO -0.067177286
## VALLE DE SAN JOSE -0.065294695
## LA BELLEZA -0.064496348
## RAGONVALIA ns -0.063666952
## CUCUTILLA ns -0.063556526
## PALMAR -0.060404594
## AGUADA -0.059010365
## EL ZULIA ns -0.058669878
## JESUS MARIA -0.058631801
## EL TARRA ns -0.057947375
## CONTRATACION -0.054418452
## FLORIAN -0.051542349
## PAMPLONITA ns -0.050505194
## SARDINATA ns -0.046610263
## HACARI ns -0.046447263
## COROMORO -0.045779801
## ARBOLEDAS ns -0.041039183
## SAN ANDRES -0.040451574
## CHIMA -0.038610939
## RIONEGRO -0.033080527
## GUACA -0.028896176
## CABRERA -0.028403911
## VETAS -0.026200480
## VILLA CARO ns -0.021910212
## SAN CAYETANO ns -0.019969780
## ENCINO -0.019487943
## BARBOSA -0.018734788
## BETULIA -0.018384645
## TIBU ns -0.017771095
## CACOTA ns -0.012475189
## VILLA DEL ROSARIO ns -0.012369960
## CALIFORNIA -0.007844470
## SANTA BARBARA -0.005721343
## ONZAGA -0.001537523
## MATANZA 0.009789123
## SAN VICENTE DE CHUCURI 0.010671893
## CONVENCION ns 0.011807360
## CONFINES 0.014578666
## OCAMONTE 0.020755290
## SURATA 0.021138505
## PUERTO SANTANDER ns 0.029075575
## MALAGA 0.030506310
## HATO 0.030693358
## GALAN 0.033606637
## OCAÑA ns 0.034943877
## BUCARASICA ns 0.036694945
## EL CARMEN ns 0.042714836
## CERRITO 0.047723747
## CARCASI 0.054635625
## EL CARMEN DE CHUCURI 0.061491417
## SOCORRO 0.066146727
## CONCEPCION 0.078500390
## PIEDECUESTA 0.079613240
## ABREGO ns 0.081791186
## ENCISO 0.085942133
## SUAITA 0.086112578
## EL PLAYON 0.087768011
## MOGOTES 0.091918960
## LA PAZ 0.103853569
## SAN MIGUEL 0.111390429
## PINCHOTE 0.112664252
## LA PLAYA ns 0.151807399
## ZAPATOCA 0.156313870
## TEORAMA ns 0.161992024
## SAN CALIXTO ns 0.164283135
## GÜEPSA 0.164612214
## SAN JOSE DE MIRANDA 0.172888044
## FLORIDABLANCA 0.177898595
## PARAMO 0.179580974
## SAN JOAQUIN 0.184094566
## MACARAVITA 0.197426535
## SAN GIL 0.198764782
## BARICHARA 0.257187510
## CURITI 0.285640400
## LOS PATIOS ns 0.356144436
## MOLAGAVITA 0.371676689
## LOS SANTOS 0.404024401
## BUCARAMANGA 0.446953913
## ARATOCA 0.599728954
## VILLANUEVA 0.636822760
## CAPITANEJO 0.819082925
## JORDAN 0.879260557
## CEPITA 0.989640072
## LEBRIJA 2.048702988
## SAN BENITO 2.088409223
## GIRON 3.765671666
##
## $`Dim 1`$col
## coord
## Bufalos -0.54799904
## Remolacha -0.31883761
## Cebolla.De.Rama -0.28702970
## Lechuga -0.27319799
## Zanahoria -0.27289311
## Papa -0.23403356
## Haba -0.21115743
## Repollo -0.12632471
## Durazno -0.10550021
## Bovinos -0.04842915
## Plátano -0.04769896
## Mora -0.03251157
## Avena -0.01821067
## Arveja 0.04929681
## Cacao 0.15913996
## Ovinos 0.17066259
## Maíz 0.17681104
## Trigo 0.17690449
## Cebolla.De.Bulbo 0.21205166
## Aguacate 0.21476631
## Frijol 0.24407754
## Yuca 0.36669939
## Aptitud 0.42211753
## Tomate 0.58231783
## Caprinos 0.60658930
## Habichuela 0.68984525
## Tabaco 3.14699623
## Limón 4.28089289
## Piña 5.23189708
##
##
## $`Dim 2`
## $`Dim 2`$row
## coord
## LA PAZ -0.517863799
## PUERTO WILCHES -0.394850547
## EL CARMEN DE CHUCURI -0.329681023
## SAN VICENTE DE CHUCURI -0.312675933
## TEORAMA ns -0.277956380
## GIRON -0.258963016
## BARRANCABERMEJA -0.212843551
## LEBRIJA -0.201608413
## BUCARASICA ns -0.192451556
## EL PLAYON -0.181165745
## SABANA DE TORRES -0.146008660
## SARDINATA ns -0.133146947
## RIONEGRO -0.127061462
## LA BELLEZA -0.126672348
## FLORIAN -0.107779983
## EL TARRA ns -0.106558241
## SUCRE -0.099878219
## SAN BENITO -0.098582664
## CONVENCION ns -0.085841093
## SANTA HELENA DEL OPON -0.083549035
## CIMITARRA -0.077783862
## SOCORRO -0.074660318
## EL GUACAMAYO -0.073916139
## PALMAS DEL SOCORRO -0.069916564
## CHIMA -0.069459327
## LOURDES ns -0.068687551
## VALLE DE SAN JOSE -0.067441120
## TIBU ns -0.065645555
## BOLIVAR -0.064833156
## SIMACOTA -0.059623012
## EL PEÑON -0.059010342
## PUERTO PARRA -0.058893460
## GUAPOTA -0.058046851
## LA ESPERANZA ns -0.054456901
## VELEZ -0.049235083
## LANDAZURI -0.049209909
## CUCUTILLA ns -0.049180311
## BUCARAMANGA -0.047523314
## GUADALUPE -0.044040526
## HACARI ns -0.039044595
## EL ZULIA ns -0.037501912
## SUAITA -0.036929924
## CUCUTA ns -0.033175792
## OIBA -0.032190077
## CACOTA ns -0.029772576
## DURANIA ns -0.027129222
## GRAMALOTE ns -0.025851623
## BOCHALEMA ns -0.023835489
## SANTIAGO ns -0.023016397
## GUAVATA -0.022809675
## GAMBITA -0.019569732
## CHARALA -0.017893263
## SALAZAR ns -0.017332682
## PALMAR -0.017193220
## SAN CAYETANO ns -0.016240741
## PUENTE NACIONAL -0.014034854
## TOLEDO ns -0.012860607
## BETULIA -0.006161315
## CHINACOTA ns -0.004388167
## CONTRATACION 0.001738680
## JESUS MARIA 0.003161343
## ARBOLEDAS ns 0.007099756
## OCAMONTE 0.012126372
## LABATECA ns 0.015842076
## HATO 0.017256843
## HERRAN ns 0.027849668
## MATANZA 0.031475182
## BARBOSA 0.032447763
## CHARTA 0.045498961
## AGUADA 0.054154312
## COROMORO 0.056260460
## CABRERA 0.057008347
## CHIPATA 0.057779319
## VILLA DEL ROSARIO ns 0.065975848
## SAN ANDRES 0.072830993
## RAGONVALIA ns 0.085262620
## PAMPLONITA ns 0.087911956
## MOGOTES 0.091224196
## PARAMO 0.092086468
## EL CARMEN ns 0.093083913
## VILLA CARO ns 0.094911797
## CONFINES 0.113881491
## GALAN 0.119048629
## PUERTO SANTANDER ns 0.122163485
## ALBANIA 0.129481736
## PIEDECUESTA 0.141082037
## SURATA 0.156308267
## MALAGA 0.164171316
## ENCINO 0.171051752
## PINCHOTE 0.173195380
## CACHIRA ns 0.183230168
## SANTA BARBARA 0.185674213
## SAN GIL 0.203190414
## CALIFORNIA 0.203699432
## ONZAGA 0.209471379
## ENCISO 0.232365143
## GÜEPSA 0.247604884
## ZAPATOCA 0.261717915
## LOS SANTOS 0.278337069
## SAN MIGUEL 0.307325815
## CONCEPCION 0.316562500
## MACARAVITA 0.330179585
## SAN JOAQUIN 0.334417008
## SAN JOSE DE MIRANDA 0.338888797
## VETAS 0.357389865
## FLORIDABLANCA 0.360282891
## SAN CALIXTO ns 0.361679818
## ABREGO ns 0.364127889
## BARICHARA 0.380167487
## CERRITO 0.387576854
## VILLANUEVA 0.392099078
## CURITI 0.443808022
## GUACA 0.459706967
## MUTISCUA ns 0.464024209
## LOS PATIOS ns 0.465279875
## ARATOCA 0.492037053
## CARCASI 0.503039060
## PAMPLONA ns 0.558260951
## JORDAN 0.560978147
## MOLAGAVITA 0.592669023
## LA PLAYA ns 0.614817904
## CHITAGA ns 0.715421235
## CEPITA 0.769073198
## CAPITANEJO 0.999115777
## OCAÑA ns 1.245269887
## SILOS ns 1.920428802
## TONA 4.621608944
##
## $`Dim 2`$col
## coord
## Cacao -0.56573917
## Bufalos -0.54044120
## Piña -0.44569112
## Limón -0.27701327
## Plátano -0.15846741
## Aptitud -0.03599419
## Yuca -0.01999220
## Bovinos -0.01959077
## Aguacate 0.19514746
## Mora 0.24643051
## Maíz 0.30432064
## Ovinos 0.35234511
## Tabaco 0.46895976
## Avena 0.49603872
## Habichuela 0.52660651
## Trigo 0.67807762
## Tomate 0.70695290
## Caprinos 0.78511867
## Frijol 0.85234715
## Arveja 0.87568024
## Haba 1.09445812
## Cebolla.De.Bulbo 1.10417979
## Lechuga 1.12188986
## Remolacha 1.27472563
## Durazno 1.32090404
## Zanahoria 1.46433739
## Repollo 1.81288665
## Papa 2.59700332
## Cebolla.De.Rama 8.65938257
## CA plantas
row.names(dat.CA1)<-dat$Municipio
res<-CA(dat.CA1, ncp=3, row.sup=NULL, col.sup=NULL, graph = T)
get_eigenvalue(res)
## eigenvalue variance.percent cumulative.variance.percent
## Dim.1 0.7558015961 19.57240577 19.57241
## Dim.2 0.6094010012 15.78118349 35.35359
## Dim.3 0.5285142740 13.68652287 49.04011
## Dim.4 0.3519951913 9.11534555 58.15546
## Dim.5 0.2969107203 7.68886587 65.84432
## Dim.6 0.2481422081 6.42594566 72.27027
## Dim.7 0.1913562922 4.95540499 77.22567
## Dim.8 0.1594588306 4.12938125 81.35506
## Dim.9 0.1314853705 3.40497432 84.76003
## Dim.10 0.1153322119 2.98666854 87.74670
## Dim.11 0.0765693565 1.98285704 89.72956
## Dim.12 0.0702440668 1.81905593 91.54861
## Dim.13 0.0632314944 1.63745680 93.18607
## Dim.14 0.0602219547 1.55952110 94.74559
## Dim.15 0.0436709497 1.13091260 95.87650
## Dim.16 0.0420617102 1.08923938 96.96574
## Dim.17 0.0352293866 0.91230801 97.87805
## Dim.18 0.0254082138 0.65797674 98.53603
## Dim.19 0.0245871375 0.63671397 99.17274
## Dim.20 0.0219835416 0.56929067 99.74203
## Dim.21 0.0056264844 0.14570469 99.88774
## Dim.22 0.0039084734 0.10121469 99.98895
## Dim.23 0.0004267064 0.01105008 100.00000
fviz_eig(res, addlabels = TRUE, ylim = c(0, 20))
get_ca_row(res)
## Correspondence Analysis - Results for rows
## ===================================================
## Name Description
## 1 "$coord" "Coordinates for the rows"
## 2 "$cos2" "Cos2 for the rows"
## 3 "$contrib" "contributions of the rows"
## 4 "$inertia" "Inertia of the rows"
get_ca_col(res)
## Correspondence Analysis - Results for columns
## ===================================================
## Name Description
## 1 "$coord" "Coordinates for the columns"
## 2 "$cos2" "Cos2 for the columns"
## 3 "$contrib" "contributions of the columns"
## 4 "$inertia" "Inertia of the columns"
fviz_ca_row(res, repel=T)
## Warning: ggrepel: 102 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
fviz_ca_col(res, repel =T)
fviz_ca_biplot(res, repel = T)
## Warning: ggrepel: 103 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
res.hcpc<-HCPC(res ,nb.clust=-1,consol=TRUE,min=3,max=10,graph=TRUE)
res.plantas<-res.hcpc$data.clust[,ncol(res.hcpc$data.clust),drop=F]
res.plantas
## clust
## SAN VICENTE DE CHUCURI 1
## LA PAZ 1
## LEBRIJA 2
## EL CARMEN DE CHUCURI 1
## GIRON 2
## GUAPOTA 1
## RIONEGRO 1
## CIMITARRA 1
## EL PLAYON 1
## LANDAZURI 1
## LA BELLEZA 1
## SANTA HELENA DEL OPON 1
## SUCRE 1
## BETULIA 1
## FLORIAN 1
## SARDINATA ns 1
## EL TARRA ns 1
## SIMACOTA 1
## GALAN 1
## TEORAMA ns 1
## BUCARAMANGA 1
## CUCUTA ns 1
## BUCARASICA ns 1
## VELEZ 1
## GUADALUPE 1
## LA ESPERANZA ns 1
## CONTRATACION 1
## CHIMA 1
## EL GUACAMAYO 1
## SAN BENITO 2
## LOURDES ns 1
## PALMAS DEL SOCORRO 1
## BARRANCABERMEJA 3
## CONVENCION ns 1
## SOCORRO 3
## ARATOCA 2
## OCAMONTE 1
## EL ZULIA ns 3
## EL PEÑON 1
## VALLE DE SAN JOSE 3
## PUERTO PARRA 3
## HACARI ns 1
## PUERTO SANTANDER ns 3
## TIBU ns 3
## LOS SANTOS 2
## SAN CAYETANO ns 3
## TOLEDO ns 3
## CHARALA 3
## SUAITA 3
## ARBOLEDAS ns 3
## DURANIA ns 3
## PALMAR 3
## CUCUTILLA ns 3
## SANTIAGO ns 3
## SABANA DE TORRES 3
## SALAZAR ns 3
## PINCHOTE 3
## LOS PATIOS ns 3
## EL CARMEN ns 3
## BOCHALEMA ns 3
## MOGOTES 3
## JORDAN 2
## OIBA 3
## MATANZA 3
## COROMORO 3
## CAPITANEJO 3
## PUERTO WILCHES 3
## CEPITA 2
## PARAMO 3
## GÜEPSA 3
## VILLA DEL ROSARIO ns 3
## BARBOSA 3
## GUAVATA 3
## CABRERA 3
## VILLANUEVA 3
## CHINACOTA ns 3
## CACOTA ns 1
## BOLIVAR 1
## JESUS MARIA 3
## GRAMALOTE ns 3
## ZAPATOCA 1
## SAN GIL 3
## HATO 1
## LABATECA ns 3
## MACARAVITA 3
## VILLA CARO ns 3
## CURITI 3
## SAN CALIXTO ns 3
## CONFINES 3
## SAN JOSE DE MIRANDA 3
## FLORIDABLANCA 1
## BARICHARA 3
## PIEDECUESTA 3
## OCAÑA ns 3
## SAN MIGUEL 3
## ENCISO 3
## SURATA 3
## SAN JOAQUIN 3
## ABREGO ns 3
## CHIPATA 3
## LA PLAYA ns 3
## MALAGA 3
## RAGONVALIA ns 3
## AGUADA 3
## SANTA BARBARA 3
## PUENTE NACIONAL 3
## MOLAGAVITA 3
## CHARTA 3
## SAN ANDRES 3
## GAMBITA 3
## CALIFORNIA 3
## PAMPLONITA ns 3
## HERRAN ns 3
## ONZAGA 4
## TONA 4
## GUACA 4
## ENCINO 4
## CONCEPCION 4
## PAMPLONA ns 4
## CARCASI 4
## ALBANIA 4
## CACHIRA ns 4
## CHITAGA ns 4
## CERRITO 4
## SILOS ns 4
## MUTISCUA ns 4
## VETAS 4
write.csv(res.plantas, "~/Library/CloudStorage/GoogleDrive-icarounam@gmail.com/Mi unidad/Agrosavia/Ingreso/Evas/resplantas.csv")
res.hcpc$desc.var
## $`1`
## Intern % glob % Intern freq Glob freq p.value
## Cacao 69.068072368 37.37985313 63678 70196 0.000000e+00
## Aguacate 5.756215020 4.12479831 5307 7746 2.272271e-272
## Haba 0.000000000 0.01863774 0 35 1.087218e-10
## Aptitud 0.482667361 0.61877300 445 1162 1.416976e-13
## Trigo 0.000000000 0.02609284 0 49 8.504839e-15
## Remolacha 0.000000000 0.02822286 0 53 5.705018e-16
## Repollo 0.009761812 0.05804325 9 109 7.221886e-20
## Lechuga 0.009761812 0.05910826 9 111 2.211686e-20
## Yuca 6.331077270 7.45616137 5837 14002 1.573819e-74
## Mora 1.520673348 2.14813276 1402 4034 6.580583e-77
## Habichuela 0.154019697 0.47765867 142 897 2.848402e-97
## Zanahoria 0.000000000 0.19436501 0 365 1.318638e-107
## Durazno 0.101956701 0.67095867 94 1260 8.939397e-229
## Arveja 0.104125992 0.74604214 96 1401 9.566011e-264
## Tabaco 0.000000000 0.55913223 0 1050 1.459749e-309
## Tomate 0.769013840 2.14067767 709 4020 0.000000e+00
## Plátano 9.057876697 12.41646298 8351 23317 0.000000e+00
## Piña 1.185517810 5.14561401 1093 9663 0.000000e+00
## Papa 0.080263786 4.01936195 74 7548 0.000000e+00
## Maíz 3.363486485 9.06805970 3101 17029 0.000000e+00
## Limón 0.812399670 4.08965286 749 7680 0.000000e+00
## Frijol 1.046683153 6.51042915 965 12226 0.000000e+00
## Cebolla.De.Bulbo 0.146427177 2.04269640 135 3836 0.000000e+00
## v.test
## Cacao Inf
## Aguacate 35.261678
## Haba -6.454298
## Aptitud -7.394728
## Trigo -7.759822
## Remolacha -8.095456
## Repollo -9.124274
## Lechuga -9.251594
## Yuca -18.264941
## Mora -18.561528
## Habichuela -20.930110
## Zanahoria -22.035404
## Durazno -32.292412
## Arveja -34.694410
## Tabaco -37.610285
## Tomate -Inf
## Plátano -Inf
## Piña -Inf
## Papa -Inf
## Maíz -Inf
## Limón -Inf
## Frijol -Inf
## Cebolla.De.Bulbo -Inf
##
## $`2`
## Intern % glob % Intern freq Glob freq p.value
## Tabaco 4.17767782 0.55913223 743 1050 0.000000e+00
## Piña 45.86449255 5.14561401 8157 9663 0.000000e+00
## Limón 28.20354231 4.08965286 5016 7680 0.000000e+00
## Habichuela 0.70846219 0.47765867 126 897 1.218619e-05
## Trigo 0.00000000 0.02609284 0 49 1.525561e-02
## Remolacha 0.00000000 0.02822286 0 53 1.024563e-02
## Tomate 1.79364633 2.14067767 319 4020 6.487693e-04
## Repollo 0.00000000 0.05804325 0 109 3.887390e-05
## Lechuga 0.00000000 0.05910826 0 111 3.185547e-05
## Zanahoria 0.00000000 0.19436501 0 365 3.259076e-16
## Arveja 0.06747259 0.74604214 12 1401 7.837665e-44
## Durazno 0.00000000 0.67095867 0 1260 4.602189e-55
## Yuca 4.53190891 7.45616137 806 14002 5.926582e-62
## Aguacate 1.32133821 4.12479831 235 7746 3.388254e-113
## Cebolla.De.Bulbo 0.00000000 2.04269640 0 3836 5.495098e-168
## Mora 0.01124543 2.14813276 2 4034 9.202720e-172
## Maíz 2.92381220 9.06805970 520 17029 8.636334e-255
## Plátano 0.82091650 12.41646298 146 23317 0.000000e+00
## Papa 0.00000000 4.01936195 0 7548 0.000000e+00
## Frijol 0.52853528 6.51042915 94 12226 0.000000e+00
## Cacao 8.47343267 37.37985313 1507 70196 0.000000e+00
## v.test
## Tabaco Inf
## Piña Inf
## Limón Inf
## Habichuela 4.374230
## Trigo -2.426253
## Remolacha -2.567427
## Tomate -3.410361
## Repollo -4.114073
## Lechuga -4.159780
## Zanahoria -8.163332
## Arveja -13.884747
## Durazno -15.629246
## Yuca -16.609723
## Aguacate -22.610866
## Cebolla.De.Bulbo -27.625517
## Mora -27.938077
## Maíz -34.095378
## Plátano -Inf
## Papa -Inf
## Frijol -Inf
## Cacao -Inf
##
## $`3`
## Intern % glob % Intern freq Glob freq p.value
## Yuca 10.972999419 7.45616137 7177 14002 0.000000e+00
## Tomate 4.407852491 2.14067767 2883 4020 0.000000e+00
## Plátano 22.565208085 12.41646298 14759 23317 0.000000e+00
## Maíz 19.557838730 9.06805970 12792 17029 0.000000e+00
## Frijol 15.565850228 6.51042915 10181 12226 0.000000e+00
## Cebolla.De.Bulbo 5.123383176 2.04269640 3351 3836 0.000000e+00
## Mora 3.511910222 2.14813276 2297 4034 5.628680e-184
## Habichuela 0.952512002 0.47765867 623 897 6.464830e-99
## Arveja 1.220071553 0.74604214 798 1401 3.456236e-64
## Aptitud 0.937222885 0.61877300 613 1162 7.699424e-36
## Trigo 0.059627557 0.02609284 39 49 3.554521e-10
## Tabaco 0.469375898 0.55913223 307 1050 1.282372e-04
## Haba 0.001528912 0.01863774 1 35 1.221391e-05
## Remolacha 0.000000000 0.02822286 0 53 2.780722e-10
## Lechuga 0.000000000 0.05910826 0 111 4.502048e-21
## Durazno 0.344005137 0.67095867 225 1260 1.009277e-40
## Zanahoria 0.003057823 0.19436501 2 365 4.299921e-64
## Limón 2.878940770 4.08965286 1883 7680 4.893457e-88
## Aguacate 2.771916950 4.12479831 1813 7746 2.218722e-109
## Piña 0.631440541 5.14561401 413 9663 0.000000e+00
## Papa 0.314955814 4.01936195 206 7548 0.000000e+00
## Cacao 7.661376632 37.37985313 5011 70196 0.000000e+00
## v.test
## Yuca Inf
## Tomate Inf
## Plátano Inf
## Maíz Inf
## Frijol Inf
## Cebolla.De.Bulbo Inf
## Mora 28.925977
## Habichuela 21.109798
## Arveja 16.915542
## Aptitud 12.497543
## Trigo 6.272432
## Tabaco -3.829820
## Haba -4.373734
## Remolacha -6.310535
## Lechuga -9.420206
## Durazno -13.361919
## Zanahoria -16.902669
## Limón -19.890753
## Aguacate -22.219634
## Piña -Inf
## Papa -Inf
## Cacao -Inf
##
## $`4`
## Intern % glob % Intern freq Glob freq p.value
## Zanahoria 2.92647533 0.19436501 363 365 0.000000e+00
## Papa 58.59400193 4.01936195 7268 7548 0.000000e+00
## Durazno 7.58626250 0.67095867 941 1260 0.000000e+00
## Arveja 3.99064818 0.74604214 495 1401 5.667426e-221
## Lechuga 0.82231538 0.05910826 102 111 1.591533e-108
## Remolacha 0.42728152 0.02822286 53 53 5.127573e-63
## Repollo 0.56433409 0.05804325 70 109 1.783971e-54
## Haba 0.27410513 0.01863774 34 35 4.722465e-39
## Frijol 7.94904869 6.51042915 986 12226 7.454794e-11
## Cebolla.De.Bulbo 2.82167043 2.04269640 350 3836 1.882605e-09
## Mora 2.68461787 2.14813276 333 4034 4.189222e-05
## Trigo 0.08061916 0.02609284 10 49 2.383733e-03
## Aguacate 3.15220896 4.12479831 391 7746 5.728786e-09
## Habichuela 0.04837149 0.47765867 6 897 4.161201e-19
## Tomate 0.87874879 2.14067767 109 4020 2.550425e-29
## Aptitud 0.01612383 0.61877300 2 1162 1.776218e-31
## Tabaco 0.00000000 0.55913223 0 1050 1.120370e-31
## Maíz 4.96613995 9.06805970 616 17029 1.522609e-70
## Limón 0.25798130 4.08965286 32 7680 1.298555e-180
## Yuca 1.46726862 7.45616137 182 14002 6.008282e-216
## Piña 0.00000000 5.14561401 0 9663 4.081506e-295
## Plátano 0.49177685 12.41646298 61 23317 0.000000e+00
## Cacao 0.00000000 37.37985313 0 70196 0.000000e+00
## v.test
## Zanahoria Inf
## Papa Inf
## Durazno Inf
## Arveja 31.731765
## Lechuga 22.130960
## Remolacha 16.755904
## Repollo 15.542669
## Haba 13.072598
## Frijol 6.511210
## Cebolla.De.Bulbo 6.007625
## Mora 4.096788
## Trigo 3.037723
## Aguacate -5.824488
## Habichuela -8.932598
## Tomate -11.241489
## Aptitud -11.671801
## Tabaco -11.710935
## Maíz -17.756960
## Limón -28.657343
## Yuca -31.365355
## Piña -36.715838
## Plátano -Inf
## Cacao -Inf
##
## attr(,"class")
## [1] "descfreq" "list"
res.hcpc$desc.axes
##
## Link between the cluster variable and the quantitative variables
## ================================================================
## Eta2 P-value
## Dim.1 0.9190926 5.926476e-67
## Dim.3 0.8003517 7.385224e-43
## Dim.2 0.7943667 4.525677e-42
##
## Description of each cluster by quantitative variables
## =====================================================
## $`1`
## v.test Mean in category Overall mean sd in category Overall sd p.value
## Dim.1 -139.8659 -0.23470147 0.38433939 0.1367092 1.1172090 0
## Dim.3 -164.8737 -0.02643592 0.40864921 0.2212057 0.6661157 0
## Dim.2 -194.6509 -0.30205326 0.05340772 0.1302027 0.4609591 0
##
## $`2`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.2 357.01610 1.6475501 0.05340772 0.3818306 0.4609591
## Dim.3 -27.91265 0.2285434 0.40864921 0.3046975 0.6661157
## Dim.1 -54.27193 -0.2029965 0.38433939 0.0956347 1.1172090
## p.value
## Dim.2 0.000000e+00
## Dim.3 1.873469e-171
## Dim.1 0.000000e+00
##
## $`3`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.3 351.58062 0.87560720 0.40864921 0.2865648 0.6661157
## Dim.2 16.91997 0.06895897 0.05340772 0.2344284 0.4609591
## Dim.1 -122.83564 0.11071045 0.38433939 0.3078932 1.1172090
## p.value
## Dim.3 0.000000e+00
## Dim.2 3.205829e-64
## Dim.1 0.000000e+00
##
## $`4`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.1 411.23755 3.24982947 0.38433939 0.58695114 1.1172090
## Dim.2 -12.88936 0.01635112 0.05340772 0.02938392 0.4609591
## Dim.3 -294.76996 -0.81597995 0.40864921 0.44818255 0.6661157
## p.value
## Dim.1 0.00000e+00
## Dim.2 5.16692e-38
## Dim.3 0.00000e+00
res.hcpc$desc.ind
## $para
## Cluster: 1
## VELEZ TEORAMA ns EL PEÑON EL GUACAMAYO EL TARRA ns
## 0.07726601 0.09030418 0.10866468 0.12430141 0.13954742
## ------------------------------------------------------------
## Cluster: 2
## ARATOCA LOS SANTOS JORDAN SAN BENITO CEPITA
## 0.3629609 0.5638780 0.5997987 0.6688524 0.6963260
## ------------------------------------------------------------
## Cluster: 3
## CONFINES GRAMALOTE ns JESUS MARIA SURATA GÜEPSA
## 0.1128080 0.1255393 0.1336569 0.1620568 0.1834908
## ------------------------------------------------------------
## Cluster: 4
## CACHIRA ns ALBANIA PAMPLONA ns CARCASI CHITAGA ns
## 0.2142482 0.2689004 0.3474496 0.3526672 0.3825106
##
## $dist
## Cluster: 1
## SAN VICENTE DE CHUCURI LA PAZ EL CARMEN DE CHUCURI
## 1.885169 1.804247 1.677543
## LA BELLEZA GUAPOTA
## 1.670649 1.541147
## ------------------------------------------------------------
## Cluster: 2
## GIRON LEBRIJA LOS SANTOS JORDAN ARATOCA
## 2.820866 2.668903 2.236251 1.674798 1.534392
## ------------------------------------------------------------
## Cluster: 3
## LA PLAYA ns MALAGA HERRAN ns CHIPATA ABREGO ns
## 1.842911 1.763051 1.705919 1.694812 1.674022
## ------------------------------------------------------------
## Cluster: 4
## VETAS MUTISCUA ns SILOS ns CERRITO CHITAGA ns
## 4.577471 4.463326 4.394968 4.225269 3.515282
summary(res, nb.dec = 3, nbelements=10, nbind = 10, ncp = 3, file="")
##
## Call:
## CA(X = dat.CA1, ncp = 3, row.sup = NULL, col.sup = NULL, graph = T)
##
## The chi square of independence between the two variables is equal to 725167.6 (p-value = 0 ).
##
## Eigenvalues
## Dim.1 Dim.2 Dim.3 Dim.4 Dim.5 Dim.6 Dim.7
## Variance 0.756 0.609 0.529 0.352 0.297 0.248 0.191
## % of var. 19.572 15.781 13.687 9.115 7.689 6.426 4.955
## Cumulative % of var. 19.572 35.354 49.040 58.155 65.844 72.270 77.226
## Dim.8 Dim.9 Dim.10 Dim.11 Dim.12 Dim.13 Dim.14
## Variance 0.159 0.131 0.115 0.077 0.070 0.063 0.060
## % of var. 4.129 3.405 2.987 1.983 1.819 1.637 1.560
## Cumulative % of var. 81.355 84.760 87.747 89.730 91.549 93.186 94.746
## Dim.15 Dim.16 Dim.17 Dim.18 Dim.19 Dim.20 Dim.21
## Variance 0.044 0.042 0.035 0.025 0.025 0.022 0.006
## % of var. 1.131 1.089 0.912 0.658 0.637 0.569 0.146
## Cumulative % of var. 95.877 96.966 97.878 98.536 99.173 99.742 99.888
## Dim.22 Dim.23
## Variance 0.004 0.000
## % of var. 0.101 0.011
## Cumulative % of var. 99.989 100.000
##
## Rows (the 10 first)
## Iner*1000 Dim.1 ctr cos2 Dim.2 ctr
## AGUADA | 6.166 | 0.421 0.019 0.023 | 0.037 0.000
## ALBANIA | 15.615 | 2.997 1.480 0.717 | 0.047 0.000
## ARATOCA | 8.758 | -0.193 0.006 0.005 | 1.452 0.446
## BARBOSA | 17.752 | 0.005 0.000 0.000 | 0.228 0.008
## BARICHARA | 20.826 | 0.244 0.034 0.012 | 0.407 0.117
## BARRANCABERMEJA | 14.696 | -0.212 0.044 0.023 | 0.472 0.270
## BETULIA | 2.848 | -0.322 0.061 0.163 | -0.314 0.072
## BOLIVAR | 42.243 | 0.059 0.004 0.001 | -0.408 0.248
## BUCARAMANGA | 4.356 | -0.295 0.055 0.096 | 0.063 0.003
## CABRERA | 14.051 | 0.034 0.000 0.000 | 0.093 0.002
## cos2 Dim.3 ctr cos2
## AGUADA 0.000 | 1.087 0.176 0.151 |
## ALBANIA 0.000 | -0.896 0.189 0.064 |
## ARATOCA 0.310 | 0.273 0.018 0.011 |
## BARBOSA 0.003 | 0.988 0.179 0.053 |
## BARICHARA 0.034 | 1.190 1.150 0.292 |
## BARRANCABERMEJA 0.112 | 0.338 0.159 0.057 |
## BETULIA 0.155 | -0.422 0.150 0.279 |
## BOLIVAR 0.036 | -0.110 0.021 0.003 |
## BUCARAMANGA 0.004 | -0.087 0.007 0.008 |
## CABRERA 0.001 | 1.039 0.343 0.129 |
##
## Columns (the 10 first)
## Iner*1000 Dim.1 ctr cos2 Dim.2 ctr
## Aguacate | 51.602 | -0.181 0.180 0.026 | -0.240 0.388
## Arveja | 65.067 | 1.282 1.621 0.188 | -0.023 0.001
## Avena | 52.621 | 1.246 0.002 0.000 | -0.089 0.000
## Cacao | 326.577 | -0.441 9.601 0.222 | -0.514 16.191
## Cebolla.De.Bulbo | 206.101 | 0.547 0.808 0.030 | 0.019 0.001
## Durazno | 177.075 | 2.962 7.790 0.332 | 0.010 0.000
## Frijol | 213.445 | 0.438 1.653 0.059 | 0.106 0.121
## Haba | 9.914 | 3.246 0.260 0.198 | -0.028 0.000
## Habichuela | 54.667 | 0.029 0.001 0.000 | 0.393 0.121
## Lechuga | 136.843 | 4.179 1.366 0.075 | -0.050 0.000
## cos2 Dim.3 ctr cos2
## Aguacate 0.046 | -0.237 0.437 0.045 |
## Arveja 0.000 | 0.512 0.370 0.030 |
## Avena 0.000 | 0.811 0.001 0.000 |
## Cacao 0.302 | -0.604 25.831 0.418 |
## Cebolla.De.Bulbo 0.000 | 1.576 9.599 0.246 |
## Durazno 0.000 | -0.566 0.406 0.012 |
## Frijol 0.003 | 1.244 19.072 0.472 |
## Haba 0.000 | -0.666 0.016 0.008 |
## Habichuela 0.014 | 0.655 0.388 0.038 |
## Lechuga 0.000 | -1.717 0.330 0.013 |
res$eig
## eigenvalue percentage of variance cumulative percentage of variance
## dim 1 0.7558015961 19.57240577 19.57241
## dim 2 0.6094010012 15.78118349 35.35359
## dim 3 0.5285142740 13.68652287 49.04011
## dim 4 0.3519951913 9.11534555 58.15546
## dim 5 0.2969107203 7.68886587 65.84432
## dim 6 0.2481422081 6.42594566 72.27027
## dim 7 0.1913562922 4.95540499 77.22567
## dim 8 0.1594588306 4.12938125 81.35506
## dim 9 0.1314853705 3.40497432 84.76003
## dim 10 0.1153322119 2.98666854 87.74670
## dim 11 0.0765693565 1.98285704 89.72956
## dim 12 0.0702440668 1.81905593 91.54861
## dim 13 0.0632314944 1.63745680 93.18607
## dim 14 0.0602219547 1.55952110 94.74559
## dim 15 0.0436709497 1.13091260 95.87650
## dim 16 0.0420617102 1.08923938 96.96574
## dim 17 0.0352293866 0.91230801 97.87805
## dim 18 0.0254082138 0.65797674 98.53603
## dim 19 0.0245871375 0.63671397 99.17274
## dim 20 0.0219835416 0.56929067 99.74203
## dim 21 0.0056264844 0.14570469 99.88774
## dim 22 0.0039084734 0.10121469 99.98895
## dim 23 0.0004267064 0.01105008 100.00000
res$var
## NULL
res$ind
## NULL
res$ind.sup
## NULL
dimdesc(res, axes=1:3)
## $`Dim 1`
## $`Dim 1`$row
## coord
## SAN VICENTE DE CHUCURI -0.4936493662
## LA PAZ -0.4715538473
## LEBRIJA -0.4700898464
## EL CARMEN DE CHUCURI -0.4480808976
## GIRON -0.4449635889
## GUAPOTA -0.4129129952
## RIONEGRO -0.3809641072
## CIMITARRA -0.3742263517
## EL PLAYON -0.3702175338
## LANDAZURI -0.3597572923
## LA BELLEZA -0.3547430587
## SANTA HELENA DEL OPON -0.3475540413
## SUCRE -0.3323033377
## BETULIA -0.3222532346
## FLORIAN -0.3220588409
## SARDINATA ns -0.3157751276
## EL TARRA ns -0.3141610503
## SIMACOTA -0.3085567735
## GALAN -0.3031418220
## TEORAMA ns -0.3027816367
## BUCARAMANGA -0.2951037682
## CUCUTA ns -0.2872592816
## BUCARASICA ns -0.2659127451
## VELEZ -0.2625722530
## GUADALUPE -0.2622714585
## LA ESPERANZA ns -0.2537174813
## CONTRATACION -0.2515514774
## CHIMA -0.2452892860
## EL GUACAMAYO -0.2447754001
## SAN BENITO -0.2356601214
## LOURDES ns -0.2292855769
## PALMAS DEL SOCORRO -0.2210202068
## BARRANCABERMEJA -0.2117645843
## CONVENCION ns -0.2111371060
## SOCORRO -0.1951607090
## ARATOCA -0.1926022741
## OCAMONTE -0.1873521440
## EL ZULIA ns -0.1848933882
## EL PEÑON -0.1797541967
## VALLE DE SAN JOSE -0.1750153127
## PUERTO PARRA -0.1726331152
## HACARI ns -0.1707102754
## PUERTO SANTANDER ns -0.1699311273
## TIBU ns -0.1661656726
## LOS SANTOS -0.1655565450
## SAN CAYETANO ns -0.1564378832
## TOLEDO ns -0.1493735605
## CHARALA -0.1469180971
## SUAITA -0.1446521169
## ARBOLEDAS ns -0.1418156651
## DURANIA ns -0.1267506007
## PALMAR -0.1249218698
## CUCUTILLA ns -0.1214521911
## SANTIAGO ns -0.1155849986
## SABANA DE TORRES -0.1107433352
## SALAZAR ns -0.1007290509
## PINCHOTE -0.0834217168
## LOS PATIOS ns -0.0774197882
## EL CARMEN ns -0.0753117827
## BOCHALEMA ns -0.0735829979
## MOGOTES -0.0713601178
## JORDAN -0.0603053513
## OIBA -0.0598006582
## MATANZA -0.0534618250
## COROMORO -0.0503404932
## CAPITANEJO -0.0487425933
## PUERTO WILCHES -0.0261594515
## CEPITA -0.0195543516
## PARAMO -0.0102165561
## GÜEPSA -0.0074054834
## VILLA DEL ROSARIO ns 0.0006951219
## BARBOSA 0.0052561339
## GUAVATA 0.0268949852
## CABRERA 0.0341649857
## VILLANUEVA 0.0391198225
## CHINACOTA ns 0.0452821548
## CACOTA ns 0.0516655403
## BOLIVAR 0.0586602749
## JESUS MARIA 0.0590021179
## GRAMALOTE ns 0.0778581526
## ZAPATOCA 0.0781751195
## SAN GIL 0.1012955842
## HATO 0.1021514490
## LABATECA ns 0.1076349441
## MACARAVITA 0.1079187903
## VILLA CARO ns 0.1755547407
## CURITI 0.1972713954
## SAN CALIXTO ns 0.1979531626
## CONFINES 0.2115298069
## SAN JOSE DE MIRANDA 0.2181164652
## FLORIDABLANCA 0.2343026388
## BARICHARA 0.2435601432
## PIEDECUESTA 0.2499901397
## OCAÑA ns 0.2551875006
## SAN MIGUEL 0.2598706027
## ENCISO 0.2631251568
## SURATA 0.2836304695
## SAN JOAQUIN 0.3165365478
## ABREGO ns 0.3513637903
## CHIPATA 0.3587397394
## LA PLAYA ns 0.3695875330
## MALAGA 0.3717854696
## RAGONVALIA ns 0.4180761213
## AGUADA 0.4213717301
## SANTA BARBARA 0.4447863645
## PUENTE NACIONAL 0.4642884999
## MOLAGAVITA 0.4692006236
## CHARTA 0.4748319828
## SAN ANDRES 0.7034069537
## GAMBITA 0.8406076574
## CALIFORNIA 1.0828604910
## PAMPLONITA ns 1.1881396688
## HERRAN ns 1.5154562687
## ONZAGA 1.5382132741
## TONA 1.7350142894
## GUACA 2.1521689747
## ENCINO 2.2387655261
## CONCEPCION 2.3513009428
## PAMPLONA ns 2.6554354414
## CARCASI 2.7499573046
## ALBANIA 2.9966523590
## CACHIRA ns 3.1188894854
## CHITAGA ns 3.2730488979
## CERRITO 3.8436708824
## SILOS ns 3.9724216754
## MUTISCUA ns 4.0127575730
## VETAS 4.1215548177
##
## $`Dim 1`$col
## coord
## Piña -0.48941427
## Cacao -0.44059620
## Limón -0.35315570
## Aptitud -0.18852039
## Aguacate -0.18149069
## Plátano -0.16095670
## Yuca -0.15516311
## Tabaco -0.13067228
## Habichuela 0.02921809
## Maíz 0.07580197
## Tomate 0.15106205
## Mora 0.39754682
## Frijol 0.43800945
## Cebolla.De.Bulbo 0.54690598
## Trigo 1.00662105
## Avena 1.24557134
## Arveja 1.28151387
## Repollo 2.47577427
## Durazno 2.96217693
## Haba 3.24638286
## Papa 3.58314995
## Zanahoria 3.77432449
## Lechuga 4.17892710
## Remolacha 4.26520028
##
##
## $`Dim 2`
## $`Dim 2`$row
## coord
## SAN VICENTE DE CHUCURI -0.641243974
## LA PAZ -0.605789033
## LA BELLEZA -0.584154886
## EL CARMEN DE CHUCURI -0.562724873
## SUCRE -0.559960206
## LANDAZURI -0.529325157
## GUAPOTA -0.524779917
## FLORIAN -0.513490619
## CIMITARRA -0.459998153
## SANTA HELENA DEL OPON -0.453427598
## EL PEÑON -0.422169032
## EL PLAYON -0.415734353
## VELEZ -0.412758856
## BOLIVAR -0.408394718
## GALAN -0.388013175
## BUCARASICA ns -0.386091923
## SARDINATA ns -0.378231307
## EL TARRA ns -0.353661236
## EL GUACAMAYO -0.341778719
## CACOTA ns -0.331618076
## HACARI ns -0.320598076
## GUADALUPE -0.319956117
## TEORAMA ns -0.314095756
## BETULIA -0.313904809
## CONTRATACION -0.310666279
## LOURDES ns -0.305910571
## PALMAS DEL SOCORRO -0.299497933
## TOLEDO ns -0.293567153
## CHIMA -0.289483791
## CONVENCION ns -0.270011194
## VALLE DE SAN JOSE -0.262934986
## LA ESPERANZA ns -0.261898011
## OCAMONTE -0.248301235
## CHARTA -0.247294588
## EL ZULIA ns -0.237435917
## ZAPATOCA -0.229580722
## FLORIDABLANCA -0.226294112
## CUCUTILLA ns -0.195461210
## EL CARMEN ns -0.177917714
## LABATECA ns -0.170400196
## TIBU ns -0.168071115
## PUERTO WILCHES -0.161345004
## DURANIA ns -0.158405838
## RAGONVALIA ns -0.144582021
## ARBOLEDAS ns -0.137822315
## CHARALA -0.130130107
## SAN CAYETANO ns -0.124802728
## TONA -0.122425993
## CHINACOTA ns -0.116589902
## SOCORRO -0.111112555
## MATANZA -0.110979093
## RIONEGRO -0.102120474
## SANTA BARBARA -0.101101880
## SANTIAGO ns -0.097736956
## PUERTO SANTANDER ns -0.088181323
## VILLA CARO ns -0.071321733
## CUCUTA ns -0.070351617
## CALIFORNIA -0.069843335
## PARAMO -0.063892023
## SIMACOTA -0.061740505
## COROMORO -0.059246989
## SAN CALIXTO ns -0.055102141
## OIBA -0.047049678
## JESUS MARIA -0.039163544
## PAMPLONA ns -0.035394697
## BOCHALEMA ns -0.032349886
## HERRAN ns -0.032147068
## MUTISCUA ns -0.021640173
## GUACA -0.021138262
## PALMAR -0.017284984
## GAMBITA -0.011127948
## GUAVATA -0.010203067
## VETAS -0.008769813
## CONFINES -0.007838066
## ENCINO -0.006768738
## SILOS ns -0.006523175
## CERRITO -0.001695313
## SURATA 0.001265788
## CACHIRA ns 0.002847406
## CHITAGA ns 0.009805008
## PIEDECUESTA 0.013676821
## PINCHOTE 0.016176330
## PUENTE NACIONAL 0.017385677
## CARCASI 0.030037840
## SAN ANDRES 0.036583380
## AGUADA 0.036923215
## GÜEPSA 0.042005443
## ALBANIA 0.047233836
## ABREGO ns 0.049353715
## CONCEPCION 0.051121933
## SALAZAR ns 0.060471480
## BUCARAMANGA 0.062850960
## SABANA DE TORRES 0.069893231
## CHIPATA 0.074339458
## LA PLAYA ns 0.085823819
## PAMPLONITA ns 0.090191752
## CABRERA 0.093457135
## OCAÑA ns 0.098373575
## VILLA DEL ROSARIO ns 0.100929910
## SAN MIGUEL 0.106640556
## ONZAGA 0.126377718
## CURITI 0.144479066
## PUERTO PARRA 0.152776694
## GRAMALOTE ns 0.168255049
## SAN JOSE DE MIRANDA 0.168919089
## HATO 0.184771891
## ENCISO 0.188343392
## MALAGA 0.222713261
## BARBOSA 0.228403212
## MOLAGAVITA 0.236741550
## SAN JOAQUIN 0.237341173
## SAN GIL 0.383453496
## BARICHARA 0.407033709
## SUAITA 0.442758969
## BARRANCABERMEJA 0.471718998
## MOGOTES 0.531276059
## LOS PATIOS ns 0.547000051
## CAPITANEJO 0.562771191
## VILLANUEVA 0.873081888
## MACARAVITA 0.899576836
## SAN BENITO 1.113631888
## CEPITA 1.188746878
## ARATOCA 1.452090068
## JORDAN 1.716332740
## LOS SANTOS 2.237055194
## LEBRIJA 2.254022056
## GIRON 2.436701572
##
## $`Dim 2`$col
## coord
## Cacao -0.513772304
## Aguacate -0.239558284
## Mora -0.229639997
## Plátano -0.228588207
## Avena -0.089469165
## Lechuga -0.049929742
## Haba -0.028009665
## Arveja -0.022857627
## Yuca -0.020249830
## Remolacha -0.012433191
## Papa -0.006846079
## Zanahoria -0.004401635
## Repollo 0.007603457
## Durazno 0.009756036
## Cebolla.De.Bulbo 0.018806623
## Maíz 0.038645718
## Aptitud 0.058290850
## Frijol 0.106337605
## Trigo 0.220008037
## Tomate 0.246871925
## Habichuela 0.393295499
## Limón 1.899828719
## Tabaco 2.032334895
## Piña 2.521073594
## Animales
row.names(dat.CA2)<-dat$Municipio
res<-CA(dat.CA2, ncp=3, row.sup=NULL, col.sup=NULL, graph = T)
get_eigenvalue(res)
## eigenvalue variance.percent cumulative.variance.percent
## Dim.1 0.14353268 49.18170 49.18170
## Dim.2 0.11262515 38.59118 87.77288
## Dim.3 0.03568383 12.22712 100.00000
fviz_eig(res, addlabels = TRUE, ylim = c(0, 50))
get_ca_row(res)
## Correspondence Analysis - Results for rows
## ===================================================
## Name Description
## 1 "$coord" "Coordinates for the rows"
## 2 "$cos2" "Cos2 for the rows"
## 3 "$contrib" "contributions of the rows"
## 4 "$inertia" "Inertia of the rows"
get_ca_col(res)
## Correspondence Analysis - Results for columns
## ===================================================
## Name Description
## 1 "$coord" "Coordinates for the columns"
## 2 "$cos2" "Cos2 for the columns"
## 3 "$contrib" "contributions of the columns"
## 4 "$inertia" "Inertia of the columns"
fviz_ca_row(res, repel=T)
## Warning: ggrepel: 113 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
fviz_ca_col(res, repel =T)
fviz_ca_biplot(res, repel = T)
## Warning: ggrepel: 115 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
res.hcpc<-HCPC(res ,nb.clust=-1,consol=TRUE,min=3,max=10,graph=TRUE)
res.animales<-res.hcpc$data.clust[,ncol(res.hcpc$data.clust),drop=F]
res.animales
## clust
## PUERTO WILCHES 1
## BARRANCABERMEJA 1
## SABANA DE TORRES 1
## RIONEGRO 2
## BOLIVAR 2
## CIMITARRA 2
## TEORAMA ns 2
## PUERTO PARRA 2
## TIBU ns 2
## SIMACOTA 2
## BOCHALEMA ns 2
## EL TARRA ns 2
## SUAITA 2
## LA ESPERANZA ns 2
## GUAVATA 2
## LOURDES ns 2
## EL GUACAMAYO 2
## SARDINATA ns 2
## CONVENCION ns 2
## VELEZ 2
## SAN CALIXTO ns 2
## CACOTA ns 2
## PALMAS DEL SOCORRO 2
## LA PAZ 2
## SANTA HELENA DEL OPON 2
## GRAMALOTE ns 2
## OIBA 2
## LA PLAYA ns 2
## CUCUTILLA ns 2
## EL PEÑON 2
## LA BELLEZA 2
## GAMBITA 2
## CHITAGA ns 2
## EL CARMEN ns 2
## SUCRE 2
## HACARI ns 2
## PUENTE NACIONAL 2
## GUADALUPE 2
## CHIPATA 2
## MUTISCUA ns 2
## LANDAZURI 2
## BUCARASICA ns 2
## ALBANIA 2
## FLORIAN 2
## GUAPOTA 2
## LABATECA ns 2
## SANTIAGO ns 2
## SAN VICENTE DE CHUCURI 2
## CUCUTA ns 2
## BARBOSA 2
## OCAÑA ns 2
## SILOS ns 2
## PAMPLONITA ns 2
## SALAZAR ns 2
## HERRAN ns 2
## DURANIA ns 2
## SAN BENITO 2
## PALMAR 2
## CHARTA 2
## RAGONVALIA ns 2
## EL ZULIA ns 2
## CHIMA 2
## JESUS MARIA 2
## CHARALA 2
## SOCORRO 2
## TOLEDO ns 2
## PAMPLONA ns 2
## SAN CAYETANO ns 2
## EL CARMEN DE CHUCURI 2
## CHINACOTA ns 2
## AGUADA 2
## PARAMO 2
## LEBRIJA 2
## VALLE DE SAN JOSE 2
## BETULIA 2
## CONFINES 2
## ARBOLEDAS ns 2
## HATO 2
## CONTRATACION 2
## EL PLAYON 2
## CACHIRA ns 2
## MATANZA 2
## ABREGO ns 2
## CABRERA 2
## VILLA DEL ROSARIO ns 2
## COROMORO 2
## TONA 2
## OCAMONTE 2
## SAN ANDRES 2
## SANTA BARBARA 2
## GIRON 2
## ONZAGA 2
## MOGOTES 2
## VILLA CARO ns 2
## GUACA 2
## PIEDECUESTA 2
## BUCARAMANGA 2
## MALAGA 2
## ENCINO 2
## SAN GIL 2
## VETAS 2
## CALIFORNIA 2
## SURATA 2
## GÜEPSA 2
## PUERTO SANTANDER ns 2
## PINCHOTE 2
## GALAN 2
## CURITI 2
## SAN MIGUEL 2
## ENCISO 2
## CERRITO 2
## CARCASI 2
## CONCEPCION 2
## LOS SANTOS 3
## BARICHARA 3
## SAN JOSE DE MIRANDA 3
## ZAPATOCA 3
## FLORIDABLANCA 3
## VILLANUEVA 3
## SAN JOAQUIN 3
## MACARAVITA 3
## LOS PATIOS ns 3
## JORDAN 3
## ARATOCA 3
## MOLAGAVITA 3
## CEPITA 3
## CAPITANEJO 3
write.csv(res.animales, "~/Library/CloudStorage/GoogleDrive-icarounam@gmail.com/Mi unidad/Agrosavia/Ingreso/Evas/resanimales.csv")
res.hcpc$desc.var
## $`1`
## Intern % glob % Intern freq Glob freq p.value v.test
## Bufalos 17.0802287 3.045130 48872 70506 0.000000e+00 Inf
## Ovinos 1.8362155 2.568100 5254 59461 4.030372e-168 -27.63672
## Caprinos 0.5396111 1.558542 1544 36086 0.000000e+00 -Inf
## Bovinos 80.5439448 92.828227 230462 2149316 0.000000e+00 -Inf
##
## $`2`
## Intern % glob % Intern freq Glob freq p.value v.test
## Bovinos 95.5059157 92.828227 1856210 2149316 0.000000e+00 Inf
## Ovinos 2.4410423 2.568100 47443 59461 1.288318e-161 -27.09009
## Caprinos 0.9436831 1.558542 18341 36086 0.000000e+00 -Inf
## Bufalos 1.1093589 3.045130 21561 70506 0.000000e+00 -Inf
##
## $`3`
## Intern % glob % Intern freq Glob freq p.value v.test
## Ovinos 7.89430686 2.568100 6764 59461 0 Inf
## Caprinos 18.90828879 1.558542 16201 36086 0 Inf
## Bufalos 0.08519876 3.045130 73 70506 0 -Inf
## Bovinos 73.11220560 92.828227 62644 2149316 0 -Inf
##
## attr(,"class")
## [1] "descfreq" "list"
res.hcpc$desc.axes
##
## Link between the cluster variable and the quantitative variables
## ================================================================
## Eta2 P-value
## Dim.2 0.77294921 1.201233e-40
## Dim.1 0.76188096 2.297395e-39
## Dim.3 0.05626283 2.759088e-02
##
## Description of each cluster by quantitative variables
## =====================================================
## $`1`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.2 199.311040 0.54480302 -0.04912071 0.30260686 0.3543282
## Dim.3 -3.928568 0.01883481 0.02839950 0.01421132 0.2894964
## Dim.1 -250.130412 -0.61964167 0.17325179 0.30083392 0.3769248
## p.value
## Dim.2 0.000000e+00
## Dim.3 8.545319e-05
## Dim.1 0.000000e+00
##
## $`2`
## v.test Mean in category Overall mean sd in category Overall sd
## Dim.3 348.9233 0.04788144 0.02839950 0.2418360 0.2894964
## Dim.1 -1190.7087 0.08669142 0.17325179 0.1366078 0.3769248
## Dim.2 -1328.1649 -0.13988531 -0.04912071 0.1158318 0.3543282
## p.value
## Dim.3 0
## Dim.1 0
## Dim.2 0
##
## $`3`
## v.test Mean in category Overall mean sd in category Overall sd p.value
## Dim.2 1318.1546 1.0428503 -0.04912071 0.4615122 0.3543282 0
## Dim.1 1309.5835 1.3273080 0.17325179 0.4688796 0.3769248 0
## Dim.3 -360.7334 -0.2157572 0.02839950 0.5871477 0.2894964 0
res.hcpc$desc.ind
## $para
## Cluster: 1
## BARRANCABERMEJA SABANA DE TORRES PUERTO WILCHES
## 0.2172066 0.4686051 0.6848756
## ------------------------------------------------------------
## Cluster: 2
## EL PLAYON ARBOLEDAS ns CACHIRA ns BETULIA
## 0.01824500 0.02891758 0.03512125 0.05122155
## VALLE DE SAN JOSE
## 0.05364961
## ------------------------------------------------------------
## Cluster: 3
## ARATOCA JORDAN MACARAVITA SAN JOAQUIN ZAPATOCA
## 0.3405600 0.3655742 0.4577093 0.5253721 0.6065316
##
## $dist
## Cluster: 1
## PUERTO WILCHES BARRANCABERMEJA SABANA DE TORRES
## 1.7634754 0.8655183 0.6121993
## ------------------------------------------------------------
## Cluster: 2
## PINCHOTE GÜEPSA SAN GIL COROMORO CALIFORNIA
## 1.637935 1.464513 1.293792 1.210574 1.206176
## ------------------------------------------------------------
## Cluster: 3
## CAPITANEJO CEPITA MOLAGAVITA LOS PATIOS ns ARATOCA
## 4.139599 3.277644 2.179185 1.875940 1.832832
summary(res, nb.dec = 3, nbelements=10, nbind = 10, ncp = 3, file="")
##
## Call:
## CA(X = dat.CA2, ncp = 3, row.sup = NULL, col.sup = NULL, graph = T)
##
## The chi square of independence between the two variables is equal to 675721.1 (p-value = 0 ).
##
## Eigenvalues
## Dim.1 Dim.2 Dim.3
## Variance 0.144 0.113 0.036
## % of var. 49.182 38.591 12.227
## Cumulative % of var. 49.182 87.773 100.000
##
## Rows (the 10 first)
## Iner*1000 Dim.1 ctr cos2 Dim.2 ctr
## AGUADA | 0.079 | 0.081 0.006 0.110 | -0.172 0.035
## ALBANIA | 0.151 | 0.025 0.001 0.012 | -0.218 0.123
## ARATOCA | 6.206 | 1.444 2.616 0.605 | 1.155 2.134
## BARBOSA | 0.032 | 0.042 0.001 0.066 | -0.153 0.024
## BARICHARA | 2.615 | 0.751 1.301 0.714 | 0.471 0.651
## BARRANCABERMEJA | 16.473 | -0.536 6.687 0.583 | 0.454 6.104
## BETULIA | 0.193 | 0.109 0.050 0.371 | -0.141 0.106
## BOLIVAR | 0.881 | -0.084 0.153 0.250 | -0.131 0.479
## BUCARAMANGA | 0.206 | 0.316 0.083 0.574 | 0.046 0.002
## CABRERA | 0.159 | 0.181 0.075 0.674 | -0.071 0.015
## cos2 Dim.3 ctr cos2
## AGUADA 0.499 | 0.152 0.087 0.391 |
## ALBANIA 0.921 | -0.058 0.028 0.066 |
## ARATOCA 0.387 | 0.163 0.134 0.008 |
## BARBOSA 0.865 | 0.043 0.006 0.069 |
## BARICHARA 0.281 | 0.065 0.039 0.005 |
## BARRANCABERMEJA 0.417 | 0.004 0.002 0.000 |
## BETULIA 0.619 | 0.019 0.006 0.011 |
## BOLIVAR 0.613 | -0.062 0.339 0.137 |
## BUCARAMANGA 0.012 | 0.269 0.240 0.414 |
## CABRERA 0.103 | 0.105 0.100 0.223 |
##
## Columns
## Iner*1000 Dim.1 ctr cos2 Dim.2 ctr
## Bovinos | 7.041 | -0.004 0.010 0.002 | -0.084 5.811
## Bufalos | 123.441 | -1.478 46.323 0.539 | 1.367 50.539
## Caprinos | 116.183 | 2.066 46.327 0.572 | 1.709 40.410
## Ovinos | 45.176 | 0.641 7.340 0.233 | 0.377 3.240
## cos2 Dim.3 ctr cos2
## Bovinos 0.930 | -0.023 1.351 0.068 |
## Bufalos 0.461 | 0.033 0.093 0.000 |
## Caprinos 0.392 | -0.518 11.704 0.036 |
## Ovinos 0.081 | 1.099 86.852 0.686 |
res$eig
## eigenvalue percentage of variance cumulative percentage of variance
## dim 1 0.14353268 49.18170 49.18170
## dim 2 0.11262515 38.59118 87.77288
## dim 3 0.03568383 12.22712 100.00000
res$var
## NULL
res$ind
## NULL
res$ind.sup
## NULL
dimdesc(res, axes=1:3)
## $`Dim 1`
## $`Dim 1`$row
## coord
## PUERTO WILCHES -1.166479353
## BARRANCABERMEJA -0.536472175
## SABANA DE TORRES -0.350564426
## RIONEGRO -0.137928800
## BOLIVAR -0.083879304
## CIMITARRA -0.082618831
## TEORAMA ns -0.077900266
## PUERTO PARRA -0.049147656
## TIBU ns -0.048527838
## SIMACOTA -0.043078328
## BOCHALEMA ns -0.023294347
## EL TARRA ns -0.017640588
## SUAITA -0.013883937
## LA ESPERANZA ns -0.009416844
## GUAVATA -0.006806679
## LOURDES ns -0.003520338
## EL GUACAMAYO -0.002303267
## SARDINATA ns 0.002718703
## CONVENCION ns 0.003012685
## VELEZ 0.003153867
## SAN CALIXTO ns 0.003911734
## CACOTA ns 0.004389185
## PALMAS DEL SOCORRO 0.005221146
## LA PAZ 0.005549573
## SANTA HELENA DEL OPON 0.006225872
## GRAMALOTE ns 0.006369901
## OIBA 0.009282297
## LA PLAYA ns 0.009859854
## CUCUTILLA ns 0.011205574
## EL PEÑON 0.011558214
## LA BELLEZA 0.012947432
## GAMBITA 0.015953678
## CHITAGA ns 0.017500619
## EL CARMEN ns 0.018196756
## SUCRE 0.019126917
## HACARI ns 0.019358096
## PUENTE NACIONAL 0.019869404
## GUADALUPE 0.020415710
## CHIPATA 0.020623963
## MUTISCUA ns 0.022068762
## LANDAZURI 0.022578832
## BUCARASICA ns 0.024520769
## ALBANIA 0.025162315
## FLORIAN 0.028219878
## GUAPOTA 0.028848328
## LABATECA ns 0.028996829
## SANTIAGO ns 0.036201325
## SAN VICENTE DE CHUCURI 0.039195614
## CUCUTA ns 0.040334592
## BARBOSA 0.042328675
## OCAÑA ns 0.045807000
## SILOS ns 0.046359602
## PAMPLONITA ns 0.046442178
## SALAZAR ns 0.048854860
## HERRAN ns 0.049271811
## DURANIA ns 0.049429733
## SAN BENITO 0.051128474
## PALMAR 0.052001822
## CHARTA 0.052355612
## RAGONVALIA ns 0.054579664
## EL ZULIA ns 0.056352550
## CHIMA 0.057206748
## JESUS MARIA 0.060745050
## CHARALA 0.063195955
## SOCORRO 0.063830976
## TOLEDO ns 0.067793302
## PAMPLONA ns 0.072703329
## SAN CAYETANO ns 0.073963355
## EL CARMEN DE CHUCURI 0.074714147
## CHINACOTA ns 0.074737265
## AGUADA 0.080688754
## PARAMO 0.089977312
## LEBRIJA 0.096367907
## VALLE DE SAN JOSE 0.096634566
## BETULIA 0.109471315
## CONFINES 0.110497223
## ARBOLEDAS ns 0.116453266
## HATO 0.116864331
## CONTRATACION 0.123611068
## EL PLAYON 0.124831131
## CACHIRA ns 0.139009282
## MATANZA 0.152868521
## ABREGO ns 0.172566897
## CABRERA 0.181450396
## VILLA DEL ROSARIO ns 0.186695092
## COROMORO 0.187013725
## TONA 0.189317326
## OCAMONTE 0.198894891
## SAN ANDRES 0.219492965
## SANTA BARBARA 0.238031053
## GIRON 0.241330094
## ONZAGA 0.246029855
## MOGOTES 0.261673095
## VILLA CARO ns 0.274983496
## GUACA 0.290978753
## PIEDECUESTA 0.294753171
## BUCARAMANGA 0.316008365
## MALAGA 0.316922484
## ENCINO 0.318249476
## SAN GIL 0.338925642
## VETAS 0.349444269
## CALIFORNIA 0.350558929
## SURATA 0.383462084
## GÜEPSA 0.383884606
## PUERTO SANTANDER ns 0.419199305
## PINCHOTE 0.450979480
## GALAN 0.468232588
## CURITI 0.491629921
## SAN MIGUEL 0.496806476
## ENCISO 0.512792792
## CERRITO 0.604269260
## CARCASI 0.629748914
## CONCEPCION 0.657817612
## LOS SANTOS 0.733423574
## BARICHARA 0.750883635
## SAN JOSE DE MIRANDA 0.804786119
## ZAPATOCA 0.837944779
## FLORIDABLANCA 0.921303803
## VILLANUEVA 0.939140721
## SAN JOAQUIN 0.969049832
## MACARAVITA 0.988962088
## LOS PATIOS ns 1.200127434
## JORDAN 1.401636218
## ARATOCA 1.443809202
## MOLAGAVITA 1.555264394
## CEPITA 2.336587553
## CAPITANEJO 2.896405308
##
## $`Dim 1`$col
## coord
## Bufalos -1.477644199
## Bovinos -0.003926439
## Ovinos 0.640500612
## Caprinos 2.065541557
##
##
## $`Dim 2`
## $`Dim 2`$row
## coord
## CONVENCION ns -0.239403582
## CACOTA ns -0.238292208
## SAN CALIXTO ns -0.238259709
## GRAMALOTE ns -0.236692996
## SARDINATA ns -0.236286349
## LOURDES ns -0.233683508
## EL PEÑON -0.231594902
## LA PLAYA ns -0.231380946
## VELEZ -0.229161819
## SANTA HELENA DEL OPON -0.228184075
## EL CARMEN ns -0.226092385
## LA BELLEZA -0.225040348
## CHITAGA ns -0.224893710
## SUCRE -0.223292736
## GUADALUPE -0.223112213
## CHIPATA -0.222475364
## BUCARASICA ns -0.222038148
## OIBA -0.221998812
## MUTISCUA ns -0.221383388
## LANDAZURI -0.221061851
## ALBANIA -0.217679250
## CUCUTILLA ns -0.217286383
## GAMBITA -0.217061913
## LABATECA ns -0.216343684
## FLORIAN -0.215117632
## EL TARRA ns -0.214600588
## SANTIAGO ns -0.211669840
## EL GUACAMAYO -0.209009957
## LA PAZ -0.205525894
## SILOS ns -0.202101823
## PUENTE NACIONAL -0.200998109
## GUAVATA -0.198496065
## SAN BENITO -0.196897923
## HACARI ns -0.196090228
## OCAÑA ns -0.196054826
## CHIMA -0.194902257
## DURANIA ns -0.194204550
## RAGONVALIA ns -0.193885114
## PALMAS DEL SOCORRO -0.193635803
## GUAPOTA -0.192972904
## PALMAR -0.191539495
## LA ESPERANZA ns -0.191050689
## PUERTO PARRA -0.187367770
## HERRAN ns -0.185241250
## EL ZULIA ns -0.185204622
## SALAZAR ns -0.183956899
## TOLEDO ns -0.180715306
## SOCORRO -0.179019752
## PAMPLONITA ns -0.175267415
## TEORAMA ns -0.173714408
## CUCUTA ns -0.173331329
## SAN CAYETANO ns -0.172516317
## AGUADA -0.172155857
## EL CARMEN DE CHUCURI -0.170726946
## CHARALA -0.166293315
## PARAMO -0.158995916
## JESUS MARIA -0.157335954
## VALLE DE SAN JOSE -0.155685136
## BARBOSA -0.153163910
## SUAITA -0.152588023
## LEBRIJA -0.152230897
## CHINACOTA ns -0.150830757
## CIMITARRA -0.149821065
## PAMPLONA ns -0.149620811
## CONFINES -0.144152405
## CHARTA -0.142633550
## BETULIA -0.141459142
## SIMACOTA -0.141290343
## CONTRATACION -0.138014094
## ARBOLEDAS ns -0.136498287
## BOLIVAR -0.131328584
## BOCHALEMA ns -0.117205799
## CACHIRA ns -0.113953482
## EL PLAYON -0.112114719
## TIBU ns -0.111887063
## HATO -0.102741877
## COROMORO -0.087440434
## VILLA DEL ROSARIO ns -0.075762211
## TONA -0.073633673
## CABRERA -0.070905056
## ABREGO ns -0.067739702
## SAN VICENTE DE CHUCURI -0.063516448
## OCAMONTE -0.061176676
## MATANZA -0.047065677
## SAN ANDRES -0.033640573
## SANTA BARBARA -0.018604490
## ONZAGA -0.011903882
## RIONEGRO -0.010016882
## MOGOTES 0.003865738
## VILLA CARO ns 0.005534458
## GUACA 0.024961280
## GIRON 0.039067062
## BUCARAMANGA 0.045534604
## MALAGA 0.054575419
## SAN GIL 0.055556977
## ENCINO 0.060410818
## CALIFORNIA 0.075824040
## VETAS 0.076659908
## GÜEPSA 0.084458220
## PIEDECUESTA 0.110865933
## SURATA 0.112592704
## PUERTO SANTANDER ns 0.169666504
## PINCHOTE 0.173809381
## GALAN 0.196278572
## CURITI 0.224907457
## SAN MIGUEL 0.239455820
## ENCISO 0.251571277
## SABANA DE TORRES 0.283189453
## CERRITO 0.306406016
## CARCASI 0.347398315
## CONCEPCION 0.352526391
## BARRANCABERMEJA 0.454034117
## BARICHARA 0.470639550
## LOS SANTOS 0.489501141
## SAN JOSE DE MIRANDA 0.537941563
## ZAPATOCA 0.563435401
## FLORIDABLANCA 0.620264876
## VILLANUEVA 0.642475573
## SAN JOAQUIN 0.699413795
## MACARAVITA 0.714979052
## LOS PATIOS ns 0.852254802
## JORDAN 1.085611128
## PUERTO WILCHES 1.098288308
## ARATOCA 1.154933845
## MOLAGAVITA 1.279335385
## CEPITA 2.036786327
## CAPITANEJO 2.588228547
##
## $`Dim 2`$col
## coord
## Bovinos -0.08396755
## Ovinos 0.37692532
## Bufalos 1.36718359
## Caprinos 1.70885626