
TeorÃa
Agrupamiento o Clustering es una
técnica de aprendizaje automático no supervisado que agrupa datos en
función de su similitud.
Algunos ustos tÃpicos de esta técnica son:
- Segmentación de clientes
- Detección de anormalidades
- Categorización de documentos
Instalar paquetes y llamar
librerÃas
#install.packages("cluster") #Análisis de Agrupamiento
library(cluster)
#install.packages("ggplot2") #Graficar
library(ggplot2)
#install.packages("data.table") #Manejo de muchos datos
library(data.table)
#install.packages("factoextra") #Gráfica optimizazción de número de clusters
library(factoextra)
## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa
Paso2. Obtener los datos
df1 <- data.frame(x=c(2,2,8,5,7,6,1,4), y=c(10,5,4,8,5,4,2,9))
Paso 3. Entender los datos
summary(df1)
## x y
## Min. :1.000 Min. : 2.000
## 1st Qu.:2.000 1st Qu.: 4.000
## Median :4.500 Median : 5.000
## Mean :4.375 Mean : 5.875
## 3rd Qu.:6.250 3rd Qu.: 8.250
## Max. :8.000 Max. :10.000
str(df1)
## 'data.frame': 8 obs. of 2 variables:
## $ x: num 2 2 8 5 7 6 1 4
## $ y: num 10 5 4 8 5 4 2 9
Paso 4. Escalar los datos
# Sólo si los datos no están en la misma escala
# datos_escalados <- scale(datos_originales)
Paso 5. Determinar el número de
grupos
# Siempre es un valor inical "cualquiera", luego se optimiza.
plot(df1)

grupos1 <-3
Paso 6. Generar los grupos
set.seed(123)
clusters1 <- kmeans(df1,grupos1)
clusters1
## K-means clustering with 3 clusters of sizes 2, 3, 3
##
## Cluster means:
## x y
## 1 1.500000 3.500000
## 2 3.666667 9.000000
## 3 7.000000 4.333333
##
## Clustering vector:
## [1] 2 1 3 2 3 3 1 2
##
## Within cluster sum of squares by cluster:
## [1] 5.000000 6.666667 2.666667
## (between_SS / total_SS = 85.8 %)
##
## Available components:
##
## [1] "cluster" "centers" "totss" "withinss" "tot.withinss"
## [6] "betweenss" "size" "iter" "ifault"
Paso 7. Optimizar el número de
grupos
set.seed(123)
optimizacion <- clusGap(df1,FUN=kmeans, nstart=1,K.max=7)
#el K.max normalmente es 10, en este ejercicio al ser 8 datos se dejó en 7
plot(optimizacion, xlab="Número de Clusters k")

# Se selecciona como óptimo el primer putno más alto
Paso 8. Graficar los grupos
fviz_cluster(clusters1,data=df1)

