Actividad 4.1

Teoría

1. Crear base de datos

df <-  data.frame(x=c(2,2,8,5,7,6,1,4), y=c(10,5,4,8,5,4,2,9))

2. Determinar el número de grupos

grupos <- 3

3. Realizar la clasificación

segmentos <- kmeans(df,grupos)
#segmentos

4. Revisar la asignación de grupos

asignacion <- cbind(df, cluster=segmentos$cluster)
#asignacion

5. Graficar resultados

library(ggplot2)
library(factoextra)

fviz_cluster(segmentos, data=df, 
             palette=c("red", "blue", "green"), 
             ellipse.type="euclid", 
             star.plot= T, 
             repel= T, 
             ggtheme = theme()
             )

6. Optimizar grupos

library(cluster)
library(data.table)

set.seed(123)

optimizacion <- clusGap(df, FUN=kmeans, nstart=1, K.max=7)
plot(optimizacion, 
     xlab= "Número de clusters k")

# El punto más alto de la gráfica será el número de grupos de clusters optimo 

Actividad 4.1

Cargando librerias y dataframe

library(readr)
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:data.table':
## 
##     between, first, last

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

ventas <- read.csv("ventas.csv")
ventas <- na.omit(ventas)

Segmentación a ticket promedio

# Creando el ticket promedio (Total) por cliente (CustomerID)
agg <- aggregate(Total ~ CustomerID + BillNo, data= ventas, FUN= sum)
agg <- aggregate(Total ~ CustomerID , data= agg, FUN= mean)
#agg

# Agrupando el número de visitas por cliente
visitas <- group_by(ventas, CustomerID)
visitas <- summarize(visitas, num.visitas=n())
#visitas

# Agrupando visitas y agg
merge <- merge(visitas, agg, by= "CustomerID")
row.names(merge) <- (merge$CustomerID)
merge$CustomerID <- NULL
#merge 

Asignación de número de grupos y clasificación

# Asignando numero de grupos de clusters
grupos2<- 4

# Clsaificación
segmentos2 <- kmeans(merge, grupos2)
#segmentos2

Revisar la asignación de grupos

asignacion2 <- cbind(merge, cluster=segmentos2$cluster)
#asignacion2

Graficar resultados

fviz_cluster(segmentos2, data=merge, 
             palette=c("red", "blue", "green", "orange"), 
             ellipse.type="euclid", 
             star.plot= T, 
             repel= T, 
             ggtheme = theme()
             )

Optimizar grupos

set.seed(123)

optimizacion2 <- clusGap(merge, FUN=kmeans, nstart=1, K.max=7)
plot(optimizacion2, 
     xlab= "Número de clusters k")

# El punto más alto de la gráfica será el número de grupos de clusters optimo 

Solución de profe 4.1

ventas2 <- ventas
summary(ventas2)

##     BillNo            Itemname            Quantity            Date          
##  Length:388023      Length:388023      Min.   :    1.00   Length:388023     
##  Class :character   Class :character   1st Qu.:    2.00   Class :character  
##  Mode  :character   Mode  :character   Median :    5.00   Mode  :character  
##                                        Mean   :   12.89                     
##                                        3rd Qu.:   12.00                     
##                                        Max.   :80995.00                     
##      Hour               Price            CustomerID      Country         
##  Length:388023      Min.   :   0.000   Min.   :12346   Length:388023     
##  Class :character   1st Qu.:   1.250   1st Qu.:13950   Class :character  
##  Mode  :character   Median :   1.950   Median :15265   Mode  :character  
##                     Mean   :   3.079   Mean   :15317                     
##                     3rd Qu.:   3.750   3rd Qu.:16837                     
##                     Max.   :8142.750   Max.   :18287                     
##      Total          
##  Min.   :     0.00  
##  1st Qu.:     4.35  
##  Median :    11.40  
##  Mean   :    22.07  
##  3rd Qu.:    19.50  
##  Max.   :168469.60

#count(ventas2, BillNo, Sort= TRUE)
#count(ventas2, Country, Sort= TRUE)
#count(ventas2, Itemname, Sort= TRUE)
#count(ventas2, Date, Sort= TRUE)
#count(ventas2, Hour, Sort= TRUE)

