Note
Click here to download the full example code
2 samples permutation test on source data with spatio-temporal clustering¶
Tests if the source space data are significantly different between 2 groups of subjects (simulated here using one subject’s data). The multiple comparisons problem is addressed with a cluster-level permutation test across space and time.
# Authors: Alexandre Gramfort <alexandre.gramfort@inria.fr>
# Eric Larson <larson.eric.d@gmail.com>
# License: BSD (3-clause)
import os.path as op
import numpy as np
from scipy import stats as stats
import mne
from mne import spatial_src_adjacency
from mne.stats import spatio_temporal_cluster_test, summarize_clusters_stc
from mne.datasets import sample
print(__doc__)
Set parameters¶
data_path = sample.data_path()
stc_fname = data_path + '/MEG/sample/sample_audvis-meg-lh.stc'
subjects_dir = data_path + '/subjects'
src_fname = subjects_dir + '/fsaverage/bem/fsaverage-ico-5-src.fif'
# Load stc to in common cortical space (fsaverage)
stc = mne.read_source_estimate(stc_fname)
stc.resample(50, npad='auto')
# Read the source space we are morphing to
src = mne.read_source_spaces(src_fname)
fsave_vertices = [s['vertno'] for s in src]
morph = mne.compute_source_morph(stc, 'sample', 'fsaverage',
spacing=fsave_vertices, smooth=20,
subjects_dir=subjects_dir)
stc = morph.apply(stc)
n_vertices_fsave, n_times = stc.data.shape
tstep = stc.tstep * 1000 # convert to milliseconds
n_subjects1, n_subjects2 = 7, 9
print('Simulating data for %d and %d subjects.' % (n_subjects1, n_subjects2))
# Let's make sure our results replicate, so set the seed.
np.random.seed(0)
X1 = np.random.randn(n_vertices_fsave, n_times, n_subjects1) * 10
X2 = np.random.randn(n_vertices_fsave, n_times, n_subjects2) * 10
X1[:, :, :] += stc.data[:, :, np.newaxis]
# make the activity bigger for the second set of subjects
X2[:, :, :] += 3 * stc.data[:, :, np.newaxis]
# We want to compare the overall activity levels for each subject
X1 = np.abs(X1) # only magnitude
X2 = np.abs(X2) # only magnitude
Out:
Reading a source space...
[done]
Reading a source space...
[done]
2 source spaces read
Simulating data for 7 and 9 subjects.
Compute statistic¶
To use an algorithm optimized for spatio-temporal clustering, we just pass the spatial adjacency matrix (instead of spatio-temporal)
print('Computing adjacency.')
adjacency = spatial_src_adjacency(src)
# Note that X needs to be a list of multi-dimensional array of shape
# samples (subjects_k) x time x space, so we permute dimensions
X1 = np.transpose(X1, [2, 1, 0])
X2 = np.transpose(X2, [2, 1, 0])
X = [X1, X2]
# Now let's actually do the clustering. This can take a long time...
# Here we set the threshold quite high to reduce computation.
p_threshold = 0.0001
f_threshold = stats.distributions.f.ppf(1. - p_threshold / 2.,
n_subjects1 - 1, n_subjects2 - 1)
print('Clustering.')
T_obs, clusters, cluster_p_values, H0 = clu =\
spatio_temporal_cluster_test(X, adjacency=adjacency, n_jobs=1,
threshold=f_threshold, buffer_size=None)
# Now select the clusters that are sig. at p < 0.05 (note that this value
# is multiple-comparisons corrected).
good_cluster_inds = np.where(cluster_p_values < 0.05)[0]
Out:
Computing adjacency.
-- number of adjacent vertices : 20484
Clustering.
stat_fun(H1): min=0.000000 max=293.441888
Running initial clustering
Found 234 clusters
Permuting 1023 times...
0%| | : 0/1023 [00:00<?, ?it/s]
0%| | : 1/1023 [00:00<01:09, 14.77it/s]
0%| | : 2/1023 [00:00<00:51, 19.85it/s]
0%| | : 3/1023 [00:00<00:45, 22.42it/s]
0%| | : 4/1023 [00:00<00:42, 23.98it/s]
1%| | : 6/1023 [00:00<00:33, 30.54it/s]
1%| | : 7/1023 [00:00<00:33, 30.36it/s]
1%| | : 8/1023 [00:00<00:33, 30.25it/s]
1%| | : 9/1023 [00:00<00:33, 30.15it/s]
1%| | : 10/1023 [00:00<00:33, 30.06it/s]
1%|1 | : 11/1023 [00:00<00:33, 30.01it/s]
1%|1 | : 13/1023 [00:00<00:30, 33.16it/s]
1%|1 | : 14/1023 [00:00<00:30, 32.79it/s]
1%|1 | : 15/1023 [00:00<00:31, 32.47it/s]
2%|1 | : 16/1023 [00:00<00:31, 32.19it/s]
2%|1 | : 18/1023 [00:00<00:29, 34.58it/s]
2%|1 | : 19/1023 [00:00<00:29, 34.14it/s]
2%|1 | : 20/1023 [00:00<00:29, 33.76it/s]
2%|2 | : 21/1023 [00:00<00:29, 33.42it/s]
2%|2 | : 23/1023 [00:00<00:28, 35.42it/s]
2%|2 | : 24/1023 [00:00<00:28, 34.98it/s]
2%|2 | : 25/1023 [00:00<00:28, 34.58it/s]
3%|2 | : 26/1023 [00:00<00:29, 34.21it/s]
3%|2 | : 27/1023 [00:00<00:29, 33.89it/s]
3%|2 | : 29/1023 [00:00<00:27, 35.63it/s]
3%|2 | : 30/1023 [00:00<00:28, 35.22it/s]
3%|3 | : 31/1023 [00:00<00:28, 34.84it/s]
3%|3 | : 32/1023 [00:00<00:28, 34.49it/s]
3%|3 | : 34/1023 [00:00<00:27, 36.08it/s]
3%|3 | : 35/1023 [00:01<00:27, 35.66it/s]
4%|3 | : 36/1023 [00:01<00:27, 35.28it/s]
4%|3 | : 37/1023 [00:01<00:28, 34.93it/s]
4%|3 | : 39/1023 [00:01<00:27, 36.41it/s]
4%|3 | : 40/1023 [00:01<00:27, 35.99it/s]
4%|4 | : 41/1023 [00:01<00:27, 35.60it/s]
4%|4 | : 42/1023 [00:01<00:27, 35.23it/s]
4%|4 | : 44/1023 [00:01<00:26, 36.63it/s]
4%|4 | : 45/1023 [00:01<00:27, 36.21it/s]
4%|4 | : 46/1023 [00:01<00:27, 35.82it/s]
5%|4 | : 47/1023 [00:01<00:27, 35.46it/s]
5%|4 | : 49/1023 [00:01<00:26, 36.79it/s]
5%|4 | : 50/1023 [00:01<00:26, 36.38it/s]
5%|4 | : 51/1023 [00:01<00:27, 35.99it/s]
5%|5 | : 52/1023 [00:01<00:27, 35.62it/s]
5%|5 | : 54/1023 [00:01<00:26, 36.91it/s]
5%|5 | : 55/1023 [00:01<00:26, 36.50it/s]
5%|5 | : 56/1023 [00:01<00:26, 36.11it/s]
6%|5 | : 57/1023 [00:01<00:27, 35.76it/s]
6%|5 | : 58/1023 [00:01<00:27, 35.42it/s]
6%|5 | : 60/1023 [00:01<00:26, 36.70it/s]
6%|5 | : 61/1023 [00:01<00:26, 36.32it/s]
6%|6 | : 62/1023 [00:01<00:26, 35.95it/s]
6%|6 | : 63/1023 [00:01<00:26, 35.61it/s]
6%|6 | : 65/1023 [00:01<00:25, 36.86it/s]
6%|6 | : 66/1023 [00:01<00:26, 36.47it/s]
7%|6 | : 67/1023 [00:01<00:26, 36.11it/s]
7%|6 | : 68/1023 [00:01<00:26, 35.76it/s]
7%|6 | : 69/1023 [00:01<00:26, 35.43it/s]
7%|6 | : 71/1023 [00:01<00:25, 36.68it/s]
7%|7 | : 72/1023 [00:02<00:26, 36.31it/s]
7%|7 | : 73/1023 [00:02<00:26, 35.96it/s]
7%|7 | : 74/1023 [00:02<00:26, 35.62it/s]
7%|7 | : 75/1023 [00:02<00:26, 35.30it/s]
7%|7 | : 76/1023 [00:02<00:27, 35.00it/s]
8%|7 | : 77/1023 [00:02<00:27, 34.72it/s]
8%|7 | : 78/1023 [00:02<00:27, 34.45it/s]
8%|7 | : 79/1023 [00:02<00:27, 34.20it/s]
8%|7 | : 81/1023 [00:02<00:26, 35.48it/s]
8%|8 | : 82/1023 [00:02<00:26, 35.18it/s]
8%|8 | : 83/1023 [00:02<00:26, 34.89it/s]
8%|8 | : 84/1023 [00:02<00:27, 34.61it/s]
8%|8 | : 85/1023 [00:02<00:27, 34.35it/s]
9%|8 | : 87/1023 [00:02<00:26, 35.62it/s]
9%|8 | : 88/1023 [00:02<00:26, 35.31it/s]
9%|8 | : 89/1023 [00:02<00:26, 35.02it/s]
9%|8 | : 90/1023 [00:02<00:26, 34.74it/s]
9%|8 | : 92/1023 [00:02<00:25, 35.98it/s]
9%|9 | : 93/1023 [00:02<00:26, 35.66it/s]
9%|9 | : 94/1023 [00:02<00:26, 35.35it/s]
9%|9 | : 95/1023 [00:02<00:26, 35.05it/s]
9%|9 | : 96/1023 [00:02<00:26, 34.77it/s]
10%|9 | : 98/1023 [00:02<00:25, 35.99it/s]
10%|9 | : 99/1023 [00:02<00:25, 35.67it/s]
10%|9 | : 100/1023 [00:02<00:26, 35.35it/s]
10%|9 | : 101/1023 [00:02<00:26, 35.06it/s]
10%|# | : 103/1023 [00:02<00:25, 36.27it/s]
10%|# | : 104/1023 [00:02<00:25, 35.94it/s]
10%|# | : 105/1023 [00:02<00:25, 35.61it/s]
10%|# | : 106/1023 [00:03<00:25, 35.30it/s]
11%|# | : 108/1023 [00:03<00:25, 36.51it/s]
11%|# | : 109/1023 [00:03<00:25, 36.16it/s]
11%|# | : 110/1023 [00:03<00:25, 35.82it/s]
11%|# | : 111/1023 [00:03<00:25, 35.51it/s]
11%|#1 | : 113/1023 [00:03<00:24, 36.69it/s]
11%|#1 | : 114/1023 [00:03<00:25, 36.34it/s]
11%|#1 | : 115/1023 [00:03<00:25, 35.99it/s]
11%|#1 | : 116/1023 [00:03<00:25, 35.67it/s]
11%|#1 | : 117/1023 [00:03<00:25, 35.35it/s]
12%|#1 | : 119/1023 [00:03<00:24, 36.54it/s]
12%|#1 | : 120/1023 [00:03<00:24, 36.19it/s]
12%|#1 | : 121/1023 [00:03<00:25, 35.85it/s]
12%|#1 | : 122/1023 [00:03<00:25, 35.53it/s]
12%|#2 | : 123/1023 [00:03<00:25, 35.22it/s]
12%|#2 | : 125/1023 [00:03<00:24, 36.42it/s]
12%|#2 | : 126/1023 [00:03<00:24, 36.07it/s]
12%|#2 | : 127/1023 [00:03<00:25, 35.74it/s]
13%|#2 | : 128/1023 [00:03<00:25, 35.43it/s]
13%|#2 | : 129/1023 [00:03<00:25, 35.14it/s]
13%|#2 | : 131/1023 [00:03<00:24, 36.34it/s]
13%|#2 | : 132/1023 [00:03<00:24, 35.99it/s]
13%|#3 | : 133/1023 [00:03<00:24, 35.66it/s]
13%|#3 | : 134/1023 [00:03<00:25, 35.36it/s]
13%|#3 | : 136/1023 [00:03<00:24, 36.54it/s]
13%|#3 | : 137/1023 [00:03<00:24, 36.19it/s]
13%|#3 | : 138/1023 [00:03<00:24, 35.86it/s]
14%|#3 | : 139/1023 [00:03<00:24, 35.54it/s]
14%|#3 | : 140/1023 [00:03<00:25, 35.24it/s]
14%|#3 | : 142/1023 [00:03<00:24, 36.44it/s]
14%|#3 | : 143/1023 [00:04<00:24, 36.10it/s]
14%|#4 | : 144/1023 [00:04<00:24, 35.77it/s]
14%|#4 | : 145/1023 [00:04<00:24, 35.46it/s]
14%|#4 | : 147/1023 [00:04<00:23, 36.64it/s]
14%|#4 | : 148/1023 [00:04<00:24, 36.29it/s]
15%|#4 | : 149/1023 [00:04<00:24, 35.95it/s]
15%|#4 | : 150/1023 [00:04<00:24, 35.63it/s]
15%|#4 | : 151/1023 [00:04<00:24, 35.32it/s]
15%|#4 | : 153/1023 [00:04<00:23, 36.52it/s]
15%|#5 | : 154/1023 [00:04<00:24, 36.17it/s]
15%|#5 | : 155/1023 [00:04<00:24, 35.84it/s]
15%|#5 | : 156/1023 [00:04<00:24, 35.52it/s]
15%|#5 | : 157/1023 [00:04<00:24, 35.22it/s]
15%|#5 | : 158/1023 [00:04<00:24, 34.94it/s]
16%|#5 | : 160/1023 [00:04<00:23, 36.15it/s]
16%|#5 | : 161/1023 [00:04<00:24, 35.82it/s]
16%|#5 | : 162/1023 [00:04<00:24, 35.50it/s]
16%|#5 | : 163/1023 [00:04<00:24, 35.21it/s]
16%|#6 | : 164/1023 [00:04<00:24, 34.92it/s]
16%|#6 | : 166/1023 [00:04<00:23, 36.13it/s]
16%|#6 | : 167/1023 [00:04<00:23, 35.80it/s]
16%|#6 | : 168/1023 [00:04<00:24, 35.49it/s]
17%|#6 | : 169/1023 [00:04<00:24, 35.19it/s]
17%|#6 | : 170/1023 [00:04<00:24, 34.91it/s]
17%|#6 | : 172/1023 [00:04<00:23, 36.12it/s]
17%|#6 | : 173/1023 [00:04<00:23, 35.79it/s]
17%|#7 | : 174/1023 [00:04<00:23, 35.48it/s]
17%|#7 | : 175/1023 [00:04<00:24, 35.18it/s]
17%|#7 | : 176/1023 [00:04<00:24, 34.90it/s]
17%|#7 | : 178/1023 [00:05<00:23, 36.11it/s]
17%|#7 | : 179/1023 [00:05<00:23, 35.78it/s]
18%|#7 | : 180/1023 [00:05<00:23, 35.47it/s]
18%|#7 | : 181/1023 [00:05<00:23, 35.17it/s]
18%|#7 | : 182/1023 [00:05<00:24, 34.89it/s]
18%|#7 | : 184/1023 [00:05<00:23, 36.10it/s]
18%|#8 | : 185/1023 [00:05<00:23, 35.78it/s]
18%|#8 | : 186/1023 [00:05<00:23, 35.46it/s]
18%|#8 | : 187/1023 [00:05<00:23, 35.16it/s]
18%|#8 | : 188/1023 [00:05<00:23, 34.88it/s]
19%|#8 | : 190/1023 [00:05<00:23, 36.10it/s]
19%|#8 | : 191/1023 [00:05<00:23, 35.77it/s]
19%|#8 | : 192/1023 [00:05<00:23, 35.46it/s]
19%|#8 | : 193/1023 [00:05<00:23, 35.16it/s]
19%|#9 | : 195/1023 [00:05<00:22, 36.35it/s]
19%|#9 | : 196/1023 [00:05<00:22, 36.01it/s]
19%|#9 | : 197/1023 [00:05<00:23, 35.68it/s]
19%|#9 | : 198/1023 [00:05<00:23, 35.38it/s]
19%|#9 | : 199/1023 [00:05<00:23, 35.09it/s]
20%|#9 | : 201/1023 [00:05<00:22, 36.29it/s]
20%|#9 | : 202/1023 [00:05<00:22, 35.95it/s]
20%|#9 | : 203/1023 [00:05<00:23, 35.63it/s]
20%|#9 | : 204/1023 [00:05<00:23, 35.30it/s]
20%|## | : 205/1023 [00:05<00:23, 35.01it/s]
20%|## | : 207/1023 [00:05<00:22, 36.22it/s]
20%|## | : 208/1023 [00:05<00:22, 35.89it/s]
20%|## | : 209/1023 [00:05<00:22, 35.57it/s]
21%|## | : 210/1023 [00:05<00:23, 35.27it/s]
21%|## | : 211/1023 [00:05<00:23, 34.99it/s]
21%|## | : 213/1023 [00:05<00:22, 36.19it/s]
21%|## | : 214/1023 [00:06<00:22, 35.86it/s]
21%|##1 | : 215/1023 [00:06<00:22, 35.55it/s]
21%|##1 | : 216/1023 [00:06<00:22, 35.25it/s]
21%|##1 | : 217/1023 [00:06<00:23, 34.96it/s]
21%|##1 | : 219/1023 [00:06<00:22, 36.16it/s]
22%|##1 | : 220/1023 [00:06<00:22, 35.84it/s]
22%|##1 | : 221/1023 [00:06<00:22, 35.52it/s]
22%|##1 | : 222/1023 [00:06<00:22, 35.22it/s]
22%|##1 | : 223/1023 [00:06<00:22, 34.94it/s]
22%|##1 | : 225/1023 [00:06<00:22, 36.15it/s]
22%|##2 | : 226/1023 [00:06<00:22, 35.81it/s]
22%|##2 | : 227/1023 [00:06<00:22, 35.50it/s]
22%|##2 | : 228/1023 [00:06<00:22, 35.20it/s]
22%|##2 | : 229/1023 [00:06<00:22, 34.92it/s]
23%|##2 | : 231/1023 [00:06<00:21, 36.13it/s]
23%|##2 | : 232/1023 [00:06<00:22, 35.80it/s]
23%|##2 | : 233/1023 [00:06<00:22, 35.49it/s]
23%|##2 | : 234/1023 [00:06<00:22, 35.19it/s]
23%|##3 | : 236/1023 [00:06<00:21, 36.39it/s]
23%|##3 | : 237/1023 [00:06<00:21, 36.05it/s]
23%|##3 | : 238/1023 [00:06<00:21, 35.72it/s]
23%|##3 | : 239/1023 [00:06<00:22, 35.42it/s]
23%|##3 | : 240/1023 [00:06<00:22, 35.12it/s]
24%|##3 | : 241/1023 [00:06<00:22, 34.84it/s]
24%|##3 | : 242/1023 [00:06<00:22, 34.58it/s]
24%|##3 | : 243/1023 [00:06<00:22, 34.33it/s]
24%|##3 | : 245/1023 [00:06<00:21, 35.57it/s]
24%|##4 | : 246/1023 [00:06<00:22, 35.27it/s]
24%|##4 | : 247/1023 [00:06<00:22, 34.98it/s]
24%|##4 | : 248/1023 [00:07<00:22, 34.71it/s]
24%|##4 | : 250/1023 [00:07<00:21, 35.93it/s]
25%|##4 | : 251/1023 [00:07<00:21, 35.61it/s]
25%|##4 | : 252/1023 [00:07<00:21, 35.30it/s]
25%|##4 | : 253/1023 [00:07<00:21, 35.01it/s]
25%|##4 | : 254/1023 [00:07<00:22, 34.73it/s]
25%|##5 | : 256/1023 [00:07<00:21, 35.94it/s]
25%|##5 | : 257/1023 [00:07<00:21, 35.61it/s]
25%|##5 | : 258/1023 [00:07<00:21, 35.30it/s]
25%|##5 | : 259/1023 [00:07<00:21, 35.01it/s]
25%|##5 | : 260/1023 [00:07<00:21, 34.73it/s]
26%|##5 | : 262/1023 [00:07<00:21, 35.94it/s]
26%|##5 | : 263/1023 [00:07<00:21, 35.61it/s]
26%|##5 | : 264/1023 [00:07<00:21, 35.30it/s]
26%|##5 | : 265/1023 [00:07<00:21, 35.01it/s]
26%|##6 | : 267/1023 [00:07<00:20, 36.20it/s]
26%|##6 | : 268/1023 [00:07<00:21, 35.86it/s]
26%|##6 | : 269/1023 [00:07<00:21, 35.54it/s]
26%|##6 | : 270/1023 [00:07<00:21, 35.24it/s]
27%|##6 | : 272/1023 [00:07<00:20, 36.42it/s]
27%|##6 | : 273/1023 [00:07<00:20, 36.07it/s]
27%|##6 | : 274/1023 [00:07<00:20, 35.74it/s]
27%|##6 | : 275/1023 [00:07<00:21, 35.42it/s]
27%|##6 | : 276/1023 [00:07<00:21, 35.12it/s]
27%|##7 | : 278/1023 [00:07<00:20, 36.31it/s]
27%|##7 | : 279/1023 [00:07<00:20, 35.97it/s]
27%|##7 | : 280/1023 [00:07<00:20, 35.65it/s]
27%|##7 | : 281/1023 [00:07<00:20, 35.35it/s]
28%|##7 | : 283/1023 [00:07<00:20, 36.53it/s]
28%|##7 | : 284/1023 [00:07<00:20, 36.19it/s]
28%|##7 | : 285/1023 [00:08<00:20, 35.85it/s]
28%|##7 | : 286/1023 [00:08<00:20, 35.52it/s]
28%|##8 | : 287/1023 [00:08<00:20, 35.22it/s]
28%|##8 | : 289/1023 [00:08<00:20, 36.40it/s]
28%|##8 | : 290/1023 [00:08<00:20, 36.06it/s]
28%|##8 | : 291/1023 [00:08<00:20, 35.73it/s]
29%|##8 | : 292/1023 [00:08<00:20, 35.42it/s]
29%|##8 | : 294/1023 [00:08<00:19, 36.60it/s]
29%|##8 | : 295/1023 [00:08<00:20, 36.24it/s]
29%|##8 | : 296/1023 [00:08<00:20, 35.90it/s]
29%|##9 | : 297/1023 [00:08<00:20, 35.58it/s]
29%|##9 | : 298/1023 [00:08<00:20, 35.29it/s]
29%|##9 | : 300/1023 [00:08<00:19, 36.48it/s]
29%|##9 | : 301/1023 [00:08<00:19, 36.13it/s]
30%|##9 | : 302/1023 [00:08<00:20, 35.80it/s]
30%|##9 | : 303/1023 [00:08<00:20, 35.49it/s]
30%|##9 | : 304/1023 [00:08<00:20, 35.19it/s]
30%|##9 | : 306/1023 [00:08<00:19, 36.38it/s]
30%|### | : 307/1023 [00:08<00:19, 36.04it/s]
30%|### | : 308/1023 [00:08<00:20, 35.72it/s]
30%|### | : 309/1023 [00:08<00:20, 35.41it/s]
30%|### | : 311/1023 [00:08<00:19, 36.60it/s]
30%|### | : 312/1023 [00:08<00:19, 36.24it/s]
31%|### | : 313/1023 [00:08<00:19, 35.91it/s]
31%|### | : 314/1023 [00:08<00:19, 35.59it/s]
31%|### | : 315/1023 [00:08<00:20, 35.28it/s]
31%|### | : 316/1023 [00:08<00:20, 34.99it/s]
31%|### | : 317/1023 [00:08<00:20, 34.72it/s]
31%|###1 | : 318/1023 [00:08<00:20, 34.46it/s]
31%|###1 | : 319/1023 [00:09<00:20, 34.21it/s]
31%|###1 | : 320/1023 [00:09<00:20, 33.98it/s]
31%|###1 | : 321/1023 [00:09<00:20, 33.76it/s]
31%|###1 | : 322/1023 [00:09<00:20, 33.55it/s]
32%|###1 | : 323/1023 [00:09<00:20, 33.34it/s]
32%|###1 | : 324/1023 [00:09<00:21, 33.14it/s]
32%|###1 | : 325/1023 [00:09<00:21, 32.96it/s]
32%|###1 | : 326/1023 [00:09<00:21, 32.78it/s]
32%|###1 | : 327/1023 [00:09<00:21, 32.62it/s]
32%|###2 | : 328/1023 [00:09<00:21, 32.47it/s]
32%|###2 | : 329/1023 [00:09<00:21, 32.32it/s]
32%|###2 | : 330/1023 [00:09<00:21, 32.18it/s]
32%|###2 | : 331/1023 [00:09<00:21, 32.04it/s]
32%|###2 | : 332/1023 [00:09<00:21, 31.91it/s]
33%|###2 | : 333/1023 [00:09<00:21, 31.78it/s]
33%|###2 | : 334/1023 [00:09<00:21, 31.67it/s]
33%|###2 | : 335/1023 [00:09<00:21, 31.55it/s]
33%|###2 | : 336/1023 [00:09<00:21, 31.45it/s]
33%|###2 | : 337/1023 [00:09<00:21, 31.35it/s]
33%|###3 | : 338/1023 [00:09<00:21, 31.26it/s]
33%|###3 | : 339/1023 [00:09<00:21, 31.17it/s]
33%|###3 | : 340/1023 [00:09<00:21, 31.08it/s]
33%|###3 | : 341/1023 [00:09<00:22, 31.00it/s]
33%|###3 | : 342/1023 [00:09<00:22, 30.92it/s]
34%|###3 | : 343/1023 [00:09<00:22, 30.85it/s]
34%|###3 | : 344/1023 [00:09<00:22, 30.78it/s]
34%|###3 | : 345/1023 [00:09<00:22, 30.72it/s]
34%|###3 | : 346/1023 [00:09<00:22, 30.65it/s]
34%|###3 | : 347/1023 [00:09<00:22, 30.60it/s]
34%|###4 | : 348/1023 [00:09<00:22, 30.54it/s]
34%|###4 | : 349/1023 [00:10<00:22, 30.50it/s]
34%|###4 | : 350/1023 [00:10<00:22, 30.45it/s]
34%|###4 | : 351/1023 [00:10<00:22, 30.40it/s]
35%|###4 | : 353/1023 [00:10<00:21, 31.83it/s]
35%|###4 | : 354/1023 [00:10<00:21, 31.72it/s]
35%|###4 | : 355/1023 [00:10<00:21, 31.61it/s]
35%|###4 | : 356/1023 [00:10<00:21, 31.50it/s]
35%|###4 | : 357/1023 [00:10<00:21, 31.40it/s]
35%|###5 | : 359/1023 [00:10<00:20, 32.78it/s]
35%|###5 | : 360/1023 [00:10<00:20, 32.63it/s]
35%|###5 | : 361/1023 [00:10<00:20, 32.47it/s]
35%|###5 | : 362/1023 [00:10<00:20, 32.33it/s]
36%|###5 | : 364/1023 [00:10<00:19, 33.66it/s]
36%|###5 | : 365/1023 [00:10<00:19, 33.46it/s]
36%|###5 | : 366/1023 [00:10<00:19, 33.27it/s]
36%|###5 | : 367/1023 [00:10<00:19, 33.08it/s]
36%|###5 | : 368/1023 [00:10<00:19, 32.91it/s]
36%|###6 | : 370/1023 [00:10<00:19, 34.21it/s]
36%|###6 | : 371/1023 [00:10<00:19, 33.98it/s]
36%|###6 | : 372/1023 [00:10<00:19, 33.76it/s]
36%|###6 | : 373/1023 [00:10<00:19, 33.55it/s]
37%|###6 | : 375/1023 [00:10<00:18, 34.83it/s]
37%|###6 | : 376/1023 [00:10<00:18, 34.57it/s]
37%|###6 | : 377/1023 [00:10<00:18, 34.32it/s]
37%|###6 | : 378/1023 [00:10<00:18, 34.07it/s]
37%|###7 | : 380/1023 [00:10<00:18, 35.32it/s]
37%|###7 | : 381/1023 [00:10<00:18, 35.04it/s]
37%|###7 | : 382/1023 [00:10<00:18, 34.76it/s]
37%|###7 | : 383/1023 [00:10<00:18, 34.50it/s]
38%|###7 | : 384/1023 [00:11<00:18, 34.25it/s]
38%|###7 | : 386/1023 [00:11<00:17, 35.49it/s]
38%|###7 | : 387/1023 [00:11<00:18, 35.20it/s]
38%|###7 | : 388/1023 [00:11<00:18, 34.91it/s]
38%|###8 | : 389/1023 [00:11<00:18, 34.64it/s]
38%|###8 | : 391/1023 [00:11<00:17, 35.86it/s]
38%|###8 | : 392/1023 [00:11<00:17, 35.54it/s]
38%|###8 | : 393/1023 [00:11<00:17, 35.25it/s]
39%|###8 | : 394/1023 [00:11<00:17, 34.96it/s]
39%|###8 | : 396/1023 [00:11<00:17, 36.17it/s]
39%|###8 | : 397/1023 [00:11<00:17, 35.84it/s]
39%|###8 | : 398/1023 [00:11<00:17, 35.53it/s]
39%|###9 | : 399/1023 [00:11<00:17, 35.23it/s]
39%|###9 | : 401/1023 [00:11<00:17, 36.42it/s]
39%|###9 | : 402/1023 [00:11<00:17, 36.08it/s]
39%|###9 | : 403/1023 [00:11<00:17, 35.74it/s]
39%|###9 | : 404/1023 [00:11<00:17, 35.43it/s]
40%|###9 | : 405/1023 [00:11<00:17, 35.13it/s]
40%|###9 | : 407/1023 [00:11<00:16, 36.33it/s]
40%|###9 | : 408/1023 [00:11<00:17, 35.99it/s]
40%|###9 | : 409/1023 [00:11<00:17, 35.66it/s]
40%|#### | : 410/1023 [00:11<00:17, 35.35it/s]
40%|#### | : 412/1023 [00:11<00:16, 36.53it/s]
40%|#### | : 413/1023 [00:11<00:16, 36.18it/s]
40%|#### | : 414/1023 [00:11<00:16, 35.84it/s]
41%|#### | : 415/1023 [00:11<00:17, 35.52it/s]
41%|#### | : 416/1023 [00:11<00:17, 35.22it/s]
41%|#### | : 418/1023 [00:11<00:16, 36.41it/s]
41%|#### | : 419/1023 [00:11<00:16, 36.07it/s]
41%|####1 | : 420/1023 [00:11<00:16, 35.74it/s]
41%|####1 | : 421/1023 [00:12<00:16, 35.43it/s]
41%|####1 | : 423/1023 [00:12<00:16, 36.61it/s]
41%|####1 | : 424/1023 [00:12<00:16, 36.26it/s]
42%|####1 | : 425/1023 [00:12<00:16, 35.92it/s]
42%|####1 | : 426/1023 [00:12<00:16, 35.60it/s]
42%|####1 | : 427/1023 [00:12<00:16, 35.30it/s]
42%|####1 | : 429/1023 [00:12<00:16, 36.49it/s]
42%|####2 | : 430/1023 [00:12<00:16, 36.14it/s]
42%|####2 | : 431/1023 [00:12<00:16, 35.81it/s]
42%|####2 | : 432/1023 [00:12<00:16, 35.50it/s]
42%|####2 | : 433/1023 [00:12<00:16, 35.20it/s]
43%|####2 | : 435/1023 [00:12<00:16, 36.39it/s]
43%|####2 | : 436/1023 [00:12<00:16, 36.05it/s]
43%|####2 | : 437/1023 [00:12<00:16, 35.73it/s]
43%|####2 | : 438/1023 [00:12<00:16, 35.42it/s]
43%|####3 | : 440/1023 [00:12<00:15, 36.60it/s]
43%|####3 | : 441/1023 [00:12<00:16, 36.25it/s]
43%|####3 | : 442/1023 [00:12<00:16, 35.92it/s]
43%|####3 | : 443/1023 [00:12<00:16, 35.60it/s]
43%|####3 | : 444/1023 [00:12<00:16, 35.30it/s]
43%|####3 | : 445/1023 [00:12<00:16, 35.01it/s]
44%|####3 | : 447/1023 [00:12<00:15, 36.21it/s]
44%|####3 | : 448/1023 [00:12<00:16, 35.88it/s]
44%|####3 | : 449/1023 [00:12<00:16, 35.56it/s]
44%|####3 | : 450/1023 [00:12<00:16, 35.25it/s]
44%|####4 | : 451/1023 [00:12<00:16, 34.96it/s]
44%|####4 | : 453/1023 [00:12<00:15, 36.16it/s]
44%|####4 | : 454/1023 [00:12<00:15, 35.82it/s]
44%|####4 | : 455/1023 [00:12<00:15, 35.50it/s]
45%|####4 | : 456/1023 [00:13<00:16, 35.20it/s]
45%|####4 | : 457/1023 [00:13<00:16, 34.91it/s]
45%|####4 | : 458/1023 [00:13<00:16, 34.64it/s]
45%|####4 | : 460/1023 [00:13<00:15, 35.86it/s]
45%|####5 | : 461/1023 [00:13<00:15, 35.54it/s]
45%|####5 | : 462/1023 [00:13<00:15, 35.23it/s]
45%|####5 | : 463/1023 [00:13<00:16, 34.94it/s]
45%|####5 | : 464/1023 [00:13<00:16, 34.66it/s]
46%|####5 | : 466/1023 [00:13<00:15, 35.87it/s]
46%|####5 | : 467/1023 [00:13<00:15, 35.55it/s]
46%|####5 | : 468/1023 [00:13<00:15, 35.24it/s]
46%|####5 | : 469/1023 [00:13<00:15, 34.96it/s]
46%|####6 | : 471/1023 [00:13<00:15, 36.15it/s]
46%|####6 | : 472/1023 [00:13<00:15, 35.82it/s]
46%|####6 | : 473/1023 [00:13<00:15, 35.50it/s]
46%|####6 | : 474/1023 [00:13<00:15, 35.19it/s]
46%|####6 | : 475/1023 [00:13<00:15, 34.90it/s]
47%|####6 | : 476/1023 [00:13<00:15, 34.63it/s]
47%|####6 | : 477/1023 [00:13<00:15, 34.38it/s]
47%|####6 | : 478/1023 [00:13<00:15, 34.14it/s]
47%|####6 | : 479/1023 [00:13<00:16, 33.91it/s]
47%|####6 | : 480/1023 [00:13<00:16, 33.69it/s]
47%|####7 | : 481/1023 [00:13<00:16, 33.48it/s]
47%|####7 | : 483/1023 [00:13<00:15, 34.76it/s]
47%|####7 | : 484/1023 [00:13<00:15, 34.50it/s]
47%|####7 | : 485/1023 [00:13<00:15, 34.25it/s]
48%|####7 | : 486/1023 [00:13<00:15, 34.02it/s]
48%|####7 | : 487/1023 [00:13<00:15, 33.80it/s]
48%|####7 | : 489/1023 [00:13<00:15, 35.06it/s]
48%|####7 | : 490/1023 [00:13<00:15, 34.78it/s]
48%|####7 | : 491/1023 [00:14<00:15, 34.52it/s]
48%|####8 | : 492/1023 [00:14<00:15, 34.27it/s]
48%|####8 | : 494/1023 [00:14<00:14, 35.51it/s]
48%|####8 | : 495/1023 [00:14<00:14, 35.22it/s]
48%|####8 | : 496/1023 [00:14<00:15, 34.93it/s]
49%|####8 | : 497/1023 [00:14<00:15, 34.66it/s]
49%|####8 | : 499/1023 [00:14<00:14, 35.88it/s]
49%|####8 | : 500/1023 [00:14<00:14, 35.56it/s]
49%|####8 | : 501/1023 [00:14<00:14, 35.26it/s]
49%|####9 | : 502/1023 [00:14<00:14, 34.97it/s]
49%|####9 | : 503/1023 [00:14<00:14, 34.70it/s]
49%|####9 | : 504/1023 [00:14<00:15, 34.44it/s]
49%|####9 | : 506/1023 [00:14<00:14, 35.67it/s]
50%|####9 | : 507/1023 [00:14<00:14, 35.37it/s]
50%|####9 | : 508/1023 [00:14<00:14, 35.07it/s]
50%|####9 | : 509/1023 [00:14<00:14, 34.79it/s]
50%|####9 | : 510/1023 [00:14<00:14, 34.53it/s]
50%|##### | : 512/1023 [00:14<00:14, 35.76it/s]
50%|##### | : 513/1023 [00:14<00:14, 35.45it/s]
50%|##### | : 514/1023 [00:14<00:14, 35.15it/s]
50%|##### | : 515/1023 [00:14<00:14, 34.87it/s]
50%|##### | : 516/1023 [00:14<00:14, 34.61it/s]
51%|##### | : 517/1023 [00:14<00:14, 34.35it/s]
51%|##### | : 519/1023 [00:14<00:14, 35.59it/s]
51%|##### | : 520/1023 [00:14<00:14, 35.29it/s]
51%|##### | : 521/1023 [00:14<00:14, 35.00it/s]
51%|#####1 | : 522/1023 [00:14<00:14, 34.73it/s]
51%|#####1 | : 523/1023 [00:14<00:14, 34.47it/s]
51%|#####1 | : 525/1023 [00:14<00:13, 35.70it/s]
51%|#####1 | : 526/1023 [00:15<00:14, 35.39it/s]
52%|#####1 | : 527/1023 [00:15<00:14, 35.10it/s]
52%|#####1 | : 528/1023 [00:15<00:14, 34.82it/s]
52%|#####1 | : 529/1023 [00:15<00:14, 34.56it/s]
52%|#####1 | : 531/1023 [00:15<00:13, 35.79it/s]
52%|#####2 | : 532/1023 [00:15<00:13, 35.48it/s]
52%|#####2 | : 533/1023 [00:15<00:13, 35.18it/s]
52%|#####2 | : 534/1023 [00:15<00:14, 34.90it/s]
52%|#####2 | : 535/1023 [00:15<00:14, 34.63it/s]
52%|#####2 | : 537/1023 [00:15<00:13, 35.86it/s]
53%|#####2 | : 538/1023 [00:15<00:13, 35.54it/s]
53%|#####2 | : 539/1023 [00:15<00:13, 35.24it/s]
53%|#####2 | : 540/1023 [00:15<00:13, 34.96it/s]
53%|#####2 | : 541/1023 [00:15<00:13, 34.69it/s]
53%|#####3 | : 543/1023 [00:15<00:13, 35.90it/s]
53%|#####3 | : 544/1023 [00:15<00:13, 35.59it/s]
53%|#####3 | : 545/1023 [00:15<00:13, 35.29it/s]
53%|#####3 | : 546/1023 [00:15<00:13, 35.00it/s]
54%|#####3 | : 548/1023 [00:15<00:13, 36.21it/s]
54%|#####3 | : 549/1023 [00:15<00:13, 35.87it/s]
54%|#####3 | : 550/1023 [00:15<00:13, 35.56it/s]
54%|#####3 | : 551/1023 [00:15<00:13, 35.26it/s]
54%|#####3 | : 552/1023 [00:15<00:13, 34.98it/s]
54%|#####4 | : 553/1023 [00:15<00:13, 34.71it/s]
54%|#####4 | : 555/1023 [00:15<00:13, 35.93it/s]
54%|#####4 | : 556/1023 [00:15<00:13, 35.61it/s]
54%|#####4 | : 557/1023 [00:15<00:13, 35.31it/s]
55%|#####4 | : 558/1023 [00:15<00:13, 35.03it/s]
55%|#####4 | : 559/1023 [00:15<00:13, 34.75it/s]
55%|#####4 | : 561/1023 [00:15<00:12, 35.97it/s]
55%|#####4 | : 562/1023 [00:16<00:12, 35.65it/s]
55%|#####5 | : 563/1023 [00:16<00:13, 35.34it/s]
55%|#####5 | : 564/1023 [00:16<00:13, 35.05it/s]
55%|#####5 | : 565/1023 [00:16<00:13, 34.78it/s]
55%|#####5 | : 567/1023 [00:16<00:12, 35.99it/s]
56%|#####5 | : 568/1023 [00:16<00:12, 35.67it/s]
56%|#####5 | : 569/1023 [00:16<00:12, 35.37it/s]
56%|#####5 | : 570/1023 [00:16<00:12, 35.07it/s]
56%|#####5 | : 571/1023 [00:16<00:12, 34.79it/s]
56%|#####6 | : 573/1023 [00:16<00:12, 36.01it/s]
56%|#####6 | : 574/1023 [00:16<00:12, 35.69it/s]
56%|#####6 | : 575/1023 [00:16<00:12, 35.38it/s]
56%|#####6 | : 576/1023 [00:16<00:12, 35.09it/s]
56%|#####6 | : 577/1023 [00:16<00:12, 34.81it/s]
57%|#####6 | : 578/1023 [00:16<00:12, 34.55it/s]
57%|#####6 | : 580/1023 [00:16<00:12, 35.78it/s]
57%|#####6 | : 581/1023 [00:16<00:12, 35.46it/s]
57%|#####6 | : 582/1023 [00:16<00:12, 35.17it/s]
57%|#####6 | : 583/1023 [00:16<00:12, 34.89it/s]
57%|#####7 | : 585/1023 [00:16<00:12, 36.10it/s]
57%|#####7 | : 586/1023 [00:16<00:12, 35.76it/s]
57%|#####7 | : 587/1023 [00:16<00:12, 35.45it/s]
57%|#####7 | : 588/1023 [00:16<00:12, 35.16it/s]
58%|#####7 | : 589/1023 [00:16<00:12, 34.88it/s]
58%|#####7 | : 591/1023 [00:16<00:11, 36.09it/s]
58%|#####7 | : 592/1023 [00:16<00:12, 35.77it/s]
58%|#####7 | : 593/1023 [00:16<00:12, 35.46it/s]
58%|#####8 | : 594/1023 [00:16<00:12, 35.16it/s]
58%|#####8 | : 596/1023 [00:16<00:11, 36.36it/s]
58%|#####8 | : 597/1023 [00:17<00:11, 36.02it/s]
58%|#####8 | : 598/1023 [00:17<00:11, 35.69it/s]
59%|#####8 | : 599/1023 [00:17<00:11, 35.38it/s]
59%|#####8 | : 600/1023 [00:17<00:12, 35.09it/s]
59%|#####8 | : 602/1023 [00:17<00:11, 36.29it/s]
59%|#####8 | : 603/1023 [00:17<00:11, 35.96it/s]
59%|#####9 | : 604/1023 [00:17<00:11, 35.64it/s]
59%|#####9 | : 605/1023 [00:17<00:11, 35.33it/s]
59%|#####9 | : 606/1023 [00:17<00:11, 35.04it/s]
59%|#####9 | : 608/1023 [00:17<00:11, 36.25it/s]
60%|#####9 | : 609/1023 [00:17<00:11, 35.91it/s]
60%|#####9 | : 610/1023 [00:17<00:11, 35.59it/s]
60%|#####9 | : 611/1023 [00:17<00:11, 35.29it/s]
60%|#####9 | : 612/1023 [00:17<00:11, 35.01it/s]
60%|#####9 | : 613/1023 [00:17<00:11, 34.73it/s]
60%|###### | : 614/1023 [00:17<00:11, 34.47it/s]
60%|###### | : 615/1023 [00:17<00:11, 34.23it/s]
60%|###### | : 616/1023 [00:17<00:11, 34.00it/s]
60%|###### | : 617/1023 [00:17<00:12, 33.77it/s]
61%|###### | : 619/1023 [00:17<00:11, 35.03it/s]
61%|###### | : 620/1023 [00:17<00:11, 34.75it/s]
61%|###### | : 621/1023 [00:17<00:11, 34.50it/s]
61%|###### | : 622/1023 [00:17<00:11, 34.25it/s]
61%|###### | : 623/1023 [00:17<00:11, 34.01it/s]
61%|######1 | : 625/1023 [00:17<00:11, 35.27it/s]
61%|######1 | : 626/1023 [00:17<00:11, 34.99it/s]
61%|######1 | : 627/1023 [00:17<00:11, 34.71it/s]
61%|######1 | : 628/1023 [00:17<00:11, 34.45it/s]
61%|######1 | : 629/1023 [00:17<00:11, 34.20it/s]
62%|######1 | : 630/1023 [00:17<00:11, 33.97it/s]
62%|######1 | : 632/1023 [00:18<00:11, 35.23it/s]
62%|######1 | : 633/1023 [00:18<00:11, 34.94it/s]
62%|######1 | : 634/1023 [00:18<00:11, 34.67it/s]
62%|######2 | : 635/1023 [00:18<00:11, 34.41it/s]
62%|######2 | : 636/1023 [00:18<00:11, 34.17it/s]
62%|######2 | : 637/1023 [00:18<00:11, 33.94it/s]
62%|######2 | : 638/1023 [00:18<00:11, 33.71it/s]
63%|######2 | : 640/1023 [00:18<00:10, 34.98it/s]
63%|######2 | : 641/1023 [00:18<00:11, 34.71it/s]
63%|######2 | : 642/1023 [00:18<00:11, 34.44it/s]
63%|######2 | : 643/1023 [00:18<00:11, 34.19it/s]
63%|######2 | : 644/1023 [00:18<00:11, 33.96it/s]
63%|######3 | : 646/1023 [00:18<00:10, 35.21it/s]
63%|######3 | : 647/1023 [00:18<00:10, 34.93it/s]
63%|######3 | : 648/1023 [00:18<00:10, 34.66it/s]
63%|######3 | : 649/1023 [00:18<00:10, 34.41it/s]
64%|######3 | : 650/1023 [00:18<00:10, 34.17it/s]
64%|######3 | : 651/1023 [00:18<00:10, 33.94it/s]
64%|######3 | : 653/1023 [00:18<00:10, 35.20it/s]
64%|######3 | : 654/1023 [00:18<00:10, 34.92it/s]
64%|######4 | : 655/1023 [00:18<00:10, 34.65it/s]
64%|######4 | : 656/1023 [00:18<00:10, 34.39it/s]
64%|######4 | : 657/1023 [00:18<00:10, 34.14it/s]
64%|######4 | : 659/1023 [00:18<00:10, 35.39it/s]
65%|######4 | : 660/1023 [00:18<00:10, 35.10it/s]
65%|######4 | : 661/1023 [00:18<00:10, 34.82it/s]
65%|######4 | : 662/1023 [00:18<00:10, 34.55it/s]
65%|######4 | : 663/1023 [00:18<00:10, 34.31it/s]
65%|######4 | : 664/1023 [00:18<00:10, 34.07it/s]
65%|######5 | : 665/1023 [00:19<00:10, 33.84it/s]
65%|######5 | : 667/1023 [00:19<00:10, 35.10it/s]
65%|######5 | : 668/1023 [00:19<00:10, 34.82it/s]
65%|######5 | : 669/1023 [00:19<00:10, 34.56it/s]
65%|######5 | : 670/1023 [00:19<00:10, 34.30it/s]
66%|######5 | : 671/1023 [00:19<00:10, 34.06it/s]
66%|######5 | : 672/1023 [00:19<00:10, 33.84it/s]
66%|######5 | : 673/1023 [00:19<00:10, 33.62it/s]
66%|######5 | : 675/1023 [00:19<00:09, 34.90it/s]
66%|######6 | : 676/1023 [00:19<00:10, 34.63it/s]
66%|######6 | : 677/1023 [00:19<00:10, 34.37it/s]
66%|######6 | : 678/1023 [00:19<00:10, 34.13it/s]
66%|######6 | : 679/1023 [00:19<00:10, 33.89it/s]
67%|######6 | : 681/1023 [00:19<00:09, 35.15it/s]
67%|######6 | : 682/1023 [00:19<00:09, 34.87it/s]
67%|######6 | : 683/1023 [00:19<00:09, 34.61it/s]
67%|######6 | : 684/1023 [00:19<00:09, 34.35it/s]
67%|######6 | : 685/1023 [00:19<00:09, 34.11it/s]
67%|######7 | : 687/1023 [00:19<00:09, 35.35it/s]
67%|######7 | : 688/1023 [00:19<00:09, 35.06it/s]
67%|######7 | : 689/1023 [00:19<00:09, 34.78it/s]
67%|######7 | : 690/1023 [00:19<00:09, 34.52it/s]
68%|######7 | : 691/1023 [00:19<00:09, 34.27it/s]
68%|######7 | : 693/1023 [00:19<00:09, 35.51it/s]
68%|######7 | : 694/1023 [00:19<00:09, 35.22it/s]
68%|######7 | : 695/1023 [00:19<00:09, 34.93it/s]
68%|######8 | : 696/1023 [00:19<00:09, 34.66it/s]
68%|######8 | : 697/1023 [00:19<00:09, 34.41it/s]
68%|######8 | : 698/1023 [00:19<00:09, 34.16it/s]
68%|######8 | : 700/1023 [00:19<00:09, 35.41it/s]
69%|######8 | : 701/1023 [00:20<00:09, 35.11it/s]
69%|######8 | : 702/1023 [00:20<00:09, 34.84it/s]
69%|######8 | : 703/1023 [00:20<00:09, 34.57it/s]
69%|######8 | : 704/1023 [00:20<00:09, 34.31it/s]
69%|######8 | : 705/1023 [00:20<00:09, 34.07it/s]
69%|######9 | : 706/1023 [00:20<00:09, 33.85it/s]
69%|######9 | : 708/1023 [00:20<00:08, 35.11it/s]
69%|######9 | : 709/1023 [00:20<00:09, 34.83it/s]
69%|######9 | : 710/1023 [00:20<00:09, 34.57it/s]
70%|######9 | : 711/1023 [00:20<00:09, 34.32it/s]
70%|######9 | : 713/1023 [00:20<00:08, 35.56it/s]
70%|######9 | : 714/1023 [00:20<00:08, 35.26it/s]
70%|######9 | : 715/1023 [00:20<00:08, 34.97it/s]
70%|######9 | : 716/1023 [00:20<00:08, 34.70it/s]
70%|####### | : 717/1023 [00:20<00:08, 34.45it/s]
70%|####### | : 718/1023 [00:20<00:08, 34.20it/s]
70%|####### | : 720/1023 [00:20<00:08, 35.45it/s]
70%|####### | : 721/1023 [00:20<00:08, 35.15it/s]
71%|####### | : 722/1023 [00:20<00:08, 34.87it/s]
71%|####### | : 723/1023 [00:20<00:08, 34.61it/s]
71%|####### | : 724/1023 [00:20<00:08, 34.35it/s]
71%|####### | : 725/1023 [00:20<00:08, 34.11it/s]
71%|####### | : 726/1023 [00:20<00:08, 33.88it/s]
71%|#######1 | : 728/1023 [00:20<00:08, 35.14it/s]
71%|#######1 | : 729/1023 [00:20<00:08, 34.86it/s]
71%|#######1 | : 730/1023 [00:20<00:08, 34.59it/s]
71%|#######1 | : 731/1023 [00:20<00:08, 34.34it/s]
72%|#######1 | : 732/1023 [00:20<00:08, 34.09it/s]
72%|#######1 | : 733/1023 [00:20<00:08, 33.86it/s]
72%|#######1 | : 734/1023 [00:21<00:08, 33.63it/s]
72%|#######1 | : 735/1023 [00:21<00:08, 33.42it/s]
72%|#######1 | : 736/1023 [00:21<00:08, 33.22it/s]
72%|#######2 | : 737/1023 [00:21<00:08, 33.04it/s]
72%|#######2 | : 739/1023 [00:21<00:08, 34.34it/s]
72%|#######2 | : 740/1023 [00:21<00:08, 34.07it/s]
72%|#######2 | : 741/1023 [00:21<00:08, 33.84it/s]
73%|#######2 | : 742/1023 [00:21<00:08, 33.63it/s]
73%|#######2 | : 743/1023 [00:21<00:08, 33.42it/s]
73%|#######2 | : 744/1023 [00:21<00:08, 33.23it/s]
73%|#######2 | : 745/1023 [00:21<00:08, 33.04it/s]
73%|#######2 | : 746/1023 [00:21<00:08, 32.86it/s]
73%|#######3 | : 747/1023 [00:21<00:08, 32.69it/s]
73%|#######3 | : 749/1023 [00:21<00:08, 34.01it/s]
73%|#######3 | : 750/1023 [00:21<00:08, 33.79it/s]
73%|#######3 | : 751/1023 [00:21<00:08, 33.58it/s]
74%|#######3 | : 752/1023 [00:21<00:08, 33.37it/s]
74%|#######3 | : 753/1023 [00:21<00:08, 33.18it/s]
74%|#######3 | : 754/1023 [00:21<00:08, 33.00it/s]
74%|#######3 | : 755/1023 [00:21<00:08, 32.82it/s]
74%|#######3 | : 757/1023 [00:21<00:07, 34.14it/s]
74%|#######4 | : 758/1023 [00:21<00:07, 33.90it/s]
74%|#######4 | : 759/1023 [00:21<00:07, 33.68it/s]
74%|#######4 | : 760/1023 [00:21<00:07, 33.48it/s]
74%|#######4 | : 761/1023 [00:21<00:07, 33.28it/s]
74%|#######4 | : 762/1023 [00:21<00:07, 33.09it/s]
75%|#######4 | : 763/1023 [00:21<00:07, 32.92it/s]
75%|#######4 | : 764/1023 [00:21<00:07, 32.75it/s]
75%|#######4 | : 765/1023 [00:21<00:07, 32.59it/s]
75%|#######4 | : 767/1023 [00:21<00:07, 33.92it/s]
75%|#######5 | : 768/1023 [00:22<00:07, 33.69it/s]
75%|#######5 | : 769/1023 [00:22<00:07, 33.49it/s]
75%|#######5 | : 770/1023 [00:22<00:07, 33.29it/s]
75%|#######5 | : 771/1023 [00:22<00:07, 33.10it/s]
76%|#######5 | : 773/1023 [00:22<00:07, 34.40it/s]
76%|#######5 | : 774/1023 [00:22<00:07, 34.15it/s]
76%|#######5 | : 775/1023 [00:22<00:07, 33.92it/s]
76%|#######5 | : 776/1023 [00:22<00:07, 33.70it/s]
76%|#######5 | : 777/1023 [00:22<00:07, 33.49it/s]
76%|#######6 | : 778/1023 [00:22<00:07, 33.29it/s]
76%|#######6 | : 779/1023 [00:22<00:07, 33.10it/s]
76%|#######6 | : 780/1023 [00:22<00:07, 32.91it/s]
76%|#######6 | : 781/1023 [00:22<00:07, 32.75it/s]
77%|#######6 | : 783/1023 [00:22<00:07, 34.06it/s]
77%|#######6 | : 784/1023 [00:22<00:07, 33.84it/s]
77%|#######6 | : 785/1023 [00:22<00:07, 33.62it/s]
77%|#######6 | : 786/1023 [00:22<00:07, 33.42it/s]
77%|#######7 | : 788/1023 [00:22<00:06, 34.70it/s]
77%|#######7 | : 789/1023 [00:22<00:06, 34.45it/s]
77%|#######7 | : 790/1023 [00:22<00:06, 34.21it/s]
77%|#######7 | : 791/1023 [00:22<00:06, 33.97it/s]
77%|#######7 | : 792/1023 [00:22<00:06, 33.75it/s]
78%|#######7 | : 793/1023 [00:22<00:06, 33.54it/s]
78%|#######7 | : 794/1023 [00:22<00:06, 33.34it/s]
78%|#######7 | : 795/1023 [00:22<00:06, 33.15it/s]
78%|#######7 | : 797/1023 [00:22<00:06, 34.45it/s]
78%|#######8 | : 798/1023 [00:22<00:06, 34.21it/s]
78%|#######8 | : 799/1023 [00:22<00:06, 33.97it/s]
78%|#######8 | : 800/1023 [00:22<00:06, 33.75it/s]
78%|#######8 | : 801/1023 [00:23<00:06, 33.54it/s]
78%|#######8 | : 803/1023 [00:23<00:06, 34.82it/s]
79%|#######8 | : 804/1023 [00:23<00:06, 34.56it/s]
79%|#######8 | : 805/1023 [00:23<00:06, 34.31it/s]
79%|#######8 | : 806/1023 [00:23<00:06, 34.07it/s]
79%|#######8 | : 807/1023 [00:23<00:06, 33.84it/s]
79%|#######8 | : 808/1023 [00:23<00:06, 33.63it/s]
79%|#######9 | : 810/1023 [00:23<00:06, 34.90it/s]
79%|#######9 | : 811/1023 [00:23<00:06, 34.63it/s]
79%|#######9 | : 812/1023 [00:23<00:06, 34.38it/s]
79%|#######9 | : 813/1023 [00:23<00:06, 34.13it/s]
80%|#######9 | : 814/1023 [00:23<00:06, 33.90it/s]
80%|#######9 | : 815/1023 [00:23<00:06, 33.68it/s]
80%|#######9 | : 817/1023 [00:23<00:05, 34.94it/s]
80%|#######9 | : 818/1023 [00:23<00:05, 34.68it/s]
80%|######## | : 819/1023 [00:23<00:05, 34.42it/s]
80%|######## | : 820/1023 [00:23<00:05, 34.18it/s]
80%|######## | : 821/1023 [00:23<00:05, 33.95it/s]
80%|######## | : 823/1023 [00:23<00:05, 35.20it/s]
81%|######## | : 824/1023 [00:23<00:05, 34.92it/s]
81%|######## | : 825/1023 [00:23<00:05, 34.65it/s]
81%|######## | : 826/1023 [00:23<00:05, 34.40it/s]
81%|######## | : 827/1023 [00:23<00:05, 34.16it/s]
81%|########1 | : 829/1023 [00:23<00:05, 35.41it/s]
81%|########1 | : 830/1023 [00:23<00:05, 35.12it/s]
81%|########1 | : 831/1023 [00:23<00:05, 34.84it/s]
81%|########1 | : 832/1023 [00:23<00:05, 34.57it/s]
81%|########1 | : 833/1023 [00:23<00:05, 34.32it/s]
82%|########1 | : 834/1023 [00:23<00:05, 34.08it/s]
82%|########1 | : 836/1023 [00:23<00:05, 35.25it/s]
82%|########1 | : 837/1023 [00:24<00:05, 34.96it/s]
82%|########1 | : 838/1023 [00:24<00:05, 34.68it/s]
82%|########2 | : 839/1023 [00:24<00:05, 34.42it/s]
82%|########2 | : 840/1023 [00:24<00:05, 34.17it/s]
82%|########2 | : 842/1023 [00:24<00:05, 35.41it/s]
82%|########2 | : 843/1023 [00:24<00:05, 35.11it/s]
83%|########2 | : 844/1023 [00:24<00:05, 34.83it/s]
83%|########2 | : 845/1023 [00:24<00:05, 34.56it/s]
83%|########2 | : 846/1023 [00:24<00:05, 34.30it/s]
83%|########2 | : 848/1023 [00:24<00:04, 35.53it/s]
83%|########2 | : 849/1023 [00:24<00:04, 35.23it/s]
83%|########3 | : 850/1023 [00:24<00:04, 34.94it/s]
83%|########3 | : 851/1023 [00:24<00:04, 34.67it/s]
83%|########3 | : 852/1023 [00:24<00:04, 34.41it/s]
83%|########3 | : 854/1023 [00:24<00:04, 35.63it/s]
84%|########3 | : 855/1023 [00:24<00:04, 35.33it/s]
84%|########3 | : 856/1023 [00:24<00:04, 35.03it/s]
84%|########3 | : 857/1023 [00:24<00:04, 34.75it/s]
84%|########3 | : 858/1023 [00:24<00:04, 34.49it/s]
84%|########3 | : 859/1023 [00:24<00:04, 34.24it/s]
84%|########4 | : 861/1023 [00:24<00:04, 35.48it/s]
84%|########4 | : 862/1023 [00:24<00:04, 35.17it/s]
84%|########4 | : 863/1023 [00:24<00:04, 34.89it/s]
84%|########4 | : 864/1023 [00:24<00:04, 34.62it/s]
85%|########4 | : 865/1023 [00:24<00:04, 34.36it/s]
85%|########4 | : 866/1023 [00:24<00:04, 34.11it/s]
85%|########4 | : 868/1023 [00:24<00:04, 35.36it/s]
85%|########4 | : 869/1023 [00:24<00:04, 35.07it/s]
85%|########5 | : 870/1023 [00:24<00:04, 34.79it/s]
85%|########5 | : 871/1023 [00:25<00:04, 34.53it/s]
85%|########5 | : 872/1023 [00:25<00:04, 34.28it/s]
85%|########5 | : 873/1023 [00:25<00:04, 34.05it/s]
85%|########5 | : 874/1023 [00:25<00:04, 33.82it/s]
86%|########5 | : 876/1023 [00:25<00:04, 35.09it/s]
86%|########5 | : 877/1023 [00:25<00:04, 34.81it/s]
86%|########5 | : 878/1023 [00:25<00:04, 34.55it/s]
86%|########5 | : 879/1023 [00:25<00:04, 34.30it/s]
86%|########6 | : 880/1023 [00:25<00:04, 34.06it/s]
86%|########6 | : 882/1023 [00:25<00:03, 35.31it/s]
86%|########6 | : 883/1023 [00:25<00:03, 35.02it/s]
86%|########6 | : 884/1023 [00:25<00:04, 34.75it/s]
87%|########6 | : 885/1023 [00:25<00:04, 34.49it/s]
87%|########6 | : 886/1023 [00:25<00:04, 34.24it/s]
87%|########6 | : 888/1023 [00:25<00:03, 35.48it/s]
87%|########6 | : 889/1023 [00:25<00:03, 35.19it/s]
87%|########6 | : 890/1023 [00:25<00:03, 34.91it/s]
87%|########7 | : 891/1023 [00:25<00:03, 34.64it/s]
87%|########7 | : 893/1023 [00:25<00:03, 35.86it/s]
87%|########7 | : 894/1023 [00:25<00:03, 35.55it/s]
87%|########7 | : 895/1023 [00:25<00:03, 35.25it/s]
88%|########7 | : 896/1023 [00:25<00:03, 34.96it/s]
88%|########7 | : 897/1023 [00:25<00:03, 34.69it/s]
88%|########7 | : 899/1023 [00:25<00:03, 35.91it/s]
88%|########7 | : 900/1023 [00:25<00:03, 35.59it/s]
88%|########8 | : 901/1023 [00:25<00:03, 35.29it/s]
88%|########8 | : 902/1023 [00:25<00:03, 35.01it/s]
88%|########8 | : 904/1023 [00:25<00:03, 36.21it/s]
88%|########8 | : 905/1023 [00:25<00:03, 35.87it/s]
89%|########8 | : 906/1023 [00:25<00:03, 35.54it/s]
89%|########8 | : 907/1023 [00:26<00:03, 35.24it/s]
89%|########8 | : 909/1023 [00:26<00:03, 36.44it/s]
89%|########8 | : 910/1023 [00:26<00:03, 36.09it/s]
89%|########9 | : 911/1023 [00:26<00:03, 35.76it/s]
89%|########9 | : 912/1023 [00:26<00:03, 35.45it/s]
89%|########9 | : 913/1023 [00:26<00:03, 35.16it/s]
89%|########9 | : 915/1023 [00:26<00:02, 36.36it/s]
90%|########9 | : 916/1023 [00:26<00:02, 36.02it/s]
90%|########9 | : 917/1023 [00:26<00:02, 35.70it/s]
90%|########9 | : 918/1023 [00:26<00:02, 35.39it/s]
90%|########9 | : 919/1023 [00:26<00:02, 35.10it/s]
90%|######### | : 921/1023 [00:26<00:02, 36.30it/s]
90%|######### | : 922/1023 [00:26<00:02, 35.96it/s]
90%|######### | : 923/1023 [00:26<00:02, 35.64it/s]
90%|######### | : 924/1023 [00:26<00:02, 35.33it/s]
91%|######### | : 926/1023 [00:26<00:02, 36.52it/s]
91%|######### | : 927/1023 [00:26<00:02, 36.18it/s]
91%|######### | : 928/1023 [00:26<00:02, 35.84it/s]
91%|######### | : 929/1023 [00:26<00:02, 35.53it/s]
91%|######### | : 930/1023 [00:26<00:02, 35.23it/s]
91%|#########1| : 932/1023 [00:26<00:02, 36.42it/s]
91%|#########1| : 933/1023 [00:26<00:02, 36.05it/s]
91%|#########1| : 934/1023 [00:26<00:02, 35.73it/s]
91%|#########1| : 935/1023 [00:26<00:02, 35.40it/s]
92%|#########1| : 937/1023 [00:26<00:02, 36.58it/s]
92%|#########1| : 938/1023 [00:26<00:02, 36.21it/s]
92%|#########1| : 939/1023 [00:26<00:02, 35.85it/s]
92%|#########1| : 940/1023 [00:26<00:02, 35.53it/s]
92%|#########2| : 942/1023 [00:26<00:02, 36.70it/s]
92%|#########2| : 943/1023 [00:27<00:02, 36.32it/s]
92%|#########2| : 944/1023 [00:27<00:02, 35.96it/s]
92%|#########2| : 945/1023 [00:27<00:02, 35.64it/s]
92%|#########2| : 946/1023 [00:27<00:02, 35.34it/s]
93%|#########2| : 948/1023 [00:27<00:02, 36.53it/s]
93%|#########2| : 949/1023 [00:27<00:02, 36.18it/s]
93%|#########2| : 950/1023 [00:27<00:02, 35.85it/s]
93%|#########2| : 951/1023 [00:27<00:02, 35.53it/s]
93%|#########3| : 952/1023 [00:27<00:02, 35.24it/s]
93%|#########3| : 954/1023 [00:27<00:01, 36.43it/s]
93%|#########3| : 955/1023 [00:27<00:01, 36.08it/s]
93%|#########3| : 956/1023 [00:27<00:01, 35.75it/s]
94%|#########3| : 957/1023 [00:27<00:01, 35.44it/s]
94%|#########3| : 958/1023 [00:27<00:01, 35.14it/s]
94%|#########3| : 960/1023 [00:27<00:01, 36.34it/s]
94%|#########3| : 961/1023 [00:27<00:01, 36.00it/s]
94%|#########4| : 962/1023 [00:27<00:01, 35.68it/s]
94%|#########4| : 963/1023 [00:27<00:01, 35.37it/s]
94%|#########4| : 964/1023 [00:27<00:01, 35.08it/s]
94%|#########4| : 965/1023 [00:27<00:01, 34.80it/s]
94%|#########4| : 966/1023 [00:27<00:01, 34.54it/s]
95%|#########4| : 967/1023 [00:27<00:01, 34.29it/s]
95%|#########4| : 968/1023 [00:27<00:01, 34.05it/s]
95%|#########4| : 970/1023 [00:27<00:01, 35.31it/s]
95%|#########4| : 971/1023 [00:27<00:01, 35.02it/s]
95%|#########5| : 972/1023 [00:27<00:01, 34.75it/s]
95%|#########5| : 973/1023 [00:27<00:01, 34.48it/s]
95%|#########5| : 974/1023 [00:27<00:01, 34.23it/s]
95%|#########5| : 975/1023 [00:27<00:01, 34.00it/s]
96%|#########5| : 977/1023 [00:27<00:01, 35.25it/s]
96%|#########5| : 978/1023 [00:28<00:01, 34.97it/s]
96%|#########5| : 979/1023 [00:28<00:01, 34.70it/s]
96%|#########5| : 980/1023 [00:28<00:01, 34.44it/s]
96%|#########5| : 981/1023 [00:28<00:01, 34.19it/s]
96%|#########6| : 983/1023 [00:28<00:01, 35.43it/s]
96%|#########6| : 984/1023 [00:28<00:01, 35.14it/s]
96%|#########6| : 985/1023 [00:28<00:01, 34.86it/s]
96%|#########6| : 986/1023 [00:28<00:01, 34.59it/s]
96%|#########6| : 987/1023 [00:28<00:01, 34.34it/s]
97%|#########6| : 989/1023 [00:28<00:00, 35.57it/s]
97%|#########6| : 990/1023 [00:28<00:00, 35.27it/s]
97%|#########6| : 991/1023 [00:28<00:00, 34.99it/s]
97%|#########6| : 992/1023 [00:28<00:00, 34.71it/s]
97%|#########7| : 993/1023 [00:28<00:00, 34.45it/s]
97%|#########7| : 994/1023 [00:28<00:00, 34.21it/s]
97%|#########7| : 996/1023 [00:28<00:00, 35.45it/s]
97%|#########7| : 997/1023 [00:28<00:00, 35.15it/s]
98%|#########7| : 998/1023 [00:28<00:00, 34.87it/s]
98%|#########7| : 999/1023 [00:28<00:00, 34.60it/s]
98%|#########7| : 1001/1023 [00:28<00:00, 35.82it/s]
98%|#########7| : 1002/1023 [00:28<00:00, 35.50it/s]
98%|#########8| : 1003/1023 [00:28<00:00, 35.21it/s]
98%|#########8| : 1004/1023 [00:28<00:00, 34.92it/s]
98%|#########8| : 1005/1023 [00:28<00:00, 34.65it/s]
98%|#########8| : 1007/1023 [00:28<00:00, 35.88it/s]
99%|#########8| : 1008/1023 [00:28<00:00, 35.56it/s]
99%|#########8| : 1009/1023 [00:28<00:00, 35.26it/s]
99%|#########8| : 1010/1023 [00:28<00:00, 34.98it/s]
99%|#########8| : 1011/1023 [00:28<00:00, 34.71it/s]
99%|#########9| : 1013/1023 [00:28<00:00, 35.92it/s]
99%|#########9| : 1014/1023 [00:29<00:00, 35.60it/s]
99%|#########9| : 1015/1023 [00:29<00:00, 35.30it/s]
99%|#########9| : 1016/1023 [00:29<00:00, 35.01it/s]
100%|#########9| : 1018/1023 [00:29<00:00, 36.21it/s]
100%|#########9| : 1019/1023 [00:29<00:00, 35.88it/s]
100%|#########9| : 1020/1023 [00:29<00:00, 35.55it/s]
100%|#########9| : 1021/1023 [00:29<00:00, 35.25it/s]
100%|#########9| : 1022/1023 [00:29<00:00, 34.97it/s]
100%|##########| : 1023/1023 [00:29<00:00, 34.94it/s]
Computing cluster p-values
Done.
Visualize the clusters¶
print('Visualizing clusters.')
# Now let's build a convenient representation of each cluster, where each
# cluster becomes a "time point" in the SourceEstimate
fsave_vertices = [np.arange(10242), np.arange(10242)]
stc_all_cluster_vis = summarize_clusters_stc(clu, tstep=tstep,
vertices=fsave_vertices,
subject='fsaverage')
# Let's actually plot the first "time point" in the SourceEstimate, which
# shows all the clusters, weighted by duration
subjects_dir = op.join(data_path, 'subjects')
# blue blobs are for condition A != condition B
brain = stc_all_cluster_vis.plot('fsaverage', hemi='both',
views='lateral', subjects_dir=subjects_dir,
time_label='temporal extent (ms)',
clim=dict(kind='value', lims=[0, 1, 40]))
Out:
Visualizing clusters.
Total running time of the script: ( 0 minutes 46.450 seconds)
Estimated memory usage: 112 MB