Cohort Analysis and LifeCycle Grids mixed segmentation with R(转)

This is the third post about LifeCycle Grids. You can find the first post about the sense of LifeCycle Grids and A-Z process for creating and visualizing with R programming language here. Lastly, here is the second post about adding monetary metrics (customer lifetime value – CLV – and customer acquisition cost – CAC) to the LifeCycle Grids.

Even after we added CLV and CAC to the LifeCycle Grids and obtained a representative view, there is always space for improvements. The main problem is that we combined customers with different characteristics like source of attraction, actual lifetime, and other features that we can use for better analysis. Therefore, in order to make our analysis more detailed, we will study how to combine Cohort Analysis and LifeCycle Grids in this post.

The main principle of Cohort Analysis is to combine customers through some common characteristics (e.g. registration date, first purchase date, medium/source/campaign of attraction and so on). Cohort Analysis allows us to split customers into homogeneous groups. Therefore, we can obtain benefits from combining homogeneous cohortsthrough acquisition characteristics with the homogeneous groups (grids) through lifecycle phase.

In addition, we will assume that we have not only actual CLV (CLV to date) but also predicted CLV (potential value). This can be helpful in some cases, for example when different advertisement campaigns are targeted to various customer segments that have different behavior and potential values as a result.

We will study how to combine customers with both first purchase date cohortsand campaign cohorts and distribute them into LifeCycle Grids, which interesting perspectives we have for analyzing, and how to visualize results.

Ok, let’s start by creating data sample with the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

# loading libraries

library(dplyr)

library(reshape2)

library(ggplot2)

library(googleVis)

set.seed(10)

# creating orders data sample

data <- data.frame(orderId=sample(c(1:5000), 25000, replace=TRUE),

product=sample(c(‘NULL‘,‘a‘,‘b‘,‘c‘), 25000, replace=TRUE,

prob=c(0.15, 0.65, 0.3, 0.15)))

order <- data.frame(orderId=c(1:5000),

clientId=sample(c(1:1500), 5000, replace=TRUE))

date <- data.frame(orderId=c(1:5000),

orderdate=sample((1:500), 5000, replace=TRUE))

orders <- merge(data, order, by=‘orderId‘)

orders <- merge(orders, date, by=‘orderId‘)

orders <- orders[orders$product!=‘NULL‘, ]

orders$orderdate <- as.Date(orders$orderdate, origin="2012-01-01")

rm(data, date, order)

# creating data frames with CAC, Gross margin, Campaigns and Potential CLV

gr.margin <- data.frame(product=c(‘a‘, ‘b‘, ‘c‘), grossmarg=c(1, 2, 3))

campaign <- data.frame(clientId=c(1:1500),

campaign=paste(‘campaign‘, sample(c(1:7), 1500, replace=TRUE), sep=‘ ‘))

cac <- data.frame(campaign=unique(campaign$campaign), cac=sample(c(10:15), 7, replace=TRUE))

campaign <- merge(campaign, cac, by=‘campaign‘)

potential <- data.frame(clientId=c(1:1500),

clv.p=sample(c(0:50), 1500, replace=TRUE))

rm(cac)

# reporting date

today <- as.Date(‘2013-05-16‘, format=‘%Y-%m-%d‘)

As a result, we’ve obtained the following data frames:

  • campaign, which includes campaign name and customer acquisition cost for each customer,
  • margin, which includes gross margin for each product,
  • potential, which includes potential values / predicted CLV for each client,
  • orders, which includes orders from our customers with products and order dates.

We will start by calculating necessary indexes (CLV, frequency, recency, potential value, CAC and average time lapses between purchases), adding campaigns, and defining cohort features based on the first purchase date for each customer with the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

# calculating CLV, frequency, recency, average time lapses between purchases and defining cohorts

orders <- merge(orders, gr.margin, by=‘product‘)

customers <- orders %>%

# combining products and summarising gross margin

