[R-es] MODELOS SURE- CON EFECTOS FIJOS USANDO SYSTEMFIT
Antonio Rodriguez Andres
@nton|orodr|guez@ndre@70 @end|ng |rom gm@||@com
Lun Dic 16 12:39:35 CET 2024
Hola estimados usuarios de R
Estoy usando un conjunto de datos de panel, con variables id country y
year. La idea es estimar models separadamente por OLS para cada pais,
y luego estimarlo de forma conjunta usando systemfit
Para ello hago esto
# Running OLS for each country
country_ols <- data_suicide %>%
group_by(country) %>%
group_map(~ lm(log_suicide_males ~ index + gdppercapitappp +
gdpgrowth + alcoholconsumption +
employmentrate + fertilityrate + divorcerate, data = .x))
# Extracting summaries
ols_summaries <- lapply(country_ols, summary)
No he incluido efectos fijos pero cuando uso la siguiente prueba,
siempre me da el mismo error, como poder incorporar solo efectos fijos
por pais y estimar por SURE como Zellner. Entiendo que se puede
obtener los
residuos de cada modelo por OLS y luego estimar conjuntamente las
ecuaciones (hay 20 paises)
Cualquier ayuda es bienvenida
eq1 <- log_suicide_males ~ index + gdppercapitappp + gdpgrowth > eq2
<- log_suicide_males ~ index + gdppercapitappp + gdpgrowth +
fertilityrate> eqSystem <- list( model1 = eq1, mdel2 = eq2 )> fitols
<- systemfit( eqSystem, data = data_suicide )Error in
systemfit(eqSystem, data = data_suicide) :
argument 'formula' must be an object of class 'formula' for panel-like models
Un saludo
Antonio
--
Member, Editorial Committee, *The Economic and Labour Relations Review: *
https://www.cambridge.org/core/journals/the-economic-and-labour-relations-review
https://www.researchgate.net/profile/Antonio_Andres (Research Gate profile)
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es