{CausalImpact}のcovariateの与え方について

TJO @TJO_datasci

ブログ更新しました。タイトル変更がサムネイルに反映されていませんがご容赦を / 機械学習のビジネス上の価値を「効果測定」して「数値評価」する方法 - 六本木で働くデータサイエンティストのブログ - 六本木で働くデータサイエンティストのブログ tjo.hatenablog.com/entry/2019/04/…

だみ〜 @guiltydammy

tjo.hatenablog.com/entry/2019/04/… あれ？この場合比較対象はベイズ推定した時系列なので、厳密な意味でのDIDになってないのでは？DIDは介入群と非介入群両方の時系列が取れてる状態で行うもので、それができないのでDIDの非介入群をベイズ推定しちゃうよ、それっぽいよというのがcausalimpactの気持ちでは？

TJO @TJO_datasci

あえて書かなかったことを。。。。。。笑 twitter.com/guiltydammy/st…

だみ〜 @guiltydammy

@TJO_datasci ちなみに、causalちゃんに非介入群の時系列ぶちこむoptionがあった気がするんですが、使ったことないのでわからんでした

TJO @TJO_datasci

@guiltydammy {bsts}でモデル作ってからそいつを引数で与えるというのは出来ますね

だみ〜 @guiltydammy

@TJO_datasci おおなるほど。ちょうど自分でDIDするのつれえなあ〜テスト書かないといけないしな〜って思ってたのがあったので、今度やってみます！ありがとうございます！

TJO @TJO_datasci

@guiltydammy 原論文のp.259以降にその辺の話が書いてあるんですが、あまりきちんと説明されてないのでGitHubで実装見た方が早そうです ai.google/research/pubs/…

TJO @TJO_datasci

@guiltydammy 今見たら介入(y)＆非介入ペア(x)で時系列与えた時は、y〜xなるモデルをbstsで推定してこれを反実仮想予測の算出に使ってますね github.com/google/CausalI…

だみ〜 @guiltydammy

@TJO_datasci ありがとうございます！便利っすね〜

TJO @TJO_datasci

@guiltydammy でもこれって結構強い仮定ですよね。うちでも当然使ってるんですが、時々謎の結果になることがあるのは仮定に合わない時なのかなと

だみ〜 @guiltydammy

@TJO_datasci 非介入群与えないときの話ですか？そうなると所詮状態空間モデルなので、そこから逸脱してるとハマらないっすよねえ。。。

TJO @TJO_datasci

@guiltydammy 非介入群与えていても（というか与えるのがデフォルトだと思ってました）、互いにかなり強い回帰関係がないと反実仮想予測が当たらないわけで。。。

だみ〜 @guiltydammy

@TJO_datasci あれ？まじで勘違いしてたかもしれません。最後の論文Resultって、empiricalとsimulation（ベイズ推定）の結果比べてるから、デフォルトの使い方は非介入群与えないんだと思ってました。論文読めてないかも。。。

TJO @TJO_datasci

@guiltydammy web vignetteでは非介入群与えてやってますね。振る舞いについては多分論文よりもコード見た方が早いです。ConstructModel関数を見ると1列しか与えなければご指摘のようなやり方で、2列与えればy~xをbstsで推定して反実仮想予測をかけます github.com/google/CausalI…

TJO @TJO_datasci

@guiltydammy ちなみにweb vignetteはこちら。 google.github.io/CausalImpact/C… 今論文読み返したら、empirical dataでの実験については"Pairs of regions were then randomly assigned to treatment/control"と言っているので、介入群と非介入群とにsplitさせているようです。

だみ〜 @guiltydammy

@TJO_datasci あれ、すみません、これ読むとやっぱりoptional covariates = xで、非介入群じゃなくて、共変量に見えるのですが……物分りが悪くてすみません

TJO @TJO_datasci

@guiltydammy これ僕も疑問に思っていて、vignetteを読むとこう書いてあります。なので「ココロ」としてはxは多分非介入時系列なんでしょうが（その後のsection 2でもそういう趣旨でx1からyを作っているように見える）、表現としては一貫してcovariateなんですよね。 google.github.io/CausalImpact/C… pic.twitter.com/hpl4z7IEE5

拡大

だみ〜 @guiltydammy

@TJO_datasci すみません、covariatesの係数betaが一定だという仮定に読めるのですが、理解が悪くてココロとしてなぜxが非介入時系列になるのかが全然わからんです、betaが一定だっていう仮定をおいてるからこそbiasの上昇が捉えられるよね、betaが一定じゃなかったらbiasの上昇に吸い付くよねという話に読めます

TJO @TJO_datasci

@guiltydammy あとは、うちの社内で広告効果をテストする際に{CausalImpact}を用いる際は、covariate xにDiDで得た非介入時系列を充てています（グローバルのオペレーションがそうなってる）

だみ〜 @guiltydammy

@TJO_datasci なるほど〜ありがとうございます。じゃあ私が論文読み間違ってるっぽいですね。。。反省します。。。

TJO @TJO_datasci

@guiltydammy いえいえ、僕も"covariate"扱いになってる理由を完全には理解してない（実際やろうと思えばcovariateに3つとか4つとか変数を充てられる）ので、真面目に論文読んでみます。。。

だみ〜 @guiltydammy

@TJO_datasci 一応この辺読むと、covariatesはinventionの介入を絶対受けてないものを使ってねとあって、非介入群の同時系列が得られるなら絶対それ入れるよねというのはわかったので、そういう話なのかなと思いました。議論ありがとうございました！ pic.twitter.com/oQenahc7AB

拡大