group_by(orderId, clientId, orderdate) %>%

summarise(grossmarg=sum(grossmarg)) %>%

ungroup() %>%

# calculating frequency, recency, average time lapses between purchases and defining cohorts

group_by(clientId) %>%

mutate(frequency=n(),

recency=as.numeric(today-max(orderdate)),

av.gap=round(as.numeric(max(orderdate)-min(orderdate))/frequency, 0),

cohort=format(min(orderdate), format=‘%Y-%m‘)) %>%

ungroup() %>%

# calculating CLV to date

group_by(clientId, cohort, frequency, recency, av.gap) %>%

summarise(clv=sum(grossmarg)) %>%

arrange(clientId)

# calculating potential CLV and CAC

customers <- merge(customers, campaign, by=‘clientId‘)

customers <- merge(customers, potential, by=‘clientId‘)

# leading the potential value to more or less real value

customers$clv.p <- round(customers$clv.p / sqrt(customers$recency) * customers$frequency, 2)

rm(potential, gr.margin, today)

Therefore, we’ve obtained the customers data frame that looks like:

clientId cohort frequency recency av.gap clv  campaign   cac  clv.p
1        2012-06    5       23      60   32  campaign 2   14  25.02
2        2012-01    2      426      36   20  campaign 4   10   4.65
3        2012-09    4       64      48   24  campaign 4   10  17.50
4        2012-03    2      286      66   25  campaign 2   14   0.24
5        2012-01    6       89      66   54  campaign 1   11  11.45
6        2012-04    5       85      64   27  campaign 3   12   3.25

Furthermore, we need to define segments based on frequency and recency values. We will do this with the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

# adding segments

customers <- customers %>%

 mutate(segm.freq=ifelse(between(frequency, 1, 1), ‘1‘,

 ifelse(between(frequency, 2, 2), ‘2‘,

 ifelse(between(frequency, 3, 3), ‘3‘,

 ifelse(between(frequency, 4, 4), ‘4‘,

 ifelse(between(frequency, 5, 5), ‘5‘, ‘>5‘)))))) %>%

 mutate(segm.rec=ifelse(between(recency, 0, 30), ‘0-30 days‘,

 ifelse(between(recency, 31, 60), ‘31-60 days‘,

 ifelse(between(recency, 61, 90), ‘61-90 days‘,

 ifelse(between(recency, 91, 120), ‘91-120 days‘,

 ifelse(between(recency, 121, 180), ‘121-180 days‘, ‘>180 days‘))))))

# defining order of boundaries

customers$segm.freq <- factor(customers$segm.freq, levels=c(‘>5‘, ‘5‘, ‘4‘, ‘3‘, ‘2‘, ‘1‘))

customers$segm.rec <- factor(customers$segm.rec, levels=c(‘>180 days‘, ‘121-180 days‘, ‘91-120 days‘, ‘61-90 days‘, ‘31-60 days‘, ‘0-30 days‘))

Ok, this is the time for combining Cohort Analysis and LifeCycle grids into the mixed segmentation model.

First purchase date cohorts

