Two-sample ALE meta-analysis

Meta-analytic projects often involve a number of common steps comparing two or more samples.

In this example, we replicate the ALE-based analyses from Enge et al.1.

A common project workflow with two meta-analytic samples involves the following:

1. Run a within-sample meta-analysis of the first sample.

2. Characterize/summarize the results of the first meta-analysis.

3. Run a within-sample meta-analysis of the second sample.

4. Characterize/summarize the results of the second meta-analysis.

5. Compare the two samples with a subtraction analysis.

6. Compare the two within-sample meta-analyses with a conjunction analysis.

import os

import matplotlib.pyplot as plt
from nilearn.plotting import plot_stat_map

Load Sleuth text files into Datasets

The data for this example are a subset of studies from a meta-analysis on semantic cognition in children 1. A first group of studies probed children’s semantic world knowledge (e.g., correctly naming an object after hearing its auditory description) while a second group of studies asked children to decide if two (or more) words were semantically related to one another or not.

from nimare.io import convert_sleuth_to_dataset
from nimare.utils import get_resource_path

knowledge_file = os.path.join(get_resource_path(), "semantic_knowledge_children.txt")
related_file = os.path.join(get_resource_path(), "semantic_relatedness_children.txt")

knowledge_dset = convert_sleuth_to_dataset(knowledge_file)
related_dset = convert_sleuth_to_dataset(related_file)

Individual group ALEs

Computing separate ALE analyses for each group is not strictly necessary for performing the subtraction analysis but will help the experimenter to appreciate the similarities and differences between the groups.

from nimare.correct import FWECorrector
from nimare.meta.cbma import ALE

ale = ALE(null_method="approximate")
knowledge_results = ale.fit(knowledge_dset)
related_results = ale.fit(related_dset)

corr = FWECorrector(method="montecarlo", voxel_thresh=0.001, n_iters=100, n_cores=2)
knowledge_corrected_results = corr.transform(knowledge_results)
related_corrected_results = corr.transform(related_results)

fig, axes = plt.subplots(figsize=(12, 10), nrows=2)
knowledge_img = knowledge_corrected_results.get_map(
    "z_desc-size_level-cluster_corr-FWE_method-montecarlo"
)
plot_stat_map(
    knowledge_img,
    cut_coords=4,
    display_mode="z",
    title="Semantic knowledge",
    threshold=2.326,  # cluster-level p < .01, one-tailed
    cmap="RdBu_r",
    vmax=4,
    axes=axes[0],
    figure=fig,
)

related_img = related_corrected_results.get_map(
    "z_desc-size_level-cluster_corr-FWE_method-montecarlo"
)
plot_stat_map(
    related_img,
    cut_coords=4,
    display_mode="z",
    title="Semantic relatedness",
    threshold=2.326,  # cluster-level p < .01, one-tailed
    cmap="RdBu_r",
    vmax=4,
    axes=axes[1],
    figure=fig,
)
fig.show()

Out:

  0%|          | 0/100 [00:00<?, ?it/s]
  1%|1         | 1/100 [00:00<01:30,  1.09it/s]
  3%|3         | 3/100 [00:01<00:57,  1.68it/s]
  5%|5         | 5/100 [00:02<00:49,  1.92it/s]
  7%|7         | 7/100 [00:03<00:47,  1.97it/s]
  9%|9         | 9/100 [00:04<00:43,  2.11it/s]
 10%|#         | 10/100 [00:04<00:35,  2.53it/s]
 11%|#1        | 11/100 [00:05<00:44,  2.00it/s]
 12%|#2        | 12/100 [00:05<00:36,  2.42it/s]
 13%|#3        | 13/100 [00:06<00:42,  2.04it/s]
 14%|#4        | 14/100 [00:06<00:33,  2.53it/s]
 15%|#5        | 15/100 [00:07<00:43,  1.95it/s]
 16%|#6        | 16/100 [00:07<00:36,  2.33it/s]
 17%|#7        | 17/100 [00:08<00:40,  2.05it/s]
 18%|#8        | 18/100 [00:08<00:34,  2.37it/s]
 19%|#9        | 19/100 [00:09<00:41,  1.97it/s]
 20%|##        | 20/100 [00:09<00:34,  2.29it/s]
 21%|##1       | 21/100 [00:10<00:37,  2.10it/s]
 22%|##2       | 22/100 [00:10<00:32,  2.40it/s]
 23%|##3       | 23/100 [00:10<00:38,  2.01it/s]
 24%|##4       | 24/100 [00:11<00:33,  2.28it/s]
 25%|##5       | 25/100 [00:11<00:36,  2.08it/s]
 26%|##6       | 26/100 [00:12<00:30,  2.42it/s]
 27%|##7       | 27/100 [00:12<00:37,  1.95it/s]
 28%|##8       | 28/100 [00:13<00:31,  2.32it/s]
 29%|##9       | 29/100 [00:13<00:34,  2.04it/s]
 30%|###       | 30/100 [00:13<00:28,  2.46it/s]
 31%|###1      | 31/100 [00:14<00:35,  1.92it/s]
 32%|###2      | 32/100 [00:14<00:28,  2.35it/s]
 33%|###3      | 33/100 [00:15<00:33,  2.03it/s]
 34%|###4      | 34/100 [00:15<00:25,  2.55it/s]
 35%|###5      | 35/100 [00:16<00:34,  1.90it/s]
 36%|###6      | 36/100 [00:16<00:26,  2.41it/s]
 37%|###7      | 37/100 [00:17<00:31,  1.97it/s]
 38%|###8      | 38/100 [00:17<00:23,  2.59it/s]
 39%|###9      | 39/100 [00:18<00:32,  1.86it/s]
 41%|####1     | 41/100 [00:19<00:28,  2.07it/s]
 42%|####2     | 42/100 [00:19<00:22,  2.55it/s]
 43%|####3     | 43/100 [00:20<00:29,  1.91it/s]
 45%|####5     | 45/100 [00:21<00:26,  2.08it/s]
 46%|####6     | 46/100 [00:21<00:21,  2.53it/s]
 47%|####6     | 47/100 [00:22<00:26,  1.97it/s]
 48%|####8     | 48/100 [00:22<00:21,  2.43it/s]
 49%|####9     | 49/100 [00:22<00:25,  2.00it/s]
 50%|#####     | 50/100 [00:23<00:19,  2.56it/s]
 51%|#####1    | 51/100 [00:23<00:25,  1.91it/s]
 52%|#####2    | 52/100 [00:24<00:20,  2.36it/s]
 53%|#####3    | 53/100 [00:24<00:23,  2.00it/s]
 54%|#####4    | 54/100 [00:25<00:19,  2.38it/s]
 55%|#####5    | 55/100 [00:25<00:23,  1.91it/s]
 56%|#####6    | 56/100 [00:26<00:19,  2.26it/s]
 57%|#####6    | 57/100 [00:26<00:21,  2.04it/s]
 58%|#####8    | 58/100 [00:26<00:18,  2.32it/s]
 59%|#####8    | 59/100 [00:27<00:20,  1.99it/s]
 60%|######    | 60/100 [00:27<00:17,  2.23it/s]
 61%|######1   | 61/100 [00:28<00:18,  2.12it/s]
 62%|######2   | 62/100 [00:28<00:16,  2.32it/s]
 63%|######3   | 63/100 [00:29<00:17,  2.07it/s]
 64%|######4   | 64/100 [00:29<00:16,  2.19it/s]
 65%|######5   | 65/100 [00:30<00:15,  2.21it/s]
 66%|######6   | 66/100 [00:30<00:15,  2.22it/s]
 67%|######7   | 67/100 [00:31<00:15,  2.12it/s]
 68%|######8   | 68/100 [00:31<00:15,  2.12it/s]
 69%|######9   | 69/100 [00:32<00:13,  2.27it/s]
 70%|#######   | 70/100 [00:32<00:13,  2.17it/s]
 71%|#######1  | 71/100 [00:33<00:13,  2.17it/s]
 72%|#######2  | 72/100 [00:33<00:13,  2.07it/s]
 73%|#######3  | 73/100 [00:33<00:11,  2.29it/s]
 74%|#######4  | 74/100 [00:34<00:11,  2.21it/s]
 75%|#######5  | 75/100 [00:34<00:11,  2.17it/s]
 76%|#######6  | 76/100 [00:35<00:11,  2.10it/s]
 77%|#######7  | 77/100 [00:35<00:09,  2.31it/s]
 78%|#######8  | 78/100 [00:36<00:10,  2.12it/s]
 79%|#######9  | 79/100 [00:36<00:09,  2.22it/s]
 80%|########  | 80/100 [00:37<00:10,  1.99it/s]
 81%|########1 | 81/100 [00:37<00:08,  2.35it/s]
 82%|########2 | 82/100 [00:38<00:08,  2.06it/s]
 83%|########2 | 83/100 [00:38<00:07,  2.29it/s]
 84%|########4 | 84/100 [00:39<00:08,  1.94it/s]
 85%|########5 | 85/100 [00:39<00:06,  2.45it/s]
 86%|########6 | 86/100 [00:40<00:07,  2.00it/s]
 87%|########7 | 87/100 [00:40<00:05,  2.33it/s]
 88%|########8 | 88/100 [00:41<00:06,  1.88it/s]
 90%|######### | 90/100 [00:41<00:04,  2.05it/s]
 91%|#########1| 91/100 [00:42<00:03,  2.43it/s]
 92%|#########2| 92/100 [00:43<00:04,  1.89it/s]
 94%|#########3| 94/100 [00:43<00:02,  2.05it/s]
 95%|#########5| 95/100 [00:44<00:02,  2.49it/s]
 96%|#########6| 96/100 [00:44<00:02,  1.93it/s]
 98%|#########8| 98/100 [00:45<00:00,  2.00it/s]
