Teoría

Agrupamiento o clustering es una técnica de aprendizaje automático que agrupa datos en función de similitud.

Algunos usos típicos de esta técnica son: * Segmentación de clientes * Detección de abnormalidades * Categorización de documentos

Paso 1. Instalar paquetes y llamar librerías

#install.packages("cluster") # Análisis de Agrupamiento
library(cluster)
#install.packages("ggplot2") # Paquete de gráficas
library(ggplot2)
#install.packages("data.table") # Manejo de muchos datos
library(data.table)
#install.packages("factoextra") # Gráfica optimización de número de clusters
library(factoextra)

Paso 2. 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

# Solo si los datos no estan en la misma escala.
# datos_escalados <- scale(datos_originales) (Este es un ejemplo obviamente)

Paso 5. Determinar numero de grupos

# Siempre es un valor inicial "cualquiera", luego se optimiza.
plot(df1$x, df1$y)

grupos1 <- 3

Paso 6. Generar los grupos

clusters1 <- kmeans(df1,grupos1)
clusters1
## K-means clustering with 3 clusters of sizes 3, 2, 3
## 
## Cluster means:
##          x        y
## 1 7.000000 4.333333
## 2 1.500000 3.500000
## 3 3.666667 9.000000
## 
## Clustering vector:
## [1] 3 2 1 3 1 1 2 3
## 
## Within cluster sum of squares by cluster:
## [1] 2.666667 5.000000 6.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 numero de grupos

set.seed(123)
optimizacion1 <- clusGap(df1, FUN=kmeans, nstart=1, K.max=7)
# El K.max normalmente es 10, en este ejercicio al ser 8 datos se dejo en 7.
plot (optimizacion1, xlab="Numero de clusters k")

# Se selecciona como optimo el primer punto mas alto.

Paso 8. Graficar los grupos

fviz_cluster(clusters1,data=df1)

Paso 8. Agregar clusters a la base de datos

df1_clusters <- cbind(df1, cluster = clusters1$cluster)
head(df1_clusters)
##   x  y cluster
## 1 2 10       3
## 2 2  5       2
## 3 8  4       1
## 4 5  8       3
## 5 7  5       1
## 6 6  4       1

Conclusiones

La técnica de clustering permite identificar patrones o grupos naturales en los datos sin necesidad de etiquetas previas.