We will start with a fairly common approach of combining cohorts, specifically with the first purchase date cohorts where the first purchase date is used for combining customers into groups (cohorts). Let’s take a look at this mixed segmentation from three perspectives, which I believe can be interesting:

  • We calculated average CLV and CAC in the previous post for each grid. And it is obvious that the grid consists of very different customers who are in the same lifecycle phase at the moment. Splitting customers into cohorts can be helpful in order to find how average values were obtained and whether or not there are difference between cohorts. We can obtain more accurate view and make more accurate decisions while working with groups of customers. The first prospect is a more visible structure of total averages (CLV and CAC) and comparison cohorts.
  • When we deal with cohorts that are based on a period characteristic (first purchase date in our case), we are sure that the composition of the cohort is constant. For instance, we will create monthly cohorts in our case and each new customer who makes the first purchase next month, she will be included in the next month’s cohort. Therefore, once the present month ends, we definitely know how many customers there are in the cohort and who they are. So, what is the benefit of such an obvious thing? If you like the LifeCycle Grids approach, you will certainly use it regularly, every week/month/quarter for example. It can be helpful to see paths of customer flows over the time and again, it can be easier and more helpful to deal with groups of customers who are combined into cohorts than with all of them together. Therefore, we can work with exact cohort flows from grid to grid since its “birthday”, analyze paths and proportions of the best and worst customers, direct activities on customers who are in the exact checkpoint of their path and analyze the effect of these activities, and compare the progress of different cohorts.
  • If we compare cohorts that are in the same grid, we will find the other interesting thing. For instance, we will have 14 cohorts in grid [3 purchases : 31-60 days] with customers who have made the same number of purchases (3) but with a different starting point of their lives/lifetimes with us. In other words, customers from cohort 2012-01 have been with us for 12 months more than customers from cohort 2013-01 have. What does their being located in the same grid mean? It means they had a different purchasing pace. The 2013-01 cohort clients have made their purchases faster than the 2012-01 cohort. Therefore, we can calculate purchasing pace for each cohort that allows us to answer the question “when is the best time for making an offer or contact the customer?”Note: you can find what the suitable offer is for the customer with the sequential analysis of shopping carts.

Let’s work with these prospects. We will start by combining LifeCycle Grids and first purchase date cohorts using the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

lcg.coh <- customers %>%

 group_by(cohort, segm.rec, segm.freq) %>%

 # calculating cumulative values

 summarise(quantity=n(),

 cac=sum(cac),

 clv=sum(clv),

 clv.p=sum(clv.p),

 av.gap=sum(av.gap)) %>%

 ungroup() %>%

 # calculating average values

 mutate(av.cac=round(cac/quantity, 2),

 av.clv=round(clv/quantity, 2),

 av.clv.p=round(clv.p/quantity, 2),

 av.clv.tot=av.clv+av.clv.p,

 av.gap=round(av.gap/quantity, 2),

 diff=av.clv-av.cac)

1. Structure of averages and comparison cohorts

We will start with two trivial charts:

  • differences between average CLV to date and average CAC for analyzing net customer value to date by cohort,
  • total average CLV (CLV to date + potential value) and average CAC for analyzing potential value by cohort. Note: actual CLV is dark, potential CLV is light, and CAC is the red line.

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

ggplot(lcg.coh, aes(x=cohort, fill=cohort)) +

theme_bw() +

theme(panel.grid = element_blank())+

geom_bar(aes(y=diff), stat=‘identity‘, alpha=0.5) +

geom_text(aes(y=diff, label=round(diff,0)), size=4) +

facet_grid(segm.freq ~ segm.rec) +

theme(axis.text.x=element_text(angle=90, hjust=.5, vjust=.5, face="plain")) +

ggtitle("Cohorts in LifeCycle Grids - difference between av.CLV to date and av.CAC")

ggplot(lcg.coh, aes(x=cohort, fill=cohort)) +

theme_bw() +

theme(panel.grid = element_blank())+

geom_bar(aes(y=av.clv.tot), stat=‘identity‘, alpha=0.2) +

geom_text(aes(y=av.clv.tot+10, label=round(av.clv.tot,0), color=cohort), size=4) +

geom_bar(aes(y=av.clv), stat=‘identity‘, alpha=0.7) +

geom_errorbar(aes(y=av.cac, ymax=av.cac, ymin=av.cac), color=‘red‘, size=1.2) +

geom_text(aes(y=av.cac, label=round(av.cac,0)), size=4, color=‘darkred‘, vjust=-.5) +

facet_grid(segm.freq ~ segm.rec) +

theme(axis.text.x=element_text(angle=90, hjust=.5, vjust=.5, face="plain")) +

ggtitle("Cohorts in LifeCycle Grids - total av.CLV and av.CAC")

