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 usos tícios de esta técnica son:

  • Segmentación de datos.
  • Detección de anormalidades.
  • Categorización de dpcumentos.

Paso 1. 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 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)

Paso 5. Determinar número de grupos

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

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)
optimizacion1 <-clusGap(df1, FUN=kmeans, nstart=1, K.max=7)
#El K,max normalmente es 10, en este ejercicio al se 8 datos se dejó en 7. 
plot(optimizacion1, xlab="Número de clusters k")

#se selecciona como óptimo el primer puntos 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,clusters = clusters1$cluster)
head(df1_clusters)
##   x  y clusters
## 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.

LS0tCnRpdGxlOiAiU2VnbWVudGFjaW9uIChDbHVzdGVycykiCmF1dGhvcjogIkpvc2UgUmlnb2JlcnRvIEEwMTczNTk5NCIKZGF0ZTogIjIwMjUtMDgtMTgiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVHJ1ZQogICAgdG9jX2Zsb2F0OiBUcnVlCiAgICBjb2RlX2Rvd25sb2FkOiBUcnVlCiAgICB0aGVtZTogeWV0aQotLS0KCiFbXShodHRwczovL2ZhY3VsdHlzdGFmZi5yaWNobW9uZC5lZHUvfnRtYXR0c29uL0lORk8zMDMvaW1hZ2VzL2NsdXN0ZXJfcDQuZ2lmKQoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBUZW9yw61hIDwvc3Bhbj4KKipBZ3J1cGFtaWVudG8qKiBvICpjbHVzdGVyaW5nKiBlcyB1bmEgdMOpY25pY2EgZGUgYXByZW5kaXphamUgYXV0b23DoXRpY28gbm8gc3VwZXJ2aXNhZG8gcXVlIGFncnVwYSBkYXRvcyBlbiBmdW5jacOzbiBkZSBzdSBzaW1pbGl0dWQuIAoKQWxndW5vcyB1c29zIHTDrWNpb3MgZGUgZXN0YSB0w6ljbmljYSBzb246CgoqIFNlZ21lbnRhY2nDs24gZGUgZGF0b3MuIAoqIERldGVjY2nDs24gZGUgYW5vcm1hbGlkYWRlcy4gICAKKiBDYXRlZ29yaXphY2nDs24gZGUgZHBjdW1lbnRvcy4gCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gMS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hcyA8L3NwYW4+CgpgYGB7ciBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQojaW5zdGFsbC5wYWNrYWdlcygiY2x1c3RlciIpICMiQW7DoWxpc2lzIGRlIGFncnVwYW1pZW50byIKbGlicmFyeShjbHVzdGVyKQojaW5zdGFsbC5wYWNrYWdlcygiZ2dwbG90MiIpICNHcmFmaWNhcgpsaWJyYXJ5KGdncGxvdDIpCiNpbnN0YWxsLnBhY2thZ2VzKCJkYXRhLnRhYmxlIikgI01hbmVqbyBkZSBtdWNob3MgZGF0b3MgCmxpYnJhcnkoZGF0YS50YWJsZSkKI2luc3RhbGwucGFja2FnZXMoImZhY3RvZXh0cmEiKSAjR3LDoWZpY2Egb3B0aW1pemFjacOzbiBkZSBuw7ptZXJvIGRlIGNsdXN0ZXJzCmxpYnJhcnkoZmFjdG9leHRyYSkKYGBgCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyAyLiBPYnRlbmVyIGxvcyBkYXRvcyA8L3NwYW4+CmBgYHtyfQpkZjE8LWRhdGEuZnJhbWUoeD1jKDIsMiw4LDUsNyw2LDEsNCksIHk9YygxMCw1LDQsOCw1LDQsMiw5KSkKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gMy4gRW50ZW5kZXIgbG9zIGRhdG9zIDwvc3Bhbj4KYGBge3J9CnN1bW1hcnkoZGYxKQpzdHIoZGYxKQpgYGAKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDQuIEVzY2FsYXIgbG9zIGRhdG9zIDwvc3Bhbj4KYGBge3J9CiMgU29sbyBzaSBsb3MgZGF0b3Mgbm8gZXN0YW4gZW4gbGEgbWlzbWEgZXNjYWxhLgojZGF0b3NfZXNjYWxhZG9zIDwtIHNjYWxlKGRhdG9zX29yaWdpbmFsZXMpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDUuIERldGVybWluYXIgbsO6bWVybyBkZSBncnVwb3MgPC9zcGFuPgoKYGBge3J9CiMgU2llbXByZSBlcyB1biB2YWxvciBpbmljaWFsICJDdWFscXVpZXJhIiwgbHVlZ28gc2Ugb3B0aW1pemEuCnBsb3QoZGYxJHgsZGYxJHkpCmdydXBvczEgPC0gMwpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA2LiBHZW5lcmFyIGxvcyBncnVwb3MgPC9zcGFuPgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpjbHVzdGVyczEgPC0ga21lYW5zKGRmMSxncnVwb3MxKQpjbHVzdGVyczEKYGBgCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA3LiBvcHRpbWl6YXIgZWwgbsO6bWVybyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpvcHRpbWl6YWNpb24xIDwtY2x1c0dhcChkZjEsIEZVTj1rbWVhbnMsIG5zdGFydD0xLCBLLm1heD03KQojRWwgSyxtYXggbm9ybWFsbWVudGUgZXMgMTAsIGVuIGVzdGUgZWplcmNpY2lvIGFsIHNlIDggZGF0b3Mgc2UgZGVqw7MgZW4gNy4gCnBsb3Qob3B0aW1pemFjaW9uMSwgeGxhYj0iTsO6bWVybyBkZSBjbHVzdGVycyBrIikKI3NlIHNlbGVjY2lvbmEgY29tbyDDs3B0aW1vIGVsIHByaW1lciBwdW50b3MgbcOhcyBhbHRvLgoKYGBgCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA4LiBHcmFmaWNhciBsb3MgZ3J1cG9zIDwvc3Bhbj4KYGBge3J9CmZ2aXpfY2x1c3RlcihjbHVzdGVyczEsZGF0YT1kZjEpCmBgYAojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gOS4gQWdyZWdhciBDbHVzdGVycyBhIGxhIEJhc2UgZGUgRGF0b3MgPC9zcGFuPgpgYGB7cn0KZGYxX2NsdXN0ZXJzIDwtY2JpbmQoZGYxLGNsdXN0ZXJzID0gY2x1c3RlcnMxJGNsdXN0ZXIpCmhlYWQoZGYxX2NsdXN0ZXJzKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gQ29uY2x1c2lvbmVzIDwvc3Bhbj4KTGEgdMOpY25pY2EgZGUgKkNsdXN0ZXJpbmcqIHBlcm1pdGUgaWRlbnRpZmljYXIgcGF0cm9uZXMgbyBncnVwb3MgbmF0dXJhbGVzIGVuIGxvcyBkYXRvcyBzaW4gbmVjZXNpZGFkIGRlIGV0aXF1ZXRhcyBwcmV2aWFzLiAKCgoKCgoKCgoKCgoK