In [4]:
cd
build_path niak psom
cd /home/pbellec/git/niak_tutorials/glm_connectome
Adding library niak to the search path.

Adding library psom to the search path.

GLM connectome pipeline - intra-run

This tutorial shows how to use the NIAK GLM connectome pipeline to compare two groups of subjects, 10 patients diagnosed with schizophrenia vs 10 healthy controls. We will illustrate how to include an intra-run contrast, i.e. manipulate the fMRI time series of each subject prior to generate connectivity measures. See the documentation of the pipeline for a more detailed description of the options. There is also inline documentation you can access by typing help niak_pipeline_glm_connectome. Download the tutorial as a notebook here and a matlab script here. To run this tutorial, we recommend to use jupyter from a niak docker container, as described in the NIAK installation page.

Preparing the files

Functional MRI

First download a small preprocessed fMRI dataset.

In [5]:
clear
path_data = [pwd filesep];
[status,msg,data_fmri] = niak_wget('cobre_lightweight20_nii');
--2017-02-01 12:08:21--  https://www.nitrc.org/frs/download.php/9401/cobre_lightweight20.zip
Resolving www.nitrc.org (www.nitrc.org)... 132.239.16.23
Connecting to www.nitrc.org (www.nitrc.org)|132.239.16.23|:443... connected.
HTTP request sent, awaiting response... 200 OK
Length: 89902898 (86M) [application/force-download]
Saving to: ‘/home/pbellec/git/niak_tutorials/glm_connectome/cobre_lightweight20/cobre_lightweight20.zip’

cobre_lightweight20 100%[===================>]  85.74M  2.29MB/s    in 44s     

2017-02-01 12:09:05 (1.97 MB/s) - ‘/home/pbellec/git/niak_tutorials/glm_connectome/cobre_lightweight20/cobre_lightweight20.zip’ saved [89902898/89902898]

Deleting file '/home/pbellec/git/niak_tutorials/glm_connectome/cobre_lightweight20/cobre_lightweight20.zip' 
ans = 0

When starting from fMRI data preprocessed with NIAK, it is possible to use niak_grab_fmri_preprocess on the output folder to collect the file names, as described in the pipeline documentation. In this case, we explicitely list all the files

In [6]:
file_pheno = [data_fmri.path filesep 'phenotypic_data.tsv.gz'];
tab = niak_read_csv_cell(file_pheno);
list_subject = tab(2:end,1);
files_in = struct;
for ss = 1:length(list_subject)
    files_in.fmri.(list_subject{ss}).sess1.rest = [data_fmri.path filesep 'fmri_' list_subject{ss} '.nii.gz'];
end

Brain parcellation

The second input of the pipeline is a set of brain parcels. We will just download the so-called Cambridge functional template.

In [8]:
[status,msg,data_template] = niak_wget('cambridge_template_mnc1');
--2017-02-01 12:16:23--  http://files.figshare.com/1861821/template_cambridge_basc_multiscale_mnc_sym.zip
Resolving files.figshare.com (files.figshare.com)... 54.194.135.36, 52.211.238.180
Connecting to files.figshare.com (files.figshare.com)|54.194.135.36|:80... connected.
HTTP request sent, awaiting response... 301 Moved Permanently
Location: https://ndownloader.figshare.com/files/1861821 [following]
--2017-02-01 12:16:23--  https://ndownloader.figshare.com/files/1861821
Resolving ndownloader.figshare.com (ndownloader.figshare.com)... 52.209.81.33, 52.49.140.44
Connecting to ndownloader.figshare.com (ndownloader.figshare.com)|52.209.81.33|:443... connected.
HTTP request sent, awaiting response... 302 Found
Location: https://s3-eu-west-1.amazonaws.com/pfigshare-u-files/1861821/template_cambridge_basc_multiscale_mnc_sym.zip [following]
--2017-02-01 12:16:24--  https://s3-eu-west-1.amazonaws.com/pfigshare-u-files/1861821/template_cambridge_basc_multiscale_mnc_sym.zip
Resolving s3-eu-west-1.amazonaws.com (s3-eu-west-1.amazonaws.com)... 54.231.131.139
Connecting to s3-eu-west-1.amazonaws.com (s3-eu-west-1.amazonaws.com)|54.231.131.139|:443... connected.
HTTP request sent, awaiting response... 200 OK
Length: 211786 (207K) [binary/octet-stream]
Saving to: ‘/home/pbellec/git/niak_tutorials/glm_connectome/template_cambridge_basc_multiscale_mnc_sym/template_cambridge_basc_multiscale_mnc_sym.zip’