Let’s look at grid [>5 purchases : 91-120 days]. The 2012-01 cohort has the highest actual customer´s net value and the highest total CLV. It is the oldest and had more chances to be more valuable. Compared to the 2012-02 cohort, the lifetime difference is only one month but values are significantly better. Therefore, we can distribute our limited advertisement budget more accurately than  just by knowing the grid’s total average.

2. Analyzing customer flows

Let’s study how we can visualize customers’ flows from grid to grid with the Sankey diagram. Assume we want to see the progress of the 2012-09 cohort as of the dates: 2012-10-01, 2013-01-01 and 2013-04-01. We will use the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

# customers flows analysis (FPD cohorts)

# defining cohort and reporting dates

coh <- ‘2012-09‘

report.dates <- c(‘2012-10-01‘, ‘2013-01-01‘, ‘2013-04-01‘)

report.dates <- as.Date(report.dates, format=‘%Y-%m-%d‘)

# defining segments for each cohort‘s customer for reporting dates

df.sankey <- data.frame()

for (i in 1:length(report.dates)) {

orders.cache <- orders %>%

filter(orderdate < report.dates[i])

customers.cache <- orders.cache %>%

select(-product, -grossmarg) %>%

unique() %>%

group_by(clientId) %>%

mutate(frequency=n(),

recency=as.numeric(report.dates[i] - max(orderdate)),

cohort=format(min(orderdate), format=‘%Y-%m‘)) %>%

ungroup() %>%

select(clientId, frequency, recency, cohort) %>%

unique() %>%

filter(cohort==coh) %>%

mutate(segm.freq=ifelse(between(frequency, 1, 1), ‘1 purch‘,

ifelse(between(frequency, 2, 2), ‘2 purch‘,

ifelse(between(frequency, 3, 3), ‘3 purch‘,

ifelse(between(frequency, 4, 4), ‘4 purch‘,

ifelse(between(frequency, 5, 5), ‘5 purch‘, ‘>5 purch‘)))))) %>%

mutate(segm.rec=ifelse(between(recency, 0, 30), ‘0-30 days‘,

ifelse(between(recency, 31, 60), ‘31-60 days‘,

ifelse(between(recency, 61, 90), ‘61-90 days‘,

ifelse(between(recency, 91, 120), ‘91-120 days‘,

ifelse(between(recency, 121, 180), ‘121-180 days‘, ‘>180 days‘)))))) %>%

mutate(cohort.segm=paste(cohort, segm.rec, segm.freq, sep=‘ : ‘),

report.date=report.dates[i]) %>%

select(clientId, cohort.segm, report.date)

df.sankey <- rbind(df.sankey, customers.cache)

}

# processing data for Sankey diagram format

df.sankey <- dcast(df.sankey, clientId ~ report.date, value.var=‘cohort.segm‘, fun.aggregate = NULL)

write.csv(df.sankey, ‘customers_path.csv‘, row.names=FALSE)

df.sankey <- df.sankey %>% select(-clientId)

df.sankey.plot <- data.frame()

for (i in 2:ncol(df.sankey)) {

df.sankey.cache <- df.sankey %>%

group_by(df.sankey[ , i-1], df.sankey[ , i]) %>%

summarise(n=n())

colnames(df.sankey.cache)[1:2] <- c(‘from‘, ‘to‘)

df.sankey.cache$from <- paste(df.sankey.cache$from, ‘ (‘, report.dates[i-1], ‘)‘, sep=‘‘)

df.sankey.cache$to <- paste(df.sankey.cache$to, ‘ (‘, report.dates[i], ‘)‘, sep=‘‘)

df.sankey.plot <- rbind(df.sankey.plot, df.sankey.cache)

}

# plotting

plot(gvisSankey(df.sankey.plot, from=‘from‘, to=‘to‘, weight=‘n‘,

options=list(height=900, width=1800, sankey="{link:{color:{fill:‘lightblue‘}}}")))

