前書き
デジタルテクノロジー統括部アナリストの安藤です。
最近マーケティングに関するプロジェクトに携わることになり、Marketing Mix Modeling(以下MMM)の門を叩くことになりました。
MMMでは、広告メディアがどれだけKPIに貢献したのかを回帰モデルによって推定し、それを元に予算の最適化を行います。
変数の値をそのまま使うのではなく、全体的なトレンドや季節性、繰越効果や飽和効果といった要因を扱っていることが特徴です。
MMMの代表的なOSSとしてGoogleのLightweightMMMやMetaのRobynなどが公表されていますが、2024年3月にGoogleから新たなMMMパッケージとしてMeridianが発表されました。Meridianは現在アーリーアクセス中であり、私もまだ使ったことはないのですが、一般公開とともにLightweightMMMの代替になる予定である、とのことです。
私のプロジェクトでもPython環境での使用の容易さからLightweightMMMを使用しており、Meridianへの移行を見据えて本記事ではLightweightMMMと比べた時のMeridianの新機能について掘り下げていきます。
すでに日本語でのMeridianの解説記事としてデータ分析な日々様の記事があるため、二番煎じにはなってしまいますが、可能な限りこちらで書かれていない部分について書いていきます。
全2回を予定しており、前半にあたる今回は時変ベースライン、Reach and Frequency データの取り扱い、ROIに事前分布への組み込みを取り扱います。
さらにMeridianには因果推論の考えが取り入れられており、こちらは次回触れる予定です(1-2ヶ月くらいで出せればいいな......)。
モデルの概要
全てのデータは事前にスケーリングされていることを前提としています。
各変数が指し示すものは以下の通りです。
- :地理単位のインデックス
- :時間単位のインデックス
- :制御変数(メディア変数以外の変数)のインデックス
- :メディア変数(Reach and Frequencyデータ**でないもの**)のインデックス
- :メディア変数(reach and frequencyデータ**であるもの**)のインデックス
- :時変ベースライン。詳細は後述
- :地理単位ごとのベースライン。地理単位を考慮しない場合は
- :制御変数に対する係数
- :制御変数の値
- :メディア変数(Reach and Frequencyデータ**でないもの**)に対する係数
- :メディア変数(Reach and Frequencyデータ**であるもの**)に対する係数
- :変数ベクトル を表す。は最大何期前までを考慮するかを表す整数パラメータ
- :とを組み合わせたもの。
`hill_before_adstock = False`(デフォルト)のとき
`hill_before_adstock = True`(デフォルト)のとき
- :投資効果の幾何的減衰を表す関数。
ただし
- :投資効果の飽和を表す関数。
ただしであり、
また、であり、の時凹型の、の時S字の形状をとる
- は地域・時間・チャンネルにおけるスケーリングされたReach数
- は地域・時間・チャンネルにおける平均Frequency数
- :正規分布に従う誤差項
一般的なMMMのモデル構造に馴染んだ人であれば、Reach and Frequencyデータに関する部分が追加されていること、トレンドや季節性の部分が時変ベースラインとしてまとめられていることなどが特徴的だと感じるかもしれません。
以降、これらについてさらに詳しく見ていきます。
Reach and Frequency データの取り扱い
通常のMMMでは、インプレッションに関する変数はインプレッションの総数のみを扱い、広告に接触したユーザーのユニーク数については扱いません。
Meridianでは、インプレッションを
- Reach(各期間における広告に接触したユニークユーザー数)
- Frequency(各期間におけるインプレッションの総数をリーチ数で割ったもの。平均接触回数)
に分けて扱うことが出来ます。
もう少し数式を詳しく見ていきましょう。
ある地域、チャンネルにおけるメディア広告の効果は次のように表されます。
最初に、ヒル関数をFrequency数に適用しています。ヒル関数ではS字型の飽和効果を表現可能です。広告への平均接触回数には特定の適切な頻度があり、想起の強化のために一定以上の接触が必要であるものの、多すぎても意味のない接触になったりかえって広告疲れを起こしたりするおそれがあるという想定です。
一方、Frequency数を固定したもとでは広告の効果とReach数は概ね線形の関係にあると想定し、
にアドストック関数を適用して広告の効果を推定します。実際にはReachに関しても飽和効果が見られmROIが減少する可能性があるため、観測されたデータにおけるReach数の範囲外に外挿しないことが推奨されています。
Reach and Frequencyデータを用いて予算最適化を行うケースでは、Frequencyを過去の平均に固定するオプションとFrequencyを最適化しつつ予算の最適化を行うオプションが存在します。
Optimization for media channels with reach and frequency dataのページにある通り、最適なFrequencyの推定は予算の割り当てに依存しないため、最適なFrequencyを推定した後その結果を組み込んで予算の配分を行うことが可能です。
時変ベースライン
LightweightMMMでは、時間によって変動する部分は
- トレンド項
- 季節項
の和で表していました(週次データの場合)。
一方、Meridianでは`knot_values`と呼ばれるパラメータに従って時変ベースラインを決定します。`knot`とは日本語で「結び目」を意味する単語です。
以下はSetting Knotsのページに従って解説していきます。
ドキュメントの記述と平仄を合わせるため、少しわかりにくいですが時変ベースラインパラメータのベクトルをとします。はです。
ここで
という形でを表現します。重み行列は行列、`knot_values`のベクトル (と呼ばれます)はです。
ベイズ推定の対象になるのはです。
ここは自分の解釈が入ってきてしまうのですが、`knot`は`knot_values`とは別に位置パラメータを持っていると考えるとわかりやすいです。
ドキュメントにはありませんが、議論をわかりやすくするため、`knot`の位置をと書くこととします。
各`knot`の位置はを等分するように決定されます。
、であるとすると、となります。
はととが定まった時点で各ととのL1距離によって自動的に決定されます。
ごとに以下の最大のとのL1距離、そして以上の最小のとのL1距離を取り、その距離に基づいてとの加重平均を取るようなイメージです。
例えばのような`knot`があり、に興味がある場合のベースラインは
と表現されます。
`knot`の数はバイアスとバリアンスのトレードオフを見ながら決定するべきである、と述べられています。`knot`の数が増えれば時間効果の推定値のバイアスが減少する一方でバリアンスは増加するという関係にあります。
これを踏まえた上で、
- 地理影響を考慮したモデルでは`knot`の数を`n_times`と一致させたところからスタートし、明らかなオーバーフィッティングや非現実的なメディア効果の推定が見られなくなるまで少しずつ`knot`の数を減らしていく
- 地理影響を考慮しないモデルでは`knot`の数を1からスタートし、明らかなオーバーフィッティングや非現実的なメディア効果が見られるようになるまで少しずつ`knot`の数を増やしていく
- `knot`の数を微増・微減させても(例:`knots=10`と`knots=11`)同様の結果が返ってくる可能性が高いため、`knot`の数は幅を持って試すことが有用である
また、`knot`の位置は必ずしも一定ではなく、ユーザーが自由に設定することが可能であり、例えば例えば祝日やスポーツイベントなどの影響を組み込むことが出来るようです。
ROIに関するドメイン知識や実験結果を事前分布に組み込む
すでに経験的な知識や特定のマーケティングに関する実験により、あるチャンネルの(大まかな)ROIがわかっている場合があります。
これまでのMMMツールではROIそのものに事前分布を設定することは出来ず、推定したROIとドメイン知識や実験によるROIが一致するように各種のパラメータや係数の事前分布を試行錯誤する必要がありました。
MeridianではROIを再パラメータ化すること、すなわちメディア効果の係数をROIの関数として表すことでROIに事前分布を設定し、設計者の信念をモデルに反映させることが容易になりました。
以下、ROI priors for calibrationのページに従って、をROIの関数として表すための式変形を追っていきます。
あるチャンネルによる売り上げの増分を次のように表現します。
まず、この増分をに関連する部分とそうでない部分に分離します。
増分のうちに関連しない部分を
とすると、次のように変形できます。
後の変形のため、に関する総和記号を外しておきましょう。
`media_effects_dist = Normal`の場合、を用いては次のように書くことができます。(`media_effects_dist = Normal`の場合の議論はドキュメントを参照してください)
従って、を主語にすると
となり、をROIに関する関数として書くことができました。
感想
ReachデータとFrequencyデータの活用について、観測範囲外のReach数の場合にモデルを外挿すべきではないというのは結構大きな制約ですが、それでも適切なFrequenctyを決定するための戦略決定に使えるのは嬉しいです。
また、今までトレンド+季節性で表現していた時変項目がMeridianのように特定の構造を持たないパラメータのベクトルで表現されるようになることでどう変わるのかというのは気になりますね。モデルの構造が複雑になればなるほど、解釈をするアナリストの数理的知識が試されるようになりそうです。
さらに、データ集めにも苦労することになりそうなのは言うまでもないですね......!GoogleとしてはGoogle Marketing Platformを使って欲しいという狙いがあるように見えます。こうしたプラットフォームを使わずにデータを蓄積していくのは大変ではありますが、「理論的にはこういうデータがあるとよい」というのがわかっていればその条件を満たすようにデータを集める仕組みを整えていけばよいわけですから、MMMを知ることでデータ収集やデータモデリング戦略にもプラスになると言えるでしょう。
また、これはLightweightMMMにもRobynにもないのですが、競合のデータを扱うフレームワークがあればさらに嬉しいなと感じます。例えば競合他社のTVCMのGRPや指名検索数などのデータが入手できているケースがあり、当然KPIに大きな影響をもたらしうる項目なのですが、残念ながらこれらをサポートするような機能はなさそうです(もちろんコントロール変数として組み入れることはできるのですが......)。
次回はMeridianがどのような形で因果推論の考え方を取り入れているのかについて解説していきます。
記事執筆者
安藤 有瑠聡
パーソルキャリア株式会社
デジタルテクノロジー統括部 デジタルビジネス部 アナリティクスグループ
2022年パーソルキャリア株式会社に入社。現在はアナリストとしてAzure MLを利用した社内アプリケーションのモニタリング、求人推薦アルゴリズムの開発、因果推論を活かした営業支援のための分析、社内データ基盤プロジェクトの横断調査など、幅広い分野の業務を行っている。