def METHOD_NAME(self):
    z_ints = [el**2 for el in self.y_ints]
    dataset = Dataset({'x':self.xs, 'y':self.y_ints, 'z':z_ints},
                      kdims=['x', 'y'], vdims=['z'])
    self.assertEqual(dataset.aggregate(['x'], np.mean).sort(),
                     Dataset({'x':self.xs, 'z':z_ints}, kdims=['x'], vdims=['z']).sort())