Note: if you plot this chart on your computer, it is interactive and you can highlight any paths and checkpoints.

Therefore, we can easily identify dominant paths, find the proportion of the best or worst clients, direct activities on customers who are in the exact checkpoint of their path and analyze the effect of these activities, and compare the progress of different cohorts. Lastly, we saved the path for each client in the customers_path.csv file that you can use for future work.

3. Analyzing purchasing pace

We will start by plotting the average time lapses between purchases. Actually, we already calculated this index when we created lcg.coh data frame:


1

2

3

4

5

6

7

8

ggplot(lcg.coh, aes(x=cohort, fill=cohort)) +

theme_bw() +

theme(panel.grid = element_blank())+

geom_bar(aes(y=av.gap), stat=‘identity‘, alpha=0.6) +

geom_text(aes(y=av.gap, label=round(av.gap,0)), size=4) +

facet_grid(segm.freq ~ segm.rec) +

theme(axis.text.x=element_text(angle=90, hjust=.5, vjust=.5, face="plain")) +

ggtitle("Cohorts in LifeCycle Grids - average time lapses between purchases")

This is not surprising: the time lapses increase by increasing the age of cohort. Of course, we have to take into account the boundaries on the recency axis. The larger the range of the boundary, the higher the probability to find older cohorts with a higher pace than younger cohorts. However, our main goal is to work with values that we’ve calculated.

Here are some points I want you to pay attention to:

  • We used averages. This means that we mixed customers with different purchase patterns. The averages can be used as an indicator for finding clients who have reached the mean line. However, in order to improve results, you can try to find some common features that would allow you to combine customers with similar behaviors and not only to split cohorts by behavior patterns but furthermore to understand why there are some differences. Therefore, you can make more personal offers to these groups of customers at the most suitable time. Note: this approach can be helpful for finding purchase behavior patterns as well as suitable offers.
  • Of course, we could split customers into segments based on average pace of purchases, hence avoiding cohorts. However, cohort analysis helps us to understand some roots of acquisition, e.g. time of first purchase/registration, medium/source/campaign of attraction, which can be very helpful for predicting future behavior as well as  potential CLV.
Campaign cohorts

The second example of cohort analysis is to combine customers by the advertisement campaign they were attracted by. This obviously can be helpful because we usually attract different customers with different campaigns. Therefore, we can expect that clients who were attracted by one campaign have some similarities in behavior and are sensitive to the exact same offers/communication channels. Furthermore, we would easily compare progress of campaigns in terms of monetary values (CLV and CAC).

We will use the same charts as the ones for first purchase date cohorts. Because we’ve added the campaign name to the data sample earlier, we can adapt our code by changing “cohort” value to “campaign” only:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

# campaign cohorts

lcg.camp <- customers %>%

group_by(campaign, segm.rec, segm.freq) %>%

# calculating cumulative values

summarise(quantity=n(),

cac=sum(cac),

clv=sum(clv),

clv.p=sum(clv.p),

av.gap=sum(av.gap)) %>%

ungroup() %>%

# calculating average values

mutate(av.cac=round(cac/quantity, 2),

av.clv=round(clv/quantity, 2),

av.clv.p=round(clv.p/quantity, 2),

av.clv.tot=av.clv+av.clv.p,

av.gap=round(av.gap/quantity, 2),

diff=av.clv-av.cac)

ggplot(lcg.camp, aes(x=campaign, fill=campaign)) +

theme_bw() +

theme(panel.grid = element_blank())+

geom_bar(aes(y=diff), stat=‘identity‘, alpha=0.5) +

geom_text(aes(y=diff, label=round(diff,0)), size=4) +

facet_grid(segm.freq ~ segm.rec) +

theme(axis.text.x=element_text(angle=90, hjust=.5, vjust=.5, face="plain")) +

