% Change the directory to where the main model is,
% abcdAll_noA\framework\spls_fair20-5_fair_modeA_pub_noAnom_new_common_bmi for example.
cd('/cs/research/intelsys/anomaly_detection/Code/spls_cca_pca/abcdAll_noA/framework/spls_fair20-5_fair_modeA_pub_noAnom_new_common_bmi')
load('cfg_1.mat')
res = cfg;

for i = 1:6 % Change according to the significant modes, in case they are different.
res.frwork.level = i; tic
[expVar_tr_cu(1,i), expVar_tr_cu(2,i), expVar_tr_cu(3,i)] = calc_expvar(res, 'otrid', 1);
[expVar_te_cu(1,i), expVar_te_cu(2,i), expVar_te_cu(3,i)] = calc_expvar(res, 'oteid', 1);
end

expVar_tr = expVar_tr_cu;
expVar_te = expVar_te_cu;

for i = 6:-1:2
expVar_tr(1,i) = expVar_tr_cu(1,i) - expVar_tr_cu(1,i-1);
expVar_tr(2,i) = expVar_tr_cu(2,i) - expVar_tr_cu(2,i-1);
expVar_tr(3,i) = expVar_tr_cu(3,i) - expVar_tr_cu(3,i-1);
expVar_te(1,i) = expVar_te_cu(1,i) - expVar_te_cu(1,i-1);
expVar_te(2,i) = expVar_te_cu(2,i) - expVar_te_cu(2,i-1);
expVar_te(3,i) = expVar_te_cu(3,i) - expVar_te_cu(3,i-1);
end