
# La segmentación o clusters es un conjunto
# de técnicas cuyo propósito es formar
# grupos a partir de un conjunto de elementos.
# Paso 1. Instalar paquetes y llamar librerías
#install.packages("cluster")
library(cluster)
#install.packages("ggplot2")
library(ggplot2)
#install.packages("data.table")
library(data.table)
#install.packages("factoextra")
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,4)
segmentos
## K-means clustering with 4 clusters of sizes 3, 2, 1, 2
##
## Cluster means:
## x y
## 1 7.0 4.333333
## 2 4.5 8.500000
## 3 2.0 10.000000
## 4 1.5 3.500000
##
## Clustering vector:
## [1] 3 4 1 2 1 1 4 2
##
## Within cluster sum of squares by cluster:
## [1] 2.666667 1.000000 0.000000 5.000000
## (between_SS / total_SS = 91.4 %)
##
## 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 3
## 2 2 5 4
## 3 8 4 1
## 4 5 8 2
## 5 7 5 1
## 6 6 4 1
## 7 1 2 4
## 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 sus clientes y dirigir campañas de
marketing
más enfocadas y especializadas.
LS0tCnRpdGxlOiAiQ2x1c3RlciBQdW50b3MiCmF1dGhvcjogIkNvbnN0YW50aW5vIE1pbGxldCBYYWN1ciIKZGF0ZTogIjIwMjQtMDgtMTUiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVFJVRQogICAgdG9jX2Zsb2F0OiBUUlVFCiAgICBjb2RlX2Rvd25sb2FkOiBUUlVFCiAgICB0aGVtZTogInNwYWNlbGFiIgogICAgaGlnaGxpZ2h0OiAia2F0ZSIKLS0tCgohW10oL1VzZXJzL2NvbnN0YW50aW5vbWlsbGV0eGFjdXIvRGVza3RvcC9Db25jZW50cmFjaW9uL01vZHVsbyAyL0ZvdG9zIG0yL2NsdXN0ZXIucG5nKQoKYGBge3J9CiMgTGEgc2VnbWVudGFjacOzbiBvIGNsdXN0ZXJzIGVzIHVuIGNvbmp1bnRvCiMgZGUgdMOpY25pY2FzIGN1eW8gcHJvcMOzc2l0byBlcyBmb3JtYXIKIyBncnVwb3MgYSBwYXJ0aXIgZGUgdW4gY29uanVudG8gZGUgZWxlbWVudG9zLgoKIyBQYXNvIDEuIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcsOtYXMKI2luc3RhbGwucGFja2FnZXMoImNsdXN0ZXIiKQpsaWJyYXJ5KGNsdXN0ZXIpCiNpbnN0YWxsLnBhY2thZ2VzKCJnZ3Bsb3QyIikKbGlicmFyeShnZ3Bsb3QyKQojaW5zdGFsbC5wYWNrYWdlcygiZGF0YS50YWJsZSIpCmxpYnJhcnkoZGF0YS50YWJsZSkKI2luc3RhbGwucGFja2FnZXMoImZhY3RvZXh0cmEiKQpsaWJyYXJ5KGZhY3RvZXh0cmEpCmBgYAoKYGBge3J9CiMgUGFzbyAyLiBPYnRlbmVyIGxvcyBkYXRvcwpkZiA8LSBkYXRhLmZyYW1lKHg9YygyLDIsOCw1LDcsNiwxLDQpLCB5PWMoMTAsNSw0LDgsNSw0LDIsOSkpCmBgYAoKYGBge3J9CiMgUGFzbyAzLiBDYW50aWRhZCBkZSBncnVwb3MKZ3J1cG9zIDwtIDMKYGBgCgpgYGB7cn0KIyBQYXNvIDQuIEdlbmVyYXIgbG9zIHNlZ21lbnRvcwpzZWdtZW50b3MgPC0ga21lYW5zKGRmLDQpCnNlZ21lbnRvcwpgYGAKCmBgYHtyfQojIFBhc28gNS4gQXNpZ25hciBlbCBncnVwbyBhbCBxdWUgcGVydGVuZWNlIGNhZGEgb2JzZXJ2YWNpw7NuCmFzaWduYWNpb24gPC0gY2JpbmQoZGYsIGNsdXN0ZXIgPSBzZWdtZW50b3MkY2x1c3RlcikKYXNpZ25hY2lvbgpgYGAKCmBgYHtyfQojIFBhc28gNi4gR3JhZmljYXIgbG9zIGNsdXN0ZXJzCmZ2aXpfY2x1c3RlcihzZWdtZW50b3MsIGRhdGE9ZGYpCmBgYAoKYGBge3J9CiMgUGFzbyA3LiBPcHRpbWl6YXIgbGEgY2FudGlkYWQgZGUgZ3J1cG9zCiMgTGEgY2FudGlkYWQgw7NwdGltYSBkZSBncnVwb3MgY29ycmVzcG9uZGUgYWwgcHVudG8KIyBtw6FzIGFsdG8gZGUgbGEgc2lndWllbnRlIGdyw6FmaWNhLgpzZXQuc2VlZCgxMjMpCm9wdGltaXphY2lvbiA8LSBjbHVzR2FwKGRmLCBGVU49a21lYW5zLCBuc3RhcnQ9MSwgSy5tYXggPTcpCnBsb3Qob3B0aW1pemFjaW9uLCB4bGFiPSJOw7ptZXJvIGRlIGNsdXN0ZXJzIGsiKQpgYGAKCgojIENvbmNsdXNpw7NuIAojIExhIHNlZ21lbnRhY2nDs24gbyBjbHVzdGVycyBlcyB1biBhbGdvcml0bW8gw7p0aWwgcGFyYSBsYXMgZW1wcmVzYXMKIyBxdWUgZGVzZWFuIGNsYXNpZmljYXIgc3VzIGNsaWVudGVzIHkgZGlyaWdpciBjYW1wYcOxYXMgZGUgbWFya2V0aW5nCiMgbcOhcyBlbmZvY2FkYXMgeSBlc3BlY2lhbGl6YWRhcy4KCgoK4oCv4oGg