linear_fe function

Main function for fitting the fixed effect linear model