Github: https://github.com/JaquelinMorillo/ReplicacionAbortion
En este contexto, los autores plantearon el siguiente modelo:
\[Y_st = \beta_1Repeals + \beta_2DT_t + \beta_3 Repeal_s DT_t + X_st\psi +\alpha_sDS_s+\varepsilon_{st}\] donde \(Y\) es el número logarítmico de nuevos casos de gonorrea en jóvenes de 15 a 19 años (por cada 100.000 habitantes); \(Repeal_s\) es igual a 1 si el estado legalizó el aborto antes de Roe; \(DT_t\) es una variable dummy de año; \(DS_s\) es una variable dummy de estado; \(t\) es una tendencia temporal; \(X\) es una matriz de covariables. En el documento, a veces incluía tendencias lineales específicas del estado, pero para este análisis, presento el modelo más simple. Por último, \(\varepsilon_{st}\) es un término de error estructural que se supone condicionalmente independiente de los regresores. Además, todos los errores estándar fueron agruparon a nivel estatal, lo que permite una correlación serial arbitraria.Sin embargo, identificar de forma creíble los efectos causales requiere tanto encontrar efectos como descartar explicaciones alternativas. por ello, los autores presentan pruebas de una diferencia triple en la que se utiliza una cohorte no tratada como control dentro del estado. Para llevarlo a cabo, los autores seleccionaron a los jóvenes de 25 a 29 años de los mismos estados como grupos de comparación dentro de los estados en lugar de a los jóvenes de 20 a 24 años dado que necesitaban un grupo de edad que estuviera lo suficientemente cerca como para captar las tendencias comunes, pero lo suficientemente lejos como para no violar la SUTVA. Dado que los jóvenes de 15 a 19 años tenían más probabilidades que los de 25 a 29 años de mantener relaciones sexuales con jóvenes de 20 a 24 años, eligieron el grupo ligeramente mayor como control dentro de la etapa. En este caso el modelo planteado fue el siguiente:
\[Y_st = \beta_1Repeals + \beta_2DT_t + \beta_3 Repeal_s DT_t + \delta_1 DA + \delta_2 Repeal_s DA +\delta_{3t}DA.DT_t + \delta_{4t} Repeals_s.DA.DT_t + X_{st}\xi + \alpha_{1s}DS_s + \alpha_{2s}DS_s.DA + \gamma_1t +\gamma_{2s}DS_s.t + \gamma_3DA.t+\gamma_{4s}DS_s.DA.t +\varepsilon_{st}\] Dónde el parámetro DDD que estamos estimando es \(\delta_{4t}\) (la interacción completa). en el modelo hay 7 variables ficticias separadas porque el parámetro DDD tiene las tres interacciones. Por lo tanto, como hay ocho combinaciones, los autores tuvieron que descartar una como grupo omitido y controlar por separado las otras siete| Model 1 | |
|---|---|
| (Intercept) | 8.434*** |
| repeal1 | -0.588 |
| year1986 | 0.005 |
| year1987 | -0.087 |
| year1988 | 0.072 |
| year1989 | 0.185 |
| year1990 | 0.186 |
| year1991 | 0.164 |
| year1992 | 0.017 |
| year1993 | -0.255 |
| year1994 | -0.067 |
| year1995 | -0.311 |
| year1996 | -0.369 |
| year1997 | -0.320 |
| year1998 | -0.146 |
| year1999 | -0.148 |
| year2000 | -0.188 |
| acc | -0.003 |
| ir | -0.0002 |
| alcohol | 0.213 |
| crack | -0.086 |
| poverty | -0.002 |
| income | -0.00003 |
| ur | -0.009 |
| repeal1 × year1986 | 0.174 |
| repeal1 × year1987 | 0.269+ |
| repeal1 × year1988 | 0.051 |
| repeal1 × year1989 | -0.115 |
| repeal1 × year1990 | -0.114 |
| repeal1 × year1991 | -0.150 |
| repeal1 × year1992 | -0.095 |
| repeal1 × year1993 | 0.174 |
| repeal1 × year1994 | 0.072 |
| repeal1 × year1995 | 0.241 |
| repeal1 × year1996 | 0.075 |
| repeal1 × year1997 | -0.050 |
| repeal1 × year1998 | -0.033 |
| repeal1 × year1999 | -0.027 |
| repeal1 × year2000 | 0.019 |
| Num.Obs. | 1435 |
| R2 | 0.464 |
| R2 Adj. | 0.449 |
| RMSE | 0.77 |
| Std.Errors | by: fip |