template_cambridge_ 100%[===================>] 206.82K   347KB/s    in 0.6s    

2017-02-01 12:16:25 (347 KB/s) - ‘/home/pbellec/git/niak_tutorials/glm_connectome/template_cambridge_basc_multiscale_mnc_sym/template_cambridge_basc_multiscale_mnc_sym.zip’ saved [211786/211786]

Deleting file '/home/pbellec/git/niak_tutorials/glm_connectome/template_cambridge_basc_multiscale_mnc_sym/template_cambridge_basc_multiscale_mnc_sym.zip' 
ans = 0

We are going to pick the parcellation into 36 clusters.

In [10]:
files_in.network = [data_template.path filesep 'template_cambridge_basc_multiscale_sym_scale007.mnc.gz'];

Group model

We first load the phenotypic data from the COBRE sample. This command will not work if you haven't downloaded the data yet as part of the rmap tutorial.

In [ ]:
path_cobre = [pwd filesep 'cobre_lightweight20'];
file_pheno = [path_cobre filesep 'phenotypic_data.tsv.gz'];
tab = niak_read_csv_cell(file_pheno);

Now we convert the values into a series of numerical covariates, that we save in a separate .csv file:

In [ ]:
list_subject = tab(2:end,1);
patient = strcmp(tab(2:end,5),'Patient');
age = str2double(tab(2:end,2)); 
FD = str2double(tab(2:end,9));
opt_csv.labels_x = list_subject; % Labels for the rows
opt_csv.labels_y = { 'age' , 'patient' , 'fd' };
niak_write_csv('model_patient.csv', [age patient FD] , opt_csv);

We specify to the pipeline where to find the model:

In [ ]:
files_in.model.group = [pwd filesep 'model_patient.csv'];

Connectivity maps

Now we are going to get a list of the connectivity maps associated with each subject for one network, say the DMN. Labels for each network have been specified when running the connectome pipeline. We simply grab the outputs of the connectome pipeline.

In [ ]:
path_connectome = [pwd filesep 'connectome'];
files_conn = niak_grab_connectome(path_connectome);
files_in.data = files_conn.rmap;

Brain mask

We specify the mask of brain networks to the pipeline, so that it can use it to mask the grey matter.

In [ ]:
files_in.mask = files_conn.network_rois;

Set up the options of the pipeline

First specify where to save the outputs, and how many networks to use:

In [ ]:
%% General
opt.folder_out = [pwd filesep 'subtype'];

Then specify which covariates to use as confounds before the generation of subtypes.

In [ ]:
% a list of variable names to be regressed out
% If unspecified or left empty, no confounds are regressed
opt.stack.regress_conf = {'fd'};

The options for the subtypes themselves:

In [ ]:
%% Subtyping
opt.subtype.nb_subtype = 2;        % the number of subtypes to extract
opt.subtype.sub_map_type = 'mean'; % the model for the subtype maps (options are 'mean' or 'median')

Now we add an association test between subtypes and the patient label:

In [ ]:
% scalar number for the level of acceptable false-discovery rate (FDR) for the t-maps
opt.association.patient.fdr = 0.05;                           
% turn on/off normalization of covariates in model (true: apply / false: don't apply)
opt.association.patient.normalize_x = false;                   
% turn on/off normalization of all data (true: apply / false: don't apply)
opt.association.patient.normalize_y = false;                  
% turn on/off adding a constant covariate to the model
opt.association.patient.flag_intercept = true;     
% To test a main effect of a variable
opt.association.patient.contrast.patient = 1; % scalar number for the weight of the variable in the contrast
opt.association.patient.contrast.fd = 0;      % scalar number for the weight of the variable in the contrast
opt.association.patient.contrast.age = 0;     % scalar number for the weight of the variable in the contrast
% type of data for visulization (options are 'continuous' or 'categorical')
opt.association.patient.type_visu = 'continuous';

It is also possible to add a single chi-square test on the relationship between subtypes and a categorical variable:

In [ ]:
% string name of the column in files_in.model on which the contigency table will be based
opt.chi2 = 'patient';

Run the pipeline

In [ ]:
opt.flag_test = false;  % Put this flag to true to just generate the pipeline without running it.
pipeline = niak_pipeline_subtype(files_in,opt);