Teoría

Agrupamiento o Clustering Aprendizaje automático no supervisado que agrupa datos en función de su suimilitud.

Algunos usos típicos de técnica son: Segmentación de cientes Dirección de anormalidades *Categorización de documentos

Instalar paquetes y llamar librerías

#install.packeges("cluster") #Análisis de Agrupamiento
library(cluster)
#install.packeges("ggplot2") #Graficar
library(ggplot2)
#install.packeges("data.table") #Manejo de muchos datos
library(data.table)
#install.packages("factoextra") #Grafica optimización de número de clusters
library(factoextra)
## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa

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 están 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

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

library(cluster)
library(factoextra)  # opcional para graficar bonito

set.seed(123)

optimizacion <- clusGap(
  df1,
  FUNcluster = kmeans,  # <- nombre correcto del argumento
  K.max = 7,            # <- K mayúscula
  B = 50,               # bootstraps (50 rápido; 500 más estable)
  nstart = 10
)

# Gráfica base de clusGap
plot(optimizacion, xlab = "Número de clusters k")

# (Opcional) con factoextra:
fviz_gap_stat(optimizacion)

# Paso 8: Graficar los grupos

fviz_cluster(cluster1, data=df1)

Paso 9: Agregar Clusters a la Base de Datos

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

Paso 10: Conclusiones

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

