In [7]:

    

cd /sandbox/home/Desktop
load adni_dartel_ecogpt.mat


    

In [12]:

    

R=corr(stack)
hier = niak_hierarchical_clustering(R);
order = niak_hier2order(hier);


    

    




R =

   1.00000   0.69094   0.59232   0.62096   0.61504   0.67857   0.85482
   0.69094   1.00000   0.69137   0.67308   0.68500   0.77445   0.90237
   0.59232   0.69137   1.00000   0.73532   0.65905   0.64257   0.82524
   0.62096   0.67308   0.73532   1.00000   0.75856   0.62343   0.82940
   0.61504   0.68500   0.65905   0.75856   1.00000   0.70218   0.83806
   0.67857   0.77445   0.64257   0.62343   0.70218   1.00000   0.85247
   0.85482   0.90237   0.82524   0.82940   0.83806   0.85247   1.00000

     Percentage done : 0 17 33 50 67 83 100 Done ! 

In [13]:

    

hier


    

    




hier =

    0.90237    2.00000    7.00000    8.00000
    0.78383    8.00000    6.00000    9.00000
    0.75856    4.00000    5.00000   10.00000
    0.69062    1.00000    9.00000   11.00000
    0.67673    3.00000   10.00000   12.00000
    0.54367   11.00000   12.00000   13.00000

In [15]:

    

niak_visu_matrix(R(order,order))


    

    




ans =  1






    

In [36]:

    

y = niak_normalize_tseries(stack');
D = niak_build_distance(y);
hier = niak_hierarchical_clustering(-D);
order = niak_hier2order(hier);


    

    




     Percentage done : 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 Done ! 

In [39]:

    

niak_visu_matrix(1-D(order,order)/(2*sqrt(6)))


    

    




ans =  1






    

In [40]:

    

[sil,intra,inter] = niak_build_avg_silhouette(1-D(order,order)/(2*sqrt(6)),hier);


    

    




     Percentage done : 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 Done ! 

In [41]:

    

plot(sil)


    

In [42]:

    

part = niak_threshold_hierarchy(hier,struct('thresh',10));


    

In [44]:

    

niak_visu_part(part(order))
figure
niak_visu_matrix(1-D(order,order)/(2*sqrt(6)))


    

    




ans =  2






    













    

In [55]:

    

errorbar(mean(stack(part==3,:),1)',std(stack(part==9,:),[],1)')


    

In [59]:

    

x=stack(part==7,:)';
plot(repmat((1:size(x,1))',[1 size(x,2)])+0.3*rand(size(x)),x,'.')


    

In [60]:

    

4^7


    

    




ans =  16384