def METHOD_NAME(failure_logger, mask_dataset):
    mask_value, X_np, X = mask_dataset
    mode = _masked_column_mode(X, mask_value).get()
    mask = ~sk_get_mask(X_np, value_to_mask=mask_value)
    n_columns = X.shape[1]
    for i in range(n_columns):
        column_mask = mask[:, i]
        column_mode = stats.mode(X_np[:, i][column_mask], keepdims=True)[0][0]
        assert column_mode == mode[i]