linear_cre function

Main Function for fitting correlated random effect linear model