def sample_percentiles(self):
        """Sample percentiles"""
        _check_for(self.dist, "cdf")
        if self._is_frozen:
            return self.dist.cdf(self.sorted_data)
        quantiles = (self.sorted_data - self.fit_params[-2]) / self.fit_params[
            -1
        ]
        return self.dist.cdf(quantiles)