'Compute the VECM\'s loglikelihood.'

@cache_readonly
def llf(self):

    K = self.neqs
    T = self.nobs
    r = self.coint_rank
    (s00, _, _, _, _, lambd, _) = _sij(self._delta_x, self._delta_y_1_T, self._y_lag1)
    return ((((((- K) * T) * log((2 * math.pi))) / 2) - ((T * (log(np.linalg.det(s00)) + sum(np.log((1 - lambd))[:r]))) / 2)) - ((K * T) / 2))
