linear_re function

Main Function for fitting the random effect linear model