100%|##########| 100/100 [00:46<00:00,  2.35it/s]
100%|##########| 100/100 [00:46<00:00,  2.15it/s]

  0%|          | 0/100 [00:00<?, ?it/s]
  1%|1         | 1/100 [00:00<01:21,  1.22it/s]
  3%|3         | 3/100 [00:01<00:52,  1.85it/s]
  5%|5         | 5/100 [00:02<00:45,  2.08it/s]
  7%|7         | 7/100 [00:03<00:42,  2.17it/s]
  9%|9         | 9/100 [00:04<00:38,  2.37it/s]
 11%|#1        | 11/100 [00:04<00:37,  2.37it/s]
 12%|#2        | 12/100 [00:05<00:31,  2.78it/s]
 13%|#3        | 13/100 [00:05<00:36,  2.38it/s]
 15%|#5        | 15/100 [00:06<00:35,  2.36it/s]
 17%|#7        | 17/100 [00:07<00:33,  2.48it/s]
 19%|#9        | 19/100 [00:08<00:32,  2.47it/s]
 20%|##        | 20/100 [00:08<00:28,  2.77it/s]
 21%|##1       | 21/100 [00:08<00:31,  2.50it/s]
 22%|##2       | 22/100 [00:09<00:27,  2.83it/s]
 23%|##3       | 23/100 [00:09<00:33,  2.31it/s]
 24%|##4       | 24/100 [00:09<00:28,  2.65it/s]
 25%|##5       | 25/100 [00:10<00:30,  2.46it/s]
 26%|##6       | 26/100 [00:10<00:26,  2.77it/s]
 27%|##7       | 27/100 [00:11<00:31,  2.32it/s]
 28%|##8       | 28/100 [00:11<00:27,  2.66it/s]
 29%|##9       | 29/100 [00:12<00:29,  2.41it/s]
 30%|###       | 30/100 [00:12<00:25,  2.77it/s]
 31%|###1      | 31/100 [00:12<00:30,  2.28it/s]
 32%|###2      | 32/100 [00:13<00:25,  2.63it/s]
 33%|###3      | 33/100 [00:13<00:27,  2.41it/s]
 34%|###4      | 34/100 [00:13<00:24,  2.74it/s]
 35%|###5      | 35/100 [00:14<00:28,  2.31it/s]
 36%|###6      | 36/100 [00:14<00:24,  2.61it/s]
 37%|###7      | 37/100 [00:15<00:26,  2.41it/s]
 38%|###8      | 38/100 [00:15<00:22,  2.76it/s]
 39%|###9      | 39/100 [00:16<00:26,  2.27it/s]
 40%|####      | 40/100 [00:16<00:23,  2.59it/s]
 41%|####1     | 41/100 [00:16<00:24,  2.44it/s]
 42%|####2     | 42/100 [00:17<00:21,  2.73it/s]
 43%|####3     | 43/100 [00:17<00:24,  2.31it/s]
 44%|####4     | 44/100 [00:17<00:21,  2.60it/s]
 45%|####5     | 45/100 [00:18<00:22,  2.46it/s]
 46%|####6     | 46/100 [00:18<00:19,  2.75it/s]
 47%|####6     | 47/100 [00:19<00:22,  2.32it/s]
 48%|####8     | 48/100 [00:19<00:19,  2.62it/s]
 49%|####9     | 49/100 [00:19<00:20,  2.47it/s]
 50%|#####     | 50/100 [00:20<00:18,  2.74it/s]
 51%|#####1    | 51/100 [00:20<00:21,  2.30it/s]
 52%|#####2    | 52/100 [00:21<00:18,  2.65it/s]
 53%|#####3    | 53/100 [00:21<00:19,  2.43it/s]
 54%|#####4    | 54/100 [00:21<00:16,  2.74it/s]
 55%|#####5    | 55/100 [00:22<00:19,  2.30it/s]
 56%|#####6    | 56/100 [00:22<00:16,  2.59it/s]
 57%|#####6    | 57/100 [00:23<00:17,  2.48it/s]
 58%|#####8    | 58/100 [00:23<00:15,  2.70it/s]
 59%|#####8    | 59/100 [00:23<00:17,  2.37it/s]
 60%|######    | 60/100 [00:24<00:15,  2.55it/s]
 61%|######1   | 61/100 [00:24<00:15,  2.58it/s]
 62%|######2   | 62/100 [00:25<00:14,  2.59it/s]
 63%|######3   | 63/100 [00:25<00:15,  2.46it/s]
 64%|######4   | 64/100 [00:25<00:14,  2.46it/s]
 65%|######5   | 65/100 [00:26<00:13,  2.66it/s]
 66%|######6   | 66/100 [00:26<00:13,  2.57it/s]
 67%|######7   | 67/100 [00:27<00:13,  2.45it/s]
 68%|######8   | 68/100 [00:27<00:13,  2.45it/s]
 69%|######9   | 69/100 [00:27<00:11,  2.60it/s]
 70%|#######   | 70/100 [00:28<00:11,  2.58it/s]
 71%|#######1  | 71/100 [00:28<00:12,  2.41it/s]
 72%|#######2  | 72/100 [00:29<00:11,  2.45it/s]
 73%|#######3  | 73/100 [00:29<00:10,  2.57it/s]
 74%|#######4  | 74/100 [00:29<00:10,  2.56it/s]
 75%|#######5  | 75/100 [00:30<00:10,  2.42it/s]
 76%|#######6  | 76/100 [00:30<00:09,  2.47it/s]
 77%|#######7  | 77/100 [00:31<00:09,  2.55it/s]
 78%|#######8  | 78/100 [00:31<00:08,  2.54it/s]
 79%|#######9  | 79/100 [00:31<00:08,  2.40it/s]
 80%|########  | 80/100 [00:32<00:08,  2.47it/s]
 81%|########1 | 81/100 [00:32<00:07,  2.54it/s]
 82%|########2 | 82/100 [00:33<00:07,  2.50it/s]
 83%|########2 | 83/100 [00:33<00:06,  2.46it/s]
 84%|########4 | 84/100 [00:33<00:06,  2.38it/s]
 85%|########5 | 85/100 [00:34<00:05,  2.65it/s]
 86%|########6 | 86/100 [00:34<00:05,  2.50it/s]
 87%|########7 | 87/100 [00:35<00:05,  2.47it/s]
 88%|########8 | 88/100 [00:35<00:05,  2.35it/s]
 89%|########9 | 89/100 [00:35<00:04,  2.66it/s]
 90%|######### | 90/100 [00:36<00:04,  2.42it/s]
 91%|#########1| 91/100 [00:36<00:03,  2.52it/s]
 92%|#########2| 92/100 [00:37<00:03,  2.33it/s]
 93%|#########3| 93/100 [00:37<00:02,  2.74it/s]
 94%|#########3| 94/100 [00:37<00:02,  2.46it/s]
 95%|#########5| 95/100 [00:38<00:01,  2.61it/s]
 96%|#########6| 96/100 [00:38<00:01,  2.30it/s]
 97%|#########7| 97/100 [00:38<00:01,  2.72it/s]
 98%|#########8| 98/100 [00:39<00:00,  2.42it/s]
 99%|#########9| 99/100 [00:39<00:00,  2.53it/s]
