###************************************** 新老客户统计 ***************************************###
dachu <- read.csv("D:\\Dasktop\\bigdata_game\\天池\\大厨\\qijiandiankehu.csv", header = T, encoding = "utf-8", colClasses = c("character", "Date"))
str(dachu)
head(dachu,20)
temp <- table(dachu$买家昵称)
plot(table(sort(temp))/length(temp))
#library(data.table)
#month(dachu$下单日期[nrow(dachu)])
min(dachu$下单日期)
max(dachu$下单日期)
dachu$ym <- substr(dachu$下单日期, 1,7); head(dachu)
newcusts <- c()
oldcusts <- c()
ss <- sort(unique(dachu$ym))
#新客户满足一下两个条件：1）当月购买一次；2）之前无购买记录
#老客户满足一下两个条件之一：1）当月购买两次及以上；2）当月购买一次且之前有购买记录
for(i in 1:length(ss)){
#date1 = as.Date(paste(substr(kk, 1, 6), paste(as.integer(substr(kk, 7, 7))+1,"-01", sep = ""), sep = ""))
if(i == 1){
date2 = as.Date(paste(ss[i+1], "-01", sep = ""))
now = dachu$买家昵称[dachu$下单日期 < date2]
temp = table(now)
uniq = unique(now)
newcusts = c(newcusts, sum(temp == 1))
oldcusts = c(oldcusts, sum(temp > 1))
}else if(i < length(ss)){
date1 = as.Date(paste(ss[i], "-01", sep = ""))
date2 = as.Date(paste(ss[i+1], "-01", sep = ""))
now = dachu$买家昵称[(dachu$下单日期 < date2) & (dachu$下单日期 >= date1)]
temp = table(now)
#old_now = names(temp)[temp>1]
new_now = names(temp)[temp==1]
temp2 = table(c(uniq, new_now))
newcusts = c(newcusts, (length(new_now) - sum(temp2 > 1)))
#oldcusts = c(oldcusts, (length(old_now) + sum(temp2 > 1)))
oldcusts = c(oldcusts, (length(temp) - length(new_now) + sum(temp2 > 1)))
#uniq = unique(c(uniq, old_now, new_now))
uniq = unique(c(uniq, names(temp)))
}else{
date1 = as.Date(paste(ss[i], "-01", sep = ""))
now = dachu$买家昵称[dachu$下单日期 >= date1]
temp = table(now)
#old_now = names(temp)[temp>1]
new_now = names(temp)[temp==1]
temp2 = table(c(uniq, new_now))
newcusts = c(newcusts, (length(new_now) - sum(temp2 > 1)))
#oldcusts = c(oldcusts, (length(old_now) + sum(temp2 > 1)))
oldcusts = c(oldcusts, (length(temp) - length(new_now) + sum(temp2 > 1)))
#uniq = unique(c(uniq, old_now, new_now))
uniq = unique(c(uniq, names(temp)))
}
}
newcusts
oldcusts
(newcusts1 = cbind(date=ss, newcusts))
(oldcusts1 = cbind(date=ss, oldcusts))
write.csv(newcusts1, "C:\\Users\\hasee\\Desktop\\newcusts.csv",quote = F)
write.csv(oldcusts1, "C:\\Users\\hasee\\Desktop\\oldcusts.csv",quote = F)
#library(timeSeries)
win.graph()
opar <- par(no.readonly=TRUE)
par(lty=1, pch=1) #par("cex") 查看默认值
# plot.ts(ts(newcusts+oldcusts, start = c(2014, 3), frequency = 12),main="薏凡特月度新老客户购买数量变化趋势", col=1)
# lines(ts(newcusts, start = c(2014, 3), frequency = 12), col=2)
# lines(ts(oldcusts, start = c(2014, 3), frequency = 12), col=3)
time <- seq.Date(as.Date("2014/3/1"), by = "month", length = length(ss))
plot(time, newcusts+oldcusts, xlab="月份", ylab="客户数", main="薏凡特月度新老客户购买数量变化趋势",
type = "o", col=1)
# type画点/线, "p" for points, "l" for lines, "b" for both points and lines, "c" for empty points joined by lines,
# "o" for overplotted points and lines, "s" and "S" for stair steps and "h" for histogram-like vertical lines.
# Finally, "n" does not produce any points or lines.
# pch点型,
# cex点大小：
# lty线型：0=blank, 1=solid (default), 2=dashed, 3=dotted, 4=dotdash, 5=longdash, 6=twodash)
# lwd线宽
lines(time, newcusts, type = "o", col=2)
lines(time, oldcusts, type = "o", col=3)
legend("topright", c("总体客户", "新客户", "老客户"), col=1:3, lty=1, pch=1)
# “bottomright”, “bottom”, “bottomleft”, “left”, “topleft”, “top”, “topright”, “right”, “center”
par(opar)
#par(new=TRUE)
###************************************** 当月回购率 ***************************************###
# 月初统计购买一次的客户数，月末统计这部分人回购人数。
# 当月新进的客户且购买2次以上的不计入新客户
# 新客户可直接table=1的sum，但是当月回购的客户如何计算是难点。(可以unique内连接计数)
#数据导入
dachu <- read.csv("D:\\Dasktop\\bigdata_game\\天池\\大厨\\qijiandiankehu.csv", header = T, encoding = "utf-8", colClasses = c("character", "Date"))
str(dachu)
#定义保存新客户回购数据
new_customer <- data.frame()
min(dachu$下单时间)
max(dachu$下单时间)
ss=sort(unique(substr(dachu$下单时间,1,7)))
#从第二个月开始，首月新客数和回购数均为0
for(i in seq(length(ss))[-1]){
data1 = as.Date(paste(ss[i], "-01", sep = ""))
#月初之前客户购买记录
data2 = max(i-12,1)
data2 = as.Date(paste(ss[data2], "-01", sep = ""))
temp <- table(dachu$买家昵称[(dachu$下单时间 >= data2)&(dachu$下单时间 < data1)])
#月内客户购买记录
if(i < length(ss)){
data2 = as.Date(paste(ss[i+1], "-01", sep = ""))
temp2 <- table(dachu$买家昵称[(dachu$下单时间 >= data1)&(dachu$下单时间 < data2)])
}else{
temp2 <- table(dachu$买家昵称[dachu$下单时间 >= data1])
}
#月内回购记录
temp2 = merge(data.frame(k=names(temp)[temp==1]),
data.frame(k=names(temp2)),
by=c('k'))
#保存日期、月初新客数、月内回购数
new_customer = rbind(new_customer, data.frame(date=ss[i], counts=sum(temp==1), repurchase=nrow(temp2)))
}
#计算回购率
new_customer$rate <- new_customer[[3]] / new_customer[[2]]
#colnames(new_customer) = c('date','counts','repurchase','rate')
win.graph()
opar<-par(mfrow=c(2,2))
plot(new_customer$date,new_customer$counts)
plot(new_customer$date,new_customer$repurchase);plot(new_customer$date,new_customer$rate)
par(opar)
write.csv(new_customer,"C:\\Users\\hasee\\Desktop\\new_customer.csv")
###****************************************** 季度转化率 ****************************************###
#数据导入
dachu <- read.csv("C:\\Users\\hasee\\Desktop\\qijiandiankehu.csv", header = T, encoding = "utf-8", colClasses = c("character", "Date"))
str(dachu)
#定义保存新客户回购数据
new_customer <- data.frame()
min(dachu$下单时间)
max(dachu$下单时间)
ss=sort(unique(substr(dachu$下单时间,1,7)))
#从第二个月开始，首月新客数和回购数均为0
for(i in seq(length(ss)-2)[-1]){
data1 = as.Date(paste(ss[i], "-01", sep = ""))
#季度初之前客户购买记录
temp <- table(dachu$买家昵称[dachu$下单时间 < data1])
#季度内客户购买记录
if(i < length(ss)-2){
data2 = as.Date(paste(ss[i+3], "-01", sep = ""))
temp2 <- table(dachu$买家昵称[(dachu$下单时间 >= data1)&(dachu$下单时间 < data2)])
}else{
temp2 <- table(dachu$买家昵称[dachu$下单时间 >= data1])
}
#季度内回购记录
temp2 = merge(data.frame(k=names(temp)[temp==1]),
data.frame(k=names(temp2)),
by=c('k'))
#保存日期、季度初新客数、月内回购数
new_customer = rbind(new_customer, data.frame(date=ss[i], counts=sum(temp==1), repurchase=nrow(temp2)))
}
#计算回购率
new_customer$rate <- new_customer[[3]] / new_customer[[2]]
#colnames(new_customer) = c('date','counts','repurchase','rate')
win.graph()
opar<-par(mfrow=c(2,2))
plot(new_customer$date,new_customer$counts)
plot(new_customer$date,new_customer$repurchase);plot(new_customer$date,new_customer$rate)
par(opar)
write.csv(new_customer,"C:\\Users\\hasee\\Desktop\\new_customer.csv")
###************************************ 客户连带率：该段代码貌似有问题 ***********************************###
# 只针对所有一次客户
# 月连带率=本月发生连带的客户数/本月成交总客户数
# 产品连带率=购买该产品连带的客户数/购买该产品总体客户数
# 成交总客户=1次多件客户+一次一件客户
#数据导入
library(readxl)
# dachu <- read.csv("C:\\Users\\hasee\\Desktop\\liandailv.xlsx", header = T, encoding = "utf-8", colClasses = c("character", "Date", "character"))
# read_excel(path, sheet = 1, col_names = TRUE, col_types = NULL, na = "", skip = 0)
dachu <- read_excel("C:\\Users\\hasee\\Desktop\\liandailv.xlsx", sheet = 1, col_names = TRUE, col_types = c("text", "text", "text"), na = "", skip = 0)
dachu$下单日期 <- as.Date(dachu$下单日期)
str(dachu)
unique(dachu$商品ID)
#定义保存月度连带率
min(dachu$下单日期)
max(dachu$下单日期)
month_set=sort(unique(substr(dachu$下单日期,1,7)))
#月度连带率
month_associate_rate = data.frame()
date1 = min(dachu$下单日期)
for(i in seq(length(month_set))){
if(i < length(month_set)){
date2 = as.Date(paste(month_set[i+1], "-01", sep = ""))
temp <- table(dachu$买家昵称[(dachu$下单日期 >= date1)&(dachu$下单日期 < date2)])
date1 = date2
}else{
temp = table(dachu$买家昵称[dachu$下单日期 >= date1])
}
month_associate_rate = rbind(month_associate_rate, data.frame(month=month_set[i], count = length(temp), count2= sum(temp>1), rate=(sum(temp>1)/length(temp))))
}
month_associate_rate
#产品连带率
dachu$flag <- 0
head(dachu)
temp = table(dachu$买家昵称)
# library(dplyr)
# temp2 = left_join(dachu, data.frame(x = names(temp)[temp>1], flag.y = 1), by= c("买家昵称" = "x"),suffix = c("", ".y"))
temp2 = merge(dachu, data.frame(x = names(temp)[temp>1], flag.x = 1), by.x = "买家昵称", by.y = "x", all.x = TRUE)
temp2$flag[temp2$flag.x==1] = 1
temp2$flag.x = NULL
temp2
#定义保存产品连带率
prod_set=unique(dachu$商品ID)
product_associate_rate = data.frame()
#产品连带率
for(pi in prod_set){
temp <- temp2$flag[temp2$商品ID == pi]
product_associate_rate = rbind(product_associate_rate, data.frame(product=pi, count = length(temp), count2= sum(temp==1), rate=(sum(temp==1)/length(temp))))
}
product_associate_rate = product_associate_rate[order(product_associate_rate$count, decreasing = TRUE),]
product_associate_rate$product = as.character(product_associate_rate$product)
head(product_associate_rate)
#验证
dachu[dachu$买家昵称 %in% dachu[dachu$商品ID=="42303520877",]$买家昵称,]
#产品连带率前五月度变化
#temp2为产品连带率里计算的那个
prod_set = product_associate_rate$product[1:5]
product_associate_rate_top5 = data.frame()
date1 = min(temp2$下单日期)
for(i in seq(length(month_set))){
if(i < length(month_set)){
date2 = as.Date(paste(month_set[i+1], "-01", sep = ""))
temp <- temp2[(temp2$下单日期 >= date1)&(temp2$下单日期 < date2),]
date1 = date2
}else{
temp = temp2[temp2$下单日期 >= date1,]
}
temp3 = data.frame(month=month_set[i])
for(pi in prod_set){
temp4 = temp$flag[temp$商品ID==pi]
temp3 = cbind(temp3, length(temp4), sum(temp4==1), ifelse(length(temp4)==0,0,sum(temp4==1)/length(temp4)))
}
product_associate_rate_top5 = rbind(product_associate_rate_top5, temp3)
}
colnames(product_associate_rate_top5)[-1] <- paste('top',rep(1:5,each=3),c('count','count2','rate'),sep = '')
product_associate_rate_top5
#图形展示
win.graph()
opar<-par(mfrow=c(1,2))
plot(month_associate_rate$month, month_associate_rate$rate, type="l", col = "blue", main = "月度连带率", xlab = "月份", ylab="连带率")
plot(product_associate_rate$rate, main = "产品连带率", xlab = "产品", ylab="连带率")
par(opar)
write.csv(month_associate_rate,"C:\\Users\\hasee\\Desktop\\month_associate_rate.csv")
write.csv(product_associate_rate,"C:\\Users\\hasee\\Desktop\\product_associate_rate.csv") #, quote = TRUE
write.csv(product_associate_rate_top5,"C:\\Users\\hasee\\Desktop\\product_associate_rate_top5.csv") #, quote = TRUE
# dplyr包包含了各种关联查询的函数，如inner_join,left_join,full_join,rigth_join......
library(dplyr)
library("nycflights13")
# Drop unimportant variables so it's easier to understand the join results.
flights2 <-
flights %>%
select(year:day,tailnum, carrier)
flights2 %>%
left_join(airlines,by= "carrier")
#merge(data.frame(x=1:3,y=0,z=2),data.frame(x=2:3,y=1:2),by=c("x"),all.x = T)
###******************************************* 回购率与首次消费金额关系 ********************************************###
dachu <- read.csv("D:\\Dasktop\\bigdata_game\\天池\\大厨\\suoyoukehushuju.csv", header = T, encoding = "utf-8", colClasses = c("character", "Date", "numeric"))
str(dachu)
head(dachu,20)
library(dplyr)
temp=head(dachu,20)
temp = head(arrange(dachu, 买家昵称, desc(下单时间)), 100);temp
#flights[order(flights$year, flights$month, flights$day), ]
#flights[order(desc(flights$arr_delay)), ]
#filter(group_by(temp, 买家昵称))
temp <- dachu%>%
arrange(买家昵称, 下单时间) %>%
group_by(买家昵称)%>%
mutate(count = n())%>%
slice(1)%>%
filter()
win.graph()
opar<-par(mfrow=c(1,2))
#实付金额——购买次数分布图
plot(temp$实付金额, temp$count)
#实付金额——频数(人次)分布图
plot(table(temp$实付金额))
par(opar)
#通过第一个图，暂且分组0-1000等距每200,1000-2000,2000以上
temp$group <- 0
temp[temp$实付金额 < 1000, ]$group <- temp[temp$实付金额 < 1000, ]$实付金额 %/% 100
temp[(temp$实付金额 >= 1000) & (temp$实付金额 < 2000), ]$group <- 10
temp[temp$实付金额 >= 2000, ]$group <- 11
head(temp,20)
temp2 <- temp%>%
group_by(group)%>%
summarise(n1=sum(count>1), n2=n(), rate = n1/n2)
win.graph()
#各组回购率分布图
plot(temp2$group, temp2$rate)
# i <- c("gamma","a")
# switch(i,
# beta = "You typed beta",
# alpha = "You typed alpha",
# gamma = "You typed gamma",
# delta = "You typed delta"
# )
###******************************************* 客户联带对回购的影响 *******************************************###
t0 <- Sys.time()
dachu <- read.csv("D:\\Dasktop\\bigdata_game\\天池\\大厨\\AnalysisOrderDownLoad-订单信息-子订单（全量）-10027396-8025-107.csv",
header = T, encoding = "utf-8", colClasses = c(rep("character",4), rep("Date",3), rep("character",5), "integer","numeric","character",rep("numeric",2)))
str(dachu)
dachu <- dachu[,4:5]
head(dachu)
dachu$买家昵称 <- substr(dachu$买家昵称,3,nchar(dachu$买家昵称)-1)
head(dachu,20)
library(dplyr)
#首单购买件数回购率
temp <- dachu %>%
group_by(买家昵称, 下单时间) %>%
summarise(count=n()) %>%
arrange(买家昵称, 下单时间) %>%
group_by(买家昵称) %>%
mutate(count2=n()) %>%
slice(1) %>%
group_by(count) %>%
mutate(n1 = n(), n2 = sum(count2>1), rate = n2/n1) %>%
slice(1) %>%
select(count, n1, n2, rate)
temp
win.graph()
plot(temp$count, temp$rate, main="首单购买件数与回购率", xlab = "首单购买件数",
ylab = "回购客户占比", col="red")
#按月计算新客中回购客户占比
temp <- dachu %>%
group_by(买家昵称, 下单时间) %>%
summarise(count=n()) %>% #连带件数
mutate(year=as.integer(substr(下单时间,1,4)),
month=as.integer(substr(下单时间,6,7))) %>%
arrange(买家昵称, 下单时间) %>%
group_by(买家昵称) %>%
mutate(count2=n()) %>% #回购次数
slice(1) %>% #第一次出现(前面的按时间排序不可少)即为新客
group_by(year, month) %>%
mutate(n1 = n(), n2 = sum(count>1), rate = n2/n1) %>%
slice(1) %>%
select(下单时间, year, month, n1, n2, rate)
temp
win.graph()
time <- seq.Date(as.Date(paste(substr(min(temp$下单时间),1,7), "-01", sep = "")),
by = "month", length = nrow(temp))
plot(time, temp$rate, main = "各月新客中连带客户占比", xlab = "月份",
ylab = "首单购买多件客户占比", type = "l")
#按订单统计连带率(即购买多件订单与总订单之比)
temp <- dachu %>%
group_by(买家昵称, 下单时间) %>%
summarise(count=n())
sum(temp$count>1)/nrow(temp)
Sys.time()-t0
###############################################################################################################
#setwd("H:/数据分析/内部数据/薏凡特旗舰店数据/旗舰店客户数据分析/0803")
setwd("D:\\Dasktop\\bigdata_game\\天池\\大厨")
dat <- read.csv("kehushuju.csv",header=TRUE,encoding="utf-8",colClasses=c("character","Date","integer","numeric","integer"),stringsAsFactors = F)
dat <- arrange(dat, 买家昵称, 下单日期)
head(dat)
# new_dat<-unique(dat) #数据量多时，计算量很大，而且基本不会出现重复记录，所以可以省略
# head(new_dat)
library(dplyr)
##回购次数与回购概率
###
temp <- dat %>%
group_by(买家昵称)%>%
summarise(count=n())
head(temp)
rr1 <- c()
rr2 <- c()
rate <- c()
max_count <- max(temp$count)
for (i in 1:(max_count-1)){ ###可能会出错，rate分母=0
# rr1[i] <- summarise(filter(temp,count==i+1),n())
# rr2[i] <- summarise(filter(temp,count>=i),n())
# rate[i] <- summarise(filter(temp,count==i+1),n())/summarise(filter(temp,count>=i),n())
rr1 <- c(rr1, sum(temp$count == i+1)) #效率更高
rr2 <- c(rr2, sum(temp$count >= i))
rate <- c(rate, rr1[i]/rr2[i]) #避免重复计算
}
temp2<-filter(temp,count>=2)
head(temp2)
rrr<-cbind(rr1,rr2,rate)
rrr
# write.csv(rrr,"H:/数据分析/内部数据/薏凡特旗舰店数据/旗舰店客户数据分析/0803/rrr.csv")
#计算回购周期#####
#添加购买次数列
new_dat2 <- select(dat, 买家昵称,下单日期, 下单时点)
# new_dat2<-data.frame(new_dat2) #已经是数据框结构，而且即便转换格式此处也不对，应该为：new_dat2<-as.data.frame(new_dat2)
# new_dat2<-unique(new_dat2)
# head(new_dat2)
# temp2<-group_by(new_dat2,买家昵称)
# temp2<-summarise(temp2,count=n())
# temp2 <- new_dat2 %>%
# group_by(买家昵称) %>%
# summarise(count=n())
# head(temp2)
# count2<-unique(temp2$count)
#
# new_dat2$counts=0
# for(i in count2){
# rg<-temp[temp2$count==i,]$买家昵称;
# new_dat2[new_dat2$买家昵称 %in% rg,]$counts=i
#
# }
new_dat2 <- merge(new_dat2, temp, by=c('买家昵称'))
head(new_dat2)
# old_dat<-filter(new_dat2,counts>=2)
# old_dat<-arrange(old_dat,下单日期)
# old_dat <- new_dat2 %>% ##此处太慢，后面给出改进方法
# filter(count>=2) %>%
# arrange(下单日期)
# # old_dat<-unique(old_dat)
# head(old_dat)
# #max_count2<-max(old_dat$counts)
# #num<-c(1:max_count2)
# rebuy<-c()
# redays<-c()
# # t=1
# for(i in unique(old_dat$买家昵称) ){
# rg<-filter(old_dat,old_dat$买家昵称==i)
#
# for(j in 1:(rg$count[1]-1))
# {
# #t_diff <- rg$下单日期[j+1] - rg$下单日期[j]
# t_diff <- as.integer(rg$下单日期[j+1] - rg$下单日期[j])
# # rebuy[t]=j+1
# # redays[t]=t_diff
# # t=t+1
# rebuy = c(rebuy,j+1)
# redays = c(redays,t_diff)
# }
# }
#
# head(rebuy)
# head(redays)
# mydata<-data.frame(rebuy,redays)
# #write.csv(mydata,"H:/数据分析/内部数据/薏凡特旗舰店数据/旗舰店客户数据分析/0803/mydata.csv")
# head(mydata)
###各时点回购人数占比
#不考虑时间因素时
rate <- data.frame()
for(i in sort(unique(dat$下单时点))){
temp2 = new_dat2[new_dat2$下单时点 == i,]$count
rate = rbind(rate, c(i, sum(temp2>1)/length(temp2)))
}
colnames(rate) <- c("下单时点", "rate")
rate
#考虑时间因素时
###如果考虑时间因素，则需加以下代码
new_dat3 <- arrange(new_dat2, 买家昵称, 下单日期) #最好加排序，防止出错
head(new_dat3, 50)
# for(i in temp$买家昵称){ #由于循环较大故运行时间较长
# new_dat3[new_dat3$买家昵称 == i,]$count <- 1:(temp[temp$买家昵称 == i,]$count)
# }
# head(new_dat3, 50)
#改进后,此方法必须对数据先排序！！
# t0 <- Sys.time()
# i <- 1; nmax <- nrow(new_dat3)
# repeat{
# #m = i
# n = new_dat3[i,4]
# #ss = new_dat3[i,1]
# # repeat{
# # i <- i + 1
# # if((new_dat3[i,1] != ss) | (i > nmax)){
# # new_dat3[m:(i-1),4] <- 1:new_dat3[m,4]
# # break
# # }
# # }
# new_dat3[i:(i + n - 1),4] <- 1:n
# i = i+n
# if(i > nmax) break
# }
# Sys.time()-t0
#
# t0 <- Sys.time()
# i <- 1; nmax <- nrow(new_dat3)
# while(i <= nmax){
# #m = i
# n = new_dat3[i,4]
# #ss = new_dat3[i,1]
# # repeat{
# # i <- i + 1
# # if((new_dat3[i,1] != ss) | (i > nmax)){
# # new_dat3[m:(i-1),4] <- 1:new_dat3[m,4]
# # break
# # }
# # }
# new_dat3[i:(i + n - 1),4] <- 1:n
# i = i+n
# }
# Sys.time()-t0
t0 <- Sys.time()
for(i in sort(unique(temp$count))){ #必须加sort排序
df = (new_dat3$count == i)
new_dat3[df, 4] <- rep(1:i, sum(df)/i)
}
Sys.time()-t0
head(new_dat3, 50)
tail(new_dat3,50)
#计算
rate2 <- data.frame(下单时点=c(), rate=c())
for(i in sort(unique(dat$下单时点))){
temp2 = new_dat3[new_dat3$下单时点 == i,]$count
rate2 = rbind(rate2, c(i, sum(temp2>1)/length(temp2)))
}
colnames(rate2) <- c("下单时点", "rate")
rate2
#改进方法
new_dat3$t_diff <- as.integer(new_dat3$下单日期 - c(new_dat3$下单日期[1], new_dat3$下单日期[-nrow(new_dat3)]))
head(new_dat3)
new_dat3$t_diff[new_dat3$count==1] <- 0
mydata <- new_dat3 %>%
select(count, t_diff) %>%
filter(count > 1) %>%
rename(rebuy = count, redays = t_diff)
head(mydata)
plot(mydata)
#各次购买5天内回购情况
new_dat3$m5 <- (new_dat3$t_diff <5)
new_dat3$m5[new_dat3$count == 1] <- 0
setwd("H:/数据分析/内部数据/薏凡特旗舰店数据/旗舰店客户数据分析/0803/自我研究")
dat<-read.csv("kehushuju.csv",header=T,encoding="utf-8",colClasses=c("character","Date","integer","numeric","integer"))
head(dat)
library(dplyr)
dat1<-arrange(dat,下单日期)
head(dat1)
m=5 #定义回购周期，M=5表示客户在5天内回购
counts<-c(rep(0,length(dat1[,2])))
t0<-Sys.time()
for(i in 1:length(dat1[,2])){
t_run<-dat1[,2][i]+m
goal_dat1<-filter(dat1,下单日期<=t_run)
if(length(filter(goal_dat1,goal_dat1$买家昵称==dat1[,1][i])[,1])>=2){
counts[i]<-1
}
}
tt<-Sys.time()-t0
head(counts)
end_dat5<-cbind(dat1,counts)
write.csv(end_dat5,"H:/数据分析/内部数据/薏凡特旗舰店数据/旗舰店客户数据分析/0803/自我研究/end_dat5.csv")