Estoy tratando de volver a crear un análisis realizado utilizando Stata función xtreg (aunque no tengo el código) con el R paquete plm, y estoy teniendo problemas para traducir entre los dos. Un ejemplo mínimo de la modelo que estoy tratando de cálculo sería como sigue:
$$ \text{Dependiente}_{i,t} = \beta_0 + \gamma_i + \delta_t + \beta_1*\text{Dependiente}_{i,t-1} + \beta_2*\text{Independiente}_{i,t} + \epsilon_{i,t} $$
where $\gamma_i$ is a control for the cross-sectional factor (States), $\delta_t$ is a control for each time period (Year), $\text{Dependiente}_{i,t}$ is a dependent variable in the form of a proportion, and $\text{Independiente}_{i,t}$ es un atributo de un registro que cambia con el tanto de la sección transversal y el componente de tiempo de los datos.
He intentado estimar los coeficientes para el modelo de uso de la R con el siguiente (generalizada) código:
model <- plm(Dependent_it ~ Independent_it + lag(Dependent_it,1),
data = data,
method = "within", #fixed effects model
effect = "twoway", #does the gamma_i and delta_t parts (I think)
index = c("State", "Year")
)
He tratado de volver a crear en Stata mediante la siguiente (después de ejecutar xtset State Year, yearly del curso)
xtreg Dependent Independent l.Dependent i.Year, fe
Puedo obtener resultados diferentes, así que debo estar haciendo algo mal, ya sea con xtreg o plm, o ambos. Alguien puede aclararme si he codificado el modelo de arriba adecuadamente en cualquiera de los idiomas?
