def test_MKDAChi2_fwe_1core(testdata_cbma):
"""Smoke test for MKDAChi2."""
meta = MKDAChi2()
res = meta.fit(testdata_cbma, testdata_cbma)
valid_methods = FWECorrector.inspect(res)
assert "montecarlo" in valid_methods
corr = FWECorrector(method="montecarlo", n_iters=5, n_cores=1)
cres = corr.transform(res)
assert isinstance(res, nimare.results.MetaResult)
assert isinstance(cres, nimare.results.MetaResult)
def test_MKDAChi2_fwe_1core(testdata_cbma):
"""Smoke test for MKDAChi2."""
meta = MKDAChi2()
res = meta.fit(testdata_cbma, testdata_cbma)
valid_methods = FWECorrector.inspect(res)
assert "montecarlo" in valid_methods
corr = FWECorrector(method="montecarlo", n_iters=5, n_cores=1)
cres = corr.transform(res)
assert isinstance(res, nimare.results.MetaResult)
assert isinstance(cres, nimare.results.MetaResult)
assert ("values_desc-pFgA_level-voxel_corr-fwe_method-montecarlo"
in cres.estimator.null_distributions_.keys())
assert ("values_desc-pAgF_level-voxel_corr-fwe_method-montecarlo"
in cres.estimator.null_distributions_.keys())
# Multiple comparisons correction in coordinate-based meta-analyses
# -----------------------------------------------------------------------------
# .. tip::
# For more information multiple comparisons correction and CBMA in NiMARE,
# see :ref:`multiple comparisons correction`.
from nimare.meta.cbma.ale import ALE
# First, we need to fit the Estimator to the Dataset.
meta = ALE(null_method="approximate")
results = meta.fit(dset)
# We can check which FWE correction methods are available for the ALE Estimator
# with the ``inspect`` class method.
from nimare.correct import FWECorrector
print(FWECorrector.inspect(results))
###############################################################################
# Apply the Corrector to the MetaResult
# =============================================================================
# Now that we know what FWE correction methods are available, we can use one.
#
# The "montecarlo" method is a special one that is implemented within the
# Estimator, rather than in the Corrector.
corr = FWECorrector(method="montecarlo", n_iters=50, n_cores=2)
cres = corr.transform(results)
DISTS_TO_PLOT = [
"values_desc-size_level-cluster_corr-fwe_method-montecarlo",
"values_desc-mass_level-cluster_corr-fwe_method-montecarlo",
"values_level-voxel_corr-fwe_method-montecarlo",