base = iris
names(base) = c("y", "x1", "x2", "x3", "species")
# Multiple estimation
res = feols(y ~ csw(x1, x2, x3), base, split = ~species)
# Let's the results for the setosa species
res$setosa
# now for versicolor
etable(res$versicolor)
Run the code above in your browser using DataLab