def compute_con(self, con_idx, n_epochs):
        if self.con_scores is None:
            self.con_scores = np.zeros(self.csd_shape)
        pli_mean = self._acc[con_idx] / n_epochs

        # See Vinck paper Eq. (30)
        con = (n_epochs * pli_mean ** 2 - 1) / (n_epochs - 1)

        self.con_scores[con_idx] = con
Compute final con. score for some connections
dataset/ETHPy150Open mne-tools/mne-python/mne/connectivity/spectral.py/_PLIUnbiasedEst.compute_con