Paso 1. Instalar paquetes y llamar librerias
library(cluster)
library(ggplot2)
library(data.table)
library(factoextra)
## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa
Paso 2. Obtener los datos
df <- data.frame(x=c(2,2,8,5,7,6,1,4), y=c(10,5,4,8,5,4,2,9))
Paso 3. Cantidad de grupos
grupos <- 3
Paso 4. Generar los segmentos
segmentos <- kmeans(df,grupos)
segmentos
## K-means clustering with 3 clusters of sizes 3, 3, 2
##
## Cluster means:
## x y
## 1 7.000000 4.333333
## 2 3.666667 9.000000
## 3 1.500000 3.500000
##
## Clustering vector:
## [1] 2 3 1 2 1 1 3 2
##
## Within cluster sum of squares by cluster:
## [1] 2.666667 6.666667 5.000000
## (between_SS / total_SS = 85.8 %)
##
## Available components:
##
## [1] "cluster" "centers" "totss" "withinss" "tot.withinss"
## [6] "betweenss" "size" "iter" "ifault"
Paso 5. Asignar el grupo al que pertenece cada observación
asignacion <- cbind(df, cluster = segmentos$cluster)
asignacion
## x y cluster
## 1 2 10 2
## 2 2 5 3
## 3 8 4 1
## 4 5 8 2
## 5 7 5 1
## 6 6 4 1
## 7 1 2 3
## 8 4 9 2
Paso 6. Graficar los clusters
fviz_cluster(segmentos, data = df)

Paso 7. Optimizar la cantidad de grupos
La cantidad óptima de grupos corresponde al punto más alto de la
siguiente gráfica.
set.seed(123)
optimizacion <- clusGap(df, FUN = kmeans, nstart = 1, K.max = 7)
plot(optimizacion, xlab = "Número de clusters k")

Conclusión
La segmentación o clusters es un algoritmo útil para las empresas
que desean clasificar a sus clientes y dirigir campañas de marketing más
enfocadas y especializadas.
LS0tCnRpdGxlOiAiTW9kdWxvIDIgU2VzaW9uIDEiCmF1dGhvcjogIkjDqWN0b3IgR3VhZGFsdXBlIGRlIGxhIEdhcnphIFRyZXZpw7FvIC0gQTAxMTc3OTYwIgpkYXRlOiAiMjAyNC0wMi0xOSIKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICB0b2M6IHRydWUKICAgIHRvY19mbG9hdDogdHJ1ZQogICAgY29kZV9kb3dubG9hZDogdHJ1ZQotLS0KCiMjIyMgTGEgKipzZWdtZW50YWNpw7NuKiogbyAqKmNsdXN0ZXJzKiogZXMgdW4gY29uanVudG8gZGUgdMOpY25pY2FzIGN1eW8gcHJvcMOzc2l0byBlcyBmb3JtYXIgZ3J1cG9zIGEgcGFydGlyIGRlIHVuIGNvbmp1bnRvIGRlIGVsZW1lbnRvcy4KCiFbXSgvVXNlcnMvaGVjdG9yZGVsYWdhcnphdHJldmluby9MaWJyYXJ5L0Nsb3VkU3RvcmFnZS9Hb29nbGVEcml2ZS1hMDExNzc5NjBAdGVjLm14L01pIHVuaWRhZC9MSVQvU2V4dG8gc2VtZXN0cmUvSW50ZWxpZ2VuY2lhIEFydGlmaWNpYWwgY29uIEltcGFjdG8gRW1wcmVzYXJpYWwvTW9kdWxvIDIvcG9yaWRlbnRpZGFkLnBuZykKCiMjIFBhc28gMS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyaWFzCmBgYHtyfQpsaWJyYXJ5KGNsdXN0ZXIpCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShkYXRhLnRhYmxlKQpsaWJyYXJ5KGZhY3RvZXh0cmEpCmBgYAoKIyMgUGFzbyAyLiBPYnRlbmVyIGxvcyBkYXRvcwpgYGB7cn0KZGYgPC0gZGF0YS5mcmFtZSh4PWMoMiwyLDgsNSw3LDYsMSw0KSwgeT1jKDEwLDUsNCw4LDUsNCwyLDkpKQpgYGAKCiMjIFBhc28gMy4gQ2FudGlkYWQgZGUgZ3J1cG9zCmBgYHtyfQpncnVwb3MgPC0gMwpgYGAKCiMjIFBhc28gNC4gR2VuZXJhciBsb3Mgc2VnbWVudG9zCmBgYHtyfQpzZWdtZW50b3MgPC0ga21lYW5zKGRmLGdydXBvcykKc2VnbWVudG9zCmBgYAoKIyMgUGFzbyA1LiBBc2lnbmFyIGVsIGdydXBvIGFsIHF1ZSBwZXJ0ZW5lY2UgY2FkYSBvYnNlcnZhY2nDs24KYGBge3J9CmFzaWduYWNpb24gPC0gY2JpbmQoZGYsIGNsdXN0ZXIgPSBzZWdtZW50b3MkY2x1c3RlcikKYXNpZ25hY2lvbgpgYGAKCiMjIFBhc28gNi4gR3JhZmljYXIgbG9zIGNsdXN0ZXJzCmBgYHtyfQpmdml6X2NsdXN0ZXIoc2VnbWVudG9zLCBkYXRhID0gZGYpIApgYGAKCiMjIFBhc28gNy4gT3B0aW1pemFyIGxhIGNhbnRpZGFkIGRlIGdydXBvcwpMYSBjYW50aWRhZCDDs3B0aW1hIGRlIGdydXBvcyBjb3JyZXNwb25kZSBhbCBwdW50byBtw6FzIGFsdG8gZGUgbGEgc2lndWllbnRlIGdyw6FmaWNhLgpgYGB7cn0Kc2V0LnNlZWQoMTIzKQpvcHRpbWl6YWNpb24gPC0gY2x1c0dhcChkZiwgRlVOID0ga21lYW5zLCBuc3RhcnQgPSAxLCBLLm1heCA9IDcpCnBsb3Qob3B0aW1pemFjaW9uLCB4bGFiID0gIk7Dum1lcm8gZGUgY2x1c3RlcnMgayIpCmBgYAoKIyMgQ29uY2x1c2nDs24KIyMjIyBMYSBzZWdtZW50YWNpw7NuIG8gY2x1c3RlcnMgZXMgdW4gYWxnb3JpdG1vIMO6dGlsIHBhcmEgbGFzIGVtcHJlc2FzIHF1ZSBkZXNlYW4gY2xhc2lmaWNhciBhIHN1cyBjbGllbnRlcyB5IGRpcmlnaXIgY2FtcGHDsWFzIGRlIG1hcmtldGluZyBtw6FzIGVuZm9jYWRhcyB5IGVzcGVjaWFsaXphZGFzLgoKIyMjIyBNw6FzIGluZm9ybWFjacOzbjogCltSIGZvciBEYXRhIFNjaWVuY2UgKDJlZCldKGh0dHBzOi8vcjRkcy5oYWQuY28ubnovKQ==