因果推論:H3・メッシュ単位でのSynthetic Controlの実装 8
今回はH3・メッシュ単位でのSynthetic Controlの実装をやってみたいと思います。
今日やること
- 結果の表示
- やったことの概要と流れ
所要時間は20分ほどとなっています。
それでは、さっそく始めていきましょう!
Synthetic Controlの実装でやったことの概要と流れ
①ライブラリのインポート
②データの読み込み
③データの確認
・欠損がないか
・重複行がないか
・必要な行数が正しく揃っているか
④データの前処理
・型の指定
・処置系列の作成
・ドナー系列の作成
・処置系列とドナー系列をPre/Postへ分割
⑤重みの推定
・目的関数を定義
・初期値と制約条件の作成
・Preの系列で最適化を実行し重みを推定
⑥RMSPE、Pre/Post RMSPE ratio、GAPの抽出と可視化
・全期間のドナー系列に重みを掛け反実仮想を作成
・観測値と反実仮想を使ってRMSPE、Pre/Post RMSPE ratioを計算
・GAPの抽出
・施策後平均GAPの計算
・グラフで可視化
⑦結果の出力