LS0tCnRpdGxlOiAiU2VnbWVudGFjacOzbiAoQ2x1c3RlcnMpIgphdXRob3I6ICJBbnRvbmlvIETDrWF6IEEwMDgzNzI1OSIKZGF0ZTogIjIwMjUtMDgtMTgiCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICAgIGNvZGVfZG93bmxvYWQ6IHRydWUKICAgIHRoZW1lOiB5ZXRpCi0tLQoKIVtdKGh0dHBzOi8vZGV2ZWxvcGVycy5nb29nbGUuY29tL21hY2hpbmUtbGVhcm5pbmcvY2x1c3RlcmluZy9vdmVydmlldykKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFRlb3LDrWEgPC9zcGFuPgoqKkFncnVwYW1pZW50byoqIG8gKkNsdXN0ZXJpbmcqCkFwcmVuZGl6YWplIGF1dG9tw6F0aWNvIG5vIHN1cGVydmlzYWRvIHF1ZSBhZ3J1cGEgZGF0b3MgZW4gZnVuY2nDs24gZGUgc3Ugc3VpbWlsaXR1ZC4gIAoKQWxndW5vcyB1c29zIHTDrXBpY29zIGRlIHTDqWNuaWNhIHNvbjoKKlNlZ21lbnRhY2nDs24gZGUgY2llbnRlcwoqRGlyZWNjacOzbiBkZSBhbm9ybWFsaWRhZGVzCipDYXRlZ29yaXphY2nDs24gZGUgZG9jdW1lbnRvcwoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hcyA8L3NwYW4+CgpgYGB7cn0KI2luc3RhbGwucGFja2VnZXMoImNsdXN0ZXIiKSAjQW7DoWxpc2lzIGRlIEFncnVwYW1pZW50bwpsaWJyYXJ5KGNsdXN0ZXIpCiNpbnN0YWxsLnBhY2tlZ2VzKCJnZ3Bsb3QyIikgI0dyYWZpY2FyCmxpYnJhcnkoZ2dwbG90MikKI2luc3RhbGwucGFja2VnZXMoImRhdGEudGFibGUiKSAjTWFuZWpvIGRlIG11Y2hvcyBkYXRvcwpsaWJyYXJ5KGRhdGEudGFibGUpCiNpbnN0YWxsLnBhY2thZ2VzKCJmYWN0b2V4dHJhIikgI0dyYWZpY2Egb3B0aW1pemFjacOzbiBkZSBuw7ptZXJvIGRlIGNsdXN0ZXJzCmxpYnJhcnkoZmFjdG9leHRyYSkKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gMi4gT2J0ZW5lciBsb3MgZGF0b3MgPC9zcGFuPgpgYGB7cn0KZGYxIDwtIGRhdGEuZnJhbWUoCiAgeCA9IGMoMiwyLDgsNSw3LDYsMSw0KSwKICB5ID0gYygxMCw1LDQsOCw1LDQsMiw5KQopCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDMuIEVudGVuZGVyIGxvcyBkYXRvcyA8L3NwYW4+CmBgYHtyfQpzdW1tYXJ5KGRmMSkKc3RyKGRmMSkKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gNC4gRXNjYWxhciBsb3MgZGF0b3MgPC9zcGFuPgpgYGB7cn0KIyBTb2xvIHNpIGxvcyBkYXRvcyBubyBlc3TDoW4gZW4gbGEgbWlzbWEgZXNjYWxhCiMgRGF0b3NfZXNjYWxhZG9zIDwtIHNjYWxlKGRhdG9zX29yaWdpbmFsZXMpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDUuIERldGVybWluYXIgbsO6bWVybyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0KIyBTaWVtcHJlIGVzIHVuIHZhbG9yIGluaWNpYWwgImN1YWxxdWllcmEiLCBsdWVnbyBzZSBvcHRpbWl6YS4KcGxvdChkZjEkeCxkZjEkeSkKZ3J1cG9zMSA8LSAzCgpgYGAKCgojIDxzcGFuIHN0eWxlPSJjbG9sb3I6Ymx1ZTsiPiBQYXNvIDY6IEdlbmVyYXIgbG9zIGdydXBvcyA8L3NwYW4+CmBgYHtyfQpjbHVzdGVyMSA8LSBrbWVhbnMoZGYxLGdydXBvczEpCmNsdXN0ZXIxCmBgYAoKIyA8c3BhbiBzdHlsZT0iY2xvbG9yOmJsdWU7Ij4gUGFzbyA3OiBPcHRpbWl6YXIgZWwgbsO6bWVybyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0KbGlicmFyeShjbHVzdGVyKQpsaWJyYXJ5KGZhY3RvZXh0cmEpICAjIG9wY2lvbmFsIHBhcmEgZ3JhZmljYXIgYm9uaXRvCgpzZXQuc2VlZCgxMjMpCgpvcHRpbWl6YWNpb24gPC0gY2x1c0dhcCgKICBkZjEsCiAgRlVOY2x1c3RlciA9IGttZWFucywgICMgPC0gbm9tYnJlIGNvcnJlY3RvIGRlbCBhcmd1bWVudG8KICBLLm1heCA9IDcsICAgICAgICAgICAgIyA8LSBLIG1hecO6c2N1bGEKICBCID0gNTAsICAgICAgICAgICAgICAgIyBib290c3RyYXBzICg1MCByw6FwaWRvOyA1MDAgbcOhcyBlc3RhYmxlKQogIG5zdGFydCA9IDEwCikKCiMgR3LDoWZpY2EgYmFzZSBkZSBjbHVzR2FwCnBsb3Qob3B0aW1pemFjaW9uLCB4bGFiID0gIk7Dum1lcm8gZGUgY2x1c3RlcnMgayIpCgojIChPcGNpb25hbCkgY29uIGZhY3RvZXh0cmE6CmZ2aXpfZ2FwX3N0YXQob3B0aW1pemFjaW9uKQpgYGAKIyA8c3BhbiBzdHlsZT0iY2xvbG9yOmJsdWU7Ij4gUGFzbyA4OiBHcmFmaWNhciBsb3MgZ3J1cG9zIDwvc3Bhbj4KYGBge3J9CmZ2aXpfY2x1c3RlcihjbHVzdGVyMSwgZGF0YT1kZjEpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY2xvbG9yOmJsdWU7Ij4gUGFzbyA5OiBBZ3JlZ2FyIENsdXN0ZXJzIGEgbGEgQmFzZSBkZSBEYXRvcyA8L3NwYW4+CmBgYHtyfQpkZjFfY2x1c3RlcnMgPC0gY2JpbmQoZGYxLCBjbHVzdGVyID0gY2x1c3RlcjEkY2x1c3RlcikKaGVhZChkZjFfY2x1c3RlcnMpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY2xvbG9yOmJsdWU7Ij4gUGFzbyAxMDogQ29uY2x1c2lvbmVzIDwvc3Bhbj4KTGEgdMOpY25pY28gZGUgKmNsdXN0ZXJpbmcqIHBlcm1pdGUgaWRlbnRpZmljYXIgcGF0cm9uZXMgbyBncnVwb3MgbmF0dXJhbGVzIGVuIGxvcyBkYXRvcyBzaW4gbmVjZXNpZGFkZSBkZSBldGlxdWV0YXMgcHJldmlhcy4K