
# 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
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,4)
segmentos
## K-means clustering with 4 clusters of sizes 2, 3, 1, 2
##
## Cluster means:
## x y
## 1 1.5 3.500000
## 2 7.0 4.333333
## 3 2.0 10.000000
## 4 4.5 8.500000
##
## Clustering vector:
## [1] 3 1 2 4 2 2 1 4
##
## Within cluster sum of squares by cluster:
## [1] 5.000000 2.666667 0.000000 1.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 1
## 3 8 4 2
## 4 5 8 4
## 5 7 5 2
## 6 6 4 2
## 7 1 2 1
## 8 4 9 4
# 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="Numero 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.
LS0tDQp0aXRsZTogIkNsdXN0ZXJzIg0KYXV0aG9yOiAiTHVpcyBDYXJsb3MgQm9yYm9uIE1hcnRpbmV6Ig0KZGF0ZTogIjIwMjQtMDgtMTIiDQpvdXRwdXQ6DQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiB0cnVlDQogICAgdG9jX2ZvYXQ6IHRydWUNCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlDQotLS0NCiFbXShDOlxcVXNlcnNcXGxjYm9yXFxEb2N1bWVudHNcXHNpZ25pZmljYWRvLWNsdXN0ZXIuanBnKQ0KYGBge3J9DQojIExhIHNlZ21lbnRhY2nDs24gbyBjbHVzdGVycyBlcyB1biBjb25qdW50bw0KIyBkZSB0w6ljbmljYXMgY3V5byBwcm9ww7NzaXRvIGVzIGZvcm1hcg0KIyBncnVwb3MgYSBwYXJ0aXIgZGUgdW4gY29uanVudG8gZGUgZWxlbWVudG9zLg0KDQojIFBhc28gMS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hcw0KbGlicmFyeShjbHVzdGVyKQ0KbGlicmFyeShnZ3Bsb3QyKQ0KbGlicmFyeShkYXRhLnRhYmxlKQ0KbGlicmFyeShmYWN0b2V4dHJhKQ0KDQojIFBhc28gMi4gT2J0ZW5lciBsb3MgZGF0b3MNCmRmIDwtIGRhdGEuZnJhbWUoeD1jKDIsMiw4LDUsNyw2LDEsNCksIHk9YygxMCw1LDQsOCw1LDQsMiw5KSkNCg0KIyBQYXNvIDMuIENhbnRpZGFkIGRlIGdydXBvcw0KZ3J1cG9zIDwtIDMgDQoNCiMgUGFzbyA0LiBHZW5lcmFyIGxvcyBzZWdtZW50b3MNCnNlZ21lbnRvcyA8LSBrbWVhbnMoZGYsNCkNCnNlZ21lbnRvcw0KIyBQYXNvIDUuIEFzaWduYXIgZWwgZ3J1cG8gYWwgcXVlIHBlcnRlbmVjZSBjYWRhIG9ic2VydmFjacOzbg0KYXNpZ25hY2lvbiA8LSBjYmluZChkZiwgY2x1c3RlciA9IHNlZ21lbnRvcyRjbHVzdGVyKQ0KYXNpZ25hY2lvbg0KDQojIFBhc28gNi4gR3JhZmljYXIgbG9zIGNsdXN0ZXJzDQpmdml6X2NsdXN0ZXIoc2VnbWVudG9zLCBkYXRhPWRmKQ0KDQojIFBhc28gNy4gT3B0aW1pemFyIGxhIGNhbnRpZGFkIGRlIGdydXBvcw0KIyBMYSBjYW50aWRhZCDDs3B0aW1hIGRlIGdydXBvcyBjb3JyZXNwb25kZSBhbCBwdW50bw0KIyBtw6FzIGFsdG8gZGUgbGEgc2lndWllbnRlIGdyw6FmaWNhLg0Kc2V0LnNlZWQoMTIzKQ0Kb3B0aW1pemFjaW9uIDwtIGNsdXNHYXAoZGYsIEZVTj1rbWVhbnMsIG5zdGFydD0xLCBLLm1heCA9NykNCnBsb3Qob3B0aW1pemFjaW9uLCB4bGFiPSJOdW1lcm8gZGUgY2x1c3RlcnMgayIpDQoNCiMgQ29uY2x1c2nDs24gDQojIExhIHNlZ21lbnRhY2nDs24gbyBjbHVzdGVycyBlcyB1biBhbGdvcml0bW8gw7p0aWwgcGFyYSBsYXMgZW1wcmVzYXMgcXVlIGRlc2VhbiBjbGFzaWZpY2FyIHN1cyBjbGllbnRlcyB5IGRpcmlnaXIgY2FtcGHDsWFzIGRlIG1hcmtldGluZyBtw6FzIGVuZm9jYWRhcyB5IGVzcGVjaWFsaXphZGFzLg0KYGBgDQoNCg==