class: title-slide-section, center, middle
name: logistics
# 差の差分法の実践


---
## 差の差分法の実践

- $2 \times 2$ の差の差分法
- 二方向固定効果モデル
- イベントスタディ定式化



---
## 利用するパッケージ

新たに `ggfixest` パッケージを利用する。

```{r}
library(tidyverse)
library(fixest)
library(ggfixest)
```


---
## データの読み込み

- `fixest`パッケージの`base_did` データを利用する。(疑似データ)
  - 詳細は[こちら](https://lrberge.github.io/fixest/reference/base_did.html)

```{r}
data("base_did")
head(base_did)
```

---
## 変数一覧

| name   | description            |
|--------|------------------------|
| y      | 被説明変数                  |
| x1     | 説明変数（分析では利用しない）        |
| id     | 個人の識別ID                |
| period | Time period (1-10)     |
| post   | 処置のAfterダミー。6期目以降1を取る。 |
| treat  | トリートメントグループのダミー。       |


---
## 分析１： 2*2 のDID

- 推定モデル： $y_{it} = \beta_0 + \beta_1 Post_t + \beta_2 Treat_i + \gamma Treat_i \times Post_t + \epsilon_{it}$

```{r}
result_base <- feols( y ~ post + treat + post*treat | 0 , 
       cluster = "id", data = base_did)
etable(result_base)

```

---
## 分析２：二方向固定効果モデル

- 推定モデル： $y_{it} =  \gamma Treat_i \times Post_t + \theta_i + \theta_t + \epsilon_{it}$

```{r}
result_TWFE <- feols( y ~ post*treat | id + period , 
       cluster = "id", data = base_did)
```


---
## 分析２の結果

```{r}
etable(result_TWFE)
```

---
## 分析３：イベントスタディ

- 推定モデル： $y_{it} =  \sum_{k \neq 5} \gamma_k D_{it}^k + \theta_i + \theta_t + \epsilon_{it}$
  - 処置は6期目に発生しているので、その一つ前(t=5)を基準化( $\gamma_k = 0$ )する。
  
```{r}
result_ES <- feols( y ~ i(period, treat, ref = 5) | id + period , 
       cluster = "id", data = base_did)

```

---
## 分析３の結果

```{r}
etable(result_ES)
```

---
## イベントスタディのプロット（その１）

- `fixest`パッケージの`iplot`関数

```{r}
fig1 <- iplot(result_ES) 
```


---
## イベントスタディのプロット（その２）

- `ggfixest` パッケージを使うと、ggplot ベースのイベントスタディ図を作れる。

```{r}
fig2 <- ggiplot(
  result_ES, main = 'Event Study Example') +
  labs(y = 'Effect (point estimate and 95% confidence interval', x = 'Period')+
  theme_bw()
```

---
```{r}
plot(fig2)
```

---
## おまけ：`fixest`でのStaggered DID

- Sun and Abraham (2021) の推定方法が`fixest`において実装されている。

- 詳細は[こちら](https://cran.r-project.org/web/packages/fixest/vignettes/fixest_walkthrough.html#423_Staggered_difference-in-differences_(Sun_and_Abraham,_2020)を参照。


