
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ípicos de esta técnica son:
Segmentación de clientes Detección de anormalidades
*Categorización de documentos
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)
## 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
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 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
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", main="Optimización de Clusters")

#Se selecciona como óptimo el primer punto 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 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écnica de clustering permite identificar patrones o
grupos naturales en los datos sin necesidad de etiquetas previas.
LS0tCnRpdGxlOiAiQWdydXBhbWllbnRvIChDbHVzdGVycykiCmF1dGhvcjogIk1pZ3VlbCBBbmdlbCBMb3BleiIKZGF0ZTogIjIwMjUtMDgtMTgiCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiBUUlVFCiAgICB0b2NfZmxvYXQ6IFRSVUUKICAgIGNvZGVfZG93bmxvYWQ6IFRSVUUKICAgIHRoZW1lOiB5ZXRpCi0tLQoKIVtdKGh0dHBzOi8vcmVzLmNsb3VkaW5hcnkuY29tL2RyeThyemJ5eC9pbWFnZS9mZXRjaC9zLS1maXJrbDd5LS0tL2ZfYXV0by9xX2F1dG8vY19zY2FsZSx3XzE1MzYvaHR0cHM6Ly93d3cua25pbWUuY29tL3NpdGVzL2RlZmF1bHQvZmlsZXMvcHVibGljLzIwMjQtMDUvMS13aGF0LWlzLWNsdXN0ZXItYW5hbHlzaXMucG5nKQoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBUZW9yw61hIDwvc3Bhbj4KKipBZ3J1cGFtaWVudG8qKiBvICpjbHVzdGVyaW5nKiBlcyB1bmEgdMOpY25pY2EgZGUgYXByZW5kaXphamUgYXV0b23DoXRpY28gbm8gc3VwZXJ2aXNhZG8gcXVlIGFncnVwYSBkYXRvcyBlbiBmdW5jacOzbiBkZSBzdSBzaW1pbGl0dWQuICAKCkFsZ3Vub3MgdXNvcyB0w61waWNvcyBkZSBlc3RhIHTDqWNuaWNhIHNvbjoKCipTZWdtZW50YWNpw7NuIGRlIGNsaWVudGVzCipEZXRlY2Npw7NuIGRlIGFub3JtYWxpZGFkZXMKKkNhdGVnb3JpemFjacOzbiBkZSBkb2N1bWVudG9zCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gMS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hcyA8L3NwYW4+CmBgYHtyfQojaW5zdGFsbC5wYWNrYWdlcygiY2x1c3RlciIpICNBbsOhbGlzaXMgZGUgQWdydXBhbWllbnRvCmxpYnJhcnkoY2x1c3RlcikKI2luc3RhbGwucGFja2FnZXMoImdncGxvdDIiKSAjR3JhZmljYXIKbGlicmFyeShnZ3Bsb3QyKQojaW5zdGFsbC5wYWNrYWdlcygiZGF0YS50YWJsZSIpICNNYW5lam8gZGUgbXVjaG9zIGRhdG9zCmxpYnJhcnkoZGF0YS50YWJsZSkKI2luc3RhbGwucGFja2FnZXMoImZhY3RvZXh0cmEiKSAjR3LDoWZpY2Egb3B0aW1pemFjacOzbiBkZSBuw7ptZXJvIGRlIGNsdXN0ZXJzCmxpYnJhcnkoZmFjdG9leHRyYSkKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gMi4gT2J0ZW5lciBsb3MgZGF0b3MgPC9zcGFuPgpgYGB7cn0KZGYxIDwtIGRhdGEuZnJhbWUoeD1jKDIsMiw4LDUsNyw2LDEsNCksIHk9YygxMCw1LDQsOCw1LDQsMiw5KSkKYGBgCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlOyI+IFBhc28gMy4gRW50ZW5kZXIgbG9zIGRhdG9zIDwvc3Bhbj4KYGBge3J9CnN1bW1hcnkoZGYxKQpzdHIoZGYxKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA0LiBFc2NhbGFyIGxvcyBkYXRvcyA8L3NwYW4+CmBgYHtyfQojU29sbyBzaSBsb3MgZGF0b3Mgbm8gZXN0w6FuIGVuIGxhIG1pc21hIGVzY2FsYQojZGF0b3NfZXNjYWxhZG9zIDwtIHNjYWxlKGRhdG9zX29yaWdpbmFsZXMpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZTsiPiBQYXNvIDUuIERldGVybWluYXIgbsO6bWVybyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0KI1NpZW1wcmUgZXMgdW4gdmFsb3IgaW5pY2lhbCAiY3VhbHF1aWVyYSIsIGx1ZWdvIHNlIG9wdGltaXphCnBsb3QoZGYxJHgsZGYxJHkpCmdydXBvczEgPC0gMwpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA2LiBHZW5lcmFyIGxvcyBncnVwb3MgPC9zcGFuPgpgYGB7cn0KY2x1c3RlcnMxIDwtIGttZWFucyhkZjEsIGdydXBvczEpCmNsdXN0ZXJzMQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA3LiBPcHRpbWl6YXIgZWwgbsO6bWVybyBkZSBncnVwb3MgPC9zcGFuPgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpvcHRpbWl6YWNpb24gPC0gY2x1c0dhcChkZjEsIEZVTj1rbWVhbnMsIG5zdGFydD0xLCBLLm1heD03KQojRWwgSy5tYXggbm9ybWFsbWVudGUgZXMgMTAsIGVuIGVzdGUgZWplcmNpY2lvIGFsIHNlciA4IGRhdG9zIHNlIGRlasOzIGVuIDcKcGxvdChvcHRpbWl6YWNpb24sIHhsYWI9Ik7Dum1lcm8gZGUgY2x1c3RlcnMgayIsIG1haW49Ik9wdGltaXphY2nDs24gZGUgQ2x1c3RlcnMiKQojU2Ugc2VsZWNjaW9uYSBjb21vIMOzcHRpbW8gZWwgcHJpbWVyIHB1bnRvIG3DoXMgYWx0bwpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA4LiBHcmFmaWNhciBsb3MgZ3J1cG9zICA8L3NwYW4+CmBgYHtyfQpmdml6X2NsdXN0ZXIoY2x1c3RlcnMxLCBkYXRhPWRmMSkKYGBgCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyA5LiBBZ3JlZ2FyIENsdXN0ZXJzIGEgbGEgQmFzZSBkZSBEYXRvcyA8L3NwYW4+CmBgYHtyfQpkZjFfY2x1c3RlcnMgPC0gY2JpbmQoZGYxLCBjbHVzdGVyID0gY2x1c3RlcnMxJGNsdXN0ZXIpCmhlYWQoZGYxX2NsdXN0ZXJzKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWU7Ij4gUGFzbyAxMC4gQ29uY2x1c2lvbmVzIDwvc3Bhbj4KTGEgdMOpY25pY2EgZGUgKmNsdXN0ZXJpbmcqIHBlcm1pdGUgaWRlbnRpZmljYXIgcGF0cm9uZXMgbyBncnVwb3MgbmF0dXJhbGVzIGVuIGxvcyBkYXRvcyBzaW4gbmVjZXNpZGFkIGRlIGV0aXF1ZXRhcyBwcmV2aWFzLg==