r语言商品购物篮分析

1047521767

商品购物篮分析
现代商品种类繁多，顾客往往会由于需要购买的商品众多而变得疲于选择，且顾客并不会因为商品选择丰富而选择购买更多的商品。
- q7 L& a6 m- _) q: ~' Z1 Z6 F3 ?. I
对于某些商品，顾客会选择同时购买，如面包与牛奶、薯片与可乐等，当面包与牛奶或者薯片与可乐分布在商场的两侧，且距离十分遥远时，顾客购买的欲望就会减少，在时间紧迫的情况下顾客甚至会放弃购买某些计划购买的商品。相反，把牛奶与面包摆放在相邻的位置，既给顾客提供便利，提升购物体验，又提高顾客购买的概率，达到了促销的目的。

8 n0 Y+ q) H) g  W& q0 T7 U$ }某商品零售企业共收集了9835个购物篮的数据，其中包含169个不同的商品类别，售出商品总数为43367件。其数据示例如表所示。
4 Z) `; I3 y; T# L1 ~1 o, H
. F, [( A  @4 D% }0 W4 D
针对原始数据中不同商品销量进行统计，结果如表所示。
! w3 E1 B( f, ^; p0 Y3 nØ全脂牛奶销售量最高，销量为2513件，占比5.795%。
+ w! T0 V8 g2 N' Z- j& [Ø其余热销商品还有其他蔬菜销量1903件，占比4.388%。
) ~9 q, |3 i6 a' S% |Ø面包卷销量1809件，占比4.171%。
Ø苏打销量1715件，占比3.955%。酸奶销量1372件，占比3.164%等。
& p% x/ E+ x. p. {" B8 F前20种商品销量占据的比例约为50%，基本符合“二八定律”。
) p7 N* i* W% e
* G/ ]& t" y2 A: J/ c$ o对每一类商品的热销程度进行分析，有利于商家制定商品在货架的摆放策略和位置。若是某类商品较为热销，它的摆放位置可以有如下选择。
4 J7 c1 I3 j! I+ wØ可以把此类商品摆放到商场的中心位置，方便顾客选购。
% N( J9 t' B" f. {2 L9 p! ^Ø或者把此类商品摆放到商场深处位置，使顾客在购买热销商品前经过非热销商品，增加在非热销商品处的停留时间，促进非热销产品的销量。
关联规则算法主要用于寻找数据中项集之间的关联关系，基于样本的统计规律，进行关联规则分析。根据所分析的关联关系，可从一个特征的信息来推断另一个特征的信息。当信息置信度达到某一阈值时，就可以认为规则成立。
8 p0 {( \; o* {( k1 d关于这个算法有一个非常有名的故事：“尿布和啤酒”。故事是这样的：美国的妇女们经常会嘱咐她们的丈夫下班后为孩子买尿布，而丈夫在买完尿布后又要顺 手买回自己爱喝的啤酒，因此啤酒和尿布在一起被购买的机会很多。这个举措使尿布和啤酒的销量双双增加，并一直为众商家所津津乐道。
4 O0 B5 Q$ a0 A# |* V% ?" O' e" C
Ø1. Apriori算法应用广泛，可用于消费市场价格分析，猜测顾客的消费习惯，比如较有名的“尿布和啤酒”的故事；
. C% u9 g$ a8 c: g) Z% YØ2.网络安全领域中的入侵检测技术；
Ø3.可用在用于高校管理中，根据挖掘规则可以有效地辅助学校管理部门有针对性的开展贫困助学工作；
* d$ u0 A9 P: |Ø4.也可用在移动通信领域中，指导运营商的业务运营和辅助业务提供商的决策制定。
+ Y3 s8 z" ]0 y$ lØ关联规则算法的主要应用是购物篮分析，是为了从大量的订单中发现商品潜在的关联。其中常用的一个算法叫Apriori先验算法。
& F: I5 {6 l  e5 O1 y3 d$ S, V
" @; `8 k" v0 xØ关联规则（association rule）：是形如 X → Y 的蕴含表达式，其中X和Y是不相交的项集，即：X∩Y=∅。关联规则的强度可以用它的支持度（support）和置信度（confidence）来度量。
Ø支持度：一个项集或者规则在所有事物中出现的频率，确定规则可以用于给定数据集的频繁程度。σ(X):表示项集X的支持度计数
2 i3 j  F: _# P% y1 [) z- \Ø项集X的支持度：s(X)=σ(X)/N；规则X → Y的支持度：s(X → Y) = σ(X∪Y) / N
Ø通俗解释：简单地说，X==>Y的支持度就是指物品集X和物品集Y同时出现的概率。
2 H" x! w! [9 _Ø概率描述：物品集X对物品集Y的支持度support（X==>Y）=P（X n Y）
2 `/ Z) P2 w+ hØ实例说明：某天共有1000 个顾客到商场购买物品，其中有150个顾客同时购买了圆珠笔和笔记本，那么上述的关联规则的支持度就是15％。
##商品购物篮分析 数据科学实验 20210428
& J' W/ t7 o5 _, x! M& S  }1 _
: n) P" i0 c3 o* T0 P#install.packages('arules')
' }' G8 z( N& C3 ^! L7 p4 {# 用来画很神奇的云词图！！
#install.packages('htmlwidgets')
! ]8 P; C8 O8 r: @% L4 U0 B$ y#install.packages('wordcloud2')
#install.packages('ggplot2')
6 s* y  J  }/ ?: T* p% ]#install.packages('rlang')
library(Matrix)
library(arules)
6 l9 _% G  J" k, r# h/ q& G! W6 Jlibrary(wordcloud2)
library(ggplot2)
library(rlang)
3 U) ?; F  e! l6 m3 e, w

8 }& k/ p. f' ]) o3 T4 G- ?GoodOrder <- read.csv("/Users/janine/testt/GoodsOrder.csv")
7 L% M) X% ?( C% M/ {) fGoodsTypes <- read.csv("/Users/janine/testt/GoodsTypes.csv")
" o. Y6 \- F  hhead(GoodOrder,10)

( W0 k6 s0 y) |* H& YhotGoods <- data.frame(table(GoodOrder[,2]))
names(hotGoods) <- c("Goods","Num")
+ C8 w. M7 m8 b9 A8 o' x1 n, B1 z% r2 hhotGoods$Percent <- hotGoods$Num / sum(hotGoods$Num)
/ u! `  J3 Z- }: j0 whotGoods <- hotGoods[order(hotGoods$Percent, decreasing = T),] #商品按销量降序排列
write.csv(hotGoods, "/Users/janine/testt/hotGoods.csv", row.names = F) #导出数据
set.seed(3)
1 u* I2 j  c& i' q0 g0 B' ]) T, d#hotGoods排完序
wordcloud2(hotGoods[1:60,1:2],size = 0.75) #制作销量前60的商品名称词云，size:颜色
+ ~/ X& ?4 n" l0 A8 p# G, ?9 ahead(hotGoods,10) #展示销量前10名商品名称、销量和占比
; d# k* d- d: msum(hotGoods[1:20,3]) #前20中商品销量累积占比