ggtitle("Campaigns in LifeCycle Grids - difference between av.CLV to date and av.CAC")

ggplot(lcg.camp, aes(x=campaign, fill=campaign)) +

theme_bw() +

theme(panel.grid = element_blank())+

geom_bar(aes(y=av.clv.tot), stat=‘identity‘, alpha=0.2) +

geom_text(aes(y=av.clv.tot+10, label=round(av.clv.tot,0), color=campaign), size=4) +

geom_bar(aes(y=av.clv), stat=‘identity‘, alpha=0.7) +

geom_errorbar(aes(y=av.cac, ymax=av.cac, ymin=av.cac), color=‘red‘, size=1.2) +

geom_text(aes(y=av.cac, label=round(av.cac,0)), size=4, color=‘darkred‘, vjust=-.5) +

facet_grid(segm.freq ~ segm.rec) +

theme(axis.text.x=element_text(angle=90, hjust=.5, vjust=.5, face="plain")) +

ggtitle("Campaigns in LifeCycle Grids - total av.CLV and av.CAC")

ggplot(lcg.camp, aes(x=campaign, fill=campaign)) +

theme_bw() +

theme(panel.grid = element_blank())+

geom_bar(aes(y=av.gap), stat=‘identity‘, alpha=0.6) +

geom_text(aes(y=av.gap, label=round(av.gap,0)), size=4) +

facet_grid(segm.freq ~ segm.rec) +

theme(axis.text.x=element_text(angle=90, hjust=.5, vjust=.5, face="plain")) +

ggtitle("Campaigns in LifeCycle Grids - average time lapses between purchases")

And we’ve obtained these charts:

I believe everything is clear with these charts. I don’t think that Sankey diagram can be helpful enough for campaign cohorts. If we have some general campaigns that work for a long time period, we can obtain chaotic paths. Instead, I suggest studying a more accurate and visual approach that would be used for campaigns as well as first purchase date cohorts.

Dessert part – Lifecycle phase sequential analysis

Each customer has a path of migration from one grid to another that is based on purchasing behavior and affects CLV. They all have the same initial grid [1 purchase : 0-30 days], but since maximum 30 days (in our case) they had started a journey through grids. My idea is to analyze the path patterns of each cohort and identify cohorts that attracted customers with the path we prefer or not in order to make relevant offers. We will use thelifecycle phase sequential analysis for this. Note: you can find the example of shopping cart sequential analysis in my previous posts that started here so you can obtain other benefits of the method.

Everything we need for this is to reproduce paths through grids for each customer. We will do this with the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

# lifecycle phase sequential analysis

library(TraMineR)

min.date <- min(orders$orderdate)

max.date <- max(orders$orderdate)

l <- c(seq(0,as.numeric(max.date-min.date), 10), as.numeric(max.date-min.date))

df <- data.frame()

for (i in l) {

cur.date <- min.date + i

print(cur.date)

orders.cache <- orders %>%

filter(orderdate <= cur.date)

customers.cache <- orders.cache %>%

 select(-product, -grossmarg) %>%

 unique() %>%

 group_by(clientId) %>%

 mutate(frequency=n(),

 recency=as.numeric(cur.date - max(orderdate))) %>%

 ungroup() %>%

 select(clientId, frequency, recency) %>%

 unique() %>%

 mutate(segm=

 ifelse(between(frequency, 1, 2) & between(recency, 0, 60), ‘new customer‘,

 ifelse(between(frequency, 1, 2) & between(recency, 61, 180), ‘under risk new customer‘,

 ifelse(between(frequency, 1, 2) & recency > 180, ‘1x buyer‘,

 ifelse(between(frequency, 3, 4) & between(recency, 0, 60), ‘engaged customer‘,

 ifelse(between(frequency, 3, 4) & between(recency, 61, 180), ‘under risk engaged customer‘,

 ifelse(between(frequency, 3, 4) & recency > 180, ‘former engaged customer‘,

 ifelse(frequency > 4 & between(recency, 0, 60), ‘best customer‘,

 ifelse(frequency > 4 & between(recency, 61, 180), ‘under risk best customer‘,

 ifelse(frequency > 4 & recency > 180, ‘former best customer‘, NA)))))))))) %>%

 mutate(report.date=i) %>%

 select(clientId, segm, report.date)

df <- rbind(df, customers.cache)

}

