Teoría

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

Algunos usos típicos de esta tecnica son:

  • Segmentación de clientes
  • Detección de anormalidades.
  • Categorización de documentos.

Paso 1. Instalar paquetes y llamar librerías

# install.packages("cluster")  #Analisis 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 optimización de numeros cluster
library(factoextra)

Paso 2. Obtener los datos

# 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, 7)
)

Paso 3. Entender los datos

# Paso 3. Entender los datos
summary(df1)
##        x               y         
##  Min.   :1.000   Min.   : 2.900  
##  1st Qu.:2.000   1st Qu.: 4.000  
##  Median :4.500   Median : 5.000  
##  Mean   :4.375   Mean   : 5.737  
##  3rd Qu.:6.250   3rd Qu.: 7.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 7

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 números de grupos

# Siempre es un valor inicial "cualquiera". Luego 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 2, 3, 3
## 
## Cluster means:
##          x        y
## 1 1.500000 3.950000
## 2 3.666667 8.333333
## 3 7.000000 4.333333
## 
## Clustering vector:
## [1] 2 1 3 2 3 3 1 2
## 
## Within cluster sum of squares by cluster:
## [1] 2.705000 9.333333 2.666667
##  (between_SS / total_SS =  82.9 %)
## 
## Available components:
## 
## [1] "cluster"      "centers"      "totss"        "withinss"     "tot.withinss"
## [6] "betweenss"    "size"         "iter"         "ifault"

Paso 7. Optimizar números los grupos