LS0tCnRpdGxlOiAiQWdydXBhbWllbnRvIChDbHVzdGVycykiCmF1dGhvcjogIlBhYmxvIFNhbmNobyBBMDE3MjIyMzYiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDogCiAgICB0b2M6IFRSVUUKICAgIHRvY19mbG9hdDogVFJVRQogICAgY29kZV9kb3dubG9hZDogVFJVRQogICAgdGhlbWU6IGNvc21vCmRhdGU6ICIyMDI1LTA4LTE4IgotLS0KIVtdKGh0dHBzOi8vd3d3LnF1YWx0cmljcy5jb20vbS9hc3NldHMvd3AtY29udGVudC91cGxvYWRzLzIwMjMvMDkvY2x1c3Rlci1ncmFwaC1lbi03Njh4NTEyLndlYnApCgojIDxzcGFuIHN0eWxlPSJjb2xvcjp0dXJxdW9pc2U7Ij4gVGVvcsOtYSA8L3NwYW4+CioqQWdydXBhbWllbnRvKiogbyAqY2x1c3RlcmluZyogZXMgdW5hIHTDqWNuaWNhIGRlIGFwcmVuZGl6YWplIGF1dG9tw6F0aWNvIHF1ZSBhZ3J1cGEgZGF0b3MgZW4gZnVuY2nDs24gZGUgc2ltaWxpdHVkLiAgCgpBbGd1bm9zIHVzb3MgdMOtcGljb3MgZGUgZXN0YSB0w6ljbmljYSBzb246IAoqIFNlZ21lbnRhY2nDs24gZGUgY2xpZW50ZXMKKiBEZXRlY2Npw7NuIGRlIGFibm9ybWFsaWRhZGVzCiogQ2F0ZWdvcml6YWNpw7NuIGRlIGRvY3VtZW50b3MKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyAxLiBJbnN0YWxhciBwYXF1ZXRlcyB5IGxsYW1hciBsaWJyZXLDrWFzIDwvc3Bhbj4KYGBge3IgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KI2luc3RhbGwucGFja2FnZXMoImNsdXN0ZXIiKSAjIEFuw6FsaXNpcyBkZSBBZ3J1cGFtaWVudG8KbGlicmFyeShjbHVzdGVyKQojaW5zdGFsbC5wYWNrYWdlcygiZ2dwbG90MiIpICMgUGFxdWV0ZSBkZSBncsOhZmljYXMKbGlicmFyeShnZ3Bsb3QyKQojaW5zdGFsbC5wYWNrYWdlcygiZGF0YS50YWJsZSIpICMgTWFuZWpvIGRlIG11Y2hvcyBkYXRvcwpsaWJyYXJ5KGRhdGEudGFibGUpCiNpbnN0YWxsLnBhY2thZ2VzKCJmYWN0b2V4dHJhIikgIyBHcsOhZmljYSBvcHRpbWl6YWNpw7NuIGRlIG7Dum1lcm8gZGUgY2x1c3RlcnMKbGlicmFyeShmYWN0b2V4dHJhKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyAyLiBPYnRlbmVyIGxvcyBkYXRvcyA8L3NwYW4+CmBgYHtyfQpkZjEgPC0gZGF0YS5mcmFtZSh4PWMoMiwyLDgsNSw3LDYsMSw0KSwgeT1jKDEwLDUsNCw4LDUsNCwyLDkpKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyAzLiBFbnRlbmRlciBsb3MgZGF0b3MgPC9zcGFuPgpgYGB7cn0Kc3VtbWFyeShkZjEpCnN0cihkZjEpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDQuIEVzY2FsYXIgbG9zIGRhdG9zIDwvc3Bhbj4KYGBge3J9CiMgU29sbyBzaSBsb3MgZGF0b3Mgbm8gZXN0YW4gZW4gbGEgbWlzbWEgZXNjYWxhLgojIGRhdG9zX2VzY2FsYWRvcyA8LSBzY2FsZShkYXRvc19vcmlnaW5hbGVzKSAoRXN0ZSBlcyB1biBlamVtcGxvIG9idmlhbWVudGUpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDUuIERldGVybWluYXIgbnVtZXJvIGRlIGdydXBvcyA8L3NwYW4+CmBgYHtyfQojIFNpZW1wcmUgZXMgdW4gdmFsb3IgaW5pY2lhbCAiY3VhbHF1aWVyYSIsIGx1ZWdvIHNlIG9wdGltaXphLgpwbG90KGRmMSR4LCBkZjEkeSkKZ3J1cG9zMSA8LSAzCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDYuIEdlbmVyYXIgbG9zIGdydXBvcyA8L3NwYW4+CmBgYHtyfQpjbHVzdGVyczEgPC0ga21lYW5zKGRmMSxncnVwb3MxKQpjbHVzdGVyczEKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gNy4gT3B0aW1pemFyIGVsIG51bWVybyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpvcHRpbWl6YWNpb24xIDwtIGNsdXNHYXAoZGYxLCBGVU49a21lYW5zLCBuc3RhcnQ9MSwgSy5tYXg9NykKIyBFbCBLLm1heCBub3JtYWxtZW50ZSBlcyAxMCwgZW4gZXN0ZSBlamVyY2ljaW8gYWwgc2VyIDggZGF0b3Mgc2UgZGVqbyBlbiA3LgpwbG90IChvcHRpbWl6YWNpb24xLCB4bGFiPSJOdW1lcm8gZGUgY2x1c3RlcnMgayIpCiMgU2Ugc2VsZWNjaW9uYSBjb21vIG9wdGltbyBlbCBwcmltZXIgcHVudG8gbWFzIGFsdG8uCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDguIEdyYWZpY2FyIGxvcyBncnVwb3MgPC9zcGFuPgpgYGB7cn0KZnZpel9jbHVzdGVyKGNsdXN0ZXJzMSxkYXRhPWRmMSkKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gOC4gQWdyZWdhciBjbHVzdGVycyBhIGxhIGJhc2UgZGUgZGF0b3MgPC9zcGFuPgpgYGB7cn0KZGYxX2NsdXN0ZXJzIDwtIGNiaW5kKGRmMSwgY2x1c3RlciA9IGNsdXN0ZXJzMSRjbHVzdGVyKQpoZWFkKGRmMV9jbHVzdGVycykKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjp0dXJxdW9pc2U7Ij4gQ29uY2x1c2lvbmVzIDwvc3Bhbj4KTGEgdMOpY25pY2EgZGUgKmNsdXN0ZXJpbmcqIHBlcm1pdGUgaWRlbnRpZmljYXIgcGF0cm9uZXMgbyBncnVwb3MgbmF0dXJhbGVzIGVuIGxvcyBkYXRvcyBzaW4gbmVjZXNpZGFkIGRlIGV0aXF1ZXRhcyBwcmV2aWFzLiAKCg==