diff --git a/thicket/tests/conftest.py b/thicket/tests/conftest.py
index fc717833..d7379096 100644
--- a/thicket/tests/conftest.py
+++ b/thicket/tests/conftest.py
@@ -51,6 +51,35 @@ def columnar_join_thicket(mpi_scaling_cali, rajaperf_basecuda_xl_cali):
     return thicket_list, thicket_list_cp, combined_th
 
 
+@pytest.fixture
+def stats_columnar_join_thicket(rajaperf_basecuda_xl_cali):
+    """Generator for 'columnar_join' thicket for test_stats.py.
+
+    Arguments:
+        mpi_scaling_cali (list): List of Caliper files for MPI scaling study.
+        rajaperf_basecuda_xl_cali (list): List of Caliper files for base cuda variant.
+
+    Returns:
+        list: List of original thickets, list of deepcopies of original thickets, and
+            columnar-joined thicket.
+    """
+    th_cuda128_1 = Thicket.from_caliperreader(rajaperf_basecuda_xl_cali[0:4])
+    th_cuda128_2 = Thicket.from_caliperreader(rajaperf_basecuda_xl_cali[5:9])
+
+    # To check later if modifications were unexpectedly made
+    th_cuda128_1_deep = th_cuda128_1.deepcopy()
+    th_cuda128_2_deep = th_cuda128_2.deepcopy()
+    thicket_list = [th_cuda128_1, th_cuda128_2]
+    thicket_list_cp = [th_cuda128_1_deep, th_cuda128_2_deep]
+
+    combined_th = Thicket.columnar_join(
+        thicket_list=thicket_list,
+        header_list=["Cuda 1", "Cuda 2"],
+    )
+
+    return thicket_list, thicket_list_cp, combined_th
+
+
 @pytest.fixture
 def data_dir():
     """Return path to the top-level data directory for tests."""
