feols function

Fixed-effects OLS estimation