100%|##########| 100/100 [00:40<00:00,  2.73it/s]
100%|##########| 100/100 [00:40<00:00,  2.48it/s]

Characterize the relative contributions of experiments in the ALE results

NiMARE contains two methods for this: Jackknife and FocusCounter. We will show both below, but for the sake of speed we will only apply one to each subgroup meta-analysis.

from nimare.diagnostics import FocusCounter

counter = FocusCounter(
    target_image="z_desc-size_level-cluster_corr-FWE_method-montecarlo",
    voxel_thresh=None,
)
knowledge_count_table, _ = counter.transform(knowledge_corrected_results)
knowledge_count_table.head(10)

Out:

  0%|          | 0/21 [00:00<?, ?it/s]
 10%|9         | 2/21 [00:00<00:01, 18.55it/s]
 19%|#9        | 4/21 [00:00<00:00, 18.55it/s]
 29%|##8       | 6/21 [00:00<00:00, 18.52it/s]
 38%|###8      | 8/21 [00:00<00:00, 18.44it/s]
 48%|####7     | 10/21 [00:00<00:00, 18.45it/s]
 57%|#####7    | 12/21 [00:00<00:00, 18.35it/s]
 67%|######6   | 14/21 [00:00<00:00, 18.39it/s]
 76%|#######6  | 16/21 [00:00<00:00, 18.41it/s]
 86%|########5 | 18/21 [00:00<00:00, 18.46it/s]
 95%|#########5| 20/21 [00:01<00:00, 18.43it/s]
100%|##########| 21/21 [00:01<00:00, 18.43it/s]

	1	2	3	4	5	6	7
Cluster ID
Center of Mass	(54, -28, 2)	(38, 22, -8)	(0, 20, 46)	(-6, -14, 12)	(-34, 20, -2)	(-42, 10, 28)	(-52, -38, 4)
arnoldussen2006nc-	0	0	0	0	0	1	0
arnoldussen2006rm-	0	0	0	0	0	1	0
backes2002-	0	0	1	0	0	1	0
balsamo2002-	0	0	0	0	0	0	1
balsamo2006-	0	0	1	0	0	1	0
bauer2017-	0	1	0	1	0	0	0
berl2014-	1	1	2	1	0	2	0
brauer2007-	0	0	1	0	0	0	0
gaillard2001-	0	0	0	0	0	1	0

from nimare.diagnostics import Jackknife

jackknife = Jackknife(
    target_image="z_desc-size_level-cluster_corr-FWE_method-montecarlo",
    voxel_thresh=None,
)
related_jackknife_table, _ = jackknife.transform(related_corrected_results)
related_jackknife_table.head(10)

