do_dmatinfo8

L. Bottemanne

·

Follow

3 min read

·

Mar 13, 2018

--

Share

function dmatinfo = do_dmatinfo8(dmat, ddmpath)

% size(dmat.output.Minimum, 1)% number of condition

a1 = vertcat (dmat.dm_a_1);
a2 = vertcat (dmat.dm_a_2);
a3 = vertcat (dmat.dm_a_3);
a4 = vertcat (dmat.dm_a_4);
a5 = vertcat (dmat.dm_a_5);
a6 = vertcat (dmat.dm_a_6);
a7 = vertcat (dmat.dm_a_7);
a8 = vertcat (dmat.dm_a_8);

dmatinfo.a = horzcat (a1, a2, a3, a4, a5, a6, a7, a8); 

Ter1 = vertcat (dmat.dm_Ter_1);
Ter2 = vertcat (dmat.dm_Ter_2);
Ter3 = vertcat (dmat.dm_Ter_3);
Ter4 = vertcat (dmat.dm_Ter_4);
Ter5 = vertcat (dmat.dm_Ter_5);
Ter6 = vertcat (dmat.dm_Ter_6);
Ter7 = vertcat (dmat.dm_Ter_7);
Ter8 = vertcat (dmat.dm_Ter_8);

dmatinfo.Ter = horzcat (Ter1, Ter2, Ter3, Ter4, Ter5, Ter6, Ter7, Ter8); 

eta1 = vertcat (dmat.dm_eta_1);
eta2 = vertcat (dmat.dm_eta_2);
eta3 = vertcat (dmat.dm_eta_3);
eta4 = vertcat (dmat.dm_eta_4);
eta5 = vertcat (dmat.dm_eta_5);
eta6 = vertcat (dmat.dm_eta_6);
eta7 = vertcat (dmat.dm_eta_7);
eta8 = vertcat (dmat.dm_eta_8);

dmatinfo.eta = horzcat (eta1, eta2, eta3, eta4, eta5, eta6, eta7, eta8);

z1 = vertcat (dmat.dm_z_1);
z2 = vertcat (dmat.dm_z_2);
z3 = vertcat (dmat.dm_z_3);
z4 = vertcat (dmat.dm_z_4);
z5 = vertcat (dmat.dm_z_5);
z6 = vertcat (dmat.dm_z_6);
z7 = vertcat (dmat.dm_z_7);
z8 = vertcat (dmat.dm_z_8);

dmatinfo.z = horzcat (z1, z2, z3, z4, z5, z6, z7, z8); 

sz1 = vertcat (dmat.dm_sz_1);
sz2 = vertcat (dmat.dm_sz_2);
sz3 = vertcat (dmat.dm_sz_3);
sz4 = vertcat (dmat.dm_sz_4);
sz5 = vertcat (dmat.dm_sz_5);
sz6 = vertcat (dmat.dm_sz_6);
sz7 = vertcat (dmat.dm_sz_7);
sz8 = vertcat (dmat.dm_sz_8);

dmatinfo.sz = horzcat (sz1, sz2, sz3, sz4, sz5, sz6, sz7, sz8); 

st1 = vertcat (dmat.dm_st_1);
st2 = vertcat (dmat.dm_st_2);
st3 = vertcat (dmat.dm_st_3);
st4 = vertcat (dmat.dm_st_4);
st5 = vertcat (dmat.dm_st_5);
st6 = vertcat (dmat.dm_st_6);
st7 = vertcat (dmat.dm_st_7);
st8 = vertcat (dmat.dm_st_8);

dmatinfo.st = horzcat (st1, st2, st3, st4, st5, st6, st7, st8); 

v1 = vertcat (dmat.dm_v_1);
v2 = vertcat (dmat.dm_v_2);
v3 = vertcat (dmat.dm_v_3);
v4 = vertcat (dmat.dm_v_4);
v5 = vertcat (dmat.dm_v_5);
v6 = vertcat (dmat.dm_v_6);
v7 = vertcat (dmat.dm_v_7);
v8 = vertcat (dmat.dm_v_8);

dmatinfo.v = horzcat (v1, v2, v3, v4, v5, v6, v7, v8);pi1 = vertcat (dmat.dm_pi_1);
pi2 = vertcat (dmat.dm_pi_2);
pi3 = vertcat (dmat.dm_pi_3);
pi4 = vertcat (dmat.dm_pi_4);
pi5 = vertcat (dmat.dm_pi_5);
pi6 = vertcat (dmat.dm_pi_6);
pi7 = vertcat (dmat.dm_pi_7);
pi8 = vertcat (dmat.dm_pi_8);