Observaciones: 1. Tenemos cantidades, precios y totales negativos 2. Fecha y hora no tienen formatos 3. Tenemos NA’s en CustomerID

Limpieza de datos

# Número de NA's en ventas
sum(is.na(ventas2))

## [1] 0

# Cuantos NA's por variable
sapply(ventas2, function(x) sum(is.na(x)))

##     BillNo   Itemname   Quantity       Date       Hour      Price CustomerID 
##          0          0          0          0          0          0          0 
##    Country      Total 
##          0          0

# Eliminar NA's de CustomerID
ventas2 <- na.omit(ventas2)

# Eliminar totales negativos
ventas2 <- ventas2[ventas2$Total>0,]

# Identificar Outliers
boxplot(ventas2$Total, horizontal=TRUE)

Observaciones: 4. Hay outliers en “Total”

Cantidad de visitas por cliente, ticket promedio y union de ambas tablas

# Cantidad de visitas
visitas3 <- group_by(ventas2, CustomerID) %>% summarize(Visitas = n_distinct(BillNo))

# Obtener total por ticket
ticket_promedio <- aggregate(Total ~ CustomerID + BillNo, data= ventas2, FUN= sum)
ticket_promedio <- aggregate(Total ~ CustomerID , data= ticket_promedio, FUN= mean)
#ticket_promedio

# Juntar ambas tablas
objetos <- merge(visitas3, ticket_promedio, by= "CustomerID")
row.names(objetos) <- (objetos$CustomerID)
objetos$CustomerID <- NULL
#objetos

# Eliminar datos fuera de lo normal
# Limites de los datos fuera de lo normal
# Limite inferior= Q1 - 1.5*IQR (interquartil)
# Limite superior= Q3 + 1.5*IQR (interquartil)

#Columna visitas
IQR_V <- IQR(objetos$Visitas)
#IQR_V
summary(objetos)

##     Visitas            Total         
##  Min.   :  1.000   Min.   :    3.45  
##  1st Qu.:  1.000   1st Qu.:  178.30  
##  Median :  2.000   Median :  292.00  
##  Mean   :  4.227   Mean   :  415.62  
##  3rd Qu.:  5.000   3rd Qu.:  426.63  
##  Max.   :209.000   Max.   :84236.25

LI_V <- 1-1.5*IQR_V
#LI_V

LS_V<- 5+1.5*IQR_V
#LS_V

objetos<- objetos[objetos$Visitas<=11,]

# Columna ticket promedio
colnames(objetos) <- c("Visitas", "TicketPromedio")
IQR_TP <- IQR(objetos$TicketPromedio)
#IQR_TP

LI_TP<- 178.30 -1.5*IQR_TP
#LI_TP

LS_TP <- 426.63 + 1.5*IQR_TP
#LS_TP

objetos <- objetos[objetos$TicketPromedio<= 791.69,]
summary(objetos)

##     Visitas      TicketPromedio  
##  Min.   : 1.00   Min.   :  3.45  
##  1st Qu.: 1.00   1st Qu.:168.62  
##  Median : 2.00   Median :267.00  
##  Mean   : 2.97   Mean   :293.32  
##  3rd Qu.: 4.00   3rd Qu.:383.98  
##  Max.   :11.00   Max.   :789.56

Asignación de grupos

# Asignando numero de grupos de clusters
grupos3<- 4
df <- objetos 

# Clsaificación
segmentos3 <- kmeans(df, grupos3)
#segmentos3

# Revisar la asignación de grupos
asignacion3 <- cbind(df, cluster=segmentos3$cluster)
#asignacion3

# Grafica
fviz_cluster(segmentos3, data=df, 
             palette=c("red", "blue", "green", "orange"), 
             ellipse.type="euclid", 
             star.plot= T, 
             repel= T, 
             ggtheme = theme()
             )

# Optimizar grupos
set.seed(123)

optimizacion <- clusGap(df, FUN=kmeans, nstart=1, K.max=7)

## Warning: Quick-TRANSfer stage steps exceeded maximum (= 188200)

plot(optimizacion, 
     xlab= "Número de clusters k")