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]))
30 cluster ftest spatiotemporal

Out:

Visualizing clusters.

Total running time of the script: ( 0 minutes 46.450 seconds)

Estimated memory usage: 112 MB

Gallery generated by Sphinx-Gallery