Learn R Programming

⚠️There's a newer version (0.12.1) of this package.Take me there.

fixest: Fast and user-friendly fixed-effects estimation

The fixest package offers a family of functions to perform estimations with multiple fixed-effects in both an OLS and a GLM context. Please refer to the introduction for a walk-through.

At the time of writing of this page (February 2020), fixest is the fastest existing method to perform fixed-effects estimations, often by orders of magnitude. See below for a benchmarking with the fastest alternative software.

Benchmarking

Here is a comparison of the performance of fixest functions to other state of the art methods to perform estimations with multiple fixed-effects. The results are reported in the five figures below. Package fixest (black lines) is consistently faster in all situations.

Setup

The benchmarking was performed as follows: In the OLS context, we estimate the following equation:

The same functional form (one variable, three fixed-effects) is estimated for the Poisson, the Negative Binomial and the Logit cases (with ad hoc modifications to fit each model). See Berge (2018) for more details on the setup.

For the "difficult" benchmark (OLS only), the data is generated in a way that makes the convergence of the fixed-effects slow. The phenomenon of slow convergence is frequent for real micro-level data sets involving employee and firm fixed-effects for instance.

Each estimation is replicated 10 times and the average computing time is reported in the figures.

The alternative methods used for comparison are:

  • OLS: felm (R: package lfe), reghdfe (Stata) and FixedEffectModels (Julia)
  • Poisson: glmmboot (R: package glmmML), feglm (R: package alpaca) and ppmlhdfe (Stata)
  • Negative Binomial: glm.nb (R: package MASS) and nbreg (Stata)
  • Logit: glmmboot (R: package glmmML), feglm (R: package alpaca) and logit (Stata)

All the aforementioned packages were updated at the benchmarking date: February 2020.

The code and data for the benchmarking can be found in this folder.

Acknowledgements

Of course the development of fixest has been inspired and pushed forward by (almost all) these (great) packages used in the benchmarking and I am deeply indebted to their authors. Although fixest contains many features, some are still uncovered and you should definitely have a look at these packages.

Copy Link

Version

Install

install.packages('fixest')

Monthly Downloads

24,077

Version

0.5.0

License

GPL-3

Maintainer

Laurent Berge

Last Published

June 10th, 2020

Functions in fixest (0.5.0)

deviance.fixest

Extracts the deviance of a fixest estimation
base_did

Sample data for difference in difference
BIC.fixest

Bayesian information criterion
collinearity

Collinearity diagnostics for fixest objects
demean

Centers a set of variables around a set of factors
confint.fixest

Confidence interval for parameters estimated with fixest
AIC.fixest

Aikake's an information criterion
coefplot

Plots confidence intervals and point estimates
coef.fixest

Extracts the coefficients from a fixest estimation
coeftable

Obtain various statistics from an estimation
femlm

Fixed-effects maximum likelihood model
feols

Fixed-effects OLS estimation
feNmlm

Fixed effects nonlinear maximum likelihood models
feglm

Fixed-effects GLM estimations
i

Interact variables with factors
print.fixest

A print facility for fixest objects.
model.matrix.fixest

Design matrix of a femlm model
f

Lags a variable in a fixest estimation
nobs.fixest

Extracts the number of observations form a fixest object
setFixest_dict

Sets/gets the dictionary relabeling the variables
fixef_reexported

Functions exported from nlme to implement fixest methods
fixest-deprecated

Deprecated functions in package fixest.
weights.fixest

Extracts the weights from a fixest object
setFixest_dof

Sets the default type of DoF correction in summary/vcov.fixest
predict.fixest

Predict method for fixest fits
plot.fixest.fixef

Displaying the most notable fixed-effects
fixef.fixest

Extract the Fixed-Effects from a fixest estimation.
dof

Type of degree of freedom in fixest summary
fitted.fixest

Extracts fitted values from a fixest fit
resid.fixest

Extracts residuals from a fixest object
setFixest_coefplot

Sets the defaults of coefplot
did_plot_yearly_effects

Plots the results of yearly treatment effects estimation
did_means

Treated and control sample descriptives
lag.formula

Lags a variable using a formula
logLik.fixest

Extracts the log-likelihood
r2

R2s of fixest models
sigma.fixest

Residual standard deviation of fixest estimations
terms.fixest

Extract the terms
obs2remove

Finds observations to be removed from ML estimation with fixed-effects
setFixest_nthreads

Sets/gets the number of threads to use in fixest functions
setFixest_print.type

Sets/gets what print does to fixest estimations
setFixest_notes

Sets/gets whether to display notes in fixest estimation functions
did_estimate_yearly_effects

Estimates yearly treatment effects
summary.fixest.fixef

Summary method for fixed-effects coefficients
fixest-package

Fast and User-Friendly Fixed-Effects Estimations
errbar

Plots confidence intervals
etable

Estimations table (export the results of multiples estimations to a DF or to Latex)
summary.fixest.obs2remove

Summary method for fixest.obs2remove objects
formula.fixest

Extract the formula of a fixest fit
summary.fixest

Summary of a fixest object. Computes different types of standard errors.
unpanel

Dissolves a fixest panel
[.fixest_panel

Method to subselect from a fixest_panel
panel

Constructs a fixest panel data base
xpd

Sets/gets and expands formula macros
to_integer

Fast transform of any type of vector(s) into an integer vector
update.fixest

Updates a fixest estimation
trade

Trade data sample
vcov.fixest

Extracts the variance/covariance of a femlm fit
setFixest_na_inf.rm

Sets/gets whether to remove NA/Inf values from fixest estimations