Out:

  0%|          | 0/16 [00:00<?, ?it/s]
  6%|6         | 1/16 [00:01<00:24,  1.60s/it]
 12%|#2        | 2/16 [00:03<00:22,  1.61s/it]
 19%|#8        | 3/16 [00:04<00:20,  1.60s/it]
 25%|##5       | 4/16 [00:06<00:19,  1.59s/it]
 31%|###1      | 5/16 [00:07<00:17,  1.58s/it]
 38%|###7      | 6/16 [00:09<00:15,  1.58s/it]
 44%|####3     | 7/16 [00:11<00:14,  1.58s/it]
 50%|#####     | 8/16 [00:12<00:12,  1.58s/it]
 56%|#####6    | 9/16 [00:14<00:11,  1.58s/it]
 62%|######2   | 10/16 [00:15<00:09,  1.58s/it]
 69%|######8   | 11/16 [00:17<00:07,  1.58s/it]
 75%|#######5  | 12/16 [00:18<00:06,  1.58s/it]
 81%|########1 | 13/16 [00:20<00:04,  1.58s/it]
 88%|########7 | 14/16 [00:22<00:03,  1.58s/it]
 94%|#########3| 15/16 [00:23<00:01,  1.58s/it]
100%|##########| 16/16 [00:25<00:00,  1.58s/it]
100%|##########| 16/16 [00:25<00:00,  1.58s/it]

	1	2	3	4	5
Cluster ID
Center of Mass	(38, 22, -8)	(-4, 16, 50)	(-10, 6, 4)	(-46, 22, 8)	(-56, -44, 0)
booth2001-	0.0	0.0	0.0	0.029415	0.054821
booth2003-	0.0	0.0	0.0	0.0	0.000008
booth2007-	0.021799	0.068333	0.000019	0.051598	0.000028
cao2008-	0.134026	0.06561	0.0	0.035607	0.094222
chou2006a-	0.211268	0.123099	0.0	0.000075	0.095687
chou2006b-	0.21024	0.117018	0.0	0.004724	0.117981
chou2009-	0.188491	0.1575	0.330668	0.22325	0.212472
chou2019-	0.0	0.061734	0.0	0.0	0.0
fan2020-	0.0	0.070051	0.0	0.002654	0.101595

Subtraction analysis

Typically, one would use at least 10000 iterations for a subtraction analysis. However, we have reduced this to 100 iterations for this example.

from nimare.meta.cbma import ALESubtraction

sub = ALESubtraction(n_iters=100, n_cores=1)
res_sub = sub.fit(knowledge_dset, related_dset)
img_sub = res_sub.get_map("z_desc-group1MinusGroup2")

plot_stat_map(
    img_sub,
    cut_coords=4,
    display_mode="z",
    title="Subtraction",
    cmap="RdBu_r",
    vmax=4,
)

