def score(self, p, x, y, sample_weight=None):
        """
        Overloading default regression score method

        Parameters
        ----------
        p: ndarray
            (Ns, d) array of predictor variables (Ns samples, d dimensions)
            for regression
        x: ndarray
            ndarray of (x, y) points. Needs to be a (Ns, 2) array
            corresponding to the lon/lat, for example.
            array of Points, (x, y, z) pairs of shape (N, 3) for 3d kriging
        y: ndarray
            array of targets (Ns, )
        """

        return r2_score(y_pred=self.predict(p, x),
                        y_true=y,
                        sample_weight=sample_weight)