Rで学ぶマーケットバスケット分析
Christopher Bruffaerts
Statistician
店内の商品

バスケット1: {"Bread", "Cheese"}
バスケット2: {"Bread", "Wine" , "Cheese"}
複数のバスケット
100人がスーパーに来たとき、同時に現れるアイテムの関係を見つけられますか?
例: Bread と Cheese

結果: 「もし〜なら、〜」
複数バスケットから学ぶ

多様な応用例
複数バスケットのデータセットを作成しましょう!
my_baskets = data.frame(
"Basket" = c(1,1,1,1, 2,2,2, 3,3, 4,4,4, 5,5, 6,6, 7,7),
"Product" = c("Bread", "Cheese", "Cheese", "Cheese",
"Bread", "Butter", "Wine",
"Butter", "Butter",
"Butter", "Wine", "Wine",
"Butter", "Cheese",
"Cheese", "Wine",
"Wine", "Wine")
)
バスケットの中身を確認
head(my_baskets)
Basket Product
1 1 Bread
2 1 Cheese
3 1 Cheese
4 1 Cheese
5 2 Bread
6 2 Butter
質問
n_distinct(my_baskets$Product)
[1] 4
n_distinct(my_baskets$Basket)
[1] 7
df_basket =
my_baskets %>%
group_by(Basket) %>%
summarize(
n_total = n(),
n_items = n_distinct(Product))
Basket n_total n_items
<dbl> <int> <int>
1 1 4 2
2 2 3 3
平均バスケットサイズ
basket_size %>%
summarize(
avg_total_items = mean(n_total),
avg_dist_items = mean(n_items))
# A tibble: 1 x 2
avg_total_items avg_dist_items
<dbl> <dbl>
1 2.57 1.86
バスケットサイズの分布
# 重複なしアイテム数の分布
ggplot(df_basket, aes(n_items)) +
geom_bar()

どのアイテムを見ますか?
そのアイテムは全バスケットで何回出現しますか?
そのアイテムを含むバスケットはいくつですか?
例:

Rでチーズを抽出
# チーズを含むバスケット数
my_baskets %>%
filter(Product == "Cheese") %>%
summarize(
n_tot_items = n(),
n_basket_item = n_distinct(Basket))
n_tot_items n_basket_item
1 5 3
アソシエーションルールマイニング: アイテム集合間の頻出共起関係を見つける手法。

ルール抽出の例:
この後のアジェンダ:

Rで学ぶマーケットバスケット分析