def METHOD_NAME(image, sigma_threshold=3):
    """Calculate the sigma-clipped mean and stdev of an input array
    Parameters
    ----------
    image : numpy.ndarray
        Array of which to calculate statistics
    sigma_threshold : float
        Number of sigma to use when sigma-clipping
    Returns
    -------
    mean_value : float
        Sigma-clipped mean of image
    stdev_value : float
        Sigma-clipped standard deviation of image
    """
    clipped, lower, upper = sigmaclip(image, low=sigma_threshold, high=sigma_threshold)
    mean_value = np.mean(clipped)
    stdev_value = np.std(clipped)
    return mean_value, stdev_value