def METHOD_NAME(failure_logger, sparse_random_dataset):
    X_np, _, _, X_sparse = sparse_random_dataset
    if X_sparse.format != "csc":
        pytest.skip("Skip non CSC matrices")
    means, variances = csc_mean_variance_axis0(X_sparse)
    ref_means = np.nanmean(X_np, axis=0)
    ref_variances = np.nanvar(X_np, axis=0)
    assert_allclose(means, ref_means)
    assert_allclose(variances, ref_variances)