Out:

  0%|          | 0/100 [00:00<?, ?it/s]
  1%|1         | 1/100 [00:00<00:45,  2.16it/s]
  2%|2         | 2/100 [00:00<00:43,  2.23it/s]
  3%|3         | 3/100 [00:01<00:43,  2.22it/s]
  4%|4         | 4/100 [00:01<00:43,  2.21it/s]
  5%|5         | 5/100 [00:02<00:43,  2.20it/s]
  6%|6         | 6/100 [00:02<00:42,  2.23it/s]
  7%|7         | 7/100 [00:03<00:41,  2.23it/s]
  8%|8         | 8/100 [00:03<00:41,  2.24it/s]
  9%|9         | 9/100 [00:04<00:40,  2.24it/s]
 10%|#         | 10/100 [00:04<00:40,  2.22it/s]
 11%|#1        | 11/100 [00:04<00:40,  2.22it/s]
 12%|#2        | 12/100 [00:05<00:39,  2.21it/s]
 13%|#3        | 13/100 [00:05<00:39,  2.20it/s]
 14%|#4        | 14/100 [00:06<00:38,  2.21it/s]
 15%|#5        | 15/100 [00:06<00:38,  2.22it/s]
 16%|#6        | 16/100 [00:07<00:37,  2.22it/s]
 17%|#7        | 17/100 [00:07<00:37,  2.24it/s]
 18%|#8        | 18/100 [00:08<00:36,  2.24it/s]
 19%|#9        | 19/100 [00:08<00:36,  2.22it/s]
 20%|##        | 20/100 [00:08<00:35,  2.24it/s]
 21%|##1       | 21/100 [00:09<00:35,  2.23it/s]
 22%|##2       | 22/100 [00:09<00:34,  2.24it/s]
 23%|##3       | 23/100 [00:10<00:34,  2.24it/s]
 24%|##4       | 24/100 [00:10<00:34,  2.20it/s]
 25%|##5       | 25/100 [00:11<00:33,  2.22it/s]
 26%|##6       | 26/100 [00:11<00:32,  2.24it/s]
 27%|##7       | 27/100 [00:12<00:32,  2.23it/s]
 28%|##8       | 28/100 [00:12<00:32,  2.24it/s]
 29%|##9       | 29/100 [00:13<00:31,  2.28it/s]
 30%|###       | 30/100 [00:13<00:30,  2.28it/s]
 31%|###1      | 31/100 [00:13<00:30,  2.26it/s]
 32%|###2      | 32/100 [00:14<00:30,  2.23it/s]
 33%|###3      | 33/100 [00:14<00:30,  2.22it/s]
 34%|###4      | 34/100 [00:15<00:29,  2.24it/s]
 35%|###5      | 35/100 [00:15<00:28,  2.25it/s]
 36%|###6      | 36/100 [00:16<00:28,  2.27it/s]
 37%|###7      | 37/100 [00:16<00:27,  2.29it/s]
 38%|###8      | 38/100 [00:16<00:27,  2.26it/s]
 39%|###9      | 39/100 [00:17<00:26,  2.26it/s]
 40%|####      | 40/100 [00:17<00:26,  2.25it/s]
 41%|####1     | 41/100 [00:18<00:26,  2.25it/s]
 42%|####2     | 42/100 [00:18<00:25,  2.27it/s]
 43%|####3     | 43/100 [00:19<00:24,  2.29it/s]
 44%|####4     | 44/100 [00:19<00:24,  2.30it/s]
 45%|####5     | 45/100 [00:20<00:24,  2.29it/s]
 46%|####6     | 46/100 [00:20<00:23,  2.29it/s]
 47%|####6     | 47/100 [00:20<00:23,  2.27it/s]
 48%|####8     | 48/100 [00:21<00:23,  2.25it/s]
 49%|####9     | 49/100 [00:21<00:22,  2.26it/s]
 50%|#####     | 50/100 [00:22<00:22,  2.23it/s]
 51%|#####1    | 51/100 [00:22<00:21,  2.23it/s]
 52%|#####2    | 52/100 [00:23<00:21,  2.26it/s]
 53%|#####3    | 53/100 [00:23<00:20,  2.27it/s]
 54%|#####4    | 54/100 [00:24<00:20,  2.26it/s]
 55%|#####5    | 55/100 [00:24<00:20,  2.23it/s]
 56%|#####6    | 56/100 [00:24<00:19,  2.22it/s]
 57%|#####6    | 57/100 [00:25<00:19,  2.23it/s]
 58%|#####8    | 58/100 [00:25<00:18,  2.25it/s]
 59%|#####8    | 59/100 [00:26<00:18,  2.24it/s]
 60%|######    | 60/100 [00:26<00:17,  2.24it/s]
 61%|######1   | 61/100 [00:27<00:17,  2.26it/s]
 62%|######2   | 62/100 [00:27<00:16,  2.28it/s]
 63%|######3   | 63/100 [00:28<00:16,  2.28it/s]
 64%|######4   | 64/100 [00:28<00:15,  2.26it/s]
 65%|######5   | 65/100 [00:28<00:15,  2.26it/s]
 66%|######6   | 66/100 [00:29<00:15,  2.25it/s]
 67%|######7   | 67/100 [00:29<00:14,  2.23it/s]
 68%|######8   | 68/100 [00:30<00:14,  2.24it/s]
 69%|######9   | 69/100 [00:30<00:13,  2.26it/s]
 70%|#######   | 70/100 [00:31<00:13,  2.27it/s]
 71%|#######1  | 71/100 [00:31<00:12,  2.28it/s]
 72%|#######2  | 72/100 [00:32<00:12,  2.27it/s]
 73%|#######3  | 73/100 [00:32<00:11,  2.27it/s]
 74%|#######4  | 74/100 [00:32<00:11,  2.26it/s]
 75%|#######5  | 75/100 [00:33<00:11,  2.24it/s]
 76%|#######6  | 76/100 [00:33<00:10,  2.25it/s]
 77%|#######7  | 77/100 [00:34<00:10,  2.15it/s]
 78%|#######8  | 78/100 [00:34<00:10,  2.19it/s]
 79%|#######9  | 79/100 [00:35<00:09,  2.21it/s]
 80%|########  | 80/100 [00:35<00:08,  2.23it/s]
 81%|########1 | 81/100 [00:36<00:08,  2.25it/s]
 82%|########2 | 82/100 [00:36<00:08,  2.23it/s]
 83%|########2 | 83/100 [00:37<00:07,  2.24it/s]
 84%|########4 | 84/100 [00:37<00:07,  2.26it/s]
 85%|########5 | 85/100 [00:37<00:06,  2.28it/s]
 86%|########6 | 86/100 [00:38<00:06,  2.26it/s]
 87%|########7 | 87/100 [00:38<00:05,  2.26it/s]
 88%|########8 | 88/100 [00:39<00:05,  2.28it/s]
 89%|########9 | 89/100 [00:39<00:04,  2.26it/s]
 90%|######### | 90/100 [00:40<00:04,  2.25it/s]
 91%|#########1| 91/100 [00:40<00:03,  2.27it/s]
 92%|#########2| 92/100 [00:40<00:03,  2.28it/s]
 93%|#########3| 93/100 [00:41<00:03,  2.30it/s]
 94%|#########3| 94/100 [00:41<00:02,  2.31it/s]
 95%|#########5| 95/100 [00:42<00:02,  2.30it/s]
 96%|#########6| 96/100 [00:42<00:01,  2.31it/s]
 97%|#########7| 97/100 [00:43<00:01,  2.33it/s]
 98%|#########8| 98/100 [00:43<00:00,  2.32it/s]
 99%|#########9| 99/100 [00:43<00:00,  2.31it/s]