, t5 `: c( ?, ~7 ]1 d, [% x3 w
+ Y+ H2 K5 ?% c6 D4 S% ~! J
7 [9 u; b7 o, F, L( v8 p! U1 NGoods <- merge(GoodOrder, GoodsTypes, "Goods", all.x=T, all.y=T)
$ s! I# R$ K' ~; x4 @: Hhead(Goods, 10)
4 E( |: x! g/ Z* U1 c2 m6 ChotTypes <- data.frame(table(Goods$Types))
# ]3 T  b& k/ Z5 @names(hotTypes) <- c("Types", "Num") # 重新命名
# 求出每个大类的比例值 round是保留四位小数
: T$ A! f3 M9 E. m$ ahotTypes$Percent <- round(hotTypes$Num / sum(hotTypes$Num),4)
hotTypes <- hotTypes[order(hotTypes$Percent, decreasing = T),]#数据从大到小排序
/ D0 S6 S. N% C0 L) G1 @write.csv(hotTypes,"/Users/janine/testt/hotGoods.csv",row.names = F)
6 S7 D  q5 Z! A3 P& u2 Z* _hotTypes$Types <- factor(hotTypes$Types ,levels =hotTypes$Types ,ordered = T)
head(hotTypes,10) #展示销量前10大类商品名称、销量和占比
sum(hotTypes[1:3,3])  # 证明排名前三的商品销量和超过全部销量总和的50%
$ D* }& Y$ w, w0 F, twordcloud2(hotTypes[1:10,1:2],size = 0.65) #制作销量前30的商品大类名称词云, size:颜色
: H: D4 Q: f& A% m
" H* {1 D& o" i7 K- j6 }
# 利用which函数将 这一类 全部取出来
Drink <- Goods[which(Goods$Types=="非酒精饮料"),]
hotDrink <- data.frame(table(Drink$Goods))
( _  s# v$ M! H1 v9 L& rnames(hotDrink) <- c("Goods","Num")
hotDrink$Percent <- round(hotDrink$Num/sum(hotDrink$Num),4)
$ s7 t" d4 C7 A. GhotDrink <- hotDrink[order(hotDrink$Percent, decreasing = T), ]
head(hotDrink)
/ Y9 x# ]3 }! U: o- r0 vsum(hotDrink[1:3,3]) #计算前三种饮料的占比
# 第一大类热销商品的前三销量 接近70%


par(family='STKaiti') #输出中文
/ w, U: T3 z1 Y: o4 W& q9 dhotDrink1<- hotDrink[1:5,]
. X( c8 L# K6 P! q; X' VOthers.drk <- data.frame(Goods="其他", Num=sum(hotDrink$Num)-sum(hotDrink1$Num), Percent=1-sum(hotDrink1$Percent))
7 Q" c. `; T# p5 C0 HhotDrink1 <- rbind(hotDrink1, Others.drk)
hotDrink1 <- hotDrink1[order(hotDrink1$Num, decreasing = F),]
hotDrink1$Goods <- factor(hotDrink1$Goods ,levels =hotDrink1$Goods ,ordered = T)
, K0 B8 Q0 w. s+ q$ ^myLabel = as.vector(hotDrink1$Goods)   
, s$ c. {1 d7 o( I6 umyLabel = paste(myLabel, "(", round(hotDrink1$Percent * 100, 2), "%)        ", sep = "")  
3 s' Q" X# a6 Z1 Nlibrary(RColorBrewer)
p <- ggplot(hotDrink1, aes(x="", y=Percent, fill=Goods))+
! d" g: I! }3 ]  geom_bar(stat = "identity")+
5 h& \4 ^, S  F. k* i. x. b- R  coord_polar(theta = "y")+
% ?3 h3 u& I, E, J  labs(x="", y="", title = "")+
  theme(axis.ticks = element_blank())+
  theme(legend.title = element_blank(), legend.position = "top")+
  theme(axis.text.x = element_blank())
( f, D0 m8 V- C' Qp+scale_fill_brewer(breaks = hotDrink1$Goods, labels = myLabel)+guides(fill=guide_legend(reverse = T))

: d8 y" y! G) P
/ j  p+ p& \: U0 [0 v1 O8 D6 P# 建模之前要转换数据格式list 然后才能使用apriori
# t: P5 U0 _% u- n7 tdatalist <- list()
for(i in unique(GoodOrder$ID)){
  datalist[] <- GoodOrder[which(GoodOrder$ID == i), 2]
}
. r6 G' F" K7 C3 `#datalist  # 可以显示出所有的购物篮里面的 九千多条
0 j' s9 N/ u- O& @1 D
# 以下才正式开始关联分析apriori
& I# X; h5 r: Z( m# 导入到关联分析函数中
- o# x" [' [, F2 y" ]% W6 YTransRep <- as(datalist, "transactions")
5 [# _- f) ^3 cRulesRep <- apriori(TransRep, parameter = list(support=0.02, confidence=0.25))
+ I. f4 ?  \6 _! |inspect(sort(RulesRep, by="lift")[1:25]) #按提升度从高到低查看前25条规则


% x% y; _# m$ c8 F6 }guides(fill=guide_legend(reverse = T))
通过模型的规则得出在顾客购买商品的时候会同时购买全脂牛奶。因此，商场可以根据实际情况进行布置。
! ?5 O, q% a, T9 y2 P6 e$ z
: p* V1 u3 O4 ]! F8 Q% gØ将全脂牛奶放在顾客购买商品的必经之路，或者商场显眼位置，方便顾客拿取。
Ø其他蔬菜、根茎类蔬菜、酸奶油、猪肉、黄油、本地蛋类和多种水果同时购买的概率较高，可以考虑捆绑销售，或者适当调整商场布置，将这些商品的距离尽量拉近，提升购物体验。
1 o( d: L6 H" Q! j  |2 l
. Z( `" H( Z7 o% c结论
+ Z, M  n+ D. M( ?0 H1 ]8 f% u本案例主要结合商品零售购物篮的案例，重点介绍了关联规则算法中的Apriori算法在商品零售购物篮分析案例中的应用。过程中详细的分析了商品零售的现状与问题，同时给出某商场的商品零售数据，分析了商品的热销程度，最后通过Apriori算法构建相应模型，并根据模型结果制定销售策略。


/ M" z, o3 I% J% @$ r( A
' [4 v  H+ t; {8 ~
$ O" b9 ?( m2 ^, x

' @7 v, K6 E' C! s6 |