dmatinfo.pi = horzcat (pi1, pi2, pi3, pi4, pi5, pi6, pi7, pi8);

g1 = vertcat (dmat.dm_g_1);
g2 = vertcat (dmat.dm_g_2);
g3 = vertcat (dmat.dm_g_3);
g4 = vertcat (dmat.dm_g_4);
g5 = vertcat (dmat.dm_g_5);
g6 = vertcat (dmat.dm_g_6);
g7 = vertcat (dmat.dm_g_7);
g8 = vertcat (dmat.dm_g_8);

dmatinfo.g = horzcat (g1, g2, g3, g4, g5, g6, g7, g8); 

% moyennes des paramètres estimés par condition - null
DM1 = horzcat (mean (a1), mean(Ter1), mean(eta1), mean (z1), mean (sz1), mean (st1), mean (v1), mean (pi1),mean (g1));
DM2 = horzcat (mean (a2), mean(Ter2), mean(eta2), mean (z2), mean (sz2), mean (st2), mean (v2), mean (pi2),mean (g2));
DM3 = horzcat (mean (a3), mean(Ter3), mean(eta3), mean (z3), mean (sz3), mean (st3), mean (v3), mean (pi3),mean (g3));
DM4 = horzcat (mean (a4), mean(Ter4), mean(eta4), mean (z4), mean (sz4), mean (st4), mean (v4), mean (pi4),mean (g4));
DM5 = horzcat (mean (a5), mean(Ter5), mean(eta5), mean (z5), mean (sz5), mean (st5), mean (v5), mean (pi5),mean (g5));
DM6 = horzcat (mean (a6), mean(Ter6), mean(eta6), mean (z6), mean (sz6), mean (st6), mean (v6), mean (pi6),mean (g6));
DM7 = horzcat (mean (a7), mean(Ter7), mean(eta7), mean (z7), mean (sz7), mean (st7), mean (v7), mean (pi7),mean (g7));
DM8 = horzcat (mean (a8), mean(Ter8), mean(eta8), mean (z8), mean (sz8), mean (st8), mean (v8), mean (pi8),mean (g8));

dmatinfo.DMpar = vertcat (DM1, DM2, DM3, DM4, DM5, DM6, DM7, DM8);clear DM1 DM2 DM3 DM4 DM5 DM6 DM7 DM8...
a1 a2 a3 a4 a5 a6 a7 a8 ...
Ter1 Ter2 Ter3 Ter4 Ter5 Ter6 Ter7 Ter8 ...
eta1 eta2 eta3 eta4 eta5 eta6 eta7 eta8...
z1 z2 z3 z4 z5 z6 z7 z8 ...
sz1 sz2 sz3 sz4 sz5 sz6 sz7 sz7...
st1 st2 st3 st4 st5 st6 st7 st8 ...
g1 g2 g3 g4 g5 g6 g7 g8 ...
pi1 pi2 pi3 pi4 pi5 pi6 pi7 pi8...
v1 v2 v3 v4 v5 v6 v7 v8;

% model count
cd ([ddmpath 'ddmMean\']);
filesTables = ls ;
nbM = size(filesTables,1) - 2;
disp('model number'); disp(nbM);

% save in xls files
filename = sprintf('dm_%d', nbM);
xlswrite(filename, dmatinfo.DMpar);
clear filename;
% boundary a
filename  = sprintf('dmA_%d', nbM);
xlswrite([ddmpath 'par\' filename] , dmatinfo.a); 
clear filename;
% drift v
filename  = sprintf('dmV_%d', nbM);
xlswrite([ddmpath 'par\' filename]  , dmatinfo.v);
clear filename;
% non-decisional time Ter
filename  = sprintf('dmTer_%d', nbM);
xlswrite([ddmpath 'par\' filename], dmatinfo.Ter);
clear filename;
% model fitting pi
filename  = sprintf('dmPi_%d', nbM);
xlswrite([ddmpath 'par\' filename], dmatinfo.pi);
clear filename;
% bic 

dmatinfo.bic = vertcat(dmat.BIC);
bic = dmatinfo.bic ;
filename = sprintf('BIC_mod%d', nbM);
xlswrite([ddmpath 'FitInfo\' filename], bic); % dmatinfo.bic);
clear filename;

dmatinfo.chi = vertcat(dmat.chi);

filename = sprintf('Chi2_mod%d', nbM);
xlswrite([ddmpath 'FitInfo\' filename], dmatinfo.chi  );
clear filename;

clear nbM;
%%% the end