design_test

 design_test.m

0001 % design_test.m
0002 
0003 global rootpath
0004 
0005 rootpath = '/data1/design_tests/';
0006 
0007 %
0008 %  Create data vector
0009 %
0010 
0011 means = [5 2 8 1];
0012 ncond = length(means);
0013 nruns = 2;
0014 
0015 npts_per_cond = 20;
0016 
0017 npts = ncond*nruns*npts_per_cond;
0018 
0019 data = repmat(ones(npts_per_cond,1)*means,[1 1 nruns]);
0020 data = data(:);
0021 
0022 noise_scaling = 4;
0023 noise = rand(size(data));
0024 noise = noise - mean(noise);
0025 noise = noise*noise_scaling;
0026 
0027 noisy = noise+data;
0028 
0029 plot(noisy)
0030 
0031 %
0032 %  Create images
0033 %
0034 
0035 spm_bch('design_batch');
0036 
0037 load SPM_fMRIDesMtx
0038