'coefficients for random effects/coefficients
Display (3.4) from Laird, Lange, Stram (see help(Unit))
D Z\' W r'

def _compute_b(self, D):

    self.b = np.dot(D, np.dot(np.dot(self.Z.T, self.W), self.r))
