code: def cdf(self, y, f, var):
          r
          var = self._check_param(var)
          return norm.cdf(y, loc=f, scale=np.sqrt(var))
docs: r"""
          Cumulative density function of the likelihood.
  
          Parameters
          ----------
          y: ndarray
              query quantiles, i.e.\  :math:`P(Y \leq y)`.
          f: ndarray
              latent function from the GLM prior (:math:`\mathbf{f} =
              \boldsymbol\Phi \mathbf{w}`)
          var: float, ndarray, optional
              The variance of the distribution, if not input, the initial value
              of variance is used.
  
          Returns
          -------
          cdf: ndarray
              Cumulative density function evaluated at y.