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:01<01:43, 1.05s/it]
3%|3 | 3/100 [00:02<01:05, 1.48it/s]
5%|5 | 5/100 [00:03<00:54, 1.76it/s]
6%|6 | 6/100 [00:03<00:43, 2.15it/s]
7%|7 | 7/100 [00:04<00:55, 1.69it/s]
9%|9 | 9/100 [00:05<00:49, 1.84it/s]
10%|# | 10/100 [00:05<00:40, 2.23it/s]
11%|#1 | 11/100 [00:06<00:51, 1.74it/s]
13%|#3 | 13/100 [00:07<00:46, 1.87it/s]
14%|#4 | 14/100 [00:07<00:39, 2.19it/s]
15%|#5 | 15/100 [00:08<00:48, 1.76it/s]
17%|#7 | 17/100 [00:09<00:44, 1.87it/s]
18%|#8 | 18/100 [00:09<00:37, 2.20it/s]
19%|#9 | 19/100 [00:10<00:46, 1.74it/s]
21%|##1 | 21/100 [00:11<00:41, 1.90it/s]
22%|##2 | 22/100 [00:11<00:34, 2.26it/s]
23%|##3 | 23/100 [00:12<00:43, 1.76it/s]
25%|##5 | 25/100 [00:13<00:39, 1.88it/s]
27%|##7 | 27/100 [00:14<00:39, 1.86it/s]
29%|##9 | 29/100 [00:15<00:36, 1.95it/s]
30%|### | 30/100 [00:15<00:31, 2.24it/s]
31%|###1 | 31/100 [00:16<00:37, 1.82it/s]
33%|###3 | 33/100 [00:17<00:34, 1.92it/s]
34%|###4 | 34/100 [00:17<00:29, 2.24it/s]
35%|###5 | 35/100 [00:18<00:36, 1.76it/s]
37%|###7 | 37/100 [00:19<00:33, 1.88it/s]
39%|###9 | 39/100 [00:20<00:32, 1.90it/s]
40%|#### | 40/100 [00:20<00:26, 2.26it/s]
41%|####1 | 41/100 [00:21<00:33, 1.75it/s]
43%|####3 | 43/100 [00:22<00:30, 1.87it/s]
45%|####5 | 45/100 [00:23<00:29, 1.87it/s]
47%|####6 | 47/100 [00:24<00:27, 1.90it/s]
49%|####9 | 49/100 [00:25<00:26, 1.91it/s]
51%|#####1 | 51/100 [00:26<00:25, 1.93it/s]
53%|#####3 | 53/100 [00:27<00:24, 1.93it/s]
55%|#####5 | 55/100 [00:28<00:23, 1.95it/s]
57%|#####6 | 57/100 [00:29<00:22, 1.94it/s]
59%|#####8 | 59/100 [00:30<00:20, 1.98it/s]
61%|######1 | 61/100 [00:31<00:19, 1.96it/s]
63%|######3 | 63/100 [00:32<00:18, 1.98it/s]
65%|######5 | 65/100 [00:33<00:17, 1.96it/s]
67%|######7 | 67/100 [00:34<00:16, 1.98it/s]
69%|######9 | 69/100 [00:36<00:16, 1.93it/s]
71%|#######1 | 71/100 [00:36<00:14, 1.98it/s]
72%|#######2 | 72/100 [00:37<00:12, 2.25it/s]
73%|#######3 | 73/100 [00:38<00:14, 1.82it/s]
75%|#######5 | 75/100 [00:39<00:13, 1.90it/s]
76%|#######6 | 76/100 [00:39<00:10, 2.28it/s]
77%|#######7 | 77/100 [00:40<00:12, 1.79it/s]
79%|#######9 | 79/100 [00:41<00:11, 1.84it/s]
81%|########1 | 81/100 [00:42<00:10, 1.89it/s]
83%|########2 | 83/100 [00:43<00:08, 1.92it/s]
85%|########5 | 85/100 [00:44<00:07, 1.95it/s]
87%|########7 | 87/100 [00:45<00:06, 1.93it/s]
89%|########9 | 89/100 [00:46<00:05, 1.97it/s]
90%|######### | 90/100 [00:46<00:04, 2.28it/s]
91%|#########1| 91/100 [00:47<00:04, 1.82it/s]
93%|#########3| 93/100 [00:48<00:03, 1.94it/s]
94%|#########3| 94/100 [00:48<00:02, 2.23it/s]
95%|#########5| 95/100 [00:49<00:02, 1.83it/s]
96%|#########6| 96/100 [00:49<00:01, 2.29it/s]
97%|#########7| 97/100 [00:50<00:01, 1.84it/s]
98%|#########8| 98/100 [00:50<00:00, 2.19it/s]
99%|#########9| 99/100 [00:51<00:00, 1.74it/s]
100%|##########| 100/100 [00:51<00:00, 1.94it/s]
0%| | 0/100 [00:00<?, ?it/s]
1%|1 | 1/100 [00:00<01:32, 1.07it/s]
3%|3 | 3/100 [00:01<00:56, 1.73it/s]
5%|5 | 5/100 [00:02<00:46, 2.05it/s]
7%|7 | 7/100 [00:03<00:44, 2.08it/s]
9%|9 | 9/100 [00:04<00:41, 2.20it/s]
11%|#1 | 11/100 [00:05<00:41, 2.15it/s]
13%|#3 | 13/100 [00:06<00:38, 2.25it/s]
15%|#5 | 15/100 [00:07<00:38, 2.21it/s]
17%|#7 | 17/100 [00:07<00:35, 2.32it/s]
18%|#8 | 18/100 [00:07<00:30, 2.70it/s]
19%|#9 | 19/100 [00:09<00:40, 1.99it/s]
20%|## | 20/100 [00:09<00:33, 2.41it/s]
21%|##1 | 21/100 [00:09<00:38, 2.08it/s]
22%|##2 | 22/100 [00:09<00:29, 2.61it/s]
23%|##3 | 23/100 [00:10<00:39, 1.97it/s]
24%|##4 | 24/100 [00:10<00:30, 2.53it/s]
25%|##5 | 25/100 [00:11<00:35, 2.09it/s]
26%|##6 | 26/100 [00:11<00:27, 2.68it/s]
27%|##7 | 27/100 [00:12<00:36, 2.02it/s]
28%|##8 | 28/100 [00:12<00:29, 2.45it/s]
29%|##9 | 29/100 [00:13<00:32, 2.16it/s]
30%|### | 30/100 [00:13<00:27, 2.54it/s]
31%|###1 | 31/100 [00:14<00:33, 2.08it/s]
32%|###2 | 32/100 [00:14<00:28, 2.39it/s]
33%|###3 | 33/100 [00:14<00:30, 2.20it/s]
34%|###4 | 34/100 [00:15<00:26, 2.52it/s]
35%|###5 | 35/100 [00:15<00:30, 2.10it/s]
36%|###6 | 36/100 [00:16<00:26, 2.38it/s]
37%|###7 | 37/100 [00:16<00:28, 2.24it/s]
38%|###8 | 38/100 [00:16<00:24, 2.54it/s]
39%|###9 | 39/100 [00:17<00:28, 2.12it/s]
40%|#### | 40/100 [00:17<00:25, 2.40it/s]
41%|####1 | 41/100 [00:18<00:26, 2.25it/s]
42%|####2 | 42/100 [00:18<00:22, 2.54it/s]
43%|####3 | 43/100 [00:19<00:26, 2.13it/s]
44%|####4 | 44/100 [00:19<00:23, 2.41it/s]
45%|####5 | 45/100 [00:20<00:24, 2.27it/s]
46%|####6 | 46/100 [00:20<00:21, 2.55it/s]
47%|####6 | 47/100 [00:21<00:24, 2.12it/s]
48%|####8 | 48/100 [00:21<00:21, 2.41it/s]
49%|####9 | 49/100 [00:21<00:22, 2.23it/s]
50%|##### | 50/100 [00:22<00:19, 2.53it/s]
51%|#####1 | 51/100 [00:22<00:23, 2.12it/s]
52%|#####2 | 52/100 [00:23<00:19, 2.40it/s]
53%|#####3 | 53/100 [00:23<00:20, 2.25it/s]
54%|#####4 | 54/100 [00:23<00:18, 2.54it/s]
55%|#####5 | 55/100 [00:24<00:21, 2.13it/s]
56%|#####6 | 56/100 [00:24<00:18, 2.42it/s]
57%|#####6 | 57/100 [00:25<00:19, 2.21it/s]
58%|#####8 | 58/100 [00:25<00:16, 2.61it/s]
59%|#####8 | 59/100 [00:26<00:19, 2.07it/s]
60%|###### | 60/100 [00:26<00:16, 2.44it/s]
61%|######1 | 61/100 [00:27<00:18, 2.14it/s]
62%|######2 | 62/100 [00:27<00:14, 2.59it/s]
63%|######3 | 63/100 [00:28<00:18, 2.02it/s]
64%|######4 | 64/100 [00:28<00:14, 2.49it/s]
65%|######5 | 65/100 [00:28<00:16, 2.14it/s]
66%|######6 | 66/100 [00:29<00:13, 2.60it/s]
67%|######7 | 67/100 [00:29<00:16, 2.02it/s]
68%|######8 | 68/100 [00:30<00:13, 2.45it/s]
69%|######9 | 69/100 [00:30<00:14, 2.15it/s]
70%|####### | 70/100 [00:30<00:11, 2.58it/s]
71%|#######1 | 71/100 [00:31<00:14, 2.03it/s]
72%|#######2 | 72/100 [00:31<00:11, 2.45it/s]
73%|#######3 | 73/100 [00:32<00:12, 2.13it/s]
74%|#######4 | 74/100 [00:32<00:09, 2.67it/s]
75%|#######5 | 75/100 [00:33<00:12, 2.02it/s]
76%|#######6 | 76/100 [00:33<00:09, 2.58it/s]
77%|#######7 | 77/100 [00:34<00:10, 2.12it/s]
79%|#######9 | 79/100 [00:35<00:09, 2.14it/s]
80%|######## | 80/100 [00:35<00:07, 2.66it/s]
81%|########1 | 81/100 [00:35<00:08, 2.14it/s]
83%|########2 | 83/100 [00:36<00:07, 2.13it/s]
85%|########5 | 85/100 [00:37<00:06, 2.22it/s]
87%|########7 | 87/100 [00:38<00:06, 2.14it/s]
89%|########9 | 89/100 [00:39<00:04, 2.25it/s]
91%|#########1| 91/100 [00:40<00:04, 2.23it/s]
93%|#########3| 93/100 [00:41<00:03, 2.30it/s]
95%|#########5| 95/100 [00:42<00:02, 2.25it/s]
97%|#########7| 97/100 [00:42<00:01, 2.32it/s]
99%|#########9| 99/100 [00:43<00:00, 2.28it/s]
100%|##########| 100/100 [00:43<00:00, 2.27it/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_corrected_results = counter.transform(knowledge_corrected_results)
Out:
0%| | 0/21 [00:00<?, ?it/s]
10%|9 | 2/21 [00:00<00:01, 12.38it/s]
19%|#9 | 4/21 [00:00<00:01, 10.80it/s]
29%|##8 | 6/21 [00:00<00:01, 10.58it/s]
38%|###8 | 8/21 [00:00<00:01, 10.28it/s]
48%|####7 | 10/21 [00:00<00:01, 10.16it/s]
57%|#####7 | 12/21 [00:01<00:00, 10.04it/s]
67%|######6 | 14/21 [00:01<00:00, 9.98it/s]
76%|#######6 | 16/21 [00:01<00:00, 9.93it/s]
81%|######## | 17/21 [00:01<00:00, 9.92it/s]
86%|########5 | 18/21 [00:01<00:00, 9.90it/s]
90%|######### | 19/21 [00:01<00:00, 9.88it/s]
95%|#########5| 20/21 [00:01<00:00, 9.88it/s]
100%|##########| 21/21 [00:02<00:00, 9.85it/s]
100%|##########| 21/21 [00:02<00:00, 10.10it/s]
Clusters table.
knowledge_clusters_table = knowledge_corrected_results.tables[
"z_desc-size_level-cluster_corr-FWE_method-montecarlo_tab-clust"
]
knowledge_clusters_table.head(10)
Contribution table. Here PostiveTail
refers to clusters with positive statistics.
knowledge_count_table = knowledge_corrected_results.tables[
"z_desc-size_level-cluster_corr-FWE_method-montecarlo_diag-FocusCounter_tab-counts"
]
knowledge_count_table.head(10)
from nimare.diagnostics import Jackknife
jackknife = Jackknife(
target_image="z_desc-size_level-cluster_corr-FWE_method-montecarlo",
voxel_thresh=None,
)
related_corrected_results = jackknife.transform(related_corrected_results)
related_jackknife_table = related_corrected_results.tables[
"z_desc-size_level-cluster_corr-FWE_method-montecarlo_diag-Jackknife_tab-counts"
]
related_jackknife_table.head(10)
Out:
0%| | 0/16 [00:00<?, ?it/s]
6%|6 | 1/16 [00:02<00:33, 2.26s/it]
12%|#2 | 2/16 [00:04<00:31, 2.25s/it]
19%|#8 | 3/16 [00:06<00:29, 2.24s/it]
25%|##5 | 4/16 [00:09<00:27, 2.25s/it]
31%|###1 | 5/16 [00:11<00:24, 2.24s/it]
38%|###7 | 6/16 [00:13<00:22, 2.24s/it]
44%|####3 | 7/16 [00:15<00:20, 2.22s/it]
50%|##### | 8/16 [00:17<00:17, 2.23s/it]
56%|#####6 | 9/16 [00:20<00:15, 2.24s/it]
62%|######2 | 10/16 [00:22<00:13, 2.24s/it]
69%|######8 | 11/16 [00:24<00:11, 2.24s/it]
75%|#######5 | 12/16 [00:26<00:08, 2.24s/it]
81%|########1 | 13/16 [00:29<00:06, 2.24s/it]
88%|########7 | 14/16 [00:31<00:04, 2.24s/it]
94%|#########3| 15/16 [00:33<00:02, 2.24s/it]
100%|##########| 16/16 [00:35<00:00, 2.23s/it]
100%|##########| 16/16 [00:35<00:00, 2.24s/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<01:01, 1.60it/s]
2%|2 | 2/100 [00:01<01:00, 1.61it/s]
3%|3 | 3/100 [00:01<01:00, 1.60it/s]
4%|4 | 4/100 [00:02<00:59, 1.60it/s]
5%|5 | 5/100 [00:03<00:58, 1.61it/s]
6%|6 | 6/100 [00:03<00:58, 1.62it/s]
7%|7 | 7/100 [00:04<00:57, 1.62it/s]
8%|8 | 8/100 [00:04<00:56, 1.62it/s]
9%|9 | 9/100 [00:05<00:55, 1.63it/s]
10%|# | 10/100 [00:06<00:55, 1.63it/s]
11%|#1 | 11/100 [00:06<00:54, 1.63it/s]
12%|#2 | 12/100 [00:07<00:54, 1.62it/s]
13%|#3 | 13/100 [00:08<00:54, 1.59it/s]
14%|#4 | 14/100 [00:08<00:53, 1.60it/s]
15%|#5 | 15/100 [00:09<00:52, 1.61it/s]
16%|#6 | 16/100 [00:09<00:52, 1.61it/s]
17%|#7 | 17/100 [00:10<00:51, 1.62it/s]
18%|#8 | 18/100 [00:11<00:50, 1.62it/s]
19%|#9 | 19/100 [00:11<00:49, 1.63it/s]
20%|## | 20/100 [00:12<00:49, 1.63it/s]
21%|##1 | 21/100 [00:12<00:48, 1.62it/s]
22%|##2 | 22/100 [00:13<00:47, 1.63it/s]
23%|##3 | 23/100 [00:14<00:47, 1.63it/s]
24%|##4 | 24/100 [00:14<00:47, 1.62it/s]
25%|##5 | 25/100 [00:15<00:46, 1.61it/s]
26%|##6 | 26/100 [00:16<00:45, 1.61it/s]
27%|##7 | 27/100 [00:16<00:45, 1.61it/s]
28%|##8 | 28/100 [00:17<00:44, 1.62it/s]
29%|##9 | 29/100 [00:17<00:43, 1.62it/s]
30%|### | 30/100 [00:18<00:42, 1.63it/s]
31%|###1 | 31/100 [00:19<00:42, 1.63it/s]
32%|###2 | 32/100 [00:19<00:41, 1.63it/s]
33%|###3 | 33/100 [00:20<00:41, 1.62it/s]
34%|###4 | 34/100 [00:21<00:41, 1.60it/s]
35%|###5 | 35/100 [00:21<00:40, 1.61it/s]
36%|###6 | 36/100 [00:22<00:39, 1.62it/s]
37%|###7 | 37/100 [00:22<00:39, 1.61it/s]
38%|###8 | 38/100 [00:23<00:38, 1.61it/s]
39%|###9 | 39/100 [00:24<00:37, 1.61it/s]
40%|#### | 40/100 [00:24<00:37, 1.62it/s]
41%|####1 | 41/100 [00:25<00:36, 1.62it/s]
42%|####2 | 42/100 [00:26<00:36, 1.58it/s]
43%|####3 | 43/100 [00:26<00:35, 1.60it/s]
44%|####4 | 44/100 [00:27<00:34, 1.61it/s]
45%|####5 | 45/100 [00:27<00:34, 1.62it/s]
46%|####6 | 46/100 [00:28<00:33, 1.63it/s]
47%|####6 | 47/100 [00:29<00:32, 1.63it/s]
48%|####8 | 48/100 [00:29<00:32, 1.62it/s]
49%|####9 | 49/100 [00:30<00:31, 1.62it/s]
50%|##### | 50/100 [00:30<00:30, 1.62it/s]
51%|#####1 | 51/100 [00:31<00:30, 1.62it/s]
52%|#####2 | 52/100 [00:32<00:29, 1.62it/s]
53%|#####3 | 53/100 [00:32<00:28, 1.63it/s]
54%|#####4 | 54/100 [00:33<00:28, 1.63it/s]
55%|#####5 | 55/100 [00:34<00:27, 1.63it/s]
56%|#####6 | 56/100 [00:34<00:27, 1.62it/s]
57%|#####6 | 57/100 [00:35<00:26, 1.63it/s]
58%|#####8 | 58/100 [00:35<00:25, 1.63it/s]
59%|#####8 | 59/100 [00:36<00:25, 1.63it/s]
60%|###### | 60/100 [00:37<00:24, 1.62it/s]
61%|######1 | 61/100 [00:37<00:24, 1.62it/s]
62%|######2 | 62/100 [00:38<00:23, 1.62it/s]
63%|######3 | 63/100 [00:38<00:22, 1.62it/s]
64%|######4 | 64/100 [00:39<00:22, 1.62it/s]
65%|######5 | 65/100 [00:40<00:21, 1.61it/s]
66%|######6 | 66/100 [00:40<00:21, 1.61it/s]
67%|######7 | 67/100 [00:41<00:20, 1.60it/s]
68%|######8 | 68/100 [00:42<00:19, 1.61it/s]
69%|######9 | 69/100 [00:42<00:19, 1.61it/s]
70%|####### | 70/100 [00:43<00:18, 1.61it/s]
71%|#######1 | 71/100 [00:43<00:17, 1.61it/s]
72%|#######2 | 72/100 [00:44<00:17, 1.62it/s]
73%|#######3 | 73/100 [00:45<00:16, 1.62it/s]
74%|#######4 | 74/100 [00:45<00:16, 1.62it/s]
75%|#######5 | 75/100 [00:46<00:15, 1.62it/s]
76%|#######6 | 76/100 [00:47<00:14, 1.61it/s]
77%|#######7 | 77/100 [00:47<00:14, 1.62it/s]
78%|#######8 | 78/100 [00:48<00:13, 1.61it/s]
79%|#######9 | 79/100 [00:48<00:12, 1.63it/s]
80%|######## | 80/100 [00:49<00:12, 1.64it/s]
81%|########1 | 81/100 [00:50<00:11, 1.63it/s]
82%|########2 | 82/100 [00:50<00:11, 1.63it/s]
83%|########2 | 83/100 [00:51<00:10, 1.63it/s]
84%|########4 | 84/100 [00:51<00:09, 1.63it/s]
85%|########5 | 85/100 [00:52<00:09, 1.63it/s]
86%|########6 | 86/100 [00:53<00:08, 1.63it/s]
87%|########7 | 87/100 [00:53<00:07, 1.63it/s]
88%|########8 | 88/100 [00:54<00:07, 1.64it/s]
89%|########9 | 89/100 [00:54<00:06, 1.62it/s]
90%|######### | 90/100 [00:55<00:06, 1.62it/s]
91%|#########1| 91/100 [00:56<00:05, 1.63it/s]
92%|#########2| 92/100 [00:56<00:05, 1.59it/s]
93%|#########3| 93/100 [00:57<00:04, 1.61it/s]
94%|#########3| 94/100 [00:58<00:03, 1.62it/s]
95%|#########5| 95/100 [00:58<00:03, 1.61it/s]
96%|#########6| 96/100 [00:59<00:02, 1.63it/s]
97%|#########7| 97/100 [00:59<00:01, 1.64it/s]
98%|#########8| 98/100 [01:00<00:01, 1.64it/s]
99%|#########9| 99/100 [01:01<00:00, 1.65it/s]
100%|##########| 100/100 [01:01<00:00, 1.65it/s]
100%|##########| 100/100 [01:01<00:00, 1.62it/s]
0%| | 0/228483 [00:00<?, ?it/s]
0%| | 897/228483 [00:00<00:25, 8967.38it/s]
1%| | 1798/228483 [00:00<00:25, 8991.56it/s]
1%|1 | 2698/228483 [00:00<00:25, 8964.22it/s]
2%|1 | 3609/228483 [00:00<00:24, 9020.98it/s]
2%|1 | 4517/228483 [00:00<00:24, 9041.36it/s]
2%|2 | 5424/228483 [00:00<00:24, 9048.82it/s]
3%|2 | 6329/228483 [00:00<00:24, 8896.46it/s]
3%|3 | 7220/228483 [00:00<00:25, 8850.19it/s]
4%|3 | 8116/228483 [00:00<00:24, 8881.47it/s]
4%|3 | 9005/228483 [00:01<00:24, 8863.17it/s]
4%|4 | 9900/228483 [00:01<00:24, 8888.52it/s]
5%|4 | 10796/228483 [00:01<00:24, 8907.89it/s]
5%|5 | 11704/228483 [00:01<00:24, 8957.35it/s]
6%|5 | 12617/228483 [00:01<00:23, 9007.52it/s]
6%|5 | 13524/228483 [00:01<00:23, 9024.13it/s]
6%|6 | 14436/228483 [00:01<00:23, 9051.44it/s]
7%|6 | 15350/228483 [00:01<00:23, 9076.73it/s]
7%|7 | 16266/228483 [00:01<00:23, 9099.17it/s]
8%|7 | 17176/228483 [00:01<00:23, 9059.39it/s]
8%|7 | 18092/228483 [00:02<00:23, 9089.14it/s]
8%|8 | 19002/228483 [00:02<00:23, 9092.27it/s]
9%|8 | 19912/228483 [00:02<00:22, 9070.49it/s]
9%|9 | 20832/228483 [00:02<00:22, 9107.21it/s]
10%|9 | 21743/228483 [00:02<00:22, 9103.32it/s]
10%|9 | 22663/228483 [00:02<00:22, 9131.11it/s]
10%|# | 23577/228483 [00:02<00:22, 9106.28it/s]
11%|# | 24488/228483 [00:02<00:22, 9089.39it/s]
11%|#1 | 25397/228483 [00:02<00:22, 9040.44it/s]
12%|#1 | 26302/228483 [00:02<00:22, 8995.68it/s]
12%|#1 | 27202/228483 [00:03<00:22, 8988.48it/s]
12%|#2 | 28101/228483 [00:03<00:22, 8931.38it/s]
13%|#2 | 29002/228483 [00:03<00:22, 8952.48it/s]
13%|#3 | 29898/228483 [00:03<00:22, 8946.26it/s]
13%|#3 | 30810/228483 [00:03<00:21, 8995.39it/s]
14%|#3 | 31710/228483 [00:03<00:21, 8981.91it/s]
14%|#4 | 32609/228483 [00:03<00:21, 8979.88it/s]
15%|#4 | 33508/228483 [00:03<00:21, 8967.98it/s]
15%|#5 | 34425/228483 [00:03<00:21, 9027.30it/s]
15%|#5 | 35328/228483 [00:03<00:21, 9005.18it/s]
16%|#5 | 36229/228483 [00:04<00:21, 8980.77it/s]
16%|#6 | 37138/228483 [00:04<00:21, 9010.91it/s]
17%|#6 | 38041/228483 [00:04<00:21, 9013.39it/s]
17%|#7 | 38950/228483 [00:04<00:20, 9034.64it/s]
17%|#7 | 39854/228483 [00:04<00:20, 9031.14it/s]
18%|#7 | 40758/228483 [00:04<00:32, 5786.22it/s]
18%|#8 | 41640/228483 [00:04<00:29, 6438.22it/s]
19%|#8 | 42499/228483 [00:04<00:26, 6942.16it/s]
19%|#8 | 43389/228483 [00:05<00:24, 7430.76it/s]
19%|#9 | 44280/228483 [00:05<00:23, 7816.95it/s]
20%|#9 | 45191/228483 [00:05<00:22, 8168.89it/s]
20%|## | 46091/228483 [00:05<00:21, 8401.30it/s]
21%|## | 46991/228483 [00:05<00:21, 8571.72it/s]
21%|## | 47894/228483 [00:05<00:20, 8703.89it/s]
21%|##1 | 48793/228483 [00:05<00:20, 8785.87it/s]
22%|##1 | 49697/228483 [00:05<00:20, 8859.63it/s]
22%|##2 | 50598/228483 [00:05<00:19, 8901.66it/s]
23%|##2 | 51495/228483 [00:05<00:19, 8913.28it/s]
23%|##2 | 52398/228483 [00:06<00:19, 8946.31it/s]
23%|##3 | 53296/228483 [00:06<00:19, 8947.13it/s]
24%|##3 | 54193/228483 [00:06<00:19, 8919.61it/s]
24%|##4 | 55096/228483 [00:06<00:19, 8951.07it/s]
25%|##4 | 55998/228483 [00:06<00:19, 8969.49it/s]
25%|##4 | 56905/228483 [00:06<00:19, 8997.16it/s]
25%|##5 | 57821/228483 [00:06<00:18, 9045.05it/s]
26%|##5 | 58726/228483 [00:06<00:18, 8996.67it/s]
26%|##6 | 59626/228483 [00:06<00:18, 8985.49it/s]
26%|##6 | 60525/228483 [00:06<00:18, 8924.87it/s]
27%|##6 | 61418/228483 [00:07<00:18, 8889.93it/s]
27%|##7 | 62308/228483 [00:07<00:18, 8864.48it/s]
28%|##7 | 63210/228483 [00:07<00:18, 8910.51it/s]
28%|##8 | 64113/228483 [00:07<00:18, 8943.29it/s]
28%|##8 | 65014/228483 [00:07<00:18, 8961.64it/s]
29%|##8 | 65929/228483 [00:07<00:18, 9016.76it/s]
29%|##9 | 66832/228483 [00:07<00:17, 9017.84it/s]
30%|##9 | 67741/228483 [00:07<00:17, 9039.43it/s]
30%|### | 68649/228483 [00:07<00:17, 9051.30it/s]
30%|### | 69555/228483 [00:07<00:17, 8981.86it/s]
31%|### | 70468/228483 [00:08<00:17, 9025.89it/s]
31%|###1 | 71371/228483 [00:08<00:17, 8994.64it/s]
32%|###1 | 72286/228483 [00:08<00:17, 9039.72it/s]
32%|###2 | 73199/228483 [00:08<00:17, 9065.23it/s]
32%|###2 | 74106/228483 [00:08<00:17, 9036.91it/s]
33%|###2 | 75017/228483 [00:08<00:16, 9057.92it/s]
33%|###3 | 75923/228483 [00:08<00:16, 9040.70it/s]
34%|###3 | 76830/228483 [00:08<00:16, 9047.15it/s]
34%|###4 | 77735/228483 [00:08<00:16, 8991.35it/s]
34%|###4 | 78635/228483 [00:08<00:16, 8948.63it/s]
35%|###4 | 79530/228483 [00:09<00:16, 8918.46it/s]
35%|###5 | 80426/228483 [00:09<00:16, 8928.26it/s]
36%|###5 | 81331/228483 [00:09<00:16, 8964.09it/s]
36%|###5 | 82228/228483 [00:09<00:16, 8938.25it/s]
36%|###6 | 83140/228483 [00:09<00:16, 8989.97it/s]
37%|###6 | 84040/228483 [00:09<00:16, 8987.30it/s]
37%|###7 | 84955/228483 [00:09<00:15, 9034.96it/s]
38%|###7 | 85867/228483 [00:09<00:15, 9058.64it/s]
38%|###7 | 86773/228483 [00:09<00:15, 9016.39it/s]
38%|###8 | 87675/228483 [00:09<00:15, 8996.65it/s]
39%|###8 | 88575/228483 [00:10<00:15, 8955.45it/s]
39%|###9 | 89503/228483 [00:10<00:15, 9049.34it/s]
40%|###9 | 90409/228483 [00:10<00:15, 9021.80it/s]
40%|###9 | 91314/228483 [00:10<00:15, 9029.01it/s]
40%|#### | 92217/228483 [00:10<00:24, 5590.83it/s]
41%|#### | 93107/228483 [00:10<00:21, 6281.61it/s]
41%|####1 | 94015/228483 [00:10<00:19, 6924.61it/s]
42%|####1 | 94896/228483 [00:10<00:18, 7389.92it/s]
42%|####1 | 95776/228483 [00:11<00:17, 7756.18it/s]
42%|####2 | 96655/228483 [00:11<00:16, 8035.33it/s]
43%|####2 | 97564/228483 [00:11<00:15, 8327.94it/s]
43%|####3 | 98475/228483 [00:11<00:15, 8548.59it/s]
43%|####3 | 99373/228483 [00:11<00:14, 8672.56it/s]
44%|####3 | 100276/228483 [00:11<00:14, 8775.58it/s]
44%|####4 | 101177/228483 [00:11<00:14, 8844.16it/s]
45%|####4 | 102077/228483 [00:11<00:14, 8890.30it/s]
45%|####5 | 102995/228483 [00:11<00:13, 8975.03it/s]
45%|####5 | 103898/228483 [00:11<00:13, 8977.76it/s]
46%|####5 | 104800/228483 [00:12<00:13, 8972.84it/s]
46%|####6 | 105720/228483 [00:12<00:13, 9037.95it/s]
47%|####6 | 106628/228483 [00:12<00:13, 9047.12it/s]
47%|####7 | 107539/228483 [00:12<00:13, 9064.63it/s]
47%|####7 | 108447/228483 [00:12<00:13, 9024.76it/s]
48%|####7 | 109356/228483 [00:12<00:13, 9041.89it/s]
48%|####8 | 110261/228483 [00:12<00:13, 9039.79it/s]
49%|####8 | 111167/228483 [00:12<00:12, 9045.14it/s]
49%|####9 | 112072/228483 [00:12<00:12, 9005.13it/s]
49%|####9 | 112973/228483 [00:12<00:12, 8918.01it/s]
50%|####9 | 113866/228483 [00:13<00:12, 8918.96it/s]
50%|##### | 114759/228483 [00:13<00:12, 8898.44it/s]
51%|##### | 115680/228483 [00:13<00:12, 8990.11it/s]
51%|#####1 | 116580/228483 [00:13<00:12, 8971.94it/s]
51%|#####1 | 117480/228483 [00:13<00:12, 8978.15it/s]
52%|#####1 | 118389/228483 [00:13<00:12, 9011.29it/s]
52%|#####2 | 119291/228483 [00:13<00:12, 9008.62it/s]
53%|#####2 | 120199/228483 [00:13<00:11, 9027.38it/s]
53%|#####3 | 121102/228483 [00:13<00:11, 8997.97it/s]
53%|#####3 | 122002/228483 [00:13<00:11, 8985.58it/s]
54%|#####3 | 122914/228483 [00:14<00:11, 9024.46it/s]
54%|#####4 | 123835/228483 [00:14<00:11, 9076.61it/s]
55%|#####4 | 124743/228483 [00:14<00:11, 9022.34it/s]
55%|#####4 | 125649/228483 [00:14<00:11, 9030.75it/s]
55%|#####5 | 126553/228483 [00:14<00:11, 9031.54it/s]
56%|#####5 | 127476/228483 [00:14<00:11, 9090.83it/s]
56%|#####6 | 128386/228483 [00:14<00:11, 9088.32it/s]
57%|#####6 | 129295/228483 [00:14<00:10, 9043.36it/s]
57%|#####6 | 130200/228483 [00:14<00:10, 9007.27it/s]
57%|#####7 | 131101/228483 [00:14<00:10, 8979.49it/s]
58%|#####7 | 131999/228483 [00:15<00:10, 8946.92it/s]
58%|#####8 | 132894/228483 [00:15<00:10, 8884.86it/s]
59%|#####8 | 133799/228483 [00:15<00:10, 8933.71it/s]
59%|#####8 | 134702/228483 [00:15<00:10, 8962.33it/s]
59%|#####9 | 135607/228483 [00:15<00:10, 8985.78it/s]
60%|#####9 | 136518/228483 [00:15<00:10, 9020.64it/s]
60%|###### | 137421/228483 [00:15<00:10, 9018.27it/s]
61%|###### | 138323/228483 [00:15<00:09, 9018.57it/s]
61%|###### | 139227/228483 [00:15<00:09, 9023.59it/s]
61%|######1 | 140130/228483 [00:15<00:09, 9000.99it/s]
62%|######1 | 141037/228483 [00:16<00:09, 9019.52it/s]
62%|######2 | 141939/228483 [00:16<00:09, 9009.63it/s]
63%|######2 | 142840/228483 [00:16<00:09, 9000.66it/s]
63%|######2 | 143752/228483 [00:16<00:09, 9035.81it/s]
63%|######3 | 144656/228483 [00:16<00:09, 9032.88it/s]
64%|######3 | 145564/228483 [00:16<00:09, 9045.73it/s]
64%|######4 | 146469/228483 [00:16<00:09, 9005.79it/s]
65%|######4 | 147373/228483 [00:16<00:08, 9013.35it/s]
65%|######4 | 148275/228483 [00:16<00:08, 8972.50it/s]
65%|######5 | 149173/228483 [00:16<00:08, 8913.32it/s]
66%|######5 | 150065/228483 [00:17<00:08, 8912.85it/s]
66%|######6 | 150957/228483 [00:17<00:08, 8868.44it/s]
66%|######6 | 151861/228483 [00:17<00:08, 8919.03it/s]
67%|######6 | 152769/228483 [00:17<00:08, 8963.41it/s]
67%|######7 | 153675/228483 [00:17<00:08, 8989.42it/s]
68%|######7 | 154580/228483 [00:17<00:08, 9005.35it/s]
68%|######8 | 155481/228483 [00:17<00:08, 8998.29it/s]
68%|######8 | 156381/228483 [00:17<00:08, 8983.73it/s]
69%|######8 | 157280/228483 [00:18<00:13, 5379.26it/s]
69%|######9 | 158177/228483 [00:18<00:11, 6108.60it/s]
70%|######9 | 159047/228483 [00:18<00:10, 6691.92it/s]
70%|####### | 159941/228483 [00:18<00:09, 7239.26it/s]
70%|####### | 160843/228483 [00:18<00:08, 7698.26it/s]
71%|####### | 161745/228483 [00:18<00:08, 8053.83it/s]
71%|#######1 | 162647/228483 [00:18<00:07, 8319.77it/s]
72%|#######1 | 163542/228483 [00:18<00:07, 8494.69it/s]
72%|#######1 | 164429/228483 [00:18<00:07, 8602.95it/s]
72%|#######2 | 165312/228483 [00:19<00:07, 8650.19it/s]
73%|#######2 | 166200/228483 [00:19<00:07, 8716.55it/s]
73%|#######3 | 167087/228483 [00:19<00:07, 8760.87it/s]
74%|#######3 | 167994/228483 [00:19<00:06, 8850.82it/s]
74%|#######3 | 168885/228483 [00:19<00:06, 8835.70it/s]
74%|#######4 | 169785/228483 [00:19<00:06, 8882.78it/s]
75%|#######4 | 170686/228483 [00:19<00:06, 8917.96it/s]
75%|#######5 | 171580/228483 [00:19<00:06, 8923.07it/s]
75%|#######5 | 172487/228483 [00:19<00:06, 8965.54it/s]
76%|#######5 | 173385/228483 [00:19<00:06, 8947.02it/s]
76%|#######6 | 174285/228483 [00:20<00:06, 8961.53it/s]
77%|#######6 | 175182/228483 [00:20<00:05, 8956.31it/s]
77%|#######7 | 176080/228483 [00:20<00:05, 8961.74it/s]
77%|#######7 | 176977/228483 [00:20<00:05, 8953.97it/s]
78%|#######7 | 177873/228483 [00:20<00:05, 8947.55it/s]
78%|#######8 | 178782/228483 [00:20<00:05, 8988.16it/s]
79%|#######8 | 179681/228483 [00:20<00:05, 8985.11it/s]
79%|#######9 | 180580/228483 [00:20<00:05, 8950.34it/s]
79%|#######9 | 181476/228483 [00:20<00:05, 8943.21it/s]
80%|#######9 | 182371/228483 [00:20<00:05, 8908.41it/s]
80%|######## | 183262/228483 [00:21<00:05, 8892.31it/s]
81%|######## | 184152/228483 [00:21<00:05, 8859.27it/s]
81%|######## | 185038/228483 [00:21<00:04, 8842.85it/s]
81%|########1 | 185932/228483 [00:21<00:04, 8871.11it/s]
82%|########1 | 186825/228483 [00:21<00:04, 8886.86it/s]
82%|########2 | 187724/228483 [00:21<00:04, 8914.88it/s]
83%|########2 | 188623/228483 [00:21<00:04, 8936.22it/s]
83%|########2 | 189520/228483 [00:21<00:04, 8945.00it/s]
83%|########3 | 190421/228483 [00:21<00:04, 8962.34it/s]
84%|########3 | 191318/228483 [00:21<00:04, 8958.12it/s]
84%|########4 | 192214/228483 [00:22<00:04, 8949.27it/s]
85%|########4 | 193109/228483 [00:22<00:03, 8939.30it/s]
85%|########4 | 194003/228483 [00:22<00:03, 8917.43it/s]
85%|########5 | 194901/228483 [00:22<00:03, 8933.94it/s]
86%|########5 | 195802/228483 [00:22<00:03, 8955.56it/s]
86%|########6 | 196701/228483 [00:22<00:03, 8965.65it/s]
86%|########6 | 197609/228483 [00:22<00:03, 8996.73it/s]
87%|########6 | 198509/228483 [00:22<00:03, 8960.14it/s]
87%|########7 | 199406/228483 [00:22<00:03, 8956.15it/s]
88%|########7 | 200302/228483 [00:22<00:03, 8907.47it/s]
88%|########8 | 201193/228483 [00:23<00:03, 8895.31it/s]
88%|########8 | 202083/228483 [00:23<00:02, 8873.33it/s]
89%|########8 | 202971/228483 [00:23<00:02, 8843.72it/s]
89%|########9 | 203869/228483 [00:23<00:02, 8884.20it/s]
90%|########9 | 204760/228483 [00:23<00:02, 8889.06it/s]
90%|######### | 205658/228483 [00:23<00:02, 8914.61it/s]
90%|######### | 206557/228483 [00:23<00:02, 8936.29it/s]
91%|######### | 207451/228483 [00:23<00:02, 8936.35it/s]
91%|#########1| 208356/228483 [00:23<00:02, 8970.39it/s]
92%|#########1| 209254/228483 [00:23<00:02, 8909.34it/s]
92%|#########1| 210150/228483 [00:24<00:02, 8924.24it/s]
92%|#########2| 211043/228483 [00:24<00:01, 8909.20it/s]
93%|#########2| 211934/228483 [00:24<00:01, 8888.26it/s]
93%|#########3| 212833/228483 [00:24<00:01, 8917.09it/s]
94%|#########3| 213733/228483 [00:24<00:01, 8940.74it/s]
94%|#########3| 214639/228483 [00:24<00:01, 8975.56it/s]
94%|#########4| 215537/228483 [00:24<00:01, 8956.66it/s]
95%|#########4| 216433/228483 [00:24<00:01, 8950.70it/s]
95%|#########5| 217329/228483 [00:24<00:01, 8894.41it/s]
96%|#########5| 218219/228483 [00:24<00:01, 8838.02it/s]
96%|#########5| 219103/228483 [00:25<00:01, 8813.31it/s]
96%|#########6| 219989/228483 [00:25<00:00, 8825.79it/s]
97%|#########6| 220872/228483 [00:25<00:00, 8816.05it/s]
97%|#########7| 221774/228483 [00:25<00:00, 8876.13it/s]
97%|#########7| 222663/228483 [00:25<00:00, 8880.00it/s]
98%|#########7| 223564/228483 [00:25<00:00, 8917.32it/s]
98%|#########8| 224475/228483 [00:25<00:00, 8973.98it/s]
99%|#########8| 225373/228483 [00:25<00:00, 8852.88it/s]
99%|#########9| 226259/228483 [00:25<00:00, 8671.12it/s]
99%|#########9| 227146/228483 [00:25<00:00, 8728.74it/s]
100%|#########9| 228046/228483 [00:26<00:00, 8806.92it/s]
100%|##########| 228483/228483 [00:33<00:00, 6839.24it/s]
<nilearn.plotting.displays._slicers.ZSlicer object at 0x7f7bd9fcd2e0>
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 0x7f7be43b8d60>
References
Total running time of the script: ( 4 minutes 9.083 seconds)