% vbm 

datadir='c:/keith/fMRI/francesco/NMTBI_BLUR10_data/';
d=dir([datadir '*_wm.mnc']);
n=size(d,1)
input_files=[datadir d(1).name];
for i=2:n
   input_files=char(input_files, [datadir d(i).name]);
end
input_files

input_files=input_files([1:25 27:36 38],:)
n=n-2

trauma=[1 1 0 0 0 0 1 0 0 1 1 1 1 0 1 1 1 0 1 0 1 0 1 0 0 0 0 0 1 0 0 0 1 1 1 0];
days_coma=[43 30 0 0 0 0 28 0 0 30 15 28 41 0 60 45 13   

X=[ones(n,1) trauma']
mvlm(input_files, X, contrast, output_file_base, which_stats)




% dbm 

datadir='c:/keith/fMRI/francesco/DXYZ_CONTROLLO/';
d=dir([datadir '*_dx.mnc']);
n_control=size(d,1)
input_files=[datadir d(1).name];
for i=2:n_control
   input_files=char(input_files, [datadir d(i).name]);
end
datadir='c:/keith/fMRI/francesco/DXYZ_TRAUMA/';
d=dir([datadir '*_dx.mnc'])
n_trauma=size(d,1)
for i=1:n_trauma
   input_files=char(input_files, [datadir d(i).name]);
end
input_files


saveas(gcf,'c:/keith/fmristat/figs/figvbmt.jpg');