Note
Click here to download the full example code
Subtraction and conjunction CBMAs
The (coordinate-based) ALE subtraction method tests at which voxels the meta-analytic results from two groups of studies differ reliably from one another. 1, 2
import os
import matplotlib.pyplot as plt
from nilearn.image import math_img
from nilearn.plotting import plot_stat_map
from nimare import io, utils
from nimare.correct import FWECorrector
from nimare.diagnostics import FocusCounter, Jackknife
from nimare.meta.cbma import ALE, ALESubtraction
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. 3 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.
knowledge_file = os.path.join(utils.get_resource_path(), "semantic_knowledge_children.txt")
related_file = os.path.join(utils.get_resource_path(), "semantic_relatedness_children.txt")
knowledge_dset = io.convert_sleuth_to_dataset(knowledge_file)
related_dset = io.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.
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=1)
knowledge_corrected_results = corr.transform(knowledge_results)
related_corrected_results = corr.transform(related_results)
fig, axes = plt.subplots(figsize=(12, 10), nrows=2)
img = knowledge_corrected_results.get_map("z_desc-size_level-cluster_corr-FWE_method-montecarlo")
plot_stat_map(
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,
)
img2 = related_corrected_results.get_map("z_desc-size_level-cluster_corr-FWE_method-montecarlo")
plot_stat_map(
img2,
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<00:37, 2.66it/s]
2%|2 | 2/100 [00:00<00:35, 2.74it/s]
3%|3 | 3/100 [00:01<00:34, 2.79it/s]
4%|4 | 4/100 [00:01<00:34, 2.81it/s]
5%|5 | 5/100 [00:01<00:33, 2.81it/s]
6%|6 | 6/100 [00:02<00:33, 2.82it/s]
7%|7 | 7/100 [00:02<00:32, 2.83it/s]
8%|8 | 8/100 [00:02<00:32, 2.83it/s]
9%|9 | 9/100 [00:03<00:32, 2.83it/s]
10%|# | 10/100 [00:03<00:31, 2.84it/s]
11%|#1 | 11/100 [00:03<00:31, 2.83it/s]
12%|#2 | 12/100 [00:04<00:31, 2.82it/s]
13%|#3 | 13/100 [00:04<00:30, 2.83it/s]
14%|#4 | 14/100 [00:04<00:30, 2.84it/s]
15%|#5 | 15/100 [00:05<00:29, 2.84it/s]
16%|#6 | 16/100 [00:05<00:29, 2.85it/s]
17%|#7 | 17/100 [00:06<00:29, 2.84it/s]
18%|#8 | 18/100 [00:06<00:28, 2.83it/s]
19%|#9 | 19/100 [00:06<00:28, 2.83it/s]
20%|## | 20/100 [00:07<00:28, 2.83it/s]
21%|##1 | 21/100 [00:07<00:27, 2.84it/s]
22%|##2 | 22/100 [00:07<00:27, 2.85it/s]
23%|##3 | 23/100 [00:08<00:26, 2.86it/s]
24%|##4 | 24/100 [00:08<00:26, 2.84it/s]
25%|##5 | 25/100 [00:08<00:26, 2.84it/s]
26%|##6 | 26/100 [00:09<00:26, 2.83it/s]
27%|##7 | 27/100 [00:09<00:25, 2.84it/s]
28%|##8 | 28/100 [00:09<00:25, 2.82it/s]
29%|##9 | 29/100 [00:10<00:25, 2.82it/s]
30%|### | 30/100 [00:10<00:24, 2.81it/s]
31%|###1 | 31/100 [00:10<00:24, 2.80it/s]
32%|###2 | 32/100 [00:11<00:24, 2.80it/s]
33%|###3 | 33/100 [00:11<00:23, 2.80it/s]
34%|###4 | 34/100 [00:12<00:23, 2.82it/s]
35%|###5 | 35/100 [00:12<00:23, 2.81it/s]
36%|###6 | 36/100 [00:12<00:22, 2.81it/s]
37%|###7 | 37/100 [00:13<00:22, 2.83it/s]
38%|###8 | 38/100 [00:13<00:21, 2.84it/s]
39%|###9 | 39/100 [00:13<00:21, 2.83it/s]
40%|#### | 40/100 [00:14<00:21, 2.83it/s]
41%|####1 | 41/100 [00:14<00:20, 2.84it/s]
42%|####2 | 42/100 [00:14<00:20, 2.83it/s]
43%|####3 | 43/100 [00:15<00:20, 2.82it/s]
44%|####4 | 44/100 [00:15<00:19, 2.82it/s]
45%|####5 | 45/100 [00:15<00:19, 2.81it/s]
46%|####6 | 46/100 [00:16<00:19, 2.81it/s]
47%|####6 | 47/100 [00:16<00:18, 2.81it/s]
48%|####8 | 48/100 [00:16<00:18, 2.81it/s]
49%|####9 | 49/100 [00:17<00:18, 2.80it/s]
50%|##### | 50/100 [00:17<00:17, 2.80it/s]
51%|#####1 | 51/100 [00:18<00:17, 2.80it/s]
52%|#####2 | 52/100 [00:18<00:17, 2.79it/s]
53%|#####3 | 53/100 [00:18<00:16, 2.79it/s]
54%|#####4 | 54/100 [00:19<00:16, 2.80it/s]
55%|#####5 | 55/100 [00:19<00:16, 2.81it/s]
56%|#####6 | 56/100 [00:19<00:15, 2.81it/s]
57%|#####6 | 57/100 [00:20<00:15, 2.81it/s]
58%|#####8 | 58/100 [00:20<00:14, 2.81it/s]
59%|#####8 | 59/100 [00:20<00:14, 2.80it/s]
60%|###### | 60/100 [00:21<00:14, 2.79it/s]
61%|######1 | 61/100 [00:21<00:13, 2.79it/s]
62%|######2 | 62/100 [00:22<00:13, 2.79it/s]
63%|######3 | 63/100 [00:22<00:13, 2.80it/s]
64%|######4 | 64/100 [00:22<00:12, 2.82it/s]
65%|######5 | 65/100 [00:23<00:12, 2.83it/s]
66%|######6 | 66/100 [00:23<00:11, 2.84it/s]
67%|######7 | 67/100 [00:23<00:11, 2.81it/s]
68%|######8 | 68/100 [00:24<00:11, 2.83it/s]
69%|######9 | 69/100 [00:24<00:10, 2.83it/s]
70%|####### | 70/100 [00:24<00:10, 2.83it/s]
71%|#######1 | 71/100 [00:25<00:10, 2.83it/s]
72%|#######2 | 72/100 [00:25<00:09, 2.84it/s]
73%|#######3 | 73/100 [00:25<00:09, 2.83it/s]
74%|#######4 | 74/100 [00:26<00:09, 2.83it/s]
75%|#######5 | 75/100 [00:26<00:08, 2.83it/s]
76%|#######6 | 76/100 [00:26<00:08, 2.82it/s]
77%|#######7 | 77/100 [00:27<00:08, 2.83it/s]
78%|#######8 | 78/100 [00:27<00:07, 2.82it/s]
79%|#######9 | 79/100 [00:28<00:07, 2.81it/s]
80%|######## | 80/100 [00:28<00:07, 2.81it/s]
81%|########1 | 81/100 [00:28<00:06, 2.81it/s]
82%|########2 | 82/100 [00:29<00:06, 2.82it/s]
83%|########2 | 83/100 [00:29<00:06, 2.83it/s]
84%|########4 | 84/100 [00:29<00:05, 2.82it/s]
85%|########5 | 85/100 [00:30<00:05, 2.82it/s]
86%|########6 | 86/100 [00:30<00:04, 2.82it/s]
87%|########7 | 87/100 [00:30<00:04, 2.81it/s]
88%|########8 | 88/100 [00:31<00:04, 2.82it/s]
89%|########9 | 89/100 [00:31<00:03, 2.83it/s]
90%|######### | 90/100 [00:31<00:03, 2.82it/s]
91%|#########1| 91/100 [00:32<00:03, 2.83it/s]
92%|#########2| 92/100 [00:32<00:02, 2.83it/s]
93%|#########3| 93/100 [00:32<00:02, 2.83it/s]
94%|#########3| 94/100 [00:33<00:02, 2.83it/s]
95%|#########5| 95/100 [00:33<00:01, 2.82it/s]
96%|#########6| 96/100 [00:34<00:01, 2.82it/s]
97%|#########7| 97/100 [00:34<00:01, 2.81it/s]
98%|#########8| 98/100 [00:34<00:00, 2.81it/s]
99%|#########9| 99/100 [00:35<00:00, 2.81it/s]
100%|##########| 100/100 [00:35<00:00, 2.82it/s]
100%|##########| 100/100 [00:35<00:00, 2.82it/s]
0%| | 0/100 [00:00<?, ?it/s]
1%|1 | 1/100 [00:00<00:33, 2.94it/s]
2%|2 | 2/100 [00:00<00:32, 3.03it/s]
3%|3 | 3/100 [00:00<00:31, 3.07it/s]
4%|4 | 4/100 [00:01<00:31, 3.08it/s]
5%|5 | 5/100 [00:01<00:30, 3.09it/s]
6%|6 | 6/100 [00:01<00:30, 3.10it/s]
7%|7 | 7/100 [00:02<00:29, 3.10it/s]
8%|8 | 8/100 [00:02<00:29, 3.10it/s]
9%|9 | 9/100 [00:02<00:29, 3.10it/s]
10%|# | 10/100 [00:03<00:28, 3.11it/s]
11%|#1 | 11/100 [00:03<00:28, 3.10it/s]
12%|#2 | 12/100 [00:03<00:28, 3.08it/s]
13%|#3 | 13/100 [00:04<00:28, 3.08it/s]
14%|#4 | 14/100 [00:04<00:27, 3.08it/s]
15%|#5 | 15/100 [00:04<00:27, 3.08it/s]
16%|#6 | 16/100 [00:05<00:27, 3.09it/s]
17%|#7 | 17/100 [00:05<00:26, 3.10it/s]
18%|#8 | 18/100 [00:05<00:26, 3.11it/s]
19%|#9 | 19/100 [00:06<00:25, 3.12it/s]
20%|## | 20/100 [00:06<00:25, 3.12it/s]
21%|##1 | 21/100 [00:06<00:25, 3.13it/s]
22%|##2 | 22/100 [00:07<00:24, 3.14it/s]
23%|##3 | 23/100 [00:07<00:24, 3.14it/s]
24%|##4 | 24/100 [00:07<00:24, 3.14it/s]
25%|##5 | 25/100 [00:08<00:23, 3.14it/s]
26%|##6 | 26/100 [00:08<00:23, 3.14it/s]
27%|##7 | 27/100 [00:08<00:23, 3.13it/s]
28%|##8 | 28/100 [00:09<00:22, 3.14it/s]
29%|##9 | 29/100 [00:09<00:22, 3.14it/s]
30%|### | 30/100 [00:09<00:22, 3.11it/s]
31%|###1 | 31/100 [00:09<00:22, 3.10it/s]
32%|###2 | 32/100 [00:10<00:22, 3.08it/s]
33%|###3 | 33/100 [00:10<00:21, 3.08it/s]
34%|###4 | 34/100 [00:10<00:21, 3.09it/s]
35%|###5 | 35/100 [00:11<00:20, 3.10it/s]
36%|###6 | 36/100 [00:11<00:20, 3.11it/s]
37%|###7 | 37/100 [00:11<00:20, 3.10it/s]
38%|###8 | 38/100 [00:12<00:20, 3.10it/s]
39%|###9 | 39/100 [00:12<00:19, 3.08it/s]
40%|#### | 40/100 [00:12<00:19, 3.09it/s]
41%|####1 | 41/100 [00:13<00:19, 3.09it/s]
42%|####2 | 42/100 [00:13<00:18, 3.09it/s]
43%|####3 | 43/100 [00:13<00:18, 3.11it/s]
44%|####4 | 44/100 [00:14<00:17, 3.12it/s]
45%|####5 | 45/100 [00:14<00:17, 3.13it/s]
46%|####6 | 46/100 [00:14<00:17, 3.13it/s]
47%|####6 | 47/100 [00:15<00:17, 3.12it/s]
48%|####8 | 48/100 [00:15<00:16, 3.12it/s]
49%|####9 | 49/100 [00:15<00:16, 3.10it/s]
50%|##### | 50/100 [00:16<00:16, 3.11it/s]
51%|#####1 | 51/100 [00:16<00:15, 3.12it/s]
52%|#####2 | 52/100 [00:16<00:15, 3.15it/s]
53%|#####3 | 53/100 [00:17<00:15, 3.13it/s]
54%|#####4 | 54/100 [00:17<00:14, 3.14it/s]
55%|#####5 | 55/100 [00:17<00:14, 3.13it/s]
56%|#####6 | 56/100 [00:18<00:14, 3.14it/s]
57%|#####6 | 57/100 [00:18<00:13, 3.15it/s]
58%|#####8 | 58/100 [00:18<00:13, 3.13it/s]
59%|#####8 | 59/100 [00:18<00:13, 3.14it/s]
60%|###### | 60/100 [00:19<00:12, 3.13it/s]
61%|######1 | 61/100 [00:19<00:12, 3.13it/s]
62%|######2 | 62/100 [00:19<00:12, 3.12it/s]
63%|######3 | 63/100 [00:20<00:11, 3.11it/s]
64%|######4 | 64/100 [00:20<00:11, 3.11it/s]
65%|######5 | 65/100 [00:20<00:11, 3.10it/s]
66%|######6 | 66/100 [00:21<00:10, 3.12it/s]
67%|######7 | 67/100 [00:21<00:10, 3.14it/s]
68%|######8 | 68/100 [00:21<00:10, 3.14it/s]
69%|######9 | 69/100 [00:22<00:09, 3.14it/s]
70%|####### | 70/100 [00:22<00:09, 3.16it/s]
71%|#######1 | 71/100 [00:22<00:09, 3.18it/s]
72%|#######2 | 72/100 [00:23<00:08, 3.18it/s]
73%|#######3 | 73/100 [00:23<00:08, 3.18it/s]
74%|#######4 | 74/100 [00:23<00:08, 3.18it/s]
75%|#######5 | 75/100 [00:24<00:07, 3.19it/s]
76%|#######6 | 76/100 [00:24<00:07, 3.18it/s]
77%|#######7 | 77/100 [00:24<00:07, 3.16it/s]
78%|#######8 | 78/100 [00:24<00:06, 3.17it/s]
79%|#######9 | 79/100 [00:25<00:06, 3.17it/s]
80%|######## | 80/100 [00:25<00:06, 3.16it/s]
81%|########1 | 81/100 [00:25<00:06, 3.16it/s]
82%|########2 | 82/100 [00:26<00:05, 3.16it/s]
83%|########2 | 83/100 [00:26<00:05, 3.16it/s]
84%|########4 | 84/100 [00:26<00:05, 3.17it/s]
85%|########5 | 85/100 [00:27<00:04, 3.17it/s]
86%|########6 | 86/100 [00:27<00:04, 3.16it/s]
87%|########7 | 87/100 [00:27<00:04, 3.16it/s]
88%|########8 | 88/100 [00:28<00:03, 3.16it/s]
89%|########9 | 89/100 [00:28<00:03, 3.17it/s]
90%|######### | 90/100 [00:28<00:03, 3.16it/s]
91%|#########1| 91/100 [00:29<00:02, 3.18it/s]
92%|#########2| 92/100 [00:29<00:02, 3.18it/s]
93%|#########3| 93/100 [00:29<00:02, 3.16it/s]
94%|#########3| 94/100 [00:30<00:01, 3.16it/s]
95%|#########5| 95/100 [00:30<00:01, 3.16it/s]
96%|#########6| 96/100 [00:30<00:01, 3.16it/s]
97%|#########7| 97/100 [00:30<00:00, 3.17it/s]
98%|#########8| 98/100 [00:31<00:00, 3.17it/s]
99%|#########9| 99/100 [00:31<00:00, 3.15it/s]
100%|##########| 100/100 [00:31<00:00, 3.16it/s]
100%|##########| 100/100 [00:31<00:00, 3.13it/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.
counter = FocusCounter(
target_image="z_desc-size_level-cluster_corr-FWE_method-montecarlo",
voxel_thresh=None,
)
knowledge_count_table, knowledge_cluster_img = counter.transform(knowledge_corrected_results)
related_count_table, related_cluster_img = counter.transform(related_corrected_results)
Out:
0%| | 0/21 [00:00<?, ?it/s]
10%|9 | 2/21 [00:00<00:01, 18.45it/s]
19%|#9 | 4/21 [00:00<00:00, 18.26it/s]
29%|##8 | 6/21 [00:00<00:00, 18.22it/s]
38%|###8 | 8/21 [00:00<00:00, 18.17it/s]
48%|####7 | 10/21 [00:00<00:00, 18.19it/s]
57%|#####7 | 12/21 [00:00<00:00, 18.16it/s]
67%|######6 | 14/21 [00:00<00:00, 18.19it/s]
76%|#######6 | 16/21 [00:00<00:00, 18.22it/s]
86%|########5 | 18/21 [00:00<00:00, 18.29it/s]
95%|#########5| 20/21 [00:01<00:00, 18.38it/s]
100%|##########| 21/21 [00:01<00:00, 18.29it/s]
0%| | 0/16 [00:00<?, ?it/s]
19%|#8 | 3/16 [00:00<00:00, 22.28it/s]
38%|###7 | 6/16 [00:00<00:00, 22.14it/s]
56%|#####6 | 9/16 [00:00<00:00, 22.13it/s]
75%|#######5 | 12/16 [00:00<00:00, 22.16it/s]
94%|#########3| 15/16 [00:00<00:00, 22.00it/s]
100%|##########| 16/16 [00:00<00:00, 22.05it/s]
jackknife = Jackknife(
target_image="z_desc-size_level-cluster_corr-FWE_method-montecarlo",
voxel_thresh=None,
)
knowledge_jackknife_table, _ = jackknife.transform(knowledge_corrected_results)
related_jackknife_table, _ = jackknife.transform(related_corrected_results)
Out:
0%| | 0/21 [00:00<?, ?it/s]
5%|4 | 1/21 [00:01<00:29, 1.49s/it]
10%|9 | 2/21 [00:02<00:28, 1.48s/it]
14%|#4 | 3/21 [00:04<00:26, 1.46s/it]
19%|#9 | 4/21 [00:05<00:24, 1.45s/it]
24%|##3 | 5/21 [00:07<00:23, 1.44s/it]
29%|##8 | 6/21 [00:08<00:21, 1.45s/it]
33%|###3 | 7/21 [00:10<00:20, 1.46s/it]
38%|###8 | 8/21 [00:11<00:18, 1.45s/it]
43%|####2 | 9/21 [00:13<00:17, 1.48s/it]
48%|####7 | 10/21 [00:14<00:16, 1.47s/it]
52%|#####2 | 11/21 [00:16<00:14, 1.47s/it]
57%|#####7 | 12/21 [00:17<00:13, 1.47s/it]
62%|######1 | 13/21 [00:19<00:11, 1.46s/it]
67%|######6 | 14/21 [00:20<00:10, 1.46s/it]
71%|#######1 | 15/21 [00:21<00:08, 1.46s/it]
76%|#######6 | 16/21 [00:23<00:07, 1.46s/it]
81%|######## | 17/21 [00:24<00:05, 1.47s/it]
86%|########5 | 18/21 [00:26<00:04, 1.48s/it]
90%|######### | 19/21 [00:27<00:02, 1.47s/it]
95%|#########5| 20/21 [00:29<00:01, 1.47s/it]
100%|##########| 21/21 [00:30<00:00, 1.47s/it]
100%|##########| 21/21 [00:30<00:00, 1.47s/it]
0%| | 0/16 [00:00<?, ?it/s]
6%|6 | 1/16 [00:01<00:18, 1.24s/it]
12%|#2 | 2/16 [00:02<00:17, 1.22s/it]
19%|#8 | 3/16 [00:03<00:15, 1.21s/it]
25%|##5 | 4/16 [00:04<00:14, 1.22s/it]
31%|###1 | 5/16 [00:06<00:13, 1.23s/it]
38%|###7 | 6/16 [00:07<00:12, 1.23s/it]
44%|####3 | 7/16 [00:08<00:11, 1.23s/it]
50%|##### | 8/16 [00:09<00:09, 1.23s/it]
56%|#####6 | 9/16 [00:11<00:08, 1.23s/it]
62%|######2 | 10/16 [00:12<00:07, 1.23s/it]
69%|######8 | 11/16 [00:13<00:06, 1.22s/it]
75%|#######5 | 12/16 [00:14<00:04, 1.22s/it]
81%|########1 | 13/16 [00:15<00:03, 1.22s/it]
88%|########7 | 14/16 [00:17<00:02, 1.22s/it]
94%|#########3| 15/16 [00:18<00:01, 1.22s/it]
100%|##########| 16/16 [00:19<00:00, 1.23s/it]
100%|##########| 16/16 [00:19<00:00, 1.23s/it]
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.
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]
2%|2 | 2/100 [00:00<00:05, 19.25it/s]
5%|5 | 5/100 [00:00<00:04, 21.45it/s]
8%|8 | 8/100 [00:00<00:04, 22.01it/s]
11%|#1 | 11/100 [00:00<00:03, 22.26it/s]
14%|#4 | 14/100 [00:00<00:03, 22.41it/s]
17%|#7 | 17/100 [00:00<00:03, 22.61it/s]
20%|## | 20/100 [00:00<00:03, 22.65it/s]
23%|##3 | 23/100 [00:01<00:03, 22.68it/s]
26%|##6 | 26/100 [00:01<00:03, 22.78it/s]
29%|##9 | 29/100 [00:01<00:03, 22.84it/s]
32%|###2 | 32/100 [00:01<00:02, 22.69it/s]
35%|###5 | 35/100 [00:01<00:02, 22.69it/s]
38%|###8 | 38/100 [00:01<00:02, 22.78it/s]
41%|####1 | 41/100 [00:01<00:02, 22.92it/s]
44%|####4 | 44/100 [00:01<00:02, 22.90it/s]
47%|####6 | 47/100 [00:02<00:02, 22.99it/s]
50%|##### | 50/100 [00:02<00:02, 23.06it/s]
53%|#####3 | 53/100 [00:02<00:02, 23.10it/s]
56%|#####6 | 56/100 [00:02<00:01, 23.06it/s]
59%|#####8 | 59/100 [00:02<00:01, 22.94it/s]
62%|######2 | 62/100 [00:02<00:01, 23.04it/s]
65%|######5 | 65/100 [00:02<00:01, 23.08it/s]
68%|######8 | 68/100 [00:02<00:01, 23.04it/s]
71%|#######1 | 71/100 [00:03<00:01, 22.92it/s]
74%|#######4 | 74/100 [00:03<00:01, 22.90it/s]
77%|#######7 | 77/100 [00:03<00:01, 22.87it/s]
80%|######## | 80/100 [00:03<00:00, 22.82it/s]
83%|########2 | 83/100 [00:03<00:00, 22.78it/s]
86%|########6 | 86/100 [00:03<00:00, 22.84it/s]
89%|########9 | 89/100 [00:03<00:00, 22.78it/s]
92%|#########2| 92/100 [00:04<00:00, 22.86it/s]
95%|#########5| 95/100 [00:04<00:00, 22.79it/s]
98%|#########8| 98/100 [00:04<00:00, 22.74it/s]
100%|##########| 100/100 [00:04<00:00, 22.76it/s]
0%| | 0/228483 [00:00<?, ?it/s]
0%| | 945/228483 [00:00<00:24, 9447.04it/s]
1%| | 1890/228483 [00:00<00:24, 9344.01it/s]
1%|1 | 2837/228483 [00:00<00:24, 9400.95it/s]
2%|1 | 3795/228483 [00:00<00:23, 9468.89it/s]
2%|2 | 4746/228483 [00:00<00:23, 9481.31it/s]
2%|2 | 5695/228483 [00:00<00:23, 9292.70it/s]
3%|2 | 6646/228483 [00:00<00:23, 9360.27it/s]
3%|3 | 7583/228483 [00:00<00:23, 9303.03it/s]
4%|3 | 8515/228483 [00:00<00:23, 9308.01it/s]
4%|4 | 9452/228483 [00:01<00:23, 9324.64it/s]
5%|4 | 10385/228483 [00:01<00:23, 9294.58it/s]
5%|4 | 11315/228483 [00:01<00:23, 9283.70it/s]
5%|5 | 12247/228483 [00:01<00:23, 9292.80it/s]
6%|5 | 13195/228483 [00:01<00:23, 9348.66it/s]
6%|6 | 14130/228483 [00:01<00:22, 9341.54it/s]
7%|6 | 15082/228483 [00:01<00:22, 9393.62it/s]
7%|7 | 16036/228483 [00:01<00:22, 9435.85it/s]
7%|7 | 16980/228483 [00:01<00:22, 9414.10it/s]
8%|7 | 17925/228483 [00:01<00:22, 9423.82it/s]
8%|8 | 18868/228483 [00:02<00:22, 9387.12it/s]
9%|8 | 19815/228483 [00:02<00:22, 9410.74it/s]
9%|9 | 20757/228483 [00:02<00:22, 9363.82it/s]
9%|9 | 21697/228483 [00:02<00:22, 9372.88it/s]
10%|9 | 22640/228483 [00:02<00:21, 9387.56it/s]
10%|# | 23579/228483 [00:02<00:27, 7319.34it/s]
11%|# | 24535/228483 [00:02<00:25, 7877.43it/s]
11%|#1 | 25478/228483 [00:02<00:24, 8285.73it/s]
12%|#1 | 26412/228483 [00:02<00:23, 8571.93it/s]
12%|#1 | 27345/228483 [00:03<00:22, 8781.52it/s]
12%|#2 | 28297/228483 [00:03<00:22, 8992.46it/s]
13%|#2 | 29238/228483 [00:03<00:21, 9111.19it/s]
13%|#3 | 30188/228483 [00:03<00:21, 9222.49it/s]
14%|#3 | 31121/228483 [00:03<00:23, 8474.91it/s]
14%|#4 | 31988/228483 [00:03<00:24, 8179.44it/s]
14%|#4 | 32940/228483 [00:03<00:22, 8548.59it/s]
15%|#4 | 33903/228483 [00:03<00:21, 8852.39it/s]
15%|#5 | 34841/228483 [00:03<00:21, 9003.08it/s]
16%|#5 | 35778/228483 [00:03<00:21, 9107.47it/s]
16%|#6 | 36722/228483 [00:04<00:20, 9204.00it/s]
16%|#6 | 37652/228483 [00:04<00:20, 9232.31it/s]
17%|#6 | 38589/228483 [00:04<00:20, 9272.31it/s]
17%|#7 | 39547/228483 [00:04<00:20, 9361.87it/s]
18%|#7 | 40507/228483 [00:04<00:19, 9431.95it/s]
18%|#8 | 41460/228483 [00:04<00:19, 9458.59it/s]
19%|#8 | 42416/228483 [00:04<00:19, 9487.96it/s]
19%|#8 | 43378/228483 [00:04<00:19, 9526.49it/s]
19%|#9 | 44332/228483 [00:04<00:19, 9470.87it/s]
20%|#9 | 45286/228483 [00:04<00:19, 9489.32it/s]
20%|## | 46236/228483 [00:05<00:19, 9468.76it/s]
21%|## | 47184/228483 [00:05<00:19, 9419.71it/s]
21%|##1 | 48137/228483 [00:05<00:19, 9450.67it/s]
21%|##1 | 49098/228483 [00:05<00:18, 9495.64it/s]
22%|##1 | 50048/228483 [00:05<00:18, 9490.08it/s]
22%|##2 | 51002/228483 [00:05<00:18, 9503.06it/s]
23%|##2 | 51953/228483 [00:05<00:24, 7256.93it/s]
23%|##3 | 52902/228483 [00:05<00:22, 7805.76it/s]
24%|##3 | 53853/228483 [00:05<00:21, 8246.82it/s]
24%|##3 | 54806/228483 [00:06<00:20, 8593.71it/s]
24%|##4 | 55729/228483 [00:06<00:19, 8769.76it/s]
25%|##4 | 56671/228483 [00:06<00:19, 8954.36it/s]
25%|##5 | 57629/228483 [00:06<00:18, 9133.46it/s]
26%|##5 | 58562/228483 [00:06<00:18, 9188.20it/s]
26%|##6 | 59523/228483 [00:06<00:18, 9311.74it/s]
26%|##6 | 60477/228483 [00:06<00:17, 9379.15it/s]
27%|##6 | 61438/228483 [00:06<00:17, 9447.13it/s]
27%|##7 | 62387/228483 [00:06<00:17, 9428.68it/s]
28%|##7 | 63351/228483 [00:06<00:17, 9490.75it/s]
28%|##8 | 64303/228483 [00:07<00:17, 9481.33it/s]
29%|##8 | 65253/228483 [00:07<00:17, 9436.84it/s]
29%|##8 | 66222/228483 [00:07<00:17, 9511.13it/s]
29%|##9 | 67174/228483 [00:07<00:16, 9496.97it/s]
30%|##9 | 68152/228483 [00:07<00:16, 9581.38it/s]
30%|### | 69111/228483 [00:07<00:16, 9517.91it/s]
31%|### | 70072/228483 [00:07<00:16, 9544.92it/s]
31%|###1 | 71027/228483 [00:07<00:16, 9521.31it/s]
32%|###1 | 71988/228483 [00:07<00:16, 9546.33it/s]
32%|###1 | 72943/228483 [00:07<00:16, 9545.45it/s]
32%|###2 | 73898/228483 [00:08<00:16, 9530.50it/s]
33%|###2 | 74852/228483 [00:08<00:16, 9497.39it/s]
33%|###3 | 75812/228483 [00:08<00:16, 9525.42it/s]
34%|###3 | 76786/228483 [00:08<00:15, 9587.48it/s]
34%|###4 | 77748/228483 [00:08<00:15, 9596.52it/s]
34%|###4 | 78708/228483 [00:08<00:15, 9587.26it/s]
35%|###4 | 79667/228483 [00:08<00:15, 9549.04it/s]
35%|###5 | 80622/228483 [00:08<00:16, 8845.76it/s]
36%|###5 | 81517/228483 [00:08<00:18, 7950.74it/s]
36%|###6 | 82356/228483 [00:09<00:18, 8066.23it/s]
36%|###6 | 83312/228483 [00:09<00:17, 8476.74it/s]
37%|###6 | 84269/228483 [00:09<00:16, 8784.04it/s]
37%|###7 | 85228/228483 [00:09<00:15, 9013.04it/s]
38%|###7 | 86192/228483 [00:09<00:15, 9193.03it/s]
38%|###8 | 87133/228483 [00:09<00:15, 9255.22it/s]
39%|###8 | 88064/228483 [00:09<00:20, 6870.22it/s]
39%|###8 | 89012/228483 [00:09<00:18, 7491.37it/s]
39%|###9 | 89965/228483 [00:09<00:17, 8007.13it/s]
40%|###9 | 90915/228483 [00:10<00:16, 8402.81it/s]
40%|#### | 91855/228483 [00:10<00:15, 8675.81it/s]
41%|#### | 92801/228483 [00:10<00:15, 8895.25it/s]
41%|####1 | 93762/228483 [00:10<00:14, 9099.99it/s]
41%|####1 | 94710/228483 [00:10<00:14, 9208.52it/s]
42%|####1 | 95654/228483 [00:10<00:14, 9273.49it/s]
42%|####2 | 96618/228483 [00:10<00:14, 9379.91it/s]
43%|####2 | 97576/228483 [00:10<00:13, 9437.11it/s]
43%|####3 | 98526/228483 [00:10<00:13, 9348.97it/s]
44%|####3 | 99473/228483 [00:10<00:13, 9382.58it/s]
44%|####3 | 100415/228483 [00:11<00:13, 9363.75it/s]
44%|####4 | 101369/228483 [00:11<00:13, 9414.89it/s]
45%|####4 | 102322/228483 [00:11<00:13, 9446.73it/s]
45%|####5 | 103276/228483 [00:11<00:13, 9471.92it/s]
46%|####5 | 104224/228483 [00:11<00:13, 9464.22it/s]
46%|####6 | 105175/228483 [00:11<00:13, 9476.33it/s]
46%|####6 | 106131/228483 [00:11<00:12, 9500.51it/s]
47%|####6 | 107107/228483 [00:11<00:12, 9575.64it/s]
47%|####7 | 108065/228483 [00:11<00:12, 9553.61it/s]
48%|####7 | 109021/228483 [00:11<00:12, 9464.02it/s]
48%|####8 | 109978/228483 [00:12<00:12, 9492.78it/s]
49%|####8 | 110928/228483 [00:12<00:12, 9486.08it/s]
49%|####8 | 111886/228483 [00:12<00:12, 9512.98it/s]
49%|####9 | 112838/228483 [00:12<00:12, 9469.35it/s]
50%|####9 | 113802/228483 [00:12<00:12, 9518.19it/s]
50%|##### | 114757/228483 [00:12<00:11, 9526.19it/s]
51%|##### | 115729/228483 [00:12<00:11, 9581.93it/s]
51%|#####1 | 116688/228483 [00:12<00:11, 9538.02it/s]
51%|#####1 | 117642/228483 [00:12<00:11, 9459.35it/s]
52%|#####1 | 118597/228483 [00:12<00:11, 9485.96it/s]
52%|#####2 | 119562/228483 [00:13<00:11, 9533.44it/s]
53%|#####2 | 120516/228483 [00:13<00:11, 9510.31it/s]
53%|#####3 | 121479/228483 [00:13<00:11, 9543.20it/s]
54%|#####3 | 122437/228483 [00:13<00:11, 9552.80it/s]
54%|#####4 | 123399/228483 [00:13<00:10, 9571.56it/s]
54%|#####4 | 124357/228483 [00:13<00:10, 9510.89it/s]
55%|#####4 | 125309/228483 [00:13<00:10, 9489.10it/s]
55%|#####5 | 126273/228483 [00:13<00:10, 9533.64it/s]
56%|#####5 | 127227/228483 [00:13<00:10, 9519.08it/s]
56%|#####6 | 128182/228483 [00:13<00:10, 9527.02it/s]
57%|#####6 | 129135/228483 [00:14<00:10, 9520.24it/s]
57%|#####6 | 130088/228483 [00:14<00:10, 9500.40it/s]
57%|#####7 | 131039/228483 [00:14<00:10, 9501.09it/s]
58%|#####7 | 131990/228483 [00:14<00:14, 6688.96it/s]
58%|#####8 | 132947/228483 [00:14<00:12, 7353.85it/s]
59%|#####8 | 133898/228483 [00:14<00:11, 7888.52it/s]
59%|#####9 | 134860/228483 [00:14<00:11, 8341.60it/s]
59%|#####9 | 135806/228483 [00:14<00:10, 8645.52it/s]
60%|#####9 | 136765/228483 [00:15<00:10, 8908.85it/s]
60%|###### | 137706/228483 [00:15<00:10, 9049.81it/s]
61%|###### | 138666/228483 [00:15<00:09, 9206.83it/s]
61%|######1 | 139618/228483 [00:15<00:09, 9296.80it/s]
62%|######1 | 140590/228483 [00:15<00:09, 9421.17it/s]
62%|######1 | 141542/228483 [00:15<00:09, 9425.29it/s]
62%|######2 | 142492/228483 [00:15<00:09, 9432.97it/s]
63%|######2 | 143453/228483 [00:15<00:08, 9483.53it/s]
63%|######3 | 144405/228483 [00:15<00:08, 9479.22it/s]
64%|######3 | 145356/228483 [00:15<00:08, 9467.64it/s]
64%|######4 | 146305/228483 [00:16<00:08, 9460.63it/s]
64%|######4 | 147253/228483 [00:16<00:08, 9460.68it/s]
65%|######4 | 148200/228483 [00:16<00:08, 9421.96it/s]
65%|######5 | 149154/228483 [00:16<00:08, 9455.44it/s]
66%|######5 | 150119/228483 [00:16<00:08, 9511.19it/s]
66%|######6 | 151080/228483 [00:16<00:08, 9537.53it/s]
67%|######6 | 152034/228483 [00:16<00:08, 9462.18it/s]
67%|######6 | 152981/228483 [00:16<00:07, 9451.87it/s]
67%|######7 | 153943/228483 [00:16<00:07, 9500.70it/s]
68%|######7 | 154894/228483 [00:16<00:07, 9496.29it/s]
68%|######8 | 155844/228483 [00:17<00:07, 9452.52it/s]
69%|######8 | 156790/228483 [00:17<00:07, 9453.98it/s]
69%|######9 | 157736/228483 [00:17<00:07, 9445.27it/s]
69%|######9 | 158681/228483 [00:17<00:07, 9428.54it/s]
70%|######9 | 159635/228483 [00:17<00:07, 9461.65it/s]
70%|####### | 160582/228483 [00:17<00:07, 9292.98it/s]
71%|####### | 161523/228483 [00:17<00:07, 9325.43it/s]
71%|#######1 | 162474/228483 [00:17<00:07, 9378.79it/s]
72%|#######1 | 163420/228483 [00:17<00:06, 9402.37it/s]
72%|#######1 | 164364/228483 [00:17<00:06, 9412.98it/s]
72%|#######2 | 165309/228483 [00:18<00:06, 9421.60it/s]
73%|#######2 | 166254/228483 [00:18<00:06, 9429.58it/s]
73%|#######3 | 167199/228483 [00:18<00:06, 9435.44it/s]
74%|#######3 | 168148/228483 [00:18<00:06, 9448.89it/s]
74%|#######4 | 169099/228483 [00:18<00:06, 9464.76it/s]
74%|#######4 | 170057/228483 [00:18<00:06, 9499.16it/s]
75%|#######4 | 171007/228483 [00:18<00:06, 9480.72it/s]
75%|#######5 | 171956/228483 [00:18<00:05, 9455.53it/s]
76%|#######5 | 172906/228483 [00:18<00:05, 9467.88it/s]
76%|#######6 | 173853/228483 [00:18<00:05, 9436.63it/s]
77%|#######6 | 174797/228483 [00:19<00:05, 9420.74it/s]
77%|#######6 | 175743/228483 [00:19<00:05, 9430.59it/s]
77%|#######7 | 176687/228483 [00:19<00:05, 9397.88it/s]
78%|#######7 | 177635/228483 [00:19<00:05, 9419.49it/s]
78%|#######8 | 178590/228483 [00:19<00:05, 9456.29it/s]
79%|#######8 | 179536/228483 [00:19<00:05, 9416.26it/s]
79%|#######8 | 180484/228483 [00:19<00:05, 9435.00it/s]
79%|#######9 | 181434/228483 [00:19<00:04, 9452.62it/s]
80%|#######9 | 182382/228483 [00:19<00:04, 9458.60it/s]
80%|######## | 183328/228483 [00:19<00:04, 9456.99it/s]
81%|######## | 184275/228483 [00:20<00:04, 9459.51it/s]
81%|########1 | 185221/228483 [00:20<00:04, 9386.01it/s]
81%|########1 | 186160/228483 [00:20<00:04, 9368.25it/s]
82%|########1 | 187105/228483 [00:20<00:04, 9392.37it/s]
82%|########2 | 188045/228483 [00:20<00:06, 6399.43it/s]
83%|########2 | 188991/228483 [00:20<00:05, 7087.78it/s]
83%|########3 | 189942/228483 [00:20<00:05, 7675.66it/s]
84%|########3 | 190879/228483 [00:20<00:04, 8112.24it/s]
84%|########3 | 191805/228483 [00:21<00:04, 8420.73it/s]
84%|########4 | 192748/228483 [00:21<00:04, 8699.88it/s]
85%|########4 | 193692/228483 [00:21<00:03, 8909.72it/s]
85%|########5 | 194645/228483 [00:21<00:03, 9087.25it/s]
86%|########5 | 195591/228483 [00:21<00:03, 9191.38it/s]
86%|########6 | 196526/228483 [00:21<00:03, 9214.86it/s]
86%|########6 | 197458/228483 [00:21<00:03, 9236.80it/s]
87%|########6 | 198407/228483 [00:21<00:03, 9310.39it/s]
87%|########7 | 199353/228483 [00:21<00:03, 9353.50it/s]
88%|########7 | 200304/228483 [00:21<00:02, 9399.59it/s]
88%|########8 | 201247/228483 [00:22<00:02, 9385.86it/s]
88%|########8 | 202188/228483 [00:22<00:02, 9363.18it/s]
89%|########8 | 203126/228483 [00:22<00:02, 9313.66it/s]
89%|########9 | 204059/228483 [00:22<00:02, 9309.16it/s]
90%|########9 | 205005/228483 [00:22<00:02, 9353.52it/s]
90%|######### | 205952/228483 [00:22<00:02, 9387.67it/s]
91%|######### | 206892/228483 [00:22<00:02, 9371.94it/s]
91%|######### | 207847/228483 [00:22<00:02, 9423.21it/s]
91%|#########1| 208793/228483 [00:22<00:02, 9432.81it/s]
92%|#########1| 209737/228483 [00:22<00:01, 9419.67it/s]
92%|#########2| 210680/228483 [00:23<00:01, 9376.94it/s]
93%|#########2| 211627/228483 [00:23<00:01, 9403.15it/s]
93%|#########3| 212568/228483 [00:23<00:01, 9385.17it/s]
93%|#########3| 213507/228483 [00:23<00:01, 9383.01it/s]
94%|#########3| 214448/228483 [00:23<00:01, 9390.61it/s]
94%|#########4| 215388/228483 [00:23<00:01, 9375.04it/s]
95%|#########4| 216326/228483 [00:23<00:01, 9260.49it/s]
95%|#########5| 217266/228483 [00:23<00:01, 9301.29it/s]
95%|#########5| 218201/228483 [00:23<00:01, 9314.82it/s]
96%|#########5| 219143/228483 [00:23<00:00, 9344.65it/s]
96%|#########6| 220078/228483 [00:24<00:00, 9332.65it/s]
97%|#########6| 221027/228483 [00:24<00:00, 9377.85it/s]
97%|#########7| 221965/228483 [00:24<00:00, 9327.60it/s]
98%|#########7| 222909/228483 [00:24<00:00, 9360.59it/s]
98%|#########7| 223859/228483 [00:24<00:00, 9401.06it/s]
98%|#########8| 224810/228483 [00:24<00:00, 9432.14it/s]
99%|#########8| 225756/228483 [00:24<00:00, 9439.87it/s]
99%|#########9| 226710/228483 [00:24<00:00, 9469.47it/s]
100%|#########9| 227657/228483 [00:24<00:00, 9464.82it/s]
100%|##########| 228483/228483 [00:31<00:00, 7183.66it/s]
<nilearn.plotting.displays._slicers.ZSlicer object at 0x7f3b8c269950>
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. 4 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()
.
Out:
<nilearn.plotting.displays._slicers.ZSlicer object at 0x7f3b34b101d0>
References
- 1
Laird, Angela R., et al. “ALE meta-analysis: Controlling the false discovery rate and performing statistical contrasts.” Human brain mapping 25.1 (2005): 155-164. https://doi.org/10.1002/hbm.20136
- 2
Eickhoff, Simon B., et al. “Activation likelihood estimation meta-analysis revisited.” Neuroimage 59.3 (2012): 2349-2361. https://doi.org/10.1016/j.neuroimage.2011.09.017
- 3
Enge, Alexander, et al. “A meta-analysis of fMRI studies of semantic cognition in children.” Neuroimage 241 (2021): 118436. https://doi.org/10.1016/j.neuroimage.2021.118436
- 4
Nichols, Thomas, et al. “Valid conjunction inference with the minimum statistic.” Neuroimage 25.3 (2005): 653-660. https://doi.org/10.1016/j.neuroimage.2004.12.005
Total running time of the script: ( 2 minutes 46.796 seconds)