Paso 9. Agregar clusters a la base
de datos
df1_clusters <- cbind(df1,cluster=clusters1$cluster)
head(df1_clusters)
## x y cluster
## 1 2 10 2
## 2 2 5 1
## 3 8 4 3
## 4 5 8 2
## 5 7 5 3
## 6 6 4 3
Conclusiones
La técnica de clustering permite identificar patrones o
grupos naturales en los datos sin necesidad de etiquetas previas.
LS0tCnRpdGxlOiAiU2VnbWVudGFjacOzbiAoY2x1c3RlcnMpIgphdXRob3I6ICJBMDE3Mzc2NDQiCmRhdGU6ICIyMDI1LTA4LTE4IgpvdXRwdXQ6IAogIGh0bWxfZG9jdW1lbnQ6CiAgICB0b2M6IFRSVUUKICAgIHRvY19mbG9hdDogVFJVRQogICAgY29kZV9kb3dubG9hZDogVFJVRQogICAgdGhlbWU6IHlldGkKLS0tCgohW10oaHR0cHM6Ly93d3cuY3MudXMuZXMvfmZzYW5jaG8vQmxvZy9wb3N0cy9DbHVzdGVyaW5nL2ltZy9tZWFuc2hpZnQyLmdpZikKCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVldmlvbGV0OyI+IFRlb3LDrWEgPC9zcGFuPgoKKipBZ3J1cGFtaWVudG8qKiBvICoqQ2x1c3RlcmluZyoqIGVzIHVuYSB0w6ljbmljYSBkZSBhcHJlbmRpemFqZSBhdXRvbcOhdGljbyBubyBzdXBlcnZpc2FkbyBxdWUgYWdydXBhIGRhdG9zIGVuIGZ1bmNpw7NuIGRlIHN1IHNpbWlsaXR1ZC4KCkFsZ3Vub3MgdXN0b3MgdMOtcGljb3MgZGUgZXN0YSB0w6ljbmljYSBzb246CgoqIFNlZ21lbnRhY2nDs24gZGUgY2xpZW50ZXMKKiBEZXRlY2Npw7NuIGRlIGFub3JtYWxpZGFkZXMKKiBDYXRlZ29yaXphY2nDs24gZGUgZG9jdW1lbnRvcwoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZXZpb2xldDsiPiBJbnN0YWxhciBwYXF1ZXRlcyB5IGxsYW1hciBsaWJyZXLDrWFzIDwvc3Bhbj4KCjwvc3Bhbj4KCmBgYHtyIHdhcm5pbmc9RkFMU0V9CiNpbnN0YWxsLnBhY2thZ2VzKCJjbHVzdGVyIikgICNBbsOhbGlzaXMgZGUgQWdydXBhbWllbnRvCmxpYnJhcnkoY2x1c3RlcikKI2luc3RhbGwucGFja2FnZXMoImdncGxvdDIiKSAjR3JhZmljYXIKbGlicmFyeShnZ3Bsb3QyKQojaW5zdGFsbC5wYWNrYWdlcygiZGF0YS50YWJsZSIpICNNYW5lam8gZGUgbXVjaG9zIGRhdG9zIApsaWJyYXJ5KGRhdGEudGFibGUpCiNpbnN0YWxsLnBhY2thZ2VzKCJmYWN0b2V4dHJhIikgI0dyw6FmaWNhIG9wdGltaXphemNpw7NuIGRlIG7Dum1lcm8gZGUgY2x1c3RlcnMKbGlicmFyeShmYWN0b2V4dHJhKQpgYGAKCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVldmlvbGV0OyI+IFBhc28yLiBPYnRlbmVyIGxvcyBkYXRvcyA8L3NwYW4+CmBgYHtyfQpkZjEgPC0gZGF0YS5mcmFtZSh4PWMoMiwyLDgsNSw3LDYsMSw0KSwgeT1jKDEwLDUsNCw4LDUsNCwyLDkpKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWV2aW9sZXQ7Ij4gUGFzbyAzLiBFbnRlbmRlciBsb3MgZGF0b3MgPC9zcGFuPgpgYGB7cn0Kc3VtbWFyeShkZjEpCnN0cihkZjEpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZXZpb2xldDsiPiBQYXNvIDQuIEVzY2FsYXIgbG9zIGRhdG9zIDwvc3Bhbj4KYGBge3J9CiMgU8OzbG8gc2kgbG9zIGRhdG9zIG5vIGVzdMOhbiBlbiBsYSBtaXNtYSBlc2NhbGEKIyBkYXRvc19lc2NhbGFkb3MgPC0gc2NhbGUoZGF0b3Nfb3JpZ2luYWxlcykKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVldmlvbGV0OyI+IFBhc28gNS4gRGV0ZXJtaW5hciBlbCBuw7ptZXJvIGRlIGdydXBvcyAgIDwvc3Bhbj4KYGBge3J9CiMgU2llbXByZSBlcyB1biB2YWxvciBpbmljYWwgImN1YWxxdWllcmEiLCBsdWVnbyBzZSBvcHRpbWl6YS4gCnBsb3QoZGYxKQpncnVwb3MxIDwtMwpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWV2aW9sZXQ7Ij4gUGFzbyA2LiBHZW5lcmFyIGxvcyBncnVwb3MgPC9zcGFuPgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpjbHVzdGVyczEgPC0ga21lYW5zKGRmMSxncnVwb3MxKQpjbHVzdGVyczEKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVldmlvbGV0OyI+IFBhc28gNy4gT3B0aW1pemFyIGVsIG7Dum1lcm8gZGUgZ3J1cG9zIDwvc3Bhbj4KYGBge3J9CnNldC5zZWVkKDEyMykKb3B0aW1pemFjaW9uIDwtIGNsdXNHYXAoZGYxLEZVTj1rbWVhbnMsIG5zdGFydD0xLEsubWF4PTcpCiNlbCBLLm1heCBub3JtYWxtZW50ZSBlcyAxMCwgZW4gZXN0ZSBlamVyY2ljaW8gYWwgc2VyIDggZGF0b3Mgc2UgZGVqw7MgZW4gNwpwbG90KG9wdGltaXphY2lvbiwgeGxhYj0iTsO6bWVybyBkZSBDbHVzdGVycyBrIikKIyBTZSBzZWxlY2Npb25hIGNvbW8gw7NwdGltbyBlbCBwcmltZXIgcHV0bm8gbcOhcyBhbHRvCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZXZpb2xldDsiPiBQYXNvIDguICBHcmFmaWNhciBsb3MgZ3J1cG9zIDwvc3Bhbj4KYGBge3J9CmZ2aXpfY2x1c3RlcihjbHVzdGVyczEsZGF0YT1kZjEpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZXZpb2xldDsiPiBQYXNvIDkuIEFncmVnYXIgY2x1c3RlcnMgYSBsYSBiYXNlIGRlIGRhdG9zIDwvc3Bhbj4KYGBge3J9CmRmMV9jbHVzdGVycyA8LSBjYmluZChkZjEsY2x1c3Rlcj1jbHVzdGVyczEkY2x1c3RlcikKaGVhZChkZjFfY2x1c3RlcnMpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZXZpb2xldDsiPiBDb25jbHVzaW9uZXMgPC9zcGFuPgpMYSB0w6ljbmljYSBkZSAqY2x1c3RlcmluZyogcGVybWl0ZSBpZGVudGlmaWNhciBwYXRyb25lcyBvIGdydXBvcyBuYXR1cmFsZXMgZW4gbG9zIGRhdG9zIHNpbiBuZWNlc2lkYWQgZGUgZXRpcXVldGFzIHByZXZpYXMu