set.seed(123)
optimizacion <- clusGap(df1, FUN=kmeans, nstar= 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 clusters k", main="optimización de Clusters")

# Se selecciona como óptimo de clusters

Paso 8. Graficar los grupos

fviz_cluster(clusters1, data=df1)

Paso 9. Agregar clusters a la base de dsatos

df1_clusters <- cbind(df1, clusters1$cluster)
head(df1_clusters)
##   x  y clusters1$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 identicar patrones o grupos naturales en los datso sin necesidad de etiquetas previas

LS0tCnRpdGxlOiAiQWdydXBhbWllbnRvIChDbHVzdGVycykiCmF1dGhvcjogIk1hcmlhIEpvc2UgRmxvcmVzIgpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVFJVRQogICAgdG9jX2Zsb2F0OiBUUlVFCiAgICBjb2RlX2Rvd25sb2FkOiBUUlVFCiAgICB0aGVtZTogeWV0aQotLS0KCgohW10oaHR0cHM6Ly9taXJvLm1lZGl1bS5jb20vMSpKb1NJN29GWlhxS0c4SE50aTBhc3ZnLmdpZikgCgojIDxzcGFuIHN0eWxlID0gImNvbG9yOmJsdWU7Ij4gVGVvcsOtYSA8L3NwYW4+CioqQWdydXBhbWllbnRvKiogbyAqY2x1c3RlcmluZyogZXMgdW5hIHTDqWNuaWNhIGRlIGFwcmVuZGl6YWplIGF1dG9tw6F0aWNvIG5vIHN1cGVydmlzYWRvIHF1ZSBhZ3J1cGEgZGF0b3MgZW4gZnVuY2nDs24gZGUgc2ltaWxpdHVkIAoKQWxndW5vcyB1c29zIHTDrXBpY29zIGRlIGVzdGEgdGVjbmljYSBzb246IAoKKiBTZWdtZW50YWNpw7NuIGRlIGNsaWVudGVzICAKKiBEZXRlY2Npw7NuIGRlIGFub3JtYWxpZGFkZXMuIAoqIENhdGVnb3JpemFjacOzbiBkZSBkb2N1bWVudG9zLiAKCiMgPHNwYW4gc3R5bGUgPSAiY29sb3I6Ymx1ZTsiPiBQYXNvIDEuIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcsOtYXMgPC9zcGFuPgoKYGBge3IgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0KIyBpbnN0YWxsLnBhY2thZ2VzKCJjbHVzdGVyIikgICNBbmFsaXNpcyBkZSBBZ3J1cGFtaWVudG8KbGlicmFyeShjbHVzdGVyKQojIGluc3RhbGwucGFja2FnZXMoImdncGxvdDIiKSAjR3JhZmljYXIKbGlicmFyeShnZ3Bsb3QyKQojIGluc3RhbGwucGFja2FnZXMoImRhdGEudGFibGUiKSAjTWFuZWpvIGRlIG11Y2hvcyBkYXRvcwpsaWJyYXJ5KGRhdGEudGFibGUpCiMgaW5zdGFsbC5wYWNrYWdlcygiZmFjdG9leHRyYSIpICNHcsOhZmljYSBvcHRpbWl6YWNpw7NuIGRlIG51bWVyb3MgY2x1c3RlcgpsaWJyYXJ5KGZhY3RvZXh0cmEpCgpgYGAKCiMgPHNwYW4gc3R5bGUgPSAiY29sb3I6Ymx1ZTsiPiBQYXNvIDIuIE9idGVuZXIgbG9zIGRhdG9zICA8L3NwYW4+CmBgYHtyfQojIFBhc28gMi4gT2J0ZW5lciBsb3MgZGF0b3MKZGYxIDwtIGRhdGEuZnJhbWUoCiAgeCA9IGMoMiwgMiwgOCwgNSwgNywgNiwgMSwgNCksCiAgeSA9IGMoMTAsIDUsIDQsIDgsIDUsIDQsIDIuOSwgNykKKQoKCmBgYAojIDxzcGFuIHN0eWxlID0gImNvbG9yOmJsdWU7Ij4gUGFzbyAzLiBFbnRlbmRlciBsb3MgZGF0b3MgIDwvc3Bhbj4KYGBge3J9CiMgUGFzbyAzLiBFbnRlbmRlciBsb3MgZGF0b3MKc3VtbWFyeShkZjEpCnN0cihkZjEpCmBgYAoKIyA8c3BhbiBzdHlsZSA9ICJjb2xvcjpibHVlOyI+IFBhc28gNC4gRXNjYWxhciBsb3MgZGF0b3MgIDwvc3Bhbj4KYGBge3J9CiMgU8OzbG8gc2kgbG9zIGRhdG9zIG5vIGVzdMOhbiBlbiBsYSBtaXNtYSBlc2NhbGEuCiMgZGF0b3NfZXNjYWxhZG9zIDwtIHNjYWxlKGRhdG9zX29yaWdpbmFsZXMpCmBgYAoKIyA8c3BhbiBzdHlsZSA9ICJjb2xvcjpibHVlOyI+IFBhc28gNS4gRGV0ZXJtaW5hciBuw7ptZXJvcyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0KIyBTaWVtcHJlIGVzIHVuIHZhbG9yIGluaWNpYWwgImN1YWxxdWllcmEiLiBMdWVnbyBvcHRpbWl6YS4KcGxvdChkZjEkeCxkZjEkeSkKZ3J1cG9zMSA8LSAzCgpgYGAKCiMgPHNwYW4gc3R5bGUgPSAiY29sb3I6Ymx1ZTsiPiBQYXNvIDYuIEdlbmVyYXIgbG9zIGdydXBvcyA8L3NwYW4+CmBgYHtyfQpjbHVzdGVyczEgPC0ga21lYW5zKGRmMSwgZ3J1cG9zMSkKY2x1c3RlcnMxCmBgYAoKIyA8c3BhbiBzdHlsZSA9ICJjb2xvcjpibHVlOyI+IFBhc28gNy4gT3B0aW1pemFyIG7Dum1lcm9zIGxvcyBncnVwb3MgPC9zcGFuPgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpvcHRpbWl6YWNpb24gPC0gY2x1c0dhcChkZjEsIEZVTj1rbWVhbnMsIG5zdGFyPSAxLCBLLm1heD03KQojIEVsIEsubWF4IG5vcm1hbG1lbnRlIGVzIDEwLiBlbiBlc3RlIGVqZXJjaWNpbyBhbCBzZXIgOCBkYXRvcyBzZSBkZWrDsyBlbiA3CnBsb3Qob3B0aW1pemFjaW9uLCB4bGFiPSAiTsO6bWVybyBjbHVzdGVycyBrIiwgbWFpbj0ib3B0aW1pemFjacOzbiBkZSBDbHVzdGVycyIpCiMgU2Ugc2VsZWNjaW9uYSBjb21vIMOzcHRpbW8gZGUgY2x1c3RlcnMKYGBgCgojIDxzcGFuIHN0eWxlID0gImNvbG9yOmJsdWU7Ij4gUGFzbyA4LiBHcmFmaWNhciBsb3MgZ3J1cG9zIDwvc3Bhbj4KYGBge3J9CmZ2aXpfY2x1c3RlcihjbHVzdGVyczEsIGRhdGE9ZGYxKQpgYGAKCiMgPHNwYW4gc3R5bGUgPSAiY29sb3I6Ymx1ZTsiPiBQYXNvIDkuIEFncmVnYXIgY2x1c3RlcnMgYSBsYSBiYXNlIGRlIGRzYXRvcyA8L3NwYW4+CmBgYHtyfQpkZjFfY2x1c3RlcnMgPC0gY2JpbmQoZGYxLCBjbHVzdGVyczEkY2x1c3RlcikKaGVhZChkZjFfY2x1c3RlcnMpCmBgYAojIDxzcGFuIHN0eWxlID0gImNvbG9yOmJsdWU7Ij4gQ29uY2x1c2lvbmVzIDwvc3Bhbj4KTGEgdMOpY25pY2EgZGUgKipjbHVzdGVyaW5nKiogcGVybWl0ZSBpZGVudGljYXIgcGF0cm9uZXMgbyBncnVwb3MgbmF0dXJhbGVzIGVuIGxvcyBkYXRzbyBzaW4gbmVjZXNpZGFkIGRlIGV0aXF1ZXRhcyBwcmV2aWFzIAo=