%QoMGroups2018


%  Load data
%a = mcread('A0001.tsv');
b = mcread('B0001.tsv');
c = mcread('C0001.tsv');
d = mcread('D0001.tsv');
e = mcread('E0001.tsv');
f = mcread('F0001.tsv');
g = mcread('G0001.tsv');
h = mcread('H0001.tsv');
i = mcread('I0001.tsv');
j = mcread('J0001.tsv');
k = mcread('K0001.tsv');
l = mcread('L0001.tsv');
m = mcread('M0001.tsv');

% Calculate QoM for each marker (this was used in NM!)
%a_dx = mctimeder(a,1,[1 0.99999]);
%a_dx_norm = mcnorm(a_dx);
%a_dx_norm_mean = mcmean(a_dx_norm)';
%a_dx_norm_std = mcstd(a_dx_norm)';
%a_cumdis = mccumdist(a);


b_dx = mctimeder(b,1,[1 0.99999]);
b_dx_norm = mcnorm(b_dx);
b_dx_norm_mean = mcmean(b_dx_norm)';
b_dx_norm_std = mcstd(b_dx_norm)';
b_cumdis = mccumdist(b);


c_dx = mctimeder(c,1,[1 0.99999]);
c_dx_norm = mcnorm(c_dx);
c_dx_norm_mean = mcmean(c_dx_norm)';
c_dx_norm_std = mcstd(c_dx_norm)';
c_cumdis = mccumdist(c);


d_dx = mctimeder(d,1,[1 0.99999]);
d_dx_norm = mcnorm(d_dx);
d_dx_norm_mean = mcmean(d_dx_norm)';
d_dx_norm_std = mcstd(d_dx_norm)';
d_cumdis = mccumdist(d);


e_dx = mctimeder(e,1,[1 0.99999]);
e_dx_norm = mcnorm(e_dx);
e_dx_norm_mean = mcmean(e_dx_norm)';
e_dx_norm_std = mcstd(e_dx_norm)';
e_cumdis = mccumdist(e);


f_dx = mctimeder(f,1,[1 0.99999]);
f_dx_norm = mcnorm(f_dx);
f_dx_norm_mean = mcmean(f_dx_norm)';
f_dx_norm_std = mcstd(f_dx_norm)';
f_cumdis = mccumdist(f);


g_dx = mctimeder(g,1,[1 0.99999]);
g_dx_norm = mcnorm(g_dx);
g_dx_norm_mean = mcmean(g_dx_norm)';
g_dx_norm_std = mcstd(g_dx_norm)';
g_cumdis = mccumdist(g);


h_dx = mctimeder(h,1,[1 0.99999]);
h_dx_norm = mcnorm(h_dx);
h_dx_norm_mean = mcmean(h_dx_norm)';
h_dx_norm_std = mcstd(h_dx_norm)';
h_cumdis = mccumdist(h);


i_dx = mctimeder(i,1,[1 0.99999]);
i_dx_norm = mcnorm(i_dx);
i_dx_norm_mean = mcmean(i_dx_norm)';
i_dx_norm_std = mcstd(i_dx_norm)';
i_cumdis = mccumdist(i);


j_dx = mctimeder(j,1,[1 0.99999]);
j_dx_norm = mcnorm(j_dx);
j_dx_norm_mean = mcmean(j_dx_norm)';
j_dx_norm_std = mcstd(j_dx_norm)';
j_cumdis = mccumdist(j);


k_dx = mctimeder(k,1,[1 0.99999]);
k_dx_norm = mcnorm(k_dx);
k_dx_norm_mean = mcmean(k_dx_norm)';
k_dx_norm_std = mcstd(k_dx_norm)';
k_cumdis = mccumdist(k);


l_dx = mctimeder(l,1,[1 0.99999]);
l_dx_norm = mcnorm(l_dx);
l_dx_norm_mean = mcmean(l_dx_norm)';
l_dx_norm_std = mcstd(l_dx_norm)';
l_cumdis = mccumdist(l);


m_dx = mctimeder(m,1,[1 0.99999]);
m_dx_norm = mcnorm(m_dx);
m_dx_norm_mean = mcmean(m_dx_norm)';
m_dx_norm_std = mcstd(m_dx_norm)';
m_cumdis = mccumdist(m);


%QoM without stands

%QoMA = a_dx_norm_mean(3:5);
%QoMB = b_dx_norm_mean(1:6);
%QoMC = c_dx_norm_mean(1:9);
%QoMD = d_dx_norm_mean(1:8);
%QoME = e_dx_norm_mean(1:4);
%QoMF = f_dx_norm_mean(1:13);
%QoMG = g_dx_norm_mean(1:4);
%QoMH = h_dx_norm_mean(1:3);
%QoMI = i_dx_norm_mean(1:6);
%QoMJ = j_dx_norm_mean(1:3);
%QoMK = k_dx_norm_mean(1:5);
%QoML = l_dx_norm_mean(1:3);