def METHOD_NAME():
    # This should be IID, resulting in a larger p-value
    samples = []
    labels = []
    for _ in range(num_samples):
        label = np.random.choice(num_classes)
        mean, sigma = class_stats[label]
        sample = np.random.normal(mean, sigma)
        samples.append(sample)
        labels.append(label)
    samples = np.array(samples)
    labels = np.array(labels)
    dataset = {"features": samples, "labels": labels}
    return dataset