def statistics(x, axis=(0, 1, 2)):
    """Return the mean and standard deviation of a dataset.

    Args:
        x (np.ndarray): Dataset for which statistics are computed.
        axis (None or int or tuple of ints): Axis or axes along which
            the statistics are computed.

    Returns:
        tuple: A tuple containing the mean and standard deviation.
    """
    mean = x.mean(axis)
    std = x.std(axis)
    return mean, std