{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# The Bootstrap Analysis of Stable Clusters pipeline\n", "This tutorial shows how to run a Bootstrap Analysis of Stable Clusters (BASC) pipeline using NIAK. \n", " * **More documentation**: This tutorial only shows a limited subset of available options. See the [documentation](http://niak.simexp-lab.org/pipe_basc.html) of the pipeline for a more comprehensive list of options. \n", " * **.m script**: The script of the tutorial can be downloaded [here](https://raw.githubusercontent.com/SIMEXP/niak_tutorials/master/niak_tutorial_basc.m) and the notebook is available [here](https://nbviewer.jupyter.org/github/SIMEXP/niak_tutorials/blob/master/niak_tutorial_basc.ipynb).\n", " * **Pre-requisites**: before running this tutorial, you need to complete the [fMRI preprocessing tutorial](https://nbviewer.jupyter.org/github/SIMEXP/niak_tutorials/blob/master/niak_tutorial_fmri_preprocessing.ipynb). \n", " * **Time for completion**: this tutorial will take approximetly 15 minutes to run. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Prepare input files\n", "\n", "## fMRI\n", "\n", "We will re-use the fMRI data that was generated during the fMRI preprocessing tutorial. The following instructions are going to \"grab\" these data. A few options are available to select subjects based on the level of motion, number of available volumes or run name. It is also possible to exclude subjects explicitely, e.g. if they failed visual quality control assessment. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "files_in =\n", "\n", " scalar structure containing the fields:\n", "\n", " data =\n", "\n", " scalar structure containing the fields:\n", "\n", " subject1 =\n", "\n", " scalar structure containing the fields:\n", "\n", " session1: 1x1 scalar struct\n", "\n", " subject2 =\n", "\n", " scalar structure containing the fields:\n", "\n", " session1: 1x1 scalar struct\n", "\n", "\n", " mask = /sandbox/home/git/niak_tutorials/fmri_preprocess/quality_control/group_coregistration/func_mask_group_stereonl.nii.gz\n", " areas = /sandbox/home/git/niak_tutorials/fmri_preprocess/anat/template_aal.nii.gz\n", "\n" ] } ], "source": [ "% Where to find data\n", "clear\n", "path_data = [pwd filesep];\n", "path_preprocess = [path_data 'fmri_preprocess'];\n", "\n", "% That option tells the grabber to prepare files for BASC\n", "opt_g.type_files = 'rest'; \n", "\n", "% Minimum number of volumes per run - typically aim at 2 mns\n", "% Shorther here because the test dataset is really small\n", "opt_g.min_nb_vol = 20;\n", "\n", "% Grab the \"motor\" run\n", "opt_g.filter.run = {'motor'}; % Just grab the \"motor\" runs\n", "\n", "% if uncommented, this option would exclude subject1\n", "% opt_g.exclude_subject = {'subject1'};\n", "\n", "% Grab the data\n", "files_in = niak_grab_fmri_preprocess(path_preprocess,opt_g)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Set the options of the pipeline\n", "Now we set up where to store the results:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% Where to store the results\n", "opt.folder_out = [path_data 'basc_test_niak'];" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This next parameter sets the size at which the region growing process stops (in mm$^3$). This does indirectly set the number of parcels covering the gray matter. A threshold of 1000 mm$^3$ will give about 1000 regions on the grey matter. Here we are going to reduce the dimension even more to speed up the pipeline." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% the size of the regions, when they stop growing. \n", "opt.region_growing.thre_size = 2000; " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The following parameter sets the values that will be explored regarding the number of clusters. A more exhaustive search, up to 500 clusters, could for example use an irregular grid like `[10:10:100 120:20:200 240:40:500]`. Also note that for each number of clusters, a number of solutions around those values, with different numbers being used at the individual and group levels, as well as the final consensus group partitions. " ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% Search for stable clusters in the range 10 to 30\n", "opt.grid_scales = [10:10:30]'; " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The following parameter is used to generate stability maps, and consensus partitions. Although stability is assessed over a wide range of parameters, those maps are only generated for select numbers in order to save disk space. For each set of results, three (integer) parameters actually need to be provided: the number of clusters at the individual level, at the group level and at the final consensus level. Each row will define a new set of results. We could for example use `opt.scales_maps = [10 10 10; 20 20 20];` to generate maps using the same number of clusters (5 and 10) at all three levels. For now we will leave that empty, which means we will not be generating any map. A first run of the pipeline is going to give us insights in the stable cluster solution, and provide data-driven suggestions about what number(s) of clusters to use. " ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% Scale parameters to generate stability maps and consensus clusters\n", "opt.scales_maps = [];" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The following parameters control the number of bootstrap samples used to assess the stability of the clustering both at the individual level, and the group level. Note that with 100 samples, the confidence interval on the stability measures is $\\pm 0.1$, while with 500 samples it reaches $\\pm 0.05$. Here we will use only a few (20 samples) to speed-up the pipeline. " ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% Number of bootstrap samples at the individual level.\n", "opt.stability_tseries.nb_samps = 20;\n", "% Number of bootstrap samples at the group level. \n", "opt.stability_group.nb_samps = 20; " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Also note that there is a minimal number of subjects for the group stability analysis to run. Because this demo has only two subjects, we will need to lower this parameter: " ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": true }, "outputs": [], "source": [ "opt.stability_group.min_subject = 2;" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The final set of flags tell the pipeline which level of the pipeline to run. BASC can be used at the individual level, for each subject independently, or at the group level, to find a parcellation that is the consensus of those generated across all subjects and bootstrap replications of individual fMRI time series. In addition, BASC can generate so-called \"mixed parcellation\", which individual parcellations generated using the group parcellation as an initialization. It is possible to turn on/off the generation of consensus parcellations and stability map for each level of the pipeline (individual, group or mixed): " ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": true }, "outputs": [], "source": [ "% Generate maps/time series at the individual level\n", "opt.flag_ind = false; \n", "% Generate maps/time series at the mixed level (group-level networks mixed with individual stability matrices).\n", "opt.flag_mixed = false; \n", "% Generate maps/time series at the group level\n", "opt.flag_group = true; " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Running the pipeline\n", "\n", "Now it is time to run the pipeline. As explained above, we will not at this stage generate any map or parcellation, but we will use the outputs to select the scales that we will explore in more details. The instruction to run the pipeline is: " ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Logs will be stored in /sandbox/home/git/niak_tutorials/basc_test_niak/logs/\n", "Generating dependencies ...\n", " Percentage completed : 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100- 0.39 sec\n", "Setting up the to-do list ...\n", " I found 133 job(s) to do.\n", "I could not find any log file. This pipeline has not been started (yet?). Press CTRL-C to cancel.\n", "Deamon started on 30-Apr-2017 01:22:29\n", "30-Apr-2017 01:22:29 Starting the pipeline manager...\n", "30-Apr-2017 01:22:30 Starting the garbage collector...\n", "30-Apr-2017 01:22:30 Starting worker number 1...\n", "30-Apr-2017 01:22:30 Starting worker number 2...\n", "30-Apr-2017 01:22:30 Starting worker number 3...\n", "30-Apr-2017 01:22:30 Starting worker number 4...\n", "\n", "Pipeline started on 30-Apr-2017 01:22:32\n", "user: , host: 3f86dc5d68dc, system: unix\n", "****************************************\n", "30-Apr-2017 01:22:34 mask_areas submitted (1 run | 0 fail | 0 done | 132 left)\n", "30-Apr-2017 01:22:34 mask_areas finished (0 run | 0 fail | 1 done | 132 left)\n", "30-Apr-2017 01:22:34 tseries_subject1_session1_motor submitted (1 run | 0 fail | 1 done | 131 left)\n", "30-Apr-2017 01:22:34 tseries_subject2_session1_motor submitted (2 run | 0 fail | 1 done | 130 left)\n", "30-Apr-2017 01:22:34 neighbourhood_areas submitted (3 run | 0 fail | 1 done | 129 left)\n", "30-Apr-2017 01:22:36 tseries_subject1_session1_motor finished (2 run | 0 fail | 2 done | 129 left)\n", "30-Apr-2017 01:22:37 tseries_subject2_session1_motor finished (1 run | 0 fail | 3 done | 129 left)\n", "30-Apr-2017 01:22:37 neighbourhood_areas finished (0 run | 0 fail | 4 done | 129 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2001 submitted (1 run | 0 fail | 4 done | 128 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2002 submitted (2 run | 0 fail | 4 done | 127 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2101 submitted (3 run | 0 fail | 4 done | 126 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2102 submitted (4 run | 0 fail | 4 done | 125 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2111 submitted (5 run | 0 fail | 4 done | 124 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2112 submitted (6 run | 0 fail | 4 done | 123 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2201 submitted (7 run | 0 fail | 4 done | 122 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2202 submitted (8 run | 0 fail | 4 done | 121 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2001 finished (7 run | 0 fail | 5 done | 121 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2002 finished (6 run | 0 fail | 6 done | 121 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2101 finished (5 run | 0 fail | 7 done | 121 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2211 submitted (6 run | 0 fail | 7 done | 120 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2212 submitted (7 run | 0 fail | 7 done | 119 left)\n", "30-Apr-2017 01:22:37 region_growing_area_2301 submitted (8 run | 0 fail | 7 done | 118 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2102 finished (7 run | 0 fail | 8 done | 118 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2111 finished (6 run | 0 fail | 9 done | 118 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2201 finished (5 run | 0 fail | 10 done | 118 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2301 finished (4 run | 0 fail | 11 done | 118 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2302 submitted (5 run | 0 fail | 11 done | 117 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2311 submitted (6 run | 0 fail | 11 done | 116 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2312 submitted (7 run | 0 fail | 11 done | 115 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2321 submitted (8 run | 0 fail | 11 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2202 finished (7 run | 0 fail | 12 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2211 finished (6 run | 0 fail | 13 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2112 finished (5 run | 0 fail | 14 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2212 finished (4 run | 0 fail | 15 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2302 finished (3 run | 0 fail | 16 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2311 finished (2 run | 0 fail | 17 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2312 finished (1 run | 0 fail | 18 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2321 finished (0 run | 0 fail | 19 done | 114 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2322 submitted (1 run | 0 fail | 19 done | 113 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2331 submitted (2 run | 0 fail | 19 done | 112 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2332 submitted (3 run | 0 fail | 19 done | 111 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2401 submitted (4 run | 0 fail | 19 done | 110 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2402 submitted (5 run | 0 fail | 19 done | 109 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2501 submitted (6 run | 0 fail | 19 done | 108 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2502 submitted (7 run | 0 fail | 19 done | 107 left)\n", "30-Apr-2017 01:22:38 region_growing_area_2601 submitted (8 run | 0 fail | 19 done | 106 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2322 finished (7 run | 0 fail | 20 done | 106 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2332 finished (6 run | 0 fail | 21 done | 106 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2401 finished (5 run | 0 fail | 22 done | 106 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2602 submitted (6 run | 0 fail | 22 done | 105 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2611 submitted (7 run | 0 fail | 22 done | 104 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2612 submitted (8 run | 0 fail | 22 done | 103 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2331 finished (7 run | 0 fail | 23 done | 103 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2402 finished (6 run | 0 fail | 24 done | 103 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2601 finished (5 run | 0 fail | 25 done | 103 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2701 submitted (6 run | 0 fail | 25 done | 102 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2702 submitted (7 run | 0 fail | 25 done | 101 left)\n", "30-Apr-2017 01:22:39 region_growing_area_2502 finished (6 run | 0 fail | 26 done | 101 left)\n", "30-Apr-2017 01:22:39 region_growing_area_3001 submitted (7 run | 0 fail | 26 done | 100 left)\n", "30-Apr-2017 01:22:40 region_growing_area_2501 finished (6 run | 0 fail | 27 done | 100 left)\n", "30-Apr-2017 01:22:40 region_growing_area_2602 finished (5 run | 0 fail | 28 done | 100 left)\n", "30-Apr-2017 01:22:40 region_growing_area_2611 finished (4 run | 0 fail | 29 done | 100 left)\n", "30-Apr-2017 01:22:40 region_growing_area_2612 finished (3 run | 0 fail | 30 done | 100 left)\n", "30-Apr-2017 01:22:40 region_growing_area_3002 submitted (4 run | 0 fail | 30 done | 99 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4001 submitted (5 run | 0 fail | 30 done | 98 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4002 submitted (6 run | 0 fail | 30 done | 97 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4011 submitted (7 run | 0 fail | 30 done | 96 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4012 submitted (8 run | 0 fail | 30 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_2701 finished (7 run | 0 fail | 31 done | 95 left)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "30-Apr-2017 01:22:40 region_growing_area_2702 finished (6 run | 0 fail | 32 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_3001 finished (5 run | 0 fail | 33 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_3002 finished (4 run | 0 fail | 34 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4001 finished (3 run | 0 fail | 35 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4011 finished (2 run | 0 fail | 36 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4012 finished (1 run | 0 fail | 37 done | 95 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4021 submitted (2 run | 0 fail | 37 done | 94 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4022 submitted (3 run | 0 fail | 37 done | 93 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4101 submitted (4 run | 0 fail | 37 done | 92 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4102 submitted (5 run | 0 fail | 37 done | 91 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4111 submitted (6 run | 0 fail | 37 done | 90 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4112 submitted (7 run | 0 fail | 37 done | 89 left)\n", "30-Apr-2017 01:22:40 region_growing_area_4201 submitted (8 run | 0 fail | 37 done | 88 left)\n", "30-Apr-2017 01:22:41 region_growing_area_4002 finished (7 run | 0 fail | 38 done | 88 left)\n", "30-Apr-2017 01:22:41 region_growing_area_4022 finished (6 run | 0 fail | 39 done | 88 left)\n", "30-Apr-2017 01:22:41 region_growing_area_4101 finished (5 run | 0 fail | 40 done | 88 left)\n", "30-Apr-2017 01:22:41 region_growing_area_4111 finished (4 run | 0 fail | 41 done | 88 left)\n", "30-Apr-2017 01:22:41 region_growing_area_4202 submitted (5 run | 0 fail | 41 done | 87 left)\n", "30-Apr-2017 01:22:41 region_growing_area_5001 submitted (6 run | 0 fail | 41 done | 86 left)\n", "30-Apr-2017 01:22:41 region_growing_area_5002 submitted (7 run | 0 fail | 41 done | 85 left)\n", "30-Apr-2017 01:22:41 region_growing_area_5011 submitted (8 run | 0 fail | 41 done | 84 left)\n", "30-Apr-2017 01:22:42 region_growing_area_4021 finished (7 run | 0 fail | 42 done | 84 left)\n", "30-Apr-2017 01:22:42 region_growing_area_4112 finished (6 run | 0 fail | 43 done | 84 left)\n", "30-Apr-2017 01:22:42 region_growing_area_4201 finished (5 run | 0 fail | 44 done | 84 left)\n", "30-Apr-2017 01:22:42 region_growing_area_4202 finished (4 run | 0 fail | 45 done | 84 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5012 submitted (5 run | 0 fail | 45 done | 83 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5021 submitted (6 run | 0 fail | 45 done | 82 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5022 submitted (7 run | 0 fail | 45 done | 81 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5101 submitted (8 run | 0 fail | 45 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_4102 finished (7 run | 0 fail | 46 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5001 finished (6 run | 0 fail | 47 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5011 finished (5 run | 0 fail | 48 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5002 finished (4 run | 0 fail | 49 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5012 finished (3 run | 0 fail | 50 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5021 finished (2 run | 0 fail | 51 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5022 finished (1 run | 0 fail | 52 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5101 finished (0 run | 0 fail | 53 done | 80 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5102 submitted (1 run | 0 fail | 53 done | 79 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5201 submitted (2 run | 0 fail | 53 done | 78 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5202 submitted (3 run | 0 fail | 53 done | 77 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5301 submitted (4 run | 0 fail | 53 done | 76 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5302 submitted (5 run | 0 fail | 53 done | 75 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5401 submitted (6 run | 0 fail | 53 done | 74 left)\n", "30-Apr-2017 01:22:42 region_growing_area_5402 submitted (7 run | 0 fail | 53 done | 73 left)\n", "30-Apr-2017 01:22:42 region_growing_area_6001 submitted (8 run | 0 fail | 53 done | 72 left)\n", "30-Apr-2017 01:22:43 region_growing_area_5102 finished (7 run | 0 fail | 54 done | 72 left)\n", "30-Apr-2017 01:22:43 region_growing_area_5201 finished (6 run | 0 fail | 55 done | 72 left)\n", "30-Apr-2017 01:22:43 region_growing_area_5202 finished (5 run | 0 fail | 56 done | 72 left)\n", "30-Apr-2017 01:22:43 region_growing_area_5301 finished (4 run | 0 fail | 57 done | 72 left)\n", "30-Apr-2017 01:22:43 region_growing_area_5402 finished (3 run | 0 fail | 58 done | 72 left)\n", "30-Apr-2017 01:22:43 region_growing_area_6002 submitted (4 run | 0 fail | 58 done | 71 left)\n", "30-Apr-2017 01:22:43 region_growing_area_6101 submitted (5 run | 0 fail | 58 done | 70 left)\n", "30-Apr-2017 01:22:43 region_growing_area_6102 submitted (6 run | 0 fail | 58 done | 69 left)\n", "30-Apr-2017 01:22:43 region_growing_area_6201 submitted (7 run | 0 fail | 58 done | 68 left)\n", "30-Apr-2017 01:22:43 region_growing_area_6202 submitted (8 run | 0 fail | 58 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_5302 finished (7 run | 0 fail | 59 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_5401 finished (6 run | 0 fail | 60 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6001 finished (5 run | 0 fail | 61 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6002 finished (4 run | 0 fail | 62 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6101 finished (3 run | 0 fail | 63 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6202 finished (2 run | 0 fail | 64 done | 67 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6211 submitted (3 run | 0 fail | 64 done | 66 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6212 submitted (4 run | 0 fail | 64 done | 65 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6221 submitted (5 run | 0 fail | 64 done | 64 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6222 submitted (6 run | 0 fail | 64 done | 63 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6301 submitted (7 run | 0 fail | 64 done | 62 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6302 submitted (8 run | 0 fail | 64 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6201 finished (7 run | 0 fail | 65 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6102 finished (6 run | 0 fail | 66 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6211 finished (5 run | 0 fail | 67 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6212 finished (4 run | 0 fail | 68 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6222 finished (3 run | 0 fail | 69 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6301 finished (2 run | 0 fail | 70 done | 61 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6401 submitted (3 run | 0 fail | 70 done | 60 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6402 submitted (4 run | 0 fail | 70 done | 59 left)\n", "30-Apr-2017 01:22:44 region_growing_area_7001 submitted (5 run | 0 fail | 70 done | 58 left)\n", "30-Apr-2017 01:22:44 region_growing_area_7002 submitted (6 run | 0 fail | 70 done | 57 left)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "30-Apr-2017 01:22:44 region_growing_area_7011 submitted (7 run | 0 fail | 70 done | 56 left)\n", "30-Apr-2017 01:22:44 region_growing_area_7012 submitted (8 run | 0 fail | 70 done | 55 left)\n", "30-Apr-2017 01:22:44 region_growing_area_6221 finished (7 run | 0 fail | 71 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_6302 finished (6 run | 0 fail | 72 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_6401 finished (5 run | 0 fail | 73 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_6402 finished (4 run | 0 fail | 74 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7001 finished (3 run | 0 fail | 75 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7011 finished (2 run | 0 fail | 76 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7012 finished (1 run | 0 fail | 77 done | 55 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7021 submitted (2 run | 0 fail | 77 done | 54 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7022 submitted (3 run | 0 fail | 77 done | 53 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7101 submitted (4 run | 0 fail | 77 done | 52 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7102 submitted (5 run | 0 fail | 77 done | 51 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8101 submitted (6 run | 0 fail | 77 done | 50 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8102 submitted (7 run | 0 fail | 77 done | 49 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8111 submitted (8 run | 0 fail | 77 done | 48 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7002 finished (7 run | 0 fail | 78 done | 48 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7022 finished (6 run | 0 fail | 79 done | 48 left)\n", "30-Apr-2017 01:22:45 region_growing_area_7101 finished (5 run | 0 fail | 80 done | 48 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8112 submitted (6 run | 0 fail | 80 done | 47 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8121 submitted (7 run | 0 fail | 80 done | 46 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8122 submitted (8 run | 0 fail | 80 done | 45 left)\n", "30-Apr-2017 01:22:45 region_growing_area_8111 finished (7 run | 0 fail | 81 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_7021 finished (6 run | 0 fail | 82 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_7102 finished (5 run | 0 fail | 83 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8101 finished (4 run | 0 fail | 84 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8102 finished (3 run | 0 fail | 85 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8112 finished (2 run | 0 fail | 86 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8122 finished (1 run | 0 fail | 87 done | 45 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8201 submitted (2 run | 0 fail | 87 done | 44 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8202 submitted (3 run | 0 fail | 87 done | 43 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8211 submitted (4 run | 0 fail | 87 done | 42 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8212 submitted (5 run | 0 fail | 87 done | 41 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8301 submitted (6 run | 0 fail | 87 done | 40 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8302 submitted (7 run | 0 fail | 87 done | 39 left)\n", "30-Apr-2017 01:22:46 region_growing_area_9001 submitted (8 run | 0 fail | 87 done | 38 left)\n", "30-Apr-2017 01:22:46 region_growing_area_8121 finished (7 run | 0 fail | 88 done | 38 left)\n", "30-Apr-2017 01:22:47 region_growing_area_8201 finished (6 run | 0 fail | 89 done | 38 left)\n", "30-Apr-2017 01:22:47 region_growing_area_8202 finished (5 run | 0 fail | 90 done | 38 left)\n", "30-Apr-2017 01:22:47 region_growing_area_8211 finished (4 run | 0 fail | 91 done | 38 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9002 submitted (5 run | 0 fail | 91 done | 37 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9011 submitted (6 run | 0 fail | 91 done | 36 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9012 submitted (7 run | 0 fail | 91 done | 35 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9021 submitted (8 run | 0 fail | 91 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_8302 finished (7 run | 0 fail | 92 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9001 finished (6 run | 0 fail | 93 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_8301 finished (5 run | 0 fail | 94 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_8212 finished (4 run | 0 fail | 95 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9002 finished (3 run | 0 fail | 96 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9011 finished (2 run | 0 fail | 97 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9012 finished (1 run | 0 fail | 98 done | 34 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9022 submitted (2 run | 0 fail | 98 done | 33 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9031 submitted (3 run | 0 fail | 98 done | 32 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9032 submitted (4 run | 0 fail | 98 done | 31 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9041 submitted (5 run | 0 fail | 98 done | 30 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9042 submitted (6 run | 0 fail | 98 done | 29 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9051 submitted (7 run | 0 fail | 98 done | 28 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9052 submitted (8 run | 0 fail | 98 done | 27 left)\n", "30-Apr-2017 01:22:47 region_growing_area_9021 finished (7 run | 0 fail | 99 done | 27 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9022 finished (6 run | 0 fail | 100 done | 27 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9031 finished (5 run | 0 fail | 101 done | 27 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9032 finished (4 run | 0 fail | 102 done | 27 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9042 finished (3 run | 0 fail | 103 done | 27 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9061 submitted (4 run | 0 fail | 103 done | 26 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9062 submitted (5 run | 0 fail | 103 done | 25 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9071 submitted (6 run | 0 fail | 103 done | 24 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9072 submitted (7 run | 0 fail | 103 done | 23 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9081 submitted (8 run | 0 fail | 103 done | 22 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9041 finished (7 run | 0 fail | 104 done | 22 left)\n", "30-Apr-2017 01:22:48 region_growing_area_9082 submitted (8 run | 0 fail | 104 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9051 finished (7 run | 0 fail | 105 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9052 finished (6 run | 0 fail | 106 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9061 finished (5 run | 0 fail | 107 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9062 finished (4 run | 0 fail | 108 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9071 finished (3 run | 0 fail | 109 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9072 finished (2 run | 0 fail | 110 done | 21 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9081 finished (1 run | 0 fail | 111 done | 21 left)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "30-Apr-2017 01:22:49 region_growing_area_9100 submitted (2 run | 0 fail | 111 done | 20 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9110 submitted (3 run | 0 fail | 111 done | 19 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9120 submitted (4 run | 0 fail | 111 done | 18 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9130 submitted (5 run | 0 fail | 111 done | 17 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9140 submitted (6 run | 0 fail | 111 done | 16 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9150 submitted (7 run | 0 fail | 111 done | 15 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9160 submitted (8 run | 0 fail | 111 done | 14 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9082 finished (7 run | 0 fail | 112 done | 14 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9100 finished (6 run | 0 fail | 113 done | 14 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9120 finished (5 run | 0 fail | 114 done | 14 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9170 submitted (6 run | 0 fail | 114 done | 13 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9110 finished (5 run | 0 fail | 115 done | 13 left)\n", "30-Apr-2017 01:22:49 region_growing_area_9140 finished (4 run | 0 fail | 116 done | 13 left)\n", "30-Apr-2017 01:22:50 region_growing_area_9130 finished (3 run | 0 fail | 117 done | 13 left)\n", "30-Apr-2017 01:22:50 region_growing_area_9150 finished (2 run | 0 fail | 118 done | 13 left)\n", "30-Apr-2017 01:22:50 region_growing_area_9160 finished (1 run | 0 fail | 119 done | 13 left)\n", "30-Apr-2017 01:22:50 region_growing_area_9170 finished (0 run | 0 fail | 120 done | 13 left)\n", "30-Apr-2017 01:22:50 merge_part submitted (1 run | 0 fail | 120 done | 12 left)\n", "30-Apr-2017 01:22:51 merge_part finished (0 run | 0 fail | 121 done | 12 left)\n", "30-Apr-2017 01:22:51 tseries_atoms_subject1_session1_motor submitted (1 run | 0 fail | 121 done | 11 left)\n", "30-Apr-2017 01:22:51 tseries_atoms_subject2_session1_motor submitted (2 run | 0 fail | 121 done | 10 left)\n", "30-Apr-2017 01:22:52 tseries_atoms_subject1_session1_motor finished (1 run | 0 fail | 122 done | 10 left)\n", "30-Apr-2017 01:22:52 tseries_atoms_subject2_session1_motor finished (0 run | 0 fail | 123 done | 10 left)\n", "30-Apr-2017 01:22:52 stability_ind_subject1 submitted (1 run | 0 fail | 123 done | 9 left)\n", "30-Apr-2017 01:22:52 stability_ind_subject2 submitted (2 run | 0 fail | 123 done | 8 left)\n", ".\n", "30-Apr-2017 01:23:32 stability_ind_subject1 finished (1 run | 0 fail | 124 done | 8 left)\n", "30-Apr-2017 01:23:32 msteps_ind_subject1 submitted (2 run | 0 fail | 124 done | 7 left)\n", "30-Apr-2017 01:23:33 msteps_ind_subject1 finished (1 run | 0 fail | 125 done | 7 left)\n", ".\n", "30-Apr-2017 01:23:40 stability_ind_subject2 finished (0 run | 0 fail | 126 done | 7 left)\n", "30-Apr-2017 01:23:40 stability_group_sci10 submitted (1 run | 0 fail | 126 done | 6 left)\n", "30-Apr-2017 01:23:40 stability_group_sci20 submitted (2 run | 0 fail | 126 done | 5 left)\n", "30-Apr-2017 01:23:40 stability_group_sci30 submitted (3 run | 0 fail | 126 done | 4 left)\n", "30-Apr-2017 01:23:40 summary_stability_avg_ind submitted (4 run | 0 fail | 126 done | 3 left)\n", "30-Apr-2017 01:23:40 msteps_ind_subject2 submitted (5 run | 0 fail | 126 done | 2 left)\n", "30-Apr-2017 01:23:42 msteps_ind_subject2 finished (4 run | 0 fail | 127 done | 2 left)\n", "30-Apr-2017 01:23:45 summary_stability_avg_ind finished (3 run | 0 fail | 128 done | 2 left)\n", "30-Apr-2017 01:23:45 Stopping idle worker 4 (not enough jobs left to do).\n", ".\n", "30-Apr-2017 01:24:33 stability_group_sci20 finished (2 run | 0 fail | 129 done | 2 left)\n", "30-Apr-2017 01:24:33 Stopping idle worker 2 (not enough jobs left to do).\n", "30-Apr-2017 01:24:36 stability_group_sci10 finished (1 run | 0 fail | 130 done | 2 left)\n", "30-Apr-2017 01:24:36 stability_group_sci30 finished (0 run | 0 fail | 131 done | 2 left)\n", "30-Apr-2017 01:24:36 summary_stability_group submitted (1 run | 0 fail | 131 done | 1 left)\n", "30-Apr-2017 01:24:36 msteps_group submitted (2 run | 0 fail | 131 done | 0 left)\n", ".\n", "30-Apr-2017 01:24:40 summary_stability_group finished (1 run | 0 fail | 132 done | 0 left)\n", "30-Apr-2017 01:24:40 Stopping idle worker 1 (not enough jobs left to do).\n", "Deamon terminated on 30-Apr-2017 01:24:51\n", "\n", "30-Apr-2017 01:24:51 msteps_group finished (0 run | 0 fail | 133 done | 0 left)\n", "30-Apr-2017 01:24:51 Stopping idle worker 3 (not enough jobs left to do).\n", "\n", "*******************************************\n", "Pipeline terminated on 30-Apr-2017 01:24:51\n", "All jobs have been successfully completed.\n", "\n" ] } ], "source": [ "niak_pipeline_stability_rest(files_in,opt);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now that the pipeline has successfully completed, we can check a .csv file generated as part of the outpus. This .csv file list the number of clusters identified as \"representative\" of all stable solutions across the grid of clustering parameters we selected, using a method called MSTEPS. " ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "tab =\n", "\n", " 10 7 8\n", " 10 12 11\n", " 20 16 22\n", " 20 24 27\n", " 30 30 40\n", "\n" ] } ], "source": [ "file_msteps = [opt.folder_out filesep 'stability_group' filesep 'msteps_group_table.csv'];\n", "[tab,lx,ly] = niak_read_csv(file_msteps);\n", "tab" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "So each line is a number of clusters at the individual/group/consensus levels. Here MSTEPS has selected three scales as appropriate to summarize all tested numbers of clusters. Note that, even though we selected only 10, 20 and 30 in our grid, that grid is strictly used only at the individual level. At the group and consensus levels, other parameters are explored, in the neighbourhood of the number of individual clusters. We can now set these parameters in the pipeline to generate stability maps and parcellations:" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": true }, "outputs": [], "source": [ "opt.scales_maps = tab;" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can finally restart the pipeline. It will pick up where it left, and only deal with the generation of the maps. " ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Logs will be stored in /sandbox/home/git/niak_tutorials/basc_test_niak/logs/\n", "Generating dependencies ...\n", " Percentage completed : 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100- 0.35 sec\n", "Setting up the to-do list ...\n", " I found 10 job(s) to do, and 133 job(s) already completed.\n", "\n", "Deamon started on 30-Apr-2017 01:25:00\n", "30-Apr-2017 01:25:00 Starting the pipeline manager...\n", "30-Apr-2017 01:25:00 Starting the garbage collector...\n", "30-Apr-2017 01:25:00 Starting worker number 1...\n", "30-Apr-2017 01:25:00 Starting worker number 2...\n", "30-Apr-2017 01:25:01 Starting worker number 3...\n", "30-Apr-2017 01:25:01 Starting worker number 4...\n", "\n", "Pipeline started on 30-Apr-2017 01:25:02\n", "user: , host: 3f86dc5d68dc, system: unix\n", "****************************************\n", "30-Apr-2017 01:25:05 stability_maps_group_sci10_scg7_scf8 submitted (1 run | 0 fail | 133 done | 9 left)\n", "30-Apr-2017 01:25:05 figure_stability_group_sci10_scg7_scf8 submitted (2 run | 0 fail | 133 done | 8 left)\n", "30-Apr-2017 01:25:05 stability_maps_group_sci10_scg12_scf11 submitted (3 run | 0 fail | 133 done | 7 left)\n", "30-Apr-2017 01:25:05 figure_stability_group_sci10_scg12_scf11 submitted (4 run | 0 fail | 133 done | 6 left)\n", "30-Apr-2017 01:25:05 stability_maps_group_sci20_scg16_scf22 submitted (5 run | 0 fail | 133 done | 5 left)\n", "30-Apr-2017 01:25:06 figure_stability_group_sci20_scg16_scf22 submitted (6 run | 0 fail | 133 done | 4 left)\n", "30-Apr-2017 01:25:06 stability_maps_group_sci20_scg24_scf27 submitted (7 run | 0 fail | 133 done | 3 left)\n", "30-Apr-2017 01:25:06 figure_stability_group_sci20_scg24_scf27 submitted (8 run | 0 fail | 133 done | 2 left)\n", "30-Apr-2017 01:25:08 stability_maps_group_sci10_scg7_scf8 finished (7 run | 0 fail | 134 done | 2 left)\n", "30-Apr-2017 01:25:08 stability_maps_group_sci30_scg30_scf40 submitted (8 run | 0 fail | 134 done | 1 left)\n", "30-Apr-2017 01:25:08 stability_maps_group_sci10_scg12_scf11 finished (7 run | 0 fail | 135 done | 1 left)\n", "30-Apr-2017 01:25:08 figure_stability_group_sci30_scg30_scf40 submitted (8 run | 0 fail | 135 done | 0 left)\n", "30-Apr-2017 01:25:10 stability_maps_group_sci20_scg16_scf22 finished (7 run | 0 fail | 136 done | 0 left)\n", "30-Apr-2017 01:25:10 stability_maps_group_sci20_scg24_scf27 finished (6 run | 0 fail | 137 done | 0 left)\n", "30-Apr-2017 01:25:12 stability_maps_group_sci30_scg30_scf40 finished (5 run | 0 fail | 138 done | 0 left)\n", "30-Apr-2017 01:25:12 Stopping idle worker 1 (not enough jobs left to do).\n", "30-Apr-2017 01:25:13 figure_stability_group_sci10_scg7_scf8 finished (4 run | 0 fail | 139 done | 0 left)\n", "30-Apr-2017 01:25:13 figure_stability_group_sci10_scg12_scf11 finished (3 run | 0 fail | 140 done | 0 left)\n", "30-Apr-2017 01:25:18 figure_stability_group_sci30_scg30_scf40 finished (2 run | 0 fail | 141 done | 0 left)\n", "30-Apr-2017 01:25:18 Stopping idle worker 3 (not enough jobs left to do).\n", "30-Apr-2017 01:25:18 figure_stability_group_sci20_scg16_scf22 finished (1 run | 0 fail | 142 done | 0 left)\n", "30-Apr-2017 01:25:18 Stopping idle worker 2 (not enough jobs left to do).\n", "Deamon terminated on 30-Apr-2017 01:25:19\n", "\n", "30-Apr-2017 01:25:19 figure_stability_group_sci20_scg24_scf27 finished (0 run | 0 fail | 143 done | 0 left)\n", "30-Apr-2017 01:25:19 Stopping idle worker 4 (not enough jobs left to do).\n", "\n", "*******************************************\n", "Pipeline terminated on 30-Apr-2017 01:25:19\n", "All jobs have been successfully completed.\n", "\n" ] } ], "source": [ "niak_pipeline_stability_rest(files_in,opt);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The primary output of the BASC pipeline is a set of parcels for the selected pipeline. For example, the line `[10 7 8]` in opt.scales_maps means 10 individual clusters, 7 group clusters and 8 consensus (final) group clusters. Let's load and visualize that brain parcellation: " ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAJNmlDQ1BkZWZhdWx0X3JnYi5pY2MA\nAHiclZFnUJSHFobP933bCwvssnRYepMqZQHpvUmvogJL7yxLEbEhYgQiiog0RZCggAGjUiRWRLEQ\nFBSxoFkkCCgxGEVUUPLDOxPn3vHHfX49884755yZA0ARBQBARQFSUgV8Pxd7TkhoGAe+IZKXmW7n\n4+MJ3+X9KCAAAPdWfb/zXSjRMZk8AFgGgHxeOl8AgOQCgGaOIF0AgBwFAFZUUroAADkLACx+SGgY\nAHIDAFhxX30cAFhRX30eAFj8AD8HABQHQKLFfeNR3/h/9gIAKNvxBQmxMbkc/7RYQU4kP4aT6edi\nz3FzcOD48NNiE5Jjvjn4/yp/B0FMrgAAwCEtfRM/IS5ewPmfoUYGhobw7y/e+gICAAh78L//AwDf\n9NIaAbgLANi+f7OoaoDuXQBSj//NVI8CMAoBuu7wsvjZXzMcAAAeKMAAFkiDAqiAJuiCEZiBJdiC\nE7iDNwRAKGwAHsRDCvAhB/JhBxRBCeyDg1AD9dAELdAOp6EbzsMVuA634S6MwhMQwhS8gnl4D0sI\nghAROsJEpBFFRA3RQYwQLmKNOCGeiB8SikQgcUgqkoXkIzuREqQcqUEakBbkF+QccgW5iQwjj5AJ\nZBb5G/mEYigNZaHyqDqqj3JRO9QDDUDXo3FoBpqHFqJ70Sq0ET2JdqFX0NvoKCpEX6ELGGBUjI0p\nYboYF3PAvLEwLBbjY1uxYqwSa8TasV5sALuHCbE57COOgGPiODhdnCXOFReI4+EycFtxpbga3Alc\nF64fdw83gZvHfcHT8XJ4HbwF3g0fgo/D5+CL8JX4Znwn/hp+FD+Ff08gENgEDYIZwZUQSkgkbCaU\nEg4TOgiXCcOEScICkUiUJuoQrYjexEiigFhErCaeJF4ijhCniB9IVJIiyYjkTAojpZIKSJWkVtJF\n0ghpmrREFiWrkS3I3uRo8iZyGbmJ3Eu+Q54iL1HEKBoUK0oAJZGyg1JFaadco4xT3lKpVGWqOdWX\nmkDdTq2inqLeoE5QP9LEado0B1o4LYu2l3acdpn2iPaWTqer023pYXQBfS+9hX6V/oz+QYQpoifi\nJhItsk2kVqRLZETkNYPMUGPYMTYw8hiVjDOMO4w5UbKouqiDaKToVtFa0XOiY6ILYkwxQzFvsRSx\nUrFWsZtiM+JEcXVxJ/Fo8ULxY+JXxSeZGFOF6cDkMXcym5jXmFMsAkuD5cZKZJWwfmYNseYlxCWM\nJYIkciVqJS5ICNkYW53txk5ml7FPsx+wP0nKS9pJxkjukWyXHJFclJKVspWKkSqW6pAalfokzZF2\nkk6S3i/dLf1UBiejLeMrkyNzROaazJwsS9ZSlidbLHta9rEcKqct5ye3We6Y3KDcgryCvIt8uny1\n/FX5OQW2gq1CokKFwkWFWUWmorVigmKF4iXFlxwJjh0nmVPF6efMK8kpuSplKTUoDSktKWsoByoX\nKHcoP1WhqHBVYlUqVPpU5lUVVb1U81XbVB+rkdW4avFqh9QG1BbVNdSD1Xerd6vPaEhpuGnkabRp\njGvSNW00MzQbNe9rEbS4Wklah7XuaqPaJtrx2rXad3RQHVOdBJ3DOsOr8KvMV6Wualw1pkvTtdPN\n1m3TndBj63nqFeh1673WV9UP09+vP6D/xcDEINmgyeCJobihu2GBYa/h30baRjyjWqP7q+mrnVdv\nW92z+o2xjnGM8RHjhyZMEy+T3SZ9Jp9NzUz5pu2ms2aqZhFmdWZjXBbXh1vKvWGON7c332Z+3vyj\nhamFwOK0xV+WupZJlq2WM2s01sSsaVozaaVsFWnVYCW05lhHWB+1Ftoo2UTaNNo8t1WxjbZttp22\n07JLtDtp99rewJ5v32m/6GDhsMXhsiPm6OJY7DjkJO4U6FTj9MxZ2TnOuc153sXEZbPLZVe8q4fr\nftcxN3k3nluL27y7mfsW934Pmoe/R43Hc09tT75nrxfq5e51wGt8rdra1LXd3uDt5n3A+6mPhk+G\nz6++BF8f31rfF36Gfvl+A/5M/43+rf7vA+wDygKeBGoGZgX2BTGCwoNaghaDHYPLg4Uh+iFbQm6H\nyoQmhPaEEcOCwprDFtY5rTu4bircJLwo/MF6jfW5629ukNmQvOHCRsbGyI1nIvARwRGtEcuR3pGN\nkQtRblF1UfM8B94h3qto2+iK6NkYq5jymOlYq9jy2Jk4q7gDcbPxNvGV8XMJDgk1CW8SXRPrExeT\nvJOOJ60kByd3pJBSIlLOpYqnJqX2pymk5aYNp+ukF6ULMywyDmbM8z34zZlI5vrMHgFLkC4YzNLM\n2pU1kW2dXZv9ISco50yuWG5q7uAm7U17Nk3nOef9tBm3mbe5L18pf0f+xBa7LQ1bka1RW/u2qWwr\n3Da13WX7iR2UHUk7fiswKCgveLczeGdvoXzh9sLJXS672opEivhFY7std9f/gPsh4YehPav3VO/5\nUhxdfKvEoKSyZLmUV3rrR8Mfq35c2Ru7d6jMtOzIPsK+1H0P9tvsP1EuVp5XPnnA60BXBaeiuOLd\nwY0Hb1YaV9YfohzKOiSs8qzqqVat3le9XBNfM1prX9tRJ1e3p27xcPThkSO2R9rr5etL6j8dTTj6\nsMGloatRvbHyGOFY9rEXTUFNAz9xf2pplmkuaf58PPW48ITfif4Ws5aWVrnWsja0Latt9mT4ybs/\nO/7c067b3tDB7ig5BaeyTr38JeKXB6c9Tved4Z5pP6t2tq6T2VnchXRt6prvju8W9oT2DJ9zP9fX\na9nb+aver8fPK52vvSBxoewi5WLhxZVLeZcWLqdfnrsSd2Wyb2Pfk6shV+/3+/YPXfO4duO68/Wr\nA3YDl25Y3Th/0+LmuVvcW923TW93DZoMdv5m8lvnkOlQ1x2zOz13ze/2Dq8ZvjhiM3LlnuO96/fd\n7t8eXTs6/CDwwcOx8DHhw+iHM4+SH715nP146cn2cfx48VPRp5XP5J41/q71e4fQVHhhwnFi8Ln/\n8yeTvMlXf2T+sTxV+IL+onJacbplxmjm/Kzz7N2X615OvUp/tTRX9KfYn3WvNV+f/cv2r8H5kPmp\nN/w3K3+XvpV+e/yd8bu+BZ+FZ+9T3i8tFn+Q/nDiI/fjwKfgT9NLOcvE5arPWp97v3h8GV9JWVn5\nBy6ikLxSF1/9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAAHXRFWHRTb2Z0d2FyZQBHUEwgR2hvc3Rz\nY3JpcHQgOS4wNfOvXY8AACAASURBVHic7d1Nj9vYmejxR+kbt8vdd2L6GmNpJuM7INKAV9Mw+AkC\nqPe1UT6CXFpl1aPalWpnDZDNbKTSV9DGq9mUNtk2QhgezCKYBojcNAzJE6PpxnVHfTu44SwoHh4V\nj1lSFV+OpP8PjYDFIsVHqqo8Puc5L40oigQAgLr9pO4AAAAQISEBACxBQgIAWIGEBACwAgkJAGAF\nEhIAwAokJACAFUhIAAArkJAAAFYgIQEArEBCAgBYgYQEALACCQkAYAUSEgDACiQkAIAVSEgAACuQ\nkAAAViAhAQCsQEICAFiBhAQAsAIJCQBgBRISAMAKJCQAgBVISAAAK5CQAABWICEBAKxAQgIAWIGE\nBACpMAzDMKw7igP1P+oOAABsMZ1O42zkum673a47nINDCwkAVnzf73a73W53NpvVHcshakRRVHcM\nMPjtb3/71VdflfqIb7/99sGDB6U+Ys+eUtmDLHzKYrH4zW9+U2o81/r1r3/9ww8/xMd37979+c9/\nvtXty+Xy9evX6kvP87rd7pVrptNpfOA4Di2k6tFlZ6mvvvrqn//5f5b8kNciZT9iz55S2YMsfMrr\nuvORvPvXf3W1L/+8/Sv8nXb83b/8S/aCIAjiPOT7/vYvj9siIQHAiuu6nueJSBAEdcdyiEhIAJCa\nTCYi4jhO3YEcIhISAKx0Op14lB0JqRYkpEP2hKfY+qB9esqOIRXViIQEYDf8VOTjumNAqZiHBACw\nAgkJAGAFEhIAwAokJACAFRjUAGA3PBBpFvdqfyrupVAUWkgAACuQkAAAViAhAQCsQEICAFiBhAQA\nsAKj7ADshjsid+uOAaWihQQAsAIJCQBgBRISAMAKJCQAgBUY1ABgN3wqcr+4V/umuJdCUWghAQCs\nQEICAFiBhAQAsAI1JGCXDKSXHIxu9UKjE8PJ3vhWrwncDgkJwG74pNBBDbAQXXYAACuQkAAAVqDL\nDrBeM633DGR1PFg0tJMb1ZOik176xate9oLGLetSwO3QQgIAWIGEBACwAl12AHbDA5Fm3TGgVLSQ\nAABWoIUElEUfRNAY33y8wFwbv2A02P41F8kUWG3AhIyiV/FBr/H59i8J3BYtJACAFUhIAAAr0GVX\npzAMHcepOwpgN9wRuVt3DCgVCaki0+nU9/12u91ut9WZIAgcx3FdV50E8j2N5oazjReb3KvXseZi\nmBhL6Qj1IiFVpNPpXGkMBUHQ7/dFZDgckpAAgIRUG5WfPM/Lfne5XIq81k48FrlXSVyA8lbkbXJ8\nEFt+f/HFF/HfYxiGFxcXdYdzcEhIljo6OhJ5UncUOHAPRR7WHUOlLi8v44PJZFJvJIeJhFSbMAzj\ng7i2VG8wKMNt5h7pWvqap5uVi7Z4zd3Zk28h8jI5vi/ydMvb32m3y6oTwmw2mxn7LVA2ElJFJpOJ\n7/siEgRBt9sVEdd1h8Oh4zj86gOb+IXIL29x+/312//j6OhDV/q+H9d3UTESUkXiJKTrdDpxI4mR\n3wAgJKR6kYoA29CFXiMSEgCk6EKvEQkJ2Hksior9QEICsBt+ekfusvrmXuPHCwCwAgkJAGAFuuyA\nPNF8tQhpo2WY5aqKNyJy3FhNMV2bx1oJVTpS0YoWcFph0hdX7Rk2/as+ckBHCwkAYAVaSAB2xAM2\nRNpztJAAAFaghYTa6AWYk8WqADNunmSvrLE8o+LRyzOL1urghaTRqtjU3nf6WylqodX8mlarmW7f\nN5fW1SB76UKq2WiB2tFCAgBYgYQEALACXXYAdsQdBjXsORISqpYWLbSZMK1oVfx42VhVPp5GaTnk\n5GRV/Ii0gkdRVRklOknKM9orp+UrLZ6XSXnmXM6021tSDr18le9MzrMntSDTzf0GySfJNt2wB112\nAAArkJAAAFYgIQEArEBCAgBYgUENqJo+WkFZNFb1djVhs2m6dzE2nS3aQBs6YYxWnZz30oEM49Fq\n/qkaB1EUfQ6sGnkxME1oPRct2ujqd9W9a7dU8pEW41ORv6k7BpSJFhIAwAokJACAFUhIAAArUENC\n1VS5aPPLGpWsqarmwxrLMzoV20I/m8T4IloVk9S+eQXavOSjnq5WWdWjVUEOxsUHCdwMCQnAjvhE\n5H7dMaBMdNkBAKxAQgIAWIEuO+CqgV6yyq14rW0Y2PjwdcUZ5JfTTNGadzWsJFpgK7SQAABWoIUE\nYEc8+MACHtgXtJAAAFYgIQEArEBCAgBYgYQEALACCak6YRheewYADhaj7CoyHA5d1w3D0HGcTqcj\nItPpNAgCx3Fc122323UHCFjvjsjdumNAmWghVcTzvE6n0+12fd+PzwRB0O/39TMALBEEQd0hHCJa\nSBWJ20C+76vGkOM48YHnednrl8ulyGvtxGORe2UHCax7K/I2Of6mzkAqFHdmiEj8v6gSCak6vu8H\nQRD3113r6OhI5EnZIQG5Hoo8rDuGSs1ms3a7bfw3IipAl11FstlIjWigyw6wRBAEQRBMp1P+KmtB\nC6kKvu+fnp56nuf7fhiGFxcXIuK67nA4dByHf44Bm/j6R5F3q2PnrnhbLiMU/iC+tkfhcrnMXhP/\nq9HzvOFwyB9m9UhIVfA87/Ly8srJTqcTN5JUMQlAjs8caf/jzW937q7d/vLoKHuN67pxHuKvshYk\npDrxSw9YpdPpqEENqB4JCQBWHMfp9/u+79NfVwsSEqwwkN7qYBStTvXGdQUzT4IRkQvTBWmQuvPk\nYFFF5AMtyEfRSXzwppE+Ohtk9Crdla8xzt3o77CRjepCQgKwIz4R+VndMaBMDPsGAFiBhAQAsAIJ\nCQBgBWpIqJoaMqCPFzAME2ieGG4uc7yAPpYhyzyQQTOaX7291/j8tjF92DPt+EVyoAc56DXkCu0T\njZJoGy1GN8AWtJAAAFaghQRgRzwQ2XK5oDx/Ku6lUBBaSAAAK9BCQtWa89XB4HbVi2Z0HB8sGi/y\nr9yQqmmdzdOTmwfZS65UxaRIqyoVVapRha6WaC+oqkW5D9Enw0ZnJQYJ3AwtJACAFUhIAAAr0GUH\nYEfcEblbdwwoEwkJJVILgOpVGbP8pVSTOUmqbiQi80UrPmjkV06ukw2y8VSfcmQIbHSyumVtMdMk\njGN18tp3vZkoM8NJRKR1zV1qxdWUNjHp3HQ7k5NQL7rsAABWICEBAKxAQgIAWIEaEko0aK6KMYOn\n6cnoZWaNtes054vsybFxsbvtZYPUIzRWU1TpaJBbvtIjHEWv4oOiFrhTL7j2mloprtfLe5BaCs+4\nA6GlPhW5X9yrfVPcS6EgtJAAAFYgIQEArEBCAgBYgYQEALACgxp21SB/N7nt54quzTntJXNOx7eb\nc7pYjQ7QJ8aqUr++IV6reXX7u+g4HVmwaKyufKFN9jxJCvjn0Vy7cuuFVlWQ6UxS7aNNg+yZoj3T\nAn61eh21ieC5/u3UphGqR+uTWFWQJ9r4hWMxaGV2FFQR6kGe6UNDzjcMDSgFCQnAjvik0FF2sA9d\ndgAAK5CQAABWoMtuVw20oouaN3rSuG7FzYzoRG34llYiWqOkKvO5tmV0/vqnJuOkuvPIVJ3QC0Kj\npF6itrlrvNBKIM2rl+lBqoqXbLzQql4wGydv/MwUpJrbutDffVIb0qsy6YdmCiEty21c3rvtlNXk\n56WWTF1otbGnkWHl15fb//4ABaKFBACwAi0kADviQdpWxl6ihQQAsAItpF010jZtU2trqhUzn+m1\nls2KFvoMnnThTm2SSv5inYZ7RY4bqkCVxjA3zaDSLjBUqtSkq14jfZ1RtCr56NUdVRLLn0GlF0uy\n1ZS0GvShiU1JeUYvCEXJlnnqdr0olQaprwm7yCvLPTOdVNHoQc5lo3ljzc03DDSXzoBy0UICAFiB\nhAQAsAIJqTphGIZheOVMXcEAgG2oIVVkOp2KSBiGjuN0Op34TBAEjuO4rttut+sOELDeHZG7dceA\nMpGQKhInIRE5PT2Nj4Mg6Pf7IjIcDjdPSOmIAG3b1XR1UdMggnyqDK6PNXghN9+MtafV1XujkTqr\nTqpBBIuWNoZ3s/L+2mANNcDBtOJqPn1OrqLWQr3BCq2iDygYJ7drU3vTJ24/99Q4iXXt0WoSsfaJ\nzhdb78x7myCB2yMhVWc2m81mM5WZHMeJDzzPy168XC5FXmsnHovcKz1EQPf/A5HfJ18cxI7fp6en\n8UG73abfonokpOq0223P8yaTiTEDXXF0dCTypIKogA/6yD3AX8Lnz5/XHcLhIiFVSrWKRBvR4Ps+\n/xYDrvV1kB47PxPvn7a7PfxO/H9Pv1wul9lrHMeZTCZhGLquqzozUBkSUkUmk4njOPEvenzGdd3h\ncOg4ziYNJsVYHFKFg7TYsH0VRF/Kc9DYaA7s9ZKKjr6YqSodqTVhRWT7ekda9GpqE221JVk3ndGp\nJsluXnjbMLCWNo1XzdjtbfwUVdHRJ9gaC13pEzP7HOpa2sxY9a71AtUHdhS0xWf/JO1f3vx25760\n/3f65cs/HGWvicu6IjIcDm/+JNwUCaki3W43HlOnGkmdTiduJOnNJgA2YEpGLUhI1VFtI4VUBFjl\n9PTU87wwDLfqt0BRSEgAsPL8+fMrPRmoEglpH6hygirVrNVscutJqrZxXtQqmoZd7tbopaMcA9P6\nsAPTwqx63SM6XhWjjMuMGmdxtUZJ3SWJ3BjhtZOT1Gf+VG11aLpFRfihIJXjxupHos+ayg9Dn3uk\nPrRm8zh7r/ru3LTObP7HuN+yPRmoDAkJwG74y0/v/HCX1c72GT9dAIAVSEgAACvQZbcPVGkkrV5o\nVG3DvAVfUlU509ZLe9baek28lL4qXbLP27yn1TZGhpXZVJDznmkZtfEqHuOedXr1a5AbcFo40WtR\nSelIFU5akhehHuRa4a1xtb6lB6PetfkN5lLFJNGmMenxqEfrm/+piViLcTN7i/pl0H/E0YI17FAn\nWkgAACuQkAAAVqDLDsBu+FYeLNgQaa/RQgIAWIEW0j5Qa2YaV8xU0oEMutzlNPWJqPnjBcy3m9ZM\nVVV9fXTDPCmnq9VRR/P00c/GG23fJyKDzQJbW0k2CdI42uIk+Uj1QR/qyrNROg34TTL04EyNLBin\nb2GcjDdZe8o4b5armg+rD2pQXjYMAxDWtjnMjijRRl0Yh6uo1Vfn7NCHOtBCAgBYgYQEALACXXYA\ndsOPcucHBjXsNRLSrormhmVGFWOBQV+j82RxtSxxrt2hqiB6gWmw/eqrg2S/uPEoXa5UVWVejrVH\nJoUVtbDpG61woko+erlITQLNvJXSqT0G32z87NtsA6j/4NSysBfGSzVpES5380O9snieFLoOcE1V\n2IAuOwCAFUhIAAArkJAAAFaghrSrjMtgrhUbtqTPs1Hrjd62lpCUWE4yC49evTB5Oyfbz4C5QVVm\nbU6VKkYlT9anT6mqVeNppE6qffD0D1wVvYxbHabbAPY2LTv1Gp+vHqdNIOolZbmBtomgmi6minYi\nMm8YgkyNT7YMp37v5dN38jd1R4ES0UICAFiBhAQAsAIJCQBgBRISAMAKDGrYVXolf23/04xz0/qp\n583VSbWqqT4x9gbrqOYbmxZ7HSTTS3Vq29b8VV9FW4Y13Rv1ZtTk1mR0gx7ts2RC61wb6dA0jRJQ\no0LOjcMybjF4wDxqQwtiIIaA0pmzyWiLR9rohkGy4qp+ssd8WNSKhARgN3wvn7yT+3VHgRLRZQcA\nsAIJCQBgBbrs9kFaY0jKHPrGegv5POfeRtHlIqNFw7QTnV6KSQo56VTc62suRU/pTGJYtNLAWslU\n07lps0F9BVjJlMkKL8WJaBWvURpkugWfRn26aqKuHq0qHanptyJS/EcKbIMWEgDACrSQAOyGb+XB\nQgwjM7E3aCEBAKxAC2k/lVK9KFz1O+ttyBTY2mQgU+CDrVeFvR1tYpN5olL+p5u7ax9QC1pIAAAr\nkJAAAFagy65OYRg6jlN3FMBu+FHu/CB3644CJaKFVJ0wDH/1q1/NZrP4y+l0OplMJpOJOgMAh4yE\nVJ3JZNLtdtWXQRD0+/1ut+v7fo1RAYAlSEgV8X3f8zz9jOqsu3I+tlwuRX6v/ffnKqIE1rzVfgO/\nqTuYilzpyUCVSEgVmc1m7XZ78+uPjo5Enmj/3SsvNuADHmq/gf9QdzAVudKTgSoxqKEinufNZjPf\n91XDKAzD+MD3/a1yFXCY/uvr79TxPefjx97DrW7/c/jjH/0/qS8fLA0Tx7I9GagSCaki7XY7CILp\ndKoSkuu6w+HQcRz+AIBN3Pvs75rtf1Rfvtv2fkea7b9VX3308ih7yWw26/f79NfVhYRUHdd1Ly7S\nbTw7nU7cSGLkN2CJbE8GqkQNqU6O4/B7D9ij3W67rhsEQRAEdcdyiGghAUDqSk8GqkRC2k+j6JU6\nVtuyDUZR9sro1WqVzXRnPACoAwkJwG74Xj75Tn5W1Ks9KOqFUBxqSAAAK5CQAABWoMtur+ilo6xB\nz7Qp20nyXemlV5a6v19z9cjRPH2iKnQZnSVB1ljoMn4+zehYnVw0XlQdE7BfaCEBAKxACwnAbvhe\nPn0n94t6tQfyf4t6KRSFFhIAwAokJACAFeiy2wdzVW9PRi28iE7Ud9V8WDUHVkTGo9UF56bRBFEy\n3KDRKmgQQfMk//tnme+rCEVEFqsgo7N0ZEFhsWVp0UbHhpEg49Fx9mSkjdGIlRjhB0QnaQznyQ/2\nUfLL0Gt8XnE8wFZoIQEArEBCAgBYgS47ALvhfaGj7IRRdvYhIe2qbMVC9KKFXvhoJgeGwkc6x9P4\ngqXSSxq9ZlLoMtVs1PvSgyy80JW++Lnho1hotba5rDYbbY3m6mSruTqeLwxbkRZOLxcpi7zpxWvT\neFWQ1Re6gA+hyw4AYAUSEgDACiQkAIAVqCEdkMYLbYO+8dVqg15LUNUUvWZzg2KDdvumr5MGubZU\naV5tpLggTfEky7nOxXDZvKeVi84KDsxIK3SlJ9WcrWPtx6qmduXWlUoJsiTv5dOPixzU8E1xL4Vi\n0EICAFiBhAQAsAIJCQBgBWpIO0ZVMs5z57qsbxxnWq0umcXS+HxVs5nr2/ep+kSmLrJ5hNtRM2jU\nOnL5c2qKk5011Lpuf0K1TJy+BJ9xspRym1lT5iqX6Uejr2F4nPzcz0erS6uZIAXcGAkJwG74Xj75\nqMhBDbAOXXYAACuQkAAAViAhAQCsQA1pVz3Tji+SAzVUoaUtpKoGOCwaa3NNY2tjGRJqzdCXjbQM\nrkYrXFvzV/KL6KPolTpOF1o1jmVIJn5q4aSlfv0pGwZ5w5EXCTWWQU2bNdIHnjz78GXX0t/ghemC\ngaw+Pf19qSvV7F092sJHXgC3RwsJAGAFWkgAdsN7+fSvjLLba7SQAABWoIW0Dx4l0yFbyWxJY7mo\nes3Nai1mI23Sac9QWFJ1jlsWhFSQ5wXNxDXOWc7W+eS6jyWtCGrbAJ4kZbSBqUjW1D4zNXHW+L5U\nXao5N5wEakELCQBgBRISAMAKdNnVKQxDx3HqjgLYDe/l0x8Z1LDXSEgVOT09jQ/a7Xa73RaR6XQa\nBIHjOK7rxmc2oRcJlHQZzegs57r5Ip1ypMoSi/Gq2nTLSkw+/cU3ncZkqhvpr/M0SqofjfRk+r63\nLwgZqzKGV974dfSJVsemJW43pM8GM85DyqdmTQ1yi0mir81qmJwGlI6EVJ3nz5/rXwZB0O/3RWQ4\nHG6ekACUKgxDEaHrohYkpIo4jjOZTMIwdF230+mI9hvveV72+uVyKfJaO/FY5F4VgQKptyJvk+OD\n2PB7Op1K0pce/52iSiSkisSNIREZDoebXH90dCTypMyIgGs9FHlYdwyVUkno9PSUhFQ9ElLV4g4B\n/cD3/c277Bovkv30FoZufjX9aK1m04xyXtC4IpxetNiWsdaydrLojffSYtLGkW9eEMqXP6koXaBP\n0rUFy9h3UA9jQ/mfwA3WLazA+68XIi/j44Zzv+E93e7+8N1f/Zfqq+VyabxqNpvNZjOyUS1ISBU5\nPT31PC8MQ9VB57rucDh0HMfYZQfgqs9+8ZP2L29+u3Nfv/3o5X8Yr2q3257nTSYT/jCrR0KqyPPn\nz+Mxdap01Ol0KJ8CFuJPsi4kpOq4rnvlDL/3gFUmk4njOPHgo7pjOUQkJABY6Xa7V3oyUCUS0o5R\nYxn0oQojtdmacT5jUkZfK1Anwx8aUswSpdoLp4FFx3kTLPWa/23cYAiGPrJgwwEOm7+vfJvPbN18\n+IO6Uv+t2DBI/SlvxqsvbjD9dm/QNqoRCQnAbvjxLz+VH+7WHQVKxOKqAAArkJAAAFagy25Xra2U\n2rg6e9E4nzF/hVP9jLGetOG6pXpgad3l2npI8urGCb8qNmMZ7AZB6sqYqVqeR6b1T1XJZ2D69IzU\nuz7kchFsQwsJAGAFWkgAdsR7kXd1x4Ay0UICAFiBFtKOUTUSvfKRXQrTWARqjUyrrJr2wWsmq5W2\nmum6pXO1/Gpu0UWf1jMfJ7OmjGt0jgwzgFqycZDRasHSC22DvrPsddcFmcovJmkv3TAGmblfRSgi\nT5PlVU9usXCt6Bv9Rel76I2NH9rVepL+W5H+YE3h5O9VCJSEFhIAwAokJACAFeiyA7AjGNSw72gh\nAQCsQAtpV73Qa9qZVUr12Y7PTLfPMxVvfdzBOCn6n8m5OtlobVToHutDFZI75g3Tyq09baJuZsBF\nNsIrQarY3mgXZLdtvZYaHqJq/tpIjvSjEG1EQC95ijFIdUcr3Sc2Hcuw+XgB4xAVFduJPqzFNEaj\nl4y8aJ4s4gN90MV8YRjMcJa8+OBWAy+AG6KFBACwAgkJAGAFEhIAwArUkHaMqpHMx2ldQVUqjPNP\n1aZto5P0lhcnV8sOa2WepDKiV6o2pBe0RtGr+ECvyqjihV7oGmSKMcYy2NqE34bhyg3ppaZonrcz\nYToRVdMb5dWBtKKUVoopszyjgjT+vNQGhk+j9Megfgpn2o9mwzJhbRhlt+9oIQEArEBCAgBYgYQE\nALACNaR9kJ3EMzg31BLeaOWQR7nFIVWKyM5w2oq6vddMI5wnxZ+BdjK7s9xahKZCTnrv7VYCVYUT\nVUwyTdHZolJlXAC3tf0EqfRe7Q1GpmVzFf1HLJlC10ttXVcVpO11IxwSEhKAHfE9gxr2HF12AAAr\nkJAAAFagy25XrU05UvvXmTaeU+WZ/LqR7palIwOtnGLerC9hDFLVtPRZQfmvcwPGaoqa+TQwfeB6\nDOrKG6ynt7ktSj6mIFP5uxECdaCFBACwAgkJAGAFuuwA7AiWDtp3tJAAAFaghbRfFoZSdVqNN2wm\nl+oVPUbgelq01ga5+dCJwgdZAIeGFhIAwAokpKoFQaCOwzCsMRIAsApddtUZDoeu64pI/L/T6TQI\nAsdxXNdtt9t1RwdYj0EN+46EVJHZbNZutz3PU2eCIOj3+yIyHA5JSABAQqpI3BgKgsB13TgtOY4T\nf0vPUgDqFYah+ttExUhIFQmCoNPpeJ43HA43yUDL5VLktXbisci98sIDTN6KvE2Ov6kzkKrE/epx\nTup0OnWHc3AY1FCRbMNIjWjwfT97/dHRkcgT7T+yEar3UPsN/Ie6g6mC53mdTqfb7Rr/KlE2WkgV\n6XQ6alBDzHXd4XDoOA5ddsBG/s/X8j45PnLk77f8w1mG8jpNM8vHy+wlcTXX933KurUgIVXEcZx+\nv+/7vko/nU4nbiTRYQ1s5H99Jr+4RZ44cvTbjz56abzK9/24g/3mD8JNkZAqdaUxRCoCrEI2qhcJ\nCQBERHzfPz099TzP9/0wDC8uLuqO6OCQkABARMTzvMvLy7qjOGiMsgMAWIEWEoAd8b3Id8W92oPi\nXgoFoYUEALACCQkAYAUSEgDACiQkAIAVGNQAYEd8X+h+SAxqsA8tJACAFUhIAAArkJAAAFYgIQEA\nrMCgBgA74n2hgxpgHxLS/ptLLz4wrl08kFGVwQDAh9BlBwCwAgkJAGAFEhIAwAokJACAFRjUsFdG\n0avsyRdysjpqjHNu6TU+Ly2uLQySIRiMtsBV70U+rjsGlIkWEgDACiQkAIAVSEgAACtQQ9oHxtJR\n1jiaZ0+eyXl8EM176mSjVUX9xhj2uayCHDRa6uSjaFUGs6TQBaAMJCQAO+J7kY/qjgFlossOAGAF\nEhIAwAp02e2DN8kEI1VrMXqpVWWUF6ZbVD2pjGKSmmn0ppGefJYcHMvqvTS1uNT8qeoLXYpe8aKU\nBZSBFhIAwAq0kADsiPcif607BpSJFhIAwAq0kHZVMzpWx88ayRZ8STFJXwhOFT+aeQUmM7W5n4i0\nbrG4nDHaGxhr76HUQpeiPoF0SUD9pKkCR4UJuBlaSAAAK5CQAABWICHVKQzDukMAAFtQQ6rIF198\n4XmeiIRheHFxISLT6TQIAsdxXNdtt9t1BwhY773Ij3XHgDKRkCpyeXkZH0wmk/ggCIJ+vy8iw+Hw\nBglp3ktnuS6SAzW9dGC65dywP58MxkkFXqvPH5u28rsN45zcC+34WfakFsKjZNSCHthCyqIP5VDe\nmD6TczlTx2qZWgA3Q0Kq1Gw2i9tJIuI4TnygzuiWy6XIa+3EY5F7pccHrHkr8jY5/qbOQKoynU59\n32+323Ra1IKEVCnf9+NW0bWOjo5EnpQdD5DrocjDumOoVKfTUf9SRPUY1FAbNaLB9/16IwEAG9BC\nqk7cFaC+dF13OBw6jmPssrvWQqtoqLrLM9OVqu5inBg7GF+9rEDGYky+/PVhr33KbWbvKsbilu6p\naatDlO791+lxw5HGtn84ofw1/cffcrksJioUh4RUnSuJp9PpxI0kugiAzXwmP7lNacfRbz86enn7\ngFAsElKdSEWAVSaTSdyFHgRBt9utO5yDQ0ICgBWSUL1ISHvl4vpL6nHLwE56q/pWeXOPdGfGMtZ1\nJTYVZK+IOhZwgEhIAHbEj38R+aHuIFAihn0DAKxAQgIAWIGEBACwAjWkHaO2f9U3MH2UHLxJd4xN\n54qqAQXpvrDvNwAABUlJREFUOqAa49xV88TPxoutQhVtmqq+Y6yir7iajnowr2G6YpymWshk2LUX\nHKVvXy1iqz/6WRL52mCN7TfkBaCjhQQAsAItJAC74r3Iu7pjQIloIQEArEALacf0Gqv99IxVmbPI\ncIsqLBk36DNWZUpl3Kwvu0GfMbBqZv4utGqZcSquqlrpFbhF8YvTAoeFFhIAwAokJACAFeiyA7Ar\nGNSw50hIu0qvxKhpQ8ZN9l4kW971xul8ndFJUvxIbtHnHs0XhjJPo6DpPsZJTmfJXCM1p0qf9HRu\nmkO12H5e1A3kT3IqfAoUcMjosgMAWIGEBACwAgkJAGAFaki7Sq9eREnJZ/P961RV5jxaHaxND0qq\nPI1WMTUSvd5jnEGl5E+0qqZuBFsxqGHP0UICAFiBhAQAsAIJCQBgBRISAMAKDGrYB9mhB/oGfYNk\nmICINm32ZDVkQI1lWJvjaZgXWxjjwISefL7JZQD2GAkJwK74nlF2+40uOwCAFUhIAAAr0GW3nwbX\nLfqpKjQsDwrAErSQAABWoIUEYFewdNCeo4UEALACCQkAYAUSEgDACiQkAIAVGNRQnTAMRcRxHP2M\n/iWAXAxq2HMkpIpMp9M4Ibmu22634zNBEDiOo85U7o8ij3mKlQ/ap6cAm6LLriK+73e73W63O5vN\n4jNBEPT7/W636/t+TUH9mafY+qB9esqOif/hiFrQQqqI53nT6VREVGNIddZ5npe9frlcirzWTjwW\nuVd2kMC6tyJvk+Nv6gykKhb0Wxw0ElJFgiCIf783bA8dHR2JPCk5KCDfQ5GHdcdQqbjfQkSGwyEJ\nqXokpIq4rhu3hIIgiM+ongHf9z/wq//7koP6r5Jff/+eUtmDLHzKH8uKYgv/qTXUPhb52y1v/3/6\nW/7DH5rZK/L7LVA2ElJ1JpOJaL/xrusOh0PHcYy/+l9++eWXX1YaHmC5KPq3ukNAuUhIFel0OleG\nfWfPAKjXBv0WKFEjiqK6YwAAKzCooV4kJABI0W9RIxISAMAKTIwFAFiBhHToSpqXnn3ZMh4UhuGV\nl61mmv0ePEVNPyj7QcDmGGV3QKbTaTx2SFVrSyrhDodD13XjpWM7nU5JD4pXvij7KfEjnj171u12\ny1uE8PT0ND5QP53yquvxT0dE4v+ljA+LRDgkl5eXl5eX6svnz59fOSjqKfFBv98v9UHVPOX58+f6\n51bGU9Rb0B9a+FOiKLq8vPzd735XwYOAG6CFdNBKmpeuFknacOG+G5vNZrPZLG4elfQU3/evvFoZ\nT3EcZzKZhGHoum78dkr6xOLGUBAEaukQ1iaAPaghoRS+76vl+8rTbrf7/b5aQL0Ms9msgo6seN33\nfr+vl3bKEKeiTqdT6ocG3AwJ6aDp89ILfNk4G6mGS3kPksyGh4U/xfO82WwWv6PynqKoFy/pKdmG\nUalvB9gKXXYHZDKZxP+nEwRBt9uV69bTuxnf909PTz3P830/DMOLi4uSHjSZTBzHibu54jNlPKXd\nbgdBMJ1ON1yE8GbiTywMQ/WaZTxFRDqdjhrUUOqDgBtgYuyhq2xeehkPiisi2UZS2W9n19/LlcIY\naxPAEiQkAIAVqCEBAKxAQgIAWIGEBACwAgkJAGAFEhIAwAokJACAFUhIAAArkJAAAFYgIQEArEBC\nAgBYgYQEALACCQkAYAUSEgDACiQkAIAVSEgAACuQkAAAViAhAQCsQEICAFiBhAQAsAIJCQBgBRIS\nAMAKJCQAgBVISAAAK5CQAABWICEBAKxAQgIAWIGEBACwAgkJAGAFEhIAwAokJACAFUhIAAAr/Deu\n+4A4nGehygAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "path_sc8 = [opt.folder_out filesep 'stability_group' filesep 'sci10_scg7_scf8' filesep];\n", "parcel_sc8 = [path_sc8 'brain_partition_consensus_group_sci10_scg7_scf8.nii.gz'];\n", "[hdr,sc8] = niak_read_vol(parcel_sc8);\n", "niak_montage(sc8)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "On this small dataset with big voxels, the parcellation does not look good. But that is enough to illustrate how to run the pipeline... Other (longer) tutorials will show how to explore the results on a dataset with smaller voxels and more subjects. " ] } ], "metadata": { "kernelspec": { "display_name": "Octave", "language": "octave", "name": "octave" }, "language_info": { "file_extension": ".m", "help_links": [ { "text": "GNU Octave", "url": "https://www.gnu.org/software/octave/support.html" }, { "text": "Octave Kernel", "url": "https://github.com/Calysto/octave_kernel" }, { "text": "MetaKernel Magics", "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" } ], "mimetype": "text/x-octave", "name": "octave", "version": "4.0.2" } }, "nbformat": 4, "nbformat_minor": 2 }