fixest feols demo

fixest feols demo

library(fixest)
library(broom)
library(tidyverse)
library(janitor)

options(scipen = 999)

df <- gapminder::gapminder

mlist <- list()

mlist$m1 <- feols(lifeExp ~ gdpPercap,
                  panel.id = c('country', 'year'),
                  data = df)

mlist$m2 <- feols(lifeExp ~ gdpPercap | country,
                  panel.id = c('country', 'year'),
                  data = df)

mlist$m3 <- feols(lifeExp ~ gdpPercap | country + year,
                  panel.id = c('country', 'year'),
                  data = df)

mlist$m4 <- feols(lifeExp ~ gdpPercap + pop | country + year,
                  panel.id = c('country', 'year'),
                  data = df)

mlist$m5 <- feols(lifeExp ~ gdpPercap + pop | country + year,
                  panel.id = c('country', 'year'),
                  cluster = c('continent'),
                  data = df)

etable(mlist)

coefs_exp <- coeftable(mlist$m4) %>% 
  as.data.frame() %>% 
  clean_names() %>% 
  rownames_to_column(., "predictor") %>% 
  mutate(conf_low = estimate - (1.96 * std_error),
         conf_high = estimate + (1.96 * std_error))