def test_MKDAChi2_fdr(testdata_cbma):
"""Smoke test for MKDAChi2."""
meta = MKDAChi2()
res = meta.fit(testdata_cbma, testdata_cbma)
corr = FDRCorrector(method="bh", alpha=0.001)
cres = corr.transform(res)
assert isinstance(res, nimare.results.MetaResult)
assert isinstance(cres, nimare.results.MetaResult)
def test_mkda_chi2_fdr():
"""
Smoke test for MKDAChi2
"""
meta = mkda.MKDAChi2()
res = meta.fit(pytest.cbma_testdata1, pytest.cbma_testdata1)
corr = FDRCorrector(method='fdr_bh', alpha=0.001)
cres = corr.transform(res)
assert isinstance(res, nimare.base.MetaResult)
assert isinstance(cres, nimare.base.MetaResult)
def test_Fishers(testdata_ibma):
"""
Smoke test for Fisher's.
"""
meta = ibma.Fishers()
res = meta.fit(testdata_ibma)
corr = FDRCorrector(method="indep", alpha=0.001)
cres = corr.transform(res)
assert isinstance(meta.results, nimare.results.MetaResult)
assert isinstance(res, nimare.results.MetaResult)
assert isinstance(cres, nimare.results.MetaResult)
def test_fishers():
"""
Smoke test for Fisher's.
"""
meta = ibma.Fishers()
res = meta.fit(pytest.dset_z)
corr = FDRCorrector(method='indep', alpha=0.001)
cres = corr.transform(res)
assert isinstance(meta.results, nimare.base.MetaResult)
assert isinstance(res, nimare.base.MetaResult)
assert isinstance(cres, nimare.base.MetaResult)
def analyze_and_plot(dset, ground_truth_foci=None, correct=True, return_cres=False):
meta = ALE(kernel__fwhm=10)
results = meta.fit(dset)
if correct:
corr = FDRCorrector()
cres = corr.transform(results)
else:
cres = results
# get the z coordinates
if ground_truth_foci:
stat_map_kwargs = {"cut_coords": [c[2] for c in ground_truth_foci]}
else:
stat_map_kwargs = {}
fig, ax = plt.subplots()
display = plot_stat_map(
cres.get_map("z"),
display_mode="z",
draw_cross=False,
cmap="Purples",
threshold=2.3,
symmetric_cbar=False,
figure=fig,
axes=ax,
**stat_map_kwargs,
)
if ground_truth_foci:
# place red dots indicating the ground truth foci
display.add_markers(ground_truth_foci)
if return_cres:
return fig, cres
return fig