[dev, dl, d2l] = self.loglike(tb, y, X, z, z1)