100%|##########| 100/100 [00:44<00:00,  2.31it/s]
100%|##########| 100/100 [00:44<00:00,  2.25it/s]

  0%|          | 0/228483 [00:00<?, ?it/s]
  0%|          | 945/228483 [00:00<00:24, 9444.06it/s]
  1%|          | 1892/228483 [00:00<00:23, 9455.07it/s]
  1%|1         | 2848/228483 [00:00<00:23, 9501.82it/s]
  2%|1         | 3807/228483 [00:00<00:23, 9533.68it/s]
  2%|2         | 4772/228483 [00:00<00:23, 9572.38it/s]
  3%|2         | 5730/228483 [00:00<00:23, 9550.74it/s]
  3%|2         | 6694/228483 [00:00<00:23, 9579.27it/s]
  3%|3         | 7653/228483 [00:00<00:23, 9579.88it/s]
  4%|3         | 8611/228483 [00:00<00:23, 9558.50it/s]
  4%|4         | 9573/228483 [00:01<00:22, 9574.83it/s]
  5%|4         | 10531/228483 [00:01<00:22, 9568.07it/s]
  5%|5         | 11488/228483 [00:01<00:22, 9565.14it/s]
  5%|5         | 12453/228483 [00:01<00:22, 9588.71it/s]
  6%|5         | 13412/228483 [00:01<00:22, 9559.49it/s]
  6%|6         | 14368/228483 [00:01<00:22, 9553.45it/s]
  7%|6         | 15327/228483 [00:01<00:22, 9561.40it/s]
  7%|7         | 16284/228483 [00:01<00:22, 9520.22it/s]
  8%|7         | 17242/228483 [00:01<00:22, 9536.88it/s]
  8%|7         | 18196/228483 [00:01<00:22, 9498.30it/s]
  8%|8         | 19151/228483 [00:02<00:22, 9513.40it/s]
  9%|8         | 20104/228483 [00:02<00:21, 9518.25it/s]
  9%|9         | 21056/228483 [00:02<00:21, 9513.85it/s]
 10%|9         | 22013/228483 [00:02<00:21, 9525.62it/s]
 10%|#         | 22974/228483 [00:02<00:21, 9549.81it/s]
 10%|#         | 23939/228483 [00:02<00:21, 9576.98it/s]
 11%|#         | 24899/228483 [00:02<00:21, 9582.70it/s]
 11%|#1        | 25858/228483 [00:02<00:21, 9532.48it/s]
 12%|#1        | 26815/228483 [00:02<00:21, 9542.42it/s]
 12%|#2        | 27770/228483 [00:02<00:21, 9510.29it/s]
 13%|#2        | 28727/228483 [00:03<00:20, 9526.73it/s]
 13%|#2        | 29680/228483 [00:03<00:20, 9501.04it/s]
 13%|#3        | 30636/228483 [00:03<00:20, 9518.12it/s]
 14%|#3        | 31594/228483 [00:03<00:20, 9536.48it/s]
 14%|#4        | 32548/228483 [00:03<00:29, 6660.29it/s]
 15%|#4        | 33493/228483 [00:03<00:26, 7300.18it/s]
 15%|#5        | 34433/228483 [00:03<00:24, 7817.56it/s]
 15%|#5        | 35382/228483 [00:03<00:23, 8251.87it/s]
 16%|#5        | 36312/228483 [00:03<00:22, 8535.09it/s]
 16%|#6        | 37259/228483 [00:04<00:21, 8794.01it/s]
 17%|#6        | 38199/228483 [00:04<00:21, 8966.43it/s]
 17%|#7        | 39157/228483 [00:04<00:20, 9143.00it/s]
 18%|#7        | 40107/228483 [00:04<00:20, 9246.31it/s]
 18%|#7        | 41057/228483 [00:04<00:20, 9319.86it/s]
 18%|#8        | 42023/228483 [00:04<00:19, 9418.63it/s]
 19%|#8        | 42979/228483 [00:04<00:19, 9458.36it/s]
 19%|#9        | 43933/228483 [00:04<00:19, 9482.66it/s]
 20%|#9        | 44888/228483 [00:04<00:19, 9500.44it/s]
 20%|##        | 45844/228483 [00:04<00:19, 9518.15it/s]
 20%|##        | 46806/228483 [00:05<00:19, 9547.59it/s]
 21%|##        | 47762/228483 [00:05<00:19, 9447.30it/s]
 21%|##1       | 48708/228483 [00:05<00:19, 9450.99it/s]
 22%|##1       | 49671/228483 [00:05<00:18, 9503.43it/s]
 22%|##2       | 50633/228483 [00:05<00:18, 9536.16it/s]
 23%|##2       | 51588/228483 [00:05<00:18, 9518.54it/s]
 23%|##2       | 52541/228483 [00:05<00:18, 9518.16it/s]
 23%|##3       | 53494/228483 [00:05<00:18, 9521.29it/s]
 24%|##3       | 54447/228483 [00:05<00:18, 9500.89it/s]
 24%|##4       | 55398/228483 [00:05<00:18, 9496.53it/s]
 25%|##4       | 56348/228483 [00:06<00:18, 9491.66it/s]
 25%|##5       | 57299/228483 [00:06<00:18, 9496.84it/s]
 25%|##5       | 58249/228483 [00:06<00:17, 9479.33it/s]
 26%|##5       | 59208/228483 [00:06<00:17, 9510.15it/s]
 26%|##6       | 60168/228483 [00:06<00:17, 9534.80it/s]
 27%|##6       | 61124/228483 [00:06<00:17, 9539.51it/s]
 27%|##7       | 62088/228483 [00:06<00:17, 9568.37it/s]
 28%|##7       | 63055/228483 [00:06<00:17, 9596.50it/s]
 28%|##8       | 64020/228483 [00:06<00:17, 9609.83it/s]
 28%|##8       | 64981/228483 [00:06<00:17, 9567.47it/s]
 29%|##8       | 65945/228483 [00:07<00:16, 9588.60it/s]
 29%|##9       | 66904/228483 [00:07<00:16, 9551.70it/s]
 30%|##9       | 67874/228483 [00:07<00:16, 9595.69it/s]
 30%|###       | 68834/228483 [00:07<00:16, 9587.37it/s]
 31%|###       | 69793/228483 [00:07<00:16, 9567.99it/s]
 31%|###       | 70752/228483 [00:07<00:16, 9574.08it/s]
 31%|###1      | 71710/228483 [00:07<00:16, 9561.03it/s]
 32%|###1      | 72672/228483 [00:07<00:16, 9577.85it/s]
 32%|###2      | 73630/228483 [00:07<00:16, 9523.61it/s]
 33%|###2      | 74583/228483 [00:08<00:23, 6528.77it/s]
 33%|###3      | 75529/228483 [00:08<00:21, 7190.15it/s]
 33%|###3      | 76499/228483 [00:08<00:19, 7802.72it/s]
 34%|###3      | 77462/228483 [00:08<00:18, 8273.93it/s]
 34%|###4      | 78419/228483 [00:08<00:17, 8623.15it/s]
 35%|###4      | 79387/228483 [00:08<00:16, 8916.50it/s]
 35%|###5      | 80353/228483 [00:08<00:16, 9127.42it/s]
 36%|###5      | 81330/228483 [00:08<00:15, 9311.95it/s]
 36%|###6      | 82283/228483 [00:08<00:15, 9353.19it/s]
 36%|###6      | 83258/228483 [00:09<00:15, 9469.04it/s]
 37%|###6      | 84216/228483 [00:09<00:15, 9483.04it/s]
 37%|###7      | 85192/228483 [00:09<00:14, 9562.49it/s]
 38%|###7      | 86154/228483 [00:09<00:14, 9549.70it/s]
 38%|###8      | 87113/228483 [00:09<00:14, 9540.03it/s]
 39%|###8      | 88078/228483 [00:09<00:14, 9571.65it/s]
 39%|###8      | 89038/228483 [00:09<00:14, 9573.06it/s]
 39%|###9      | 89997/228483 [00:09<00:14, 9564.54it/s]
 40%|###9      | 90955/228483 [00:09<00:14, 9543.33it/s]
 40%|####      | 91910/228483 [00:09<00:14, 9523.30it/s]
 41%|####      | 92865/228483 [00:10<00:14, 9528.64it/s]
 41%|####1     | 93826/228483 [00:10<00:14, 9551.25it/s]
 41%|####1     | 94785/228483 [00:10<00:13, 9561.90it/s]
 42%|####1     | 95742/228483 [00:10<00:13, 9554.39it/s]
 42%|####2     | 96707/228483 [00:10<00:13, 9580.87it/s]
 43%|####2     | 97674/228483 [00:10<00:13, 9605.39it/s]
 43%|####3     | 98644/228483 [00:10<00:13, 9632.96it/s]
 44%|####3     | 99608/228483 [00:10<00:13, 9617.32it/s]
 44%|####4     | 100575/228483 [00:10<00:13, 9630.58it/s]
 44%|####4     | 101539/228483 [00:10<00:13, 9595.17it/s]
 45%|####4     | 102510/228483 [00:11<00:13, 9628.70it/s]
 45%|####5     | 103473/228483 [00:11<00:13, 9610.87it/s]
 46%|####5     | 104435/228483 [00:11<00:12, 9588.68it/s]
 46%|####6     | 105394/228483 [00:11<00:12, 9588.28it/s]
 47%|####6     | 106367/228483 [00:11<00:12, 9630.45it/s]
 47%|####6     | 107338/228483 [00:11<00:12, 9651.53it/s]
 47%|####7     | 108304/228483 [00:11<00:12, 9603.33it/s]
 48%|####7     | 109265/228483 [00:11<00:12, 9547.68it/s]
 48%|####8     | 110228/228483 [00:11<00:12, 9571.28it/s]
 49%|####8     | 111188/228483 [00:11<00:12, 9579.33it/s]
 49%|####9     | 112146/228483 [00:12<00:12, 9569.93it/s]
 50%|####9     | 113104/228483 [00:12<00:12, 9522.08it/s]
 50%|####9     | 114069/228483 [00:12<00:11, 9557.67it/s]
 50%|#####     | 115044/228483 [00:12<00:11, 9612.80it/s]
 51%|#####     | 116006/228483 [00:12<00:11, 9602.15it/s]
 51%|#####1    | 116968/228483 [00:12<00:11, 9606.41it/s]
 52%|#####1    | 117929/228483 [00:12<00:11, 9602.19it/s]
 52%|#####2    | 118896/228483 [00:12<00:11, 9621.97it/s]
 52%|#####2    | 119868/228483 [00:12<00:11, 9650.57it/s]
 53%|#####2    | 120834/228483 [00:12<00:11, 9590.27it/s]
 53%|#####3    | 121795/228483 [00:13<00:11, 9595.19it/s]
 54%|#####3    | 122759/228483 [00:13<00:11, 9606.46it/s]
 54%|#####4    | 123730/228483 [00:13<00:10, 9634.52it/s]
 55%|#####4    | 124694/228483 [00:13<00:10, 9596.71it/s]
 55%|#####4    | 125654/228483 [00:13<00:10, 9597.49it/s]
 55%|#####5    | 126614/228483 [00:13<00:16, 6316.11it/s]
 56%|#####5    | 127588/228483 [00:13<00:14, 7066.84it/s]
 56%|#####6    | 128494/228483 [00:13<00:13, 7542.04it/s]
 57%|#####6    | 129444/228483 [00:14<00:12, 8038.06it/s]
 57%|#####7    | 130401/228483 [00:14<00:11, 8444.93it/s]
 57%|#####7    | 131354/228483 [00:14<00:11, 8742.41it/s]
 58%|#####7    | 132321/228483 [00:14<00:10, 9004.40it/s]
 58%|#####8    | 133281/228483 [00:14<00:10, 9173.40it/s]
 59%|#####8    | 134245/228483 [00:14<00:10, 9306.62it/s]
 59%|#####9    | 135205/228483 [00:14<00:09, 9390.72it/s]
 60%|#####9    | 136168/228483 [00:14<00:09, 9459.54it/s]
 60%|######    | 137140/228483 [00:14<00:09, 9535.74it/s]
 60%|######    | 138100/228483 [00:14<00:09, 9519.46it/s]
 61%|######    | 139067/228483 [00:15<00:09, 9563.21it/s]
 61%|######1   | 140030/228483 [00:15<00:09, 9580.22it/s]
 62%|######1   | 140998/228483 [00:15<00:09, 9608.55it/s]
 62%|######2   | 141968/228483 [00:15<00:08, 9633.04it/s]
 63%|######2   | 142933/228483 [00:15<00:08, 9602.52it/s]
 63%|######2   | 143899/228483 [00:15<00:08, 9617.48it/s]
 63%|######3   | 144862/228483 [00:15<00:08, 9612.26it/s]
 64%|######3   | 145824/228483 [00:15<00:08, 9567.81it/s]
 64%|######4   | 146782/228483 [00:15<00:08, 9542.73it/s]
 65%|######4   | 147737/228483 [00:15<00:08, 9491.98it/s]
 65%|######5   | 148691/228483 [00:16<00:08, 9506.22it/s]
 66%|######5   | 149658/228483 [00:16<00:08, 9553.52it/s]
 66%|######5   | 150614/228483 [00:16<00:08, 9514.29it/s]
 66%|######6   | 151566/228483 [00:16<00:08, 9508.38it/s]
 67%|######6   | 152532/228483 [00:16<00:07, 9553.37it/s]
 67%|######7   | 153504/228483 [00:16<00:07, 9601.44it/s]
 68%|######7   | 154470/228483 [00:16<00:07, 9618.52it/s]
 68%|######8   | 155432/228483 [00:16<00:07, 9598.66it/s]
 68%|######8   | 156392/228483 [00:16<00:07, 9598.04it/s]
 69%|######8   | 157352/228483 [00:16<00:07, 9570.66it/s]
 69%|######9   | 158319/228483 [00:17<00:07, 9597.42it/s]
 70%|######9   | 159279/228483 [00:17<00:07, 9564.81it/s]
 70%|#######   | 160236/228483 [00:17<00:07, 9562.57it/s]
 71%|#######   | 161194/228483 [00:17<00:07, 9566.43it/s]
 71%|#######   | 162151/228483 [00:17<00:06, 9550.99it/s]
 71%|#######1  | 163107/228483 [00:17<00:06, 9534.16it/s]
 72%|#######1  | 164061/228483 [00:17<00:06, 9408.95it/s]
 72%|#######2  | 165003/228483 [00:17<00:06, 9394.06it/s]
 73%|#######2  | 165950/228483 [00:17<00:06, 9416.56it/s]
 73%|#######3  | 166892/228483 [00:17<00:06, 9416.89it/s]
 73%|#######3  | 167842/228483 [00:18<00:06, 9440.23it/s]
 74%|#######3  | 168787/228483 [00:18<00:06, 9421.10it/s]
 74%|#######4  | 169737/228483 [00:18<00:06, 9442.71it/s]
 75%|#######4  | 170691/228483 [00:18<00:06, 9470.59it/s]
 75%|#######5  | 171644/228483 [00:18<00:05, 9487.77it/s]
 76%|#######5  | 172605/228483 [00:18<00:05, 9522.79it/s]
 76%|#######5  | 173558/228483 [00:18<00:05, 9402.50it/s]
 76%|#######6  | 174510/228483 [00:18<00:05, 9434.60it/s]
 77%|#######6  | 175470/228483 [00:18<00:05, 9483.03it/s]
 77%|#######7  | 176419/228483 [00:18<00:05, 9455.92it/s]
 78%|#######7  | 177378/228483 [00:19<00:05, 9495.54it/s]
 78%|#######8  | 178328/228483 [00:19<00:05, 9488.50it/s]
 78%|#######8  | 179282/228483 [00:19<00:05, 9502.01it/s]
 79%|#######8  | 180233/228483 [00:19<00:05, 9486.45it/s]
 79%|#######9  | 181188/228483 [00:19<00:04, 9504.34it/s]
 80%|#######9  | 182139/228483 [00:19<00:04, 9502.81it/s]
 80%|########  | 183090/228483 [00:19<00:04, 9498.73it/s]
 81%|########  | 184040/228483 [00:19<00:04, 9439.43it/s]
 81%|########  | 184985/228483 [00:19<00:04, 9430.40it/s]
 81%|########1 | 185929/228483 [00:19<00:04, 9398.11it/s]
 82%|########1 | 186878/228483 [00:20<00:04, 9423.71it/s]
 82%|########2 | 187821/228483 [00:20<00:04, 9424.34it/s]
 83%|########2 | 188767/228483 [00:20<00:04, 9432.57it/s]
 83%|########3 | 189721/228483 [00:20<00:04, 9464.05it/s]
 83%|########3 | 190668/228483 [00:20<00:06, 5925.10it/s]
 84%|########3 | 191466/228483 [00:20<00:05, 6358.76it/s]
 84%|########4 | 192244/228483 [00:20<00:05, 6686.98it/s]
 85%|########4 | 193182/228483 [00:20<00:04, 7355.89it/s]
 85%|########4 | 194141/228483 [00:21<00:04, 7939.58it/s]
 85%|########5 | 195081/228483 [00:21<00:04, 8337.17it/s]
 86%|########5 | 195998/228483 [00:21<00:03, 8568.51it/s]
 86%|########6 | 196940/228483 [00:21<00:03, 8810.59it/s]
 87%|########6 | 197898/228483 [00:21<00:03, 9032.87it/s]
 87%|########7 | 198845/228483 [00:21<00:03, 9157.93it/s]
 87%|########7 | 199798/228483 [00:21<00:03, 9265.21it/s]
 88%|########7 | 200735/228483 [00:21<00:02, 9288.39it/s]
 88%|########8 | 201680/228483 [00:21<00:02, 9334.16it/s]
 89%|########8 | 202620/228483 [00:21<00:02, 9351.49it/s]
 89%|########9 | 203559/228483 [00:22<00:02, 9361.98it/s]
 90%|########9 | 204502/228483 [00:22<00:02, 9378.75it/s]
 90%|########9 | 205446/228483 [00:22<00:02, 9397.01it/s]
 90%|######### | 206397/228483 [00:22<00:02, 9428.78it/s]
 91%|######### | 207343/228483 [00:22<00:02, 9435.25it/s]
 91%|#########1| 208300/228483 [00:22<00:02, 9473.99it/s]
 92%|#########1| 209248/228483 [00:22<00:02, 9468.20it/s]
 92%|#########1| 210196/228483 [00:22<00:01, 9455.43it/s]
 92%|#########2| 211146/228483 [00:22<00:01, 9468.63it/s]
 93%|#########2| 212094/228483 [00:22<00:01, 9468.19it/s]
 93%|#########3| 213041/228483 [00:23<00:01, 9463.44it/s]
 94%|#########3| 213988/228483 [00:23<00:01, 9463.61it/s]
 94%|#########4| 214941/228483 [00:23<00:01, 9482.87it/s]
 94%|#########4| 215890/228483 [00:23<00:01, 9469.38it/s]
 95%|#########4| 216844/228483 [00:23<00:01, 9489.42it/s]
 95%|#########5| 217793/228483 [00:23<00:01, 9481.98it/s]
 96%|#########5| 218743/228483 [00:23<00:01, 9484.68it/s]
 96%|#########6| 219692/228483 [00:23<00:00, 9437.01it/s]
 97%|#########6| 220636/228483 [00:23<00:00, 9424.08it/s]
 97%|#########6| 221579/228483 [00:23<00:00, 9401.70it/s]
 97%|#########7| 222520/228483 [00:24<00:00, 9402.47it/s]
 98%|#########7| 223462/228483 [00:24<00:00, 9407.10it/s]
 98%|#########8| 224411/228483 [00:24<00:00, 9431.35it/s]
 99%|#########8| 225359/228483 [00:24<00:00, 9443.56it/s]
 99%|#########9| 226315/228483 [00:24<00:00, 9476.11it/s]
 99%|#########9| 227276/228483 [00:24<00:00, 9514.73it/s]