diff --git a/thicket/tests/test_stats.py b/thicket/tests/test_stats.py
new file mode 100644
index 00000000..360eb89a
--- /dev/null
+++ b/thicket/tests/test_stats.py
@@ -0,0 +1,559 @@
+# Copyright 2022 Lawrence Livermore National Security, LLC and other
+# Thicket Project Developers. See the top-level LICENSE file for details.
+#
+# SPDX-License-Identifier: MIT
+
+import thicket as th
+
+
+def test_mean(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.mean(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_mean" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_mean"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_mean" in th_ens.statsframe.show_metric_columns()
+
+
+def test_mean_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.mean(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (idx, "Min time/rank_mean") in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_mean",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (idx, "Min time/rank_mean") in combined_th.statsframe.show_metric_columns()
+
+
+def test_median(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.median(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_median" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_median"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_median" in th_ens.statsframe.show_metric_columns()
+
+
+def test_median_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.median(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (idx, "Min time/rank_median") in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_median",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (idx, "Min time/rank_median") in combined_th.statsframe.show_metric_columns()
+
+
+def test_minimum(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.minimum(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_min" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_min"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_min" in th_ens.statsframe.show_metric_columns()
+
+
+def test_minimum_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.minimum(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (idx, "Min time/rank_min") in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_min",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (idx, "Min time/rank_min") in combined_th.statsframe.show_metric_columns()
+
+
+def test_maximum(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.maximum(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_max" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_max"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_max" in th_ens.statsframe.show_metric_columns()
+
+
+def test_maximum_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.maximum(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (idx, "Min time/rank_max") in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_max",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (idx, "Min time/rank_max") in combined_th.statsframe.show_metric_columns()
+
+
+def test_std(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.std(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_std" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_std"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_std" in th_ens.statsframe.show_metric_columns()
+
+
+def test_std_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.std(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (idx, "Min time/rank_std") in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_std",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (idx, "Min time/rank_std") in combined_th.statsframe.show_metric_columns()
+
+
+def test_percentiles(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.percentiles(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_percentiles" in th_ens.statsframe.dataframe.columns
+    assert len(th_ens.statsframe.dataframe["Min time/rank_percentiles"][0]) == 3
+
+    assert (
+        "Min time/rank_percentiles"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_percentiles" in th_ens.statsframe.show_metric_columns()
+
+
+def test_percentiles_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.percentiles(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (
+        idx,
+        "Min time/rank_percentiles",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_percentiles",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_percentiles",
+    ) in combined_th.statsframe.show_metric_columns()
+
+
+def test_variance(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.variance(th_ens, columns=["Min time/rank"])
+
+    assert "Min time/rank_var" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_var"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_var" in th_ens.statsframe.show_metric_columns()
+
+
+def test_variance_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.variance(combined_th, columns=[(idx, "Min time/rank")])
+
+    assert (idx, "Min time/rank_var") in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_var",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (idx, "Min time/rank_var") in combined_th.statsframe.show_metric_columns()
+
+
+def test_normality(rajaperf_basecuda_xl_cali):
+    th_ens = th.Thicket.from_caliperreader(rajaperf_basecuda_xl_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.check_normality(th_ens, columns=["Min time/rank"])
+
+    assert th_ens.statsframe.dataframe["Min time/rank_normality"][0] in {
+        "True",
+        "False",
+    }
+
+    assert "Min time/rank_normality" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_normality"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert "Min time/rank_normality" in th_ens.statsframe.show_metric_columns()
+
+
+def test_normality_columnar_join(columnar_join_thicket, stats_columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    sthicket_list, sthicket_list_cp, scombined_th = stats_columnar_join_thicket
+    # new data must be added before uncommenting, need 3 or more datapoints
+    # idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    for idx in ["Cuda 1", "Cuda 2"]:
+        th.check_normality(scombined_th, columns=[(idx, "Min time/rank")])
+
+        assert (
+            idx,
+            "Min time/rank_normality",
+        ) in scombined_th.statsframe.dataframe.columns
+        assert (
+            (
+                idx,
+                "Min time/rank_normality",
+            )
+            in scombined_th.statsframe.exc_metrics + scombined_th.statsframe.inc_metrics
+        )
+        assert (
+            idx,
+            "Min time/rank_normality",
+        ) in scombined_th.statsframe.show_metric_columns()
+
+
+def test_correlation(rajaperf_basecuda_xl_cali):
+    th_ens = th.Thicket.from_caliperreader(rajaperf_basecuda_xl_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.correlation_nodewise(
+        th_ens, column1="Min time/rank", column2="Max time/rank", correlation="pearson"
+    )
+
+    assert (
+        "Min time/rank_vs_Max time/rank pearson" in th_ens.statsframe.dataframe.columns
+    )
+
+
+def test_correlation_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.correlation_nodewise(
+        combined_th,
+        column1=(idx[0], "Min time/rank"),
+        column2=(idx[1], "Max time/rank"),
+        correlation="spearman",
+    )
+
+    assert (
+        "Union statistics",
+        "Min time/rank_vs_Max time/rank spearman",
+    ) in combined_th.statsframe.dataframe.columns
+
+
+def test_boxplot(example_cali):
+    th_ens = th.Thicket.from_caliperreader(example_cali)
+
+    assert sorted(th_ens.dataframe.index.get_level_values(0).unique()) == sorted(
+        th_ens.statsframe.dataframe.index.values
+    )
+
+    assert list(th_ens.statsframe.dataframe.columns) == ["name"]
+
+    th.calc_boxplot_statistics(
+        th_ens, columns=["Min time/rank"], quartiles=[0.25, 0.50, 0.75]
+    )
+
+    assert "Min time/rank_q1(0.25, 0.5, 0.75)" in th_ens.statsframe.dataframe.columns
+    assert "Min time/rank_q2(0.25, 0.5, 0.75)" in th_ens.statsframe.dataframe.columns
+    assert "Min time/rank_q3(0.25, 0.5, 0.75)" in th_ens.statsframe.dataframe.columns
+    assert "Min time/rank_iqr(0.25, 0.5, 0.75)" in th_ens.statsframe.dataframe.columns
+    assert (
+        "Min time/rank_lowerfence(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.dataframe.columns
+    )
+    assert (
+        "Min time/rank_upperfence(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.dataframe.columns
+    )
+    assert (
+        "Min time/rank_outliers(0.25, 0.5, 0.75)" in th_ens.statsframe.dataframe.columns
+    )
+
+    assert (
+        "Min time/rank_q1(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert (
+        "Min time/rank_q2(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert (
+        "Min time/rank_q3(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert (
+        "Min time/rank_iqr(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert (
+        "Min time/rank_lowerfence(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert (
+        "Min time/rank_upperfence(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+    assert (
+        "Min time/rank_outliers(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.exc_metrics + th_ens.statsframe.inc_metrics
+    )
+
+    assert (
+        "Min time/rank_q1(0.25, 0.5, 0.75)" in th_ens.statsframe.show_metric_columns()
+    )
+    assert (
+        "Min time/rank_q2(0.25, 0.5, 0.75)" in th_ens.statsframe.show_metric_columns()
+    )
+    assert (
+        "Min time/rank_q3(0.25, 0.5, 0.75)" in th_ens.statsframe.show_metric_columns()
+    )
+    assert (
+        "Min time/rank_iqr(0.25, 0.5, 0.75)" in th_ens.statsframe.show_metric_columns()
+    )
+    assert (
+        "Min time/rank_lowerfence(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.show_metric_columns()
+    )
+    assert (
+        "Min time/rank_upperfence(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.show_metric_columns()
+    )
+    assert (
+        "Min time/rank_outliers(0.25, 0.5, 0.75)"
+        in th_ens.statsframe.show_metric_columns()
+    )
+
+
+def test_boxplot_columnar_join(columnar_join_thicket):
+    thicket_list, thicket_list_cp, combined_th = columnar_join_thicket
+    idx = combined_th.dataframe.columns.levels[0][0]
+    assert sorted(combined_th.dataframe.index.get_level_values(0).unique()) == sorted(
+        combined_th.statsframe.dataframe.index.values
+    )
+
+    assert list(combined_th.statsframe.dataframe.columns) == [("name", "")]
+
+    th.calc_boxplot_statistics(
+        combined_th, columns=[(idx, "Min time/rank")], quartiles=[0.25, 0.50, 0.75]
+    )
+
+    assert (
+        idx,
+        "Min time/rank_q1(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_q2(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_q3(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_iqr(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_lowerfence(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_upperfence(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+    assert (
+        idx,
+        "Min time/rank_outliers(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.dataframe.columns
+
+    assert (
+        idx,
+        "Min time/rank_q1(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_q2(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_q3(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_iqr(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_lowerfence(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_upperfence(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+    assert (
+        idx,
+        "Min time/rank_outliers(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.exc_metrics + combined_th.statsframe.inc_metrics
+
+    assert (
+        idx,
+        "Min time/rank_q1(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()
+    assert (
+        idx,
+        "Min time/rank_q2(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()
+    assert (
+        idx,
+        "Min time/rank_q3(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()
+    assert (
+        idx,
+        "Min time/rank_iqr(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()
+    assert (
+        idx,
+        "Min time/rank_lowerfence(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()
+    assert (
+        idx,
+        "Min time/rank_upperfence(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()
+    assert (
+        idx,
+        "Min time/rank_outliers(0.25, 0.5, 0.75)",
+    ) in combined_th.statsframe.show_metric_columns()