Example #1
0
def test_stouffers_rfx():
    """
    Smoke test for Stouffer's RFX.
    """
    meta = ibma.Stouffers(inference='rfx', null='theoretical', n_iters=None)
    meta.fit(pytest.dset_z)
    assert isinstance(meta.results, nimare.base.MetaResult)
Example #2
0
def test_z_perm():
    """
    Smoke test for z permutation.
    """
    meta = ibma.Stouffers(inference='rfx', null='empirical', n_iters=10)
    meta.fit(pytest.dset_z)
    assert isinstance(meta.results, nimare.base.MetaResult)
Example #3
0
def test_Stouffers_weighted(testdata_ibma):
    """
    Smoke test for Stouffer's, weighted by sample size.
    """
    meta = ibma.Stouffers(use_sample_size=True)
    res = meta.fit(testdata_ibma)
    assert isinstance(meta.results, nimare.results.MetaResult)
    assert isinstance(res, nimare.results.MetaResult)
Example #4
0
dset = nimare.dataset.Dataset(dset_file)
dset.update_path(dset_dir)
# Calculate missing images
dset.images = nimare.transforms.transform_images(dset.images,
                                                 target='z',
                                                 masker=dset.masker,
                                                 metadata_df=dset.metadata)
dset.images = nimare.transforms.transform_images(dset.images,
                                                 target='varcope',
                                                 masker=dset.masker,
                                                 metadata_df=dset.metadata)

###############################################################################
# Stouffer's
# --------------------------------------------------
meta = ibma.Stouffers(use_sample_size=False)
meta.fit(dset)
plot_stat_map(meta.results.get_map('z'),
              cut_coords=[0, 0, -8],
              draw_cross=False,
              cmap='RdBu_r')

###############################################################################
# Stouffer's with weighting by sample size
# -----------------------------------------------------------------------------
meta = ibma.Stouffers(use_sample_size=True)
meta.fit(dset)
plot_stat_map(meta.results.get_map('z'),
              cut_coords=[0, 0, -8],
              draw_cross=False,
              cmap='RdBu_r')