df <- df %>%

mutate(grid=paste(segm.rec, segm.freq, sep=‘ : ‘)) %>%

select(clientId, grid, report.date)

We’ve checked the position of each customer in the LifeCycle Grids as of past dates. Here are two things that I want you to pay attention to: If you have quite a few customers, it would take a lot of time to reproduce grids for each past day. Therefore, we’ve used 10-day gaps in the loop that doesn´t seem detrimental in terms of accuracy as our minimal recency time lapse is 60 days.

Secondly, it would be quite tough to work with 36 grids. Therefore, we used 9 segments which you can adapt to your needs:

  • new customer (1-2 purchases and 0-60 days recency),
  • under risk new customer (1-2 purchases and 61-180 days recency),
  • 1x buyer (1-2 purchases and >180 days recency),
  • engaged customer (3-4 purchases and 0-60 days recency),
  • under risk engaged customer (3-4 purchases and 61-180 days recency),
  • former engaged customer (3-4 purchases and >180 days recency),
  • best customer (>4 purchases and 0-60 days recency),
  • under risk best customer (>4 purchases and 61-180 days recency),
  • former best customer (>4 purchases and >180 days recency).

We will do sequential analysis with the following code:


1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

# converting data to the sequince format

df <- dcast(df, clientId ~ report.date, value.var=‘grid‘, fun.aggregate = NULL)

df.seq <- seqdef(df, 2:ncol(df), left=‘DEL‘, right=‘DEL‘, xtstep=10)

# creating df with first purch.date and campaign cohort features

feat <- df %>% select(clientId)

feat <- merge(feat, campaign[, 1:2], by=‘clientId‘)

feat <- merge(feat, customers[, 1:2], by=‘clientId‘)

# plotting the 10 most frequent sequences based on campaign

seqfplot(df.seq, border=NA, group=feat$campaign)

# plotting the 10 most frequent sequences based on campaign

seqfplot(df.seq, border=NA, group=feat$campaign, cex.legend=0.9)

# plotting the 10 most frequent sequences based on first purch.date cohort

coh.list <- sort(unique(feat$cohort))

# defining cohorts for plotting

feat.coh.list <- feat[feat$cohort %in% coh.list[1:6] , ]

df.coh <- df %>% filter(clientId %in% c(feat.coh.list$clientId))

df.seq.coh <- seqdef(df.coh, 2:ncol(df.coh), left=‘DEL‘, right=‘DEL‘, xtstep=10)

seqfplot(df.seq.coh, border=NA, group=feat.coh.list$cohort, cex.legend=0.9)

What I really like about this approach is that we’ve easily led all customers to point zero of their life/lifetime with us. What I mean is that we replaced the first day of their lifetime (exact calendar date) with us with day 0 for all customers. This way, we switched from calendar dates to sequence dates. Therefore, all our sequences start from day 0. The white spaces mean that the next grid is unknown at the moment.

We can compare cohorts via share of customers with the different paths and current lifecycle phase (last color stripe). We can see that, for instance, 2012-01 cohort has brought us some part of customers who are the best now (yellow stripe), but 2012-03 cohort hasn’t.

In this way, we can identify different patterns in paths. For example, we can see the history of migrations for current best customers. Have they become the best ones avoiding under risk or former segments? Was there anything that could affect them and how we would use this for other clients?