100%|#########9| 228228/228483 [00:24<00:00, 9509.17it/s]
100%|##########| 228483/228483 [00:31<00:00, 7264.88it/s]

<nilearn.plotting.displays._slicers.ZSlicer object at 0x7f59f4281590>

Conjunction analysis

To determine the overlap of the meta-analytic results, a conjunction image can be computed by (a) identifying voxels that were statistically significant in both individual group maps and (b) selecting, for each of these voxels, the smaller of the two group-specific z values Nichols et al.2. Since this is simple arithmetic on images, conjunction is not implemented as a separate method in NiMARE but can easily be achieved with nilearn.image.math_img().

from nilearn.image import math_img

formula = "np.where(img1 * img2 > 0, np.minimum(img1, img2), 0)"
img_conj = math_img(formula, img1=knowledge_img, img2=related_img)

plot_stat_map(
    img_conj,
    cut_coords=4,
    display_mode="z",
    title="Conjunction",
    threshold=2.326,  # cluster-level p < .01, one-tailed
    cmap="RdBu_r",
    vmax=4,
)

Out:

<nilearn.plotting.displays._slicers.ZSlicer object at 0x7f59f506a0d0>

References

1(1,2)

Alexander Enge, Rasha Abdel Rahman, and Michael A Skeide. A meta-analysis of fmri studies of semantic cognition in children. NeuroImage, 241:118436, 2021. URL: https://doi.org/10.1016/j.neuroimage.2021.118436, doi:10.1016/j.neuroimage.2021.118436.

2

Thomas Nichols, Matthew Brett, Jesper Andersson, Tor Wager, and Jean-Baptiste Poline. Valid conjunction inference with the minimum statistic. Neuroimage, 25(3):653–660, 2005. URL: https://doi.org/10.1016/j.neuroimage.2004.12.005, doi:10.1016/j.neuroimage.2004.12.005.