Note
Click here to download the full example code
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:
Run a within-sample meta-analysis of the first sample.
Characterize/summarize the results of the first meta-analysis.
Run a within-sample meta-analysis of the second sample.
Characterize/summarize the results of the second meta-analysis.
Compare the two samples with a subtraction analysis.
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:32, 1.07it/s]
3%|3 | 3/100 [00:01<01:01, 1.58it/s]
5%|5 | 5/100 [00:02<00:51, 1.84it/s]
7%|7 | 7/100 [00:03<00:50, 1.86it/s]
9%|9 | 9/100 [00:04<00:45, 1.98it/s]
11%|#1 | 11/100 [00:05<00:44, 1.98it/s]
13%|#3 | 13/100 [00:06<00:42, 2.05it/s]
15%|#5 | 15/100 [00:07<00:42, 2.02it/s]
17%|#7 | 17/100 [00:08<00:39, 2.09it/s]
19%|#9 | 19/100 [00:09<00:39, 2.04it/s]
21%|##1 | 21/100 [00:10<00:37, 2.09it/s]
23%|##3 | 23/100 [00:11<00:37, 2.03it/s]
25%|##5 | 25/100 [00:12<00:35, 2.09it/s]
27%|##7 | 27/100 [00:13<00:35, 2.04it/s]
29%|##9 | 29/100 [00:14<00:34, 2.08it/s]
31%|###1 | 31/100 [00:15<00:33, 2.04it/s]
33%|###3 | 33/100 [00:16<00:31, 2.12it/s]
35%|###5 | 35/100 [00:17<00:31, 2.07it/s]
37%|###7 | 37/100 [00:18<00:29, 2.12it/s]
39%|###9 | 39/100 [00:19<00:29, 2.07it/s]
41%|####1 | 41/100 [00:20<00:28, 2.11it/s]
43%|####3 | 43/100 [00:21<00:27, 2.07it/s]
45%|####5 | 45/100 [00:22<00:25, 2.14it/s]
47%|####6 | 47/100 [00:23<00:25, 2.10it/s]
48%|####8 | 48/100 [00:23<00:21, 2.43it/s]
49%|####9 | 49/100 [00:23<00:24, 2.05it/s]
51%|#####1 | 51/100 [00:24<00:23, 2.06it/s]
52%|#####2 | 52/100 [00:25<00:19, 2.41it/s]
53%|#####3 | 53/100 [00:25<00:22, 2.07it/s]
54%|#####4 | 54/100 [00:25<00:18, 2.51it/s]
55%|#####5 | 55/100 [00:26<00:23, 1.95it/s]
56%|#####6 | 56/100 [00:26<00:18, 2.34it/s]
57%|#####6 | 57/100 [00:27<00:21, 2.00it/s]
58%|#####8 | 58/100 [00:27<00:17, 2.43it/s]
59%|#####8 | 59/100 [00:28<00:21, 1.93it/s]
60%|###### | 60/100 [00:28<00:17, 2.29it/s]
61%|######1 | 61/100 [00:29<00:19, 2.00it/s]
62%|######2 | 62/100 [00:29<00:16, 2.37it/s]
63%|######3 | 63/100 [00:30<00:19, 1.90it/s]
64%|######4 | 64/100 [00:30<00:16, 2.22it/s]
65%|######5 | 65/100 [00:31<00:17, 2.02it/s]
66%|######6 | 66/100 [00:31<00:14, 2.33it/s]
67%|######7 | 67/100 [00:32<00:17, 1.91it/s]
68%|######8 | 68/100 [00:32<00:14, 2.26it/s]
69%|######9 | 69/100 [00:33<00:15, 2.00it/s]
70%|####### | 70/100 [00:33<00:12, 2.33it/s]
71%|#######1 | 71/100 [00:34<00:15, 1.92it/s]
72%|#######2 | 72/100 [00:34<00:12, 2.24it/s]
73%|#######3 | 73/100 [00:35<00:13, 2.03it/s]
74%|#######4 | 74/100 [00:35<00:11, 2.32it/s]
75%|#######5 | 75/100 [00:36<00:12, 1.97it/s]
76%|#######6 | 76/100 [00:36<00:10, 2.19it/s]
77%|#######7 | 77/100 [00:37<00:10, 2.10it/s]
78%|#######8 | 78/100 [00:37<00:09, 2.27it/s]
79%|#######9 | 79/100 [00:38<00:10, 2.03it/s]
80%|######## | 80/100 [00:38<00:09, 2.19it/s]
81%|########1 | 81/100 [00:38<00:08, 2.14it/s]
82%|########2 | 82/100 [00:39<00:08, 2.22it/s]
83%|########2 | 83/100 [00:39<00:08, 2.06it/s]
84%|########4 | 84/100 [00:40<00:07, 2.14it/s]
85%|########5 | 85/100 [00:40<00:06, 2.18it/s]
86%|########6 | 86/100 [00:41<00:06, 2.21it/s]
87%|########7 | 87/100 [00:41<00:06, 2.10it/s]
88%|########8 | 88/100 [00:42<00:05, 2.10it/s]
89%|########9 | 89/100 [00:42<00:04, 2.28it/s]
90%|######### | 90/100 [00:43<00:04, 2.15it/s]
91%|#########1| 91/100 [00:43<00:04, 2.16it/s]
92%|#########2| 92/100 [00:44<00:03, 2.08it/s]
93%|#########3| 93/100 [00:44<00:03, 2.28it/s]
94%|#########3| 94/100 [00:44<00:02, 2.11it/s]
95%|#########5| 95/100 [00:45<00:02, 2.18it/s]
96%|#########6| 96/100 [00:45<00:01, 2.02it/s]
97%|#########7| 97/100 [00:46<00:01, 2.38it/s]
98%|#########8| 98/100 [00:46<00:00, 2.08it/s]
99%|#########9| 99/100 [00:47<00:00, 2.27it/s]
100%|##########| 100/100 [00:47<00:00, 2.36it/s]
100%|##########| 100/100 [00:47<00:00, 2.10it/s]
0%| | 0/100 [00:00<?, ?it/s]
1%|1 | 1/100 [00:00<01:18, 1.27it/s]
3%|3 | 3/100 [00:01<00:51, 1.89it/s]
5%|5 | 5/100 [00:02<00:42, 2.24it/s]
7%|7 | 7/100 [00:03<00:41, 2.26it/s]
9%|9 | 9/100 [00:04<00:37, 2.40it/s]
11%|#1 | 11/100 [00:04<00:37, 2.37it/s]
13%|#3 | 13/100 [00:05<00:35, 2.48it/s]
14%|#4 | 14/100 [00:05<00:29, 2.88it/s]
15%|#5 | 15/100 [00:06<00:36, 2.34it/s]
16%|#6 | 16/100 [00:06<00:31, 2.70it/s]
17%|#7 | 17/100 [00:07<00:34, 2.40it/s]
18%|#8 | 18/100 [00:07<00:29, 2.77it/s]
19%|#9 | 19/100 [00:08<00:35, 2.31it/s]
20%|## | 20/100 [00:08<00:30, 2.61it/s]
21%|##1 | 21/100 [00:08<00:32, 2.45it/s]
22%|##2 | 22/100 [00:09<00:28, 2.72it/s]
23%|##3 | 23/100 [00:09<00:32, 2.35it/s]
24%|##4 | 24/100 [00:09<00:29, 2.56it/s]
25%|##5 | 25/100 [00:10<00:30, 2.49it/s]
26%|##6 | 26/100 [00:10<00:27, 2.68it/s]
27%|##7 | 27/100 [00:11<00:30, 2.41it/s]
28%|##8 | 28/100 [00:11<00:29, 2.44it/s]
29%|##9 | 29/100 [00:11<00:27, 2.57it/s]
30%|### | 30/100 [00:12<00:28, 2.50it/s]
31%|###1 | 31/100 [00:12<00:28, 2.44it/s]
32%|###2 | 32/100 [00:13<00:29, 2.34it/s]
33%|###3 | 33/100 [00:13<00:25, 2.65it/s]
34%|###4 | 34/100 [00:13<00:26, 2.49it/s]
35%|###5 | 35/100 [00:14<00:26, 2.43it/s]
36%|###6 | 36/100 [00:14<00:26, 2.40it/s]
37%|###7 | 37/100 [00:15<00:23, 2.66it/s]
38%|###8 | 38/100 [00:15<00:24, 2.53it/s]
39%|###9 | 39/100 [00:15<00:24, 2.52it/s]
40%|#### | 40/100 [00:16<00:25, 2.39it/s]
41%|####1 | 41/100 [00:16<00:22, 2.67it/s]
42%|####2 | 42/100 [00:17<00:23, 2.45it/s]
43%|####3 | 43/100 [00:17<00:22, 2.52it/s]
44%|####4 | 44/100 [00:18<00:24, 2.32it/s]
45%|####5 | 45/100 [00:18<00:20, 2.68it/s]
46%|####6 | 46/100 [00:18<00:21, 2.47it/s]
47%|####6 | 47/100 [00:19<00:21, 2.50it/s]
48%|####8 | 48/100 [00:19<00:21, 2.40it/s]
49%|####9 | 49/100 [00:19<00:19, 2.67it/s]
50%|##### | 50/100 [00:20<00:20, 2.48it/s]
51%|#####1 | 51/100 [00:20<00:19, 2.50it/s]
52%|#####2 | 52/100 [00:21<00:20, 2.37it/s]
53%|#####3 | 53/100 [00:21<00:17, 2.70it/s]
54%|#####4 | 54/100 [00:21<00:18, 2.44it/s]
55%|#####5 | 55/100 [00:22<00:17, 2.54it/s]
56%|#####6 | 56/100 [00:22<00:18, 2.35it/s]
57%|#####6 | 57/100 [00:23<00:16, 2.67it/s]
58%|#####8 | 58/100 [00:23<00:17, 2.45it/s]
59%|#####8 | 59/100 [00:23<00:16, 2.54it/s]
60%|###### | 60/100 [00:24<00:17, 2.34it/s]
61%|######1 | 61/100 [00:24<00:14, 2.70it/s]
62%|######2 | 62/100 [00:25<00:15, 2.47it/s]
63%|######3 | 63/100 [00:25<00:14, 2.54it/s]
64%|######4 | 64/100 [00:25<00:15, 2.35it/s]
65%|######5 | 65/100 [00:26<00:12, 2.73it/s]
66%|######6 | 66/100 [00:26<00:13, 2.44it/s]
67%|######7 | 67/100 [00:27<00:12, 2.64it/s]
68%|######8 | 68/100 [00:27<00:13, 2.31it/s]
69%|######9 | 69/100 [00:27<00:11, 2.76it/s]
70%|####### | 70/100 [00:28<00:12, 2.46it/s]
71%|#######1 | 71/100 [00:28<00:11, 2.52it/s]
72%|#######2 | 72/100 [00:29<00:11, 2.41it/s]
73%|#######3 | 73/100 [00:29<00:10, 2.69it/s]
74%|#######4 | 74/100 [00:29<00:10, 2.52it/s]
75%|#######5 | 75/100 [00:30<00:09, 2.55it/s]
76%|#######6 | 76/100 [00:30<00:09, 2.47it/s]
77%|#######7 | 77/100 [00:30<00:08, 2.67it/s]
78%|#######8 | 78/100 [00:31<00:08, 2.57it/s]
79%|#######9 | 79/100 [00:31<00:08, 2.47it/s]
80%|######## | 80/100 [00:32<00:08, 2.49it/s]
81%|########1 | 81/100 [00:32<00:07, 2.65it/s]
82%|########2 | 82/100 [00:32<00:06, 2.61it/s]
83%|########2 | 83/100 [00:33<00:06, 2.46it/s]
84%|########4 | 84/100 [00:33<00:06, 2.49it/s]
85%|########5 | 85/100 [00:34<00:05, 2.57it/s]
86%|########6 | 86/100 [00:34<00:05, 2.55it/s]
87%|########7 | 87/100 [00:35<00:05, 2.41it/s]
88%|########8 | 88/100 [00:35<00:04, 2.46it/s]
89%|########9 | 89/100 [00:35<00:04, 2.54it/s]
90%|######### | 90/100 [00:36<00:03, 2.61it/s]
91%|#########1| 91/100 [00:36<00:03, 2.40it/s]
92%|#########2| 92/100 [00:36<00:03, 2.54it/s]
93%|#########3| 93/100 [00:37<00:02, 2.52it/s]
94%|#########3| 94/100 [00:37<00:02, 2.71it/s]
95%|#########5| 95/100 [00:38<00:02, 2.37it/s]
96%|#########6| 96/100 [00:38<00:01, 2.58it/s]
97%|#########7| 97/100 [00:38<00:01, 2.47it/s]
98%|#########8| 98/100 [00:39<00:00, 2.73it/s]
99%|#########9| 99/100 [00:39<00:00, 2.30it/s]
100%|##########| 100/100 [00:40<00:00, 2.84it/s]
100%|##########| 100/100 [00:40<00:00, 2.49it/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.50it/s]
19%|#9 | 4/21 [00:00<00:00, 18.55it/s]
29%|##8 | 6/21 [00:00<00:00, 18.43it/s]
38%|###8 | 8/21 [00:00<00:00, 18.22it/s]
48%|####7 | 10/21 [00:00<00:00, 18.22it/s]
57%|#####7 | 12/21 [00:00<00:00, 18.18it/s]
67%|######6 | 14/21 [00:00<00:00, 18.32it/s]
76%|#######6 | 16/21 [00:00<00:00, 18.25it/s]
86%|########5 | 18/21 [00:00<00:00, 18.27it/s]
95%|#########5| 20/21 [00:01<00:00, 18.39it/s]
100%|##########| 21/21 [00:01<00:00, 18.33it/s]
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:22, 1.53s/it]
12%|#2 | 2/16 [00:03<00:21, 1.52s/it]
19%|#8 | 3/16 [00:04<00:19, 1.50s/it]
25%|##5 | 4/16 [00:06<00:18, 1.51s/it]
31%|###1 | 5/16 [00:07<00:16, 1.50s/it]
38%|###7 | 6/16 [00:09<00:14, 1.50s/it]
44%|####3 | 7/16 [00:10<00:13, 1.49s/it]
50%|##### | 8/16 [00:12<00:11, 1.50s/it]
56%|#####6 | 9/16 [00:13<00:10, 1.51s/it]
62%|######2 | 10/16 [00:15<00:08, 1.50s/it]
69%|######8 | 11/16 [00:16<00:07, 1.50s/it]
75%|#######5 | 12/16 [00:18<00:06, 1.52s/it]
81%|########1 | 13/16 [00:19<00:04, 1.52s/it]
88%|########7 | 14/16 [00:21<00:03, 1.53s/it]
94%|#########3| 15/16 [00:22<00:01, 1.54s/it]
100%|##########| 16/16 [00:24<00:00, 1.54s/it]
100%|##########| 16/16 [00:24<00:00, 1.52s/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.
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:46, 2.12it/s]
2%|2 | 2/100 [00:00<00:45, 2.17it/s]
3%|3 | 3/100 [00:01<00:45, 2.14it/s]
4%|4 | 4/100 [00:01<00:44, 2.15it/s]
5%|5 | 5/100 [00:02<00:43, 2.16it/s]
6%|6 | 6/100 [00:02<00:42, 2.19it/s]
7%|7 | 7/100 [00:03<00:42, 2.17it/s]
8%|8 | 8/100 [00:03<00:42, 2.18it/s]
9%|9 | 9/100 [00:04<00:41, 2.19it/s]
10%|# | 10/100 [00:04<00:41, 2.19it/s]
11%|#1 | 11/100 [00:05<00:40, 2.20it/s]
12%|#2 | 12/100 [00:05<00:40, 2.19it/s]
13%|#3 | 13/100 [00:05<00:39, 2.19it/s]
14%|#4 | 14/100 [00:06<00:39, 2.20it/s]
15%|#5 | 15/100 [00:06<00:38, 2.19it/s]
16%|#6 | 16/100 [00:07<00:38, 2.16it/s]
17%|#7 | 17/100 [00:07<00:38, 2.16it/s]
18%|#8 | 18/100 [00:08<00:38, 2.15it/s]
19%|#9 | 19/100 [00:08<00:37, 2.17it/s]
20%|## | 20/100 [00:09<00:36, 2.18it/s]
21%|##1 | 21/100 [00:09<00:36, 2.18it/s]
22%|##2 | 22/100 [00:10<00:35, 2.18it/s]
23%|##3 | 23/100 [00:10<00:35, 2.17it/s]
24%|##4 | 24/100 [00:11<00:34, 2.18it/s]
25%|##5 | 25/100 [00:11<00:34, 2.17it/s]
26%|##6 | 26/100 [00:11<00:34, 2.16it/s]
27%|##7 | 27/100 [00:12<00:33, 2.15it/s]
28%|##8 | 28/100 [00:12<00:33, 2.17it/s]
29%|##9 | 29/100 [00:13<00:32, 2.17it/s]
30%|### | 30/100 [00:13<00:32, 2.17it/s]
31%|###1 | 31/100 [00:14<00:31, 2.18it/s]
32%|###2 | 32/100 [00:14<00:31, 2.19it/s]
33%|###3 | 33/100 [00:15<00:30, 2.19it/s]
34%|###4 | 34/100 [00:15<00:30, 2.20it/s]
35%|###5 | 35/100 [00:16<00:29, 2.18it/s]
36%|###6 | 36/100 [00:16<00:29, 2.17it/s]
37%|###7 | 37/100 [00:17<00:29, 2.15it/s]
38%|###8 | 38/100 [00:17<00:29, 2.13it/s]
39%|###9 | 39/100 [00:17<00:28, 2.15it/s]
40%|#### | 40/100 [00:18<00:28, 2.14it/s]
41%|####1 | 41/100 [00:18<00:27, 2.15it/s]
42%|####2 | 42/100 [00:19<00:27, 2.14it/s]
43%|####3 | 43/100 [00:19<00:26, 2.14it/s]
44%|####4 | 44/100 [00:20<00:26, 2.15it/s]
45%|####5 | 45/100 [00:20<00:25, 2.15it/s]
46%|####6 | 46/100 [00:21<00:25, 2.15it/s]
47%|####6 | 47/100 [00:21<00:24, 2.15it/s]
48%|####8 | 48/100 [00:22<00:24, 2.16it/s]
49%|####9 | 49/100 [00:22<00:23, 2.17it/s]
50%|##### | 50/100 [00:23<00:23, 2.15it/s]
51%|#####1 | 51/100 [00:23<00:22, 2.16it/s]
52%|#####2 | 52/100 [00:24<00:22, 2.15it/s]
53%|#####3 | 53/100 [00:24<00:21, 2.14it/s]
54%|#####4 | 54/100 [00:24<00:21, 2.15it/s]
55%|#####5 | 55/100 [00:25<00:20, 2.16it/s]
56%|#####6 | 56/100 [00:25<00:20, 2.17it/s]
57%|#####6 | 57/100 [00:26<00:19, 2.17it/s]
58%|#####8 | 58/100 [00:26<00:19, 2.16it/s]
59%|#####8 | 59/100 [00:27<00:19, 2.15it/s]
60%|###### | 60/100 [00:27<00:18, 2.16it/s]
61%|######1 | 61/100 [00:28<00:18, 2.16it/s]
62%|######2 | 62/100 [00:28<00:17, 2.15it/s]
63%|######3 | 63/100 [00:29<00:17, 2.16it/s]
64%|######4 | 64/100 [00:29<00:16, 2.15it/s]
65%|######5 | 65/100 [00:30<00:16, 2.17it/s]
66%|######6 | 66/100 [00:30<00:16, 2.05it/s]
67%|######7 | 67/100 [00:31<00:15, 2.08it/s]
68%|######8 | 68/100 [00:31<00:15, 2.11it/s]
69%|######9 | 69/100 [00:31<00:14, 2.12it/s]
70%|####### | 70/100 [00:32<00:14, 2.12it/s]
71%|#######1 | 71/100 [00:32<00:13, 2.13it/s]
72%|#######2 | 72/100 [00:33<00:13, 2.14it/s]
73%|#######3 | 73/100 [00:33<00:12, 2.15it/s]
74%|#######4 | 74/100 [00:34<00:12, 2.16it/s]
75%|#######5 | 75/100 [00:34<00:11, 2.18it/s]
76%|#######6 | 76/100 [00:35<00:10, 2.19it/s]
77%|#######7 | 77/100 [00:35<00:10, 2.18it/s]
78%|#######8 | 78/100 [00:36<00:10, 2.17it/s]
79%|#######9 | 79/100 [00:36<00:09, 2.17it/s]
80%|######## | 80/100 [00:37<00:09, 2.19it/s]
81%|########1 | 81/100 [00:37<00:08, 2.18it/s]
82%|########2 | 82/100 [00:37<00:08, 2.18it/s]
83%|########2 | 83/100 [00:38<00:07, 2.19it/s]
84%|########4 | 84/100 [00:38<00:07, 2.19it/s]
85%|########5 | 85/100 [00:39<00:06, 2.18it/s]
86%|########6 | 86/100 [00:39<00:06, 2.16it/s]
87%|########7 | 87/100 [00:40<00:05, 2.18it/s]
88%|########8 | 88/100 [00:40<00:05, 2.17it/s]
89%|########9 | 89/100 [00:41<00:05, 2.16it/s]
90%|######### | 90/100 [00:41<00:04, 2.17it/s]
91%|#########1| 91/100 [00:42<00:04, 2.16it/s]
92%|#########2| 92/100 [00:42<00:03, 2.16it/s]
93%|#########3| 93/100 [00:43<00:03, 2.15it/s]
94%|#########3| 94/100 [00:43<00:02, 2.16it/s]
95%|#########5| 95/100 [00:43<00:02, 2.15it/s]
96%|#########6| 96/100 [00:44<00:01, 2.15it/s]
97%|#########7| 97/100 [00:44<00:01, 2.16it/s]
98%|#########8| 98/100 [00:45<00:00, 2.17it/s]
99%|#########9| 99/100 [00:45<00:00, 2.18it/s]
100%|##########| 100/100 [00:46<00:00, 2.18it/s]
100%|##########| 100/100 [00:46<00:00, 2.16it/s]
0%| | 0/228483 [00:00<?, ?it/s]
0%| | 958/228483 [00:00<00:23, 9570.45it/s]
1%| | 1917/228483 [00:00<00:23, 9579.25it/s]
1%|1 | 2875/228483 [00:00<00:23, 9560.34it/s]
2%|1 | 3832/228483 [00:00<00:23, 9534.06it/s]
2%|2 | 4790/228483 [00:00<00:23, 9548.35it/s]
3%|2 | 5745/228483 [00:00<00:23, 9532.97it/s]
3%|2 | 6699/228483 [00:00<00:23, 9532.30it/s]
3%|3 | 7653/228483 [00:00<00:23, 9473.21it/s]
4%|3 | 8609/228483 [00:00<00:23, 9498.27it/s]
4%|4 | 9559/228483 [00:01<00:23, 9488.70it/s]
5%|4 | 10508/228483 [00:01<00:23, 9458.76it/s]
5%|5 | 11456/228483 [00:01<00:22, 9463.75it/s]
5%|5 | 12407/228483 [00:01<00:22, 9477.18it/s]
6%|5 | 13360/228483 [00:01<00:22, 9492.53it/s]
6%|6 | 14318/228483 [00:01<00:22, 9518.17it/s]
7%|6 | 15271/228483 [00:01<00:22, 9521.07it/s]
7%|7 | 16224/228483 [00:01<00:22, 9451.04it/s]
8%|7 | 17174/228483 [00:01<00:22, 9465.01it/s]
8%|7 | 18134/228483 [00:01<00:22, 9503.86it/s]
8%|8 | 19090/228483 [00:02<00:21, 9518.06it/s]
9%|8 | 20045/228483 [00:02<00:21, 9526.96it/s]
9%|9 | 20998/228483 [00:02<00:21, 9526.49it/s]
10%|9 | 21951/228483 [00:02<00:21, 9487.24it/s]
10%|# | 22900/228483 [00:02<00:21, 9462.97it/s]
10%|# | 23855/228483 [00:02<00:21, 9486.39it/s]
11%|# | 24805/228483 [00:02<00:21, 9487.76it/s]
11%|#1 | 25754/228483 [00:02<00:21, 9457.93it/s]
12%|#1 | 26700/228483 [00:02<00:21, 9436.15it/s]
12%|#2 | 27654/228483 [00:02<00:21, 9465.66it/s]
13%|#2 | 28606/228483 [00:03<00:21, 9481.17it/s]
13%|#2 | 29555/228483 [00:03<00:21, 9438.79it/s]
13%|#3 | 30503/228483 [00:03<00:20, 9449.10it/s]
14%|#3 | 31448/228483 [00:03<00:20, 9438.60it/s]
14%|#4 | 32392/228483 [00:03<00:29, 6695.98it/s]
15%|#4 | 33336/228483 [00:03<00:26, 7332.64it/s]
15%|#5 | 34274/228483 [00:03<00:24, 7842.28it/s]
15%|#5 | 35223/228483 [00:03<00:23, 8274.93it/s]
16%|#5 | 36159/228483 [00:03<00:22, 8570.14it/s]
16%|#6 | 37106/228483 [00:04<00:21, 8820.02it/s]
17%|#6 | 38061/228483 [00:04<00:21, 9028.19it/s]
17%|#7 | 39013/228483 [00:04<00:20, 9169.15it/s]
17%|#7 | 39966/228483 [00:04<00:20, 9272.53it/s]
18%|#7 | 40906/228483 [00:04<00:20, 9306.04it/s]
18%|#8 | 41865/228483 [00:04<00:19, 9389.97it/s]
19%|#8 | 42815/228483 [00:04<00:19, 9419.98it/s]
19%|#9 | 43762/228483 [00:04<00:19, 9422.93it/s]
20%|#9 | 44708/228483 [00:04<00:19, 9387.69it/s]
20%|#9 | 45669/228483 [00:04<00:19, 9452.50it/s]
20%|## | 46621/228483 [00:05<00:19, 9472.54it/s]
21%|## | 47570/228483 [00:05<00:19, 9401.14it/s]
21%|##1 | 48512/228483 [00:05<00:19, 9394.65it/s]
22%|##1 | 49453/228483 [00:05<00:19, 9398.85it/s]
22%|##2 | 50405/228483 [00:05<00:18, 9432.25it/s]
22%|##2 | 51353/228483 [00:05<00:18, 9443.73it/s]
23%|##2 | 52298/228483 [00:05<00:18, 9439.73it/s]
23%|##3 | 53246/228483 [00:05<00:18, 9451.44it/s]
24%|##3 | 54192/228483 [00:05<00:18, 9450.34it/s]
24%|##4 | 55138/228483 [00:05<00:18, 9438.75it/s]
25%|##4 | 56094/228483 [00:06<00:18, 9474.27it/s]
25%|##4 | 57042/228483 [00:06<00:18, 9466.39it/s]
25%|##5 | 57996/228483 [00:06<00:17, 9485.85it/s]
26%|##5 | 58949/228483 [00:06<00:17, 9496.42it/s]
26%|##6 | 59904/228483 [00:06<00:17, 9512.27it/s]
27%|##6 | 60859/228483 [00:06<00:17, 9522.19it/s]
27%|##7 | 61818/228483 [00:06<00:17, 9540.36it/s]
27%|##7 | 62773/228483 [00:06<00:17, 9500.08it/s]
28%|##7 | 63724/228483 [00:06<00:17, 9479.35it/s]
28%|##8 | 64679/228483 [00:06<00:17, 9498.75it/s]
29%|##8 | 65629/228483 [00:07<00:17, 9492.36it/s]
29%|##9 | 66579/228483 [00:07<00:17, 9488.86it/s]
30%|##9 | 67530/228483 [00:07<00:16, 9494.14it/s]
30%|##9 | 68480/228483 [00:07<00:16, 9491.11it/s]
30%|### | 69430/228483 [00:07<00:16, 9461.35it/s]
31%|### | 70380/228483 [00:07<00:16, 9471.03it/s]
31%|###1 | 71328/228483 [00:07<00:16, 9454.28it/s]
32%|###1 | 72274/228483 [00:07<00:23, 6535.92it/s]
32%|###2 | 73227/228483 [00:08<00:21, 7217.72it/s]
32%|###2 | 74176/228483 [00:08<00:19, 7775.10it/s]
33%|###2 | 75131/228483 [00:08<00:18, 8235.72it/s]
33%|###3 | 76089/228483 [00:08<00:17, 8597.89it/s]
34%|###3 | 77045/228483 [00:08<00:17, 8866.00it/s]
34%|###4 | 77999/228483 [00:08<00:16, 9056.66it/s]
35%|###4 | 78958/228483 [00:08<00:16, 9208.94it/s]
35%|###4 | 79901/228483 [00:08<00:16, 9273.25it/s]
35%|###5 | 80859/228483 [00:08<00:15, 9361.78it/s]
36%|###5 | 81814/228483 [00:08<00:15, 9416.55it/s]
36%|###6 | 82763/228483 [00:09<00:15, 9416.92it/s]
37%|###6 | 83717/228483 [00:09<00:15, 9452.75it/s]
37%|###7 | 84672/228483 [00:09<00:15, 9480.35it/s]
37%|###7 | 85631/228483 [00:09<00:15, 9512.88it/s]
38%|###7 | 86585/228483 [00:09<00:14, 9500.77it/s]
38%|###8 | 87538/228483 [00:09<00:14, 9506.99it/s]
39%|###8 | 88490/228483 [00:09<00:14, 9490.18it/s]
39%|###9 | 89442/228483 [00:09<00:14, 9497.43it/s]
40%|###9 | 90401/228483 [00:09<00:14, 9523.08it/s]
40%|###9 | 91357/228483 [00:09<00:14, 9531.68it/s]
40%|#### | 92313/228483 [00:10<00:14, 9540.11it/s]
41%|#### | 93271/228483 [00:10<00:14, 9549.26it/s]
41%|####1 | 94240/228483 [00:10<00:14, 9588.54it/s]
42%|####1 | 95199/228483 [00:10<00:13, 9566.48it/s]
42%|####2 | 96156/228483 [00:10<00:13, 9533.47it/s]
43%|####2 | 97114/228483 [00:10<00:13, 9546.25it/s]
43%|####2 | 98070/228483 [00:10<00:13, 9549.06it/s]
43%|####3 | 99032/228483 [00:10<00:13, 9567.86it/s]
44%|####3 | 99989/228483 [00:10<00:13, 9502.87it/s]
44%|####4 | 100942/228483 [00:10<00:13, 9508.45it/s]
45%|####4 | 101897/228483 [00:11<00:13, 9520.14it/s]
45%|####5 | 102850/228483 [00:11<00:13, 9521.78it/s]
45%|####5 | 103803/228483 [00:11<00:13, 9490.30it/s]
46%|####5 | 104753/228483 [00:11<00:13, 9456.40it/s]
46%|####6 | 105716/228483 [00:11<00:12, 9506.57it/s]
47%|####6 | 106688/228483 [00:11<00:12, 9568.18it/s]
47%|####7 | 107645/228483 [00:11<00:12, 9558.03it/s]
48%|####7 | 108601/228483 [00:11<00:12, 9439.36it/s]
48%|####7 | 109552/228483 [00:11<00:12, 9459.94it/s]
48%|####8 | 110519/228483 [00:11<00:12, 9521.23it/s]
49%|####8 | 111481/228483 [00:12<00:12, 9549.48it/s]
49%|####9 | 112437/228483 [00:12<00:12, 9511.70it/s]
50%|####9 | 113389/228483 [00:12<00:12, 9511.51it/s]
50%|##### | 114348/228483 [00:12<00:11, 9534.45it/s]
50%|##### | 115308/228483 [00:12<00:11, 9551.65it/s]
51%|##### | 116266/228483 [00:12<00:11, 9557.84it/s]
51%|#####1 | 117222/228483 [00:12<00:11, 9535.67it/s]
52%|#####1 | 118176/228483 [00:12<00:11, 9525.88it/s]
52%|#####2 | 119130/228483 [00:12<00:11, 9529.27it/s]
53%|#####2 | 120089/228483 [00:12<00:11, 9546.20it/s]
53%|#####2 | 121044/228483 [00:13<00:11, 9533.29it/s]
53%|#####3 | 121998/228483 [00:13<00:17, 6220.02it/s]
54%|#####3 | 122773/228483 [00:13<00:16, 6553.59it/s]
54%|#####4 | 123740/228483 [00:13<00:14, 7293.92it/s]
55%|#####4 | 124672/228483 [00:13<00:13, 7805.53it/s]
55%|#####4 | 125607/228483 [00:13<00:12, 8215.51it/s]
55%|#####5 | 126555/228483 [00:13<00:11, 8562.68it/s]
56%|#####5 | 127508/228483 [00:13<00:11, 8834.28it/s]
56%|#####6 | 128455/228483 [00:14<00:11, 9014.56it/s]
57%|#####6 | 129408/228483 [00:14<00:10, 9163.16it/s]
57%|#####7 | 130364/228483 [00:14<00:10, 9277.02it/s]
57%|#####7 | 131324/228483 [00:14<00:10, 9372.21it/s]
58%|#####7 | 132279/228483 [00:14<00:10, 9423.62it/s]
58%|#####8 | 133240/228483 [00:14<00:10, 9476.56it/s]
59%|#####8 | 134193/228483 [00:14<00:09, 9469.36it/s]
59%|#####9 | 135151/228483 [00:14<00:09, 9501.12it/s]
60%|#####9 | 136104/228483 [00:14<00:09, 9473.70it/s]
60%|#####9 | 137067/228483 [00:14<00:09, 9516.05it/s]
60%|###### | 138020/228483 [00:15<00:09, 9483.51it/s]
61%|###### | 138970/228483 [00:15<00:09, 9140.84it/s]
61%|######1 | 139888/228483 [00:15<00:10, 8278.00it/s]
62%|######1 | 140834/228483 [00:15<00:10, 8599.28it/s]
62%|######2 | 141794/228483 [00:15<00:09, 8877.86it/s]
62%|######2 | 142734/228483 [00:15<00:09, 9027.10it/s]
63%|######2 | 143689/228483 [00:15<00:09, 9178.49it/s]
63%|######3 | 144639/228483 [00:15<00:09, 9270.43it/s]
64%|######3 | 145598/228483 [00:15<00:08, 9363.66it/s]
64%|######4 | 146539/228483 [00:15<00:08, 9365.84it/s]
65%|######4 | 147490/228483 [00:16<00:08, 9407.58it/s]
65%|######4 | 148456/228483 [00:16<00:08, 9482.35it/s]
65%|######5 | 149407/228483 [00:16<00:08, 9488.11it/s]
66%|######5 | 150368/228483 [00:16<00:08, 9521.92it/s]
66%|######6 | 151321/228483 [00:16<00:08, 9513.02it/s]
67%|######6 | 152281/228483 [00:16<00:07, 9536.68it/s]
67%|######7 | 153238/228483 [00:16<00:07, 9545.88it/s]
67%|######7 | 154199/228483 [00:16<00:07, 9563.32it/s]
68%|######7 | 155156/228483 [00:16<00:07, 9481.07it/s]
68%|######8 | 156107/228483 [00:16<00:07, 9488.86it/s]
69%|######8 | 157066/228483 [00:17<00:07, 9517.50it/s]
69%|######9 | 158018/228483 [00:17<00:07, 9481.35it/s]
70%|######9 | 158967/228483 [00:17<00:07, 9456.95it/s]
70%|######9 | 159913/228483 [00:17<00:07, 9435.12it/s]
70%|####### | 160870/228483 [00:17<00:07, 9475.23it/s]
71%|####### | 161818/228483 [00:17<00:07, 9424.10it/s]
71%|#######1 | 162765/228483 [00:17<00:06, 9435.53it/s]
72%|#######1 | 163709/228483 [00:17<00:06, 9415.61it/s]
72%|#######2 | 164654/228483 [00:17<00:06, 9423.14it/s]
72%|#######2 | 165597/228483 [00:17<00:06, 9409.38it/s]
73%|#######2 | 166553/228483 [00:18<00:06, 9452.27it/s]
73%|#######3 | 167504/228483 [00:18<00:06, 9466.77it/s]
74%|#######3 | 168455/228483 [00:18<00:06, 9479.62it/s]
74%|#######4 | 169403/228483 [00:18<00:06, 9475.60it/s]
75%|#######4 | 170351/228483 [00:18<00:06, 9463.85it/s]
75%|#######4 | 171308/228483 [00:18<00:06, 9494.42it/s]
75%|#######5 | 172258/228483 [00:18<00:05, 9478.15it/s]
76%|#######5 | 173206/228483 [00:18<00:05, 9441.90it/s]
76%|#######6 | 174151/228483 [00:18<00:05, 9416.62it/s]
77%|#######6 | 175098/228483 [00:18<00:05, 9432.47it/s]
77%|#######7 | 176042/228483 [00:19<00:05, 9424.50it/s]
77%|#######7 | 176985/228483 [00:19<00:05, 9386.50it/s]
78%|#######7 | 177924/228483 [00:19<00:05, 9375.17it/s]
78%|#######8 | 178868/228483 [00:19<00:05, 9392.98it/s]
79%|#######8 | 179817/228483 [00:19<00:05, 9419.22it/s]
79%|#######9 | 180759/228483 [00:19<00:05, 9390.32it/s]
80%|#######9 | 181709/228483 [00:19<00:04, 9420.83it/s]
80%|#######9 | 182652/228483 [00:19<00:04, 9382.84it/s]
80%|######## | 183597/228483 [00:19<00:04, 9402.31it/s]
81%|######## | 184538/228483 [00:20<00:07, 5930.58it/s]
81%|########1 | 185484/228483 [00:20<00:06, 6678.59it/s]
82%|########1 | 186416/228483 [00:20<00:05, 7292.72it/s]
82%|########2 | 187363/228483 [00:20<00:05, 7835.38it/s]
82%|########2 | 188315/228483 [00:20<00:04, 8275.38it/s]
83%|########2 | 189252/228483 [00:20<00:04, 8571.89it/s]
83%|########3 | 190201/228483 [00:20<00:04, 8829.01it/s]
84%|########3 | 191127/228483 [00:20<00:04, 8951.13it/s]
84%|########4 | 192080/228483 [00:20<00:03, 9116.20it/s]
84%|########4 | 193017/228483 [00:21<00:03, 9188.92it/s]
85%|########4 | 193953/228483 [00:21<00:03, 9238.33it/s]
85%|########5 | 194887/228483 [00:21<00:03, 9259.83it/s]
86%|########5 | 195828/228483 [00:21<00:03, 9302.63it/s]
86%|########6 | 196771/228483 [00:21<00:03, 9340.49it/s]
87%|########6 | 197709/228483 [00:21<00:03, 9344.98it/s]
87%|########6 | 198646/228483 [00:21<00:03, 9335.32it/s]
87%|########7 | 199588/228483 [00:21<00:03, 9358.22it/s]
88%|########7 | 200532/228483 [00:21<00:02, 9381.38it/s]
88%|########8 | 201471/228483 [00:21<00:02, 9364.84it/s]
89%|########8 | 202418/228483 [00:22<00:02, 9393.58it/s]
89%|########9 | 203358/228483 [00:22<00:02, 9379.96it/s]
89%|########9 | 204304/228483 [00:22<00:02, 9401.51it/s]
90%|########9 | 205245/228483 [00:22<00:02, 9395.98it/s]
90%|######### | 206187/228483 [00:22<00:02, 9402.62it/s]
91%|######### | 207130/228483 [00:22<00:02, 9408.49it/s]
91%|#########1| 208071/228483 [00:22<00:02, 9383.08it/s]
91%|#########1| 209010/228483 [00:22<00:02, 9375.88it/s]
92%|#########1| 209948/228483 [00:22<00:01, 9337.15it/s]
92%|#########2| 210901/228483 [00:22<00:01, 9394.58it/s]
93%|#########2| 211841/228483 [00:23<00:01, 9371.41it/s]
93%|#########3| 212779/228483 [00:23<00:01, 9371.30it/s]
94%|#########3| 213717/228483 [00:23<00:01, 9347.33it/s]
94%|#########3| 214654/228483 [00:23<00:01, 9353.49it/s]
94%|#########4| 215600/228483 [00:23<00:01, 9384.79it/s]
95%|#########4| 216539/228483 [00:23<00:01, 9371.29it/s]
95%|#########5| 217479/228483 [00:23<00:01, 9379.68it/s]
96%|#########5| 218417/228483 [00:23<00:01, 9375.69it/s]
96%|#########6| 219363/228483 [00:23<00:00, 9400.01it/s]
96%|#########6| 220304/228483 [00:23<00:00, 9347.15it/s]
97%|#########6| 221258/228483 [00:24<00:00, 9402.04it/s]
97%|#########7| 222199/228483 [00:24<00:00, 9385.30it/s]
98%|#########7| 223149/228483 [00:24<00:00, 9419.51it/s]
98%|#########8| 224092/228483 [00:24<00:00, 9411.29it/s]
98%|#########8| 225042/228483 [00:24<00:00, 9435.54it/s]
99%|#########8| 225995/228483 [00:24<00:00, 9462.88it/s]
99%|#########9| 226942/228483 [00:24<00:00, 9458.02it/s]
100%|#########9| 227888/228483 [00:24<00:00, 9452.78it/s]
100%|##########| 228483/228483 [00:31<00:00, 7160.56it/s]
<nilearn.plotting.displays._slicers.ZSlicer object at 0x7fe0a3df4650>
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 0x7fe0a0641fd0>
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.
Total running time of the script: ( 3 minutes 23.964 seconds)