Conclusions. We’ve studied how Cohort analysis can help us to combine customers into groups by common characteristics and obtain more clear view for differences between customers who are in the same grid of LifeCycle Grids. Also we’ve touched sequential analysis which helped us to find some patterns in the customers’ journey through grids. And we’ve found that customers who are on the same phase of lifecycle can have significantly different purchasing behavior. Therefore, it can be the topic for the future work: to create cohorts based on purchasing behavior/patterns.

Thank you for reading this! Feel free to share your thoughts about.

转自:http://analyzecore.com/2015/04/01/cohort-analysis-and-lifecycle-grids-mixed-segmentation-with-r/

时间: 2024-10-31 16:39:27

Cohort Analysis and LifeCycle Grids mixed segmentation with R(转)的相关文章

Customer segmentation – LifeCycle Grids with R(转)

I want to share a very powerful approach for customer segmentation in this post. It is based on customer’s lifecycle, specifically on frequency and recency of purchases. The idea of using these metrics comes from the RFM analysis. Recency and frequen

Customer segmentation – LifeCycle Grids, CLV and CAC with R(转)

We studied a very powerful approach for customer segmentation in the previous post, which is based on the customer’s lifecycle. We used two metrics: frequency and recency. It is also possible and very helpful to add monetary value to our segmentation

[转] 同期群分析Cohort Analysis

一个漂亮的平均数完全是用数据创造出来的虚幻景象,会给我们的决策造成误导,因此我们需要掌握一个行之有效的方法来剖析真实的用户行为和用户价值,这个方法就是同期群分析(Cohort Analysis).事实上,数据不会说谎,只是分析数据的人没有做到精准分析而导致对数据呈现的错误解读! 国内对同期群分析相关的研究相对较少,也许不是所有的运营都知道同期群分析,但它是每个产品运营必备的分析方法.在著名的<精益数据分析>一书里面,作为测试数据分析的灵魂也提到了同期群分析的相关内容. 同期群分析最早用于医药研

CS231n笔记 Lecture 11, Detection and Segmentation

Other Computer Vision Tasks Semantic Segmentation. Pixel level, don't care about instances. Classification + Localization. Single object. Object Detection. Multiple object. Instance Segmentation. Multiple object. Semantic Segmentation Simple idea: sl

cvpr 2016 论文学习 Video object segmentation

Abstract— Video object segmentation, a binary labelling problem is vital in various applications including object tracking, action recognition, video summarization, video editing, object based encoding and video retrieval(检索). This paper presents an

Data Analysis Project

Data Analysis Project [15 marks]Due date: 5pm, October 4, 2019Requirement: Pls finish the tasks according to the requirements. Pls copy all of your R code in the end of this word file, and then submit your word file via the Turnitin link in iLearn. Y

计算机视觉与模式识别代码合集第二版two

Topic Name Reference code Image Segmentation Segmentation by Minimum Code Length AY Yang, J. Wright, S. Shankar Sastry, Y. Ma , Unsupervised Segmentation of Natural Images via Lossy Data Compression, CVIU, 2007 code Image Segmentation Normalized Cut

CVPR 2017 Paper list

CVPR2017 paper list Machine Learning 1 Spotlight 1-1A Exclusivity-Consistency Regularized Multi-View Subspace Clustering Xiaojie Guo, Xiaobo Wang, Zhen Lei, Changqing Zhang, Stan Z. Li Borrowing Treasures From the Wealthy: Deep Transfer Learning Thro

精益数据分析--测试分析

几个概念---测试分析 市场细分 市场细分(market segmentation)是指营销者通过市场调研,依据消费者的需要和欲望.购买行为和购买习惯等方面的差异,把某一产品的市场整体划分为若干消费者群的市场分类过程.每一个消费者群就是一个细分市场,每一个细分市场都是具有类似需求倾向的消费者构成的群体. 同期群分析 同期群分析(Cohort Analysis,亦称群组分析)的主要目的是分析相似群体随时间的变化(比如用户的回访). A/B测试 A/B 测试,简单来说,就是为同一个目标制定两个方案(