center_grand_mean function

Center with respect to grand mean