Part-2: Demonstration of two procedures

Contents

---

1. Load Data

2. Successive Normalization Example

3. Nuisance Correlation Example

%Add path
run('setup.m');
datadir = fullfile(basedir,'data');
filename = 'ABIDE_controlData_7Sites';
load(fullfile(datadir,filename));
load(fullfile(datadir,'cc200Map2Communities'))
max_recursion_depth(1024)

Load subject data from an ABIDE site

disp('Fields in dataset')
fieldnames(data{1})
disp('signals contains subjects x brain regions x time-series')
disp(size(data{1}.signals))

subjectno=1;
dataset=1;
[~,reorder_communities] = sort(mapToCommunities,'ascend');
X = squeeze(data{dataset}.signals(subjectno,reorder_communities,:))';
Y = mean(X,2);

Fields in dataset
ans = 
{
  [1,1] = signals
  [2,1] = subIDs
  [3,1] = age
  [4,1] = gender
  [5,1] = dataName
}
signals contains subjects x brain regions x time-series
    15   197   246

Demonstration of Successive Normalization

results = demo_successive_norm(X);

ans = 0

Demonstration of Nuisance Correlation Matrix with Global Signal Regression

results = demo_conditional_correlation(X,Y);

ans = 0
ans = tmp/22-Jun-2017-2125/sample_conditional_correlation
GPL Ghostscript 9.16: **** Could not open the file tmp/22-Jun-2017-2125/sample_conditional_correlation.pdf .