%% Voxel-based analyses pub_noAnom_new_common - Regression USING BRAIN&BEHAV SCORES - Requires Matlab 2019b or later

clear; clc

% In our case we had downloaded and structured accordingly the variables we
% wanted to predict future outcomes from. The way the data need to be
% structured can be shown by loading the file below, where you have to
% change th path accordingly. The data in it are not real, but for
% illustration purposes only. You will need to request permission from the
% ABCD to download the real ones.
load('C:\Users\Konstantinos Tsirlis\Dropbox\My Desktop\Mourao - Miranda\ABCD Preprocessing\Behavioral Stuff\Data_and_Code\Predict Future Outcomes\future_outcomes_checkpoint1_new_common.mat')

% You have to change the path accordingly.
basedir = 'C:\Users\Konstantinos Tsirlis\Dropbox\My Desktop\Mourao - Miranda\ABCD results\abcdAll_noA\framework\spls_fair20-5_fair_modeA_pub_noAnom_new_common_bmi\res\';
cd(basedir)

P_mean = table();
for i = 1:(length(dir('level*'))) % make sure it's only the number of significant modes for the X model - BE CAREFUL THE -1 IN THE FOR-LOOP!
cd(['level' num2str(i)])
load('P.mat')
P = reshape(P,[11288,2]);
P = (P-mean(P))./std(P);
P_mean.(['mode' num2str(i)]) = mean([P(:,1) P(:,2)],2);
cd(basedir)
end

%%%%%% RAW DATA

% Grades - LinRes
grades_mdl_y2 = fitlm(table2array(P_mean(patientOrderFinal_saag_2y_idx,:)),abcdsaag01mod_final_saag_2y{:,6});
grades_mdl_y3 = fitlm(table2array(P_mean(patientOrderFinal_saag_3y_idx,:)),abcdsaag01mod_final_saag_3y{:,6});

[grades_mdl_y2_Ftest(2,1),grades_mdl_y2_Ftest(1,1)] = coefTest(grades_mdl_y2);
[grades_mdl_y3_Ftest(2,1),grades_mdl_y3_Ftest(1,1)] = coefTest(grades_mdl_y3);

dummy_var = [grades_mdl_y2_Ftest; cell2mat(struct2cell(grades_mdl_y2.Rsquared))];
writematrix(grades_mdl_y2.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B2','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H2','Sheet','Grades');

dummy_var = [grades_mdl_y3_Ftest; cell2mat(struct2cell(grades_mdl_y3.Rsquared))];
writematrix(grades_mdl_y3.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B10','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H10','Sheet','Grades');


% Suicidality - Binomial LogRes
suicidality_mdl_01_y1 = fitglm(table2array(P_mean(patientOrderFinal_suicidality_1y_idx,:)),suicidality01mod_final_01_suicidality_1y{:,10},'Distribution','Binomial');
suicidality_mdl_01_y2 = fitglm(table2array(P_mean(patientOrderFinal_suicidality_2y_idx,:)),suicidality01mod_final_01_suicidality_2y{:,10},'Distribution','Binomial');
[suicidality_mdl_01_y1_Ftest(2,1),suicidality_mdl_01_y1_Ftest(1,1)] = coefTest(suicidality_mdl_01_y1);
[suicidality_mdl_01_y2_Ftest(2,1),suicidality_mdl_01_y2_Ftest(1,1)] = coefTest(suicidality_mdl_01_y2);

dummy_var = [suicidality_mdl_01_y1_Ftest; cell2mat(struct2cell(suicidality_mdl_01_y1.Rsquared))]; dummy_var(5:end)=[];
writematrix(suicidality_mdl_01_y1.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B2','Sheet','Suicidality');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H2','Sheet','Suicidality');

dummy_var = [suicidality_mdl_01_y2_Ftest; cell2mat(struct2cell(suicidality_mdl_01_y2.Rsquared))]; dummy_var(5:end)=[];
writematrix(suicidality_mdl_01_y2.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B10','Sheet','Suicidality');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H10','Sheet','Suicidality');


% Hair - Binomial LogRes

abcdyhr_mdl_01_y1_etg = fitglm(table2array(P_mean(patientOrderFinal_hair_1y_idx,:)),abcdyhr01mod_final_01_hair_1y{:,9},'Distribution','Binomial');
abcdyhr_mdl_01_y1_cotinine = fitglm(table2array(P_mean(patientOrderFinal_hair_1y_idx,:)),abcdyhr01mod_final_01_hair_1y{:,10},'Distribution','Binomial');
abcdyhr_mdl_01_y1_sum = fitglm(table2array(P_mean(patientOrderFinal_hair_1y_idx,:)),abcdyhr01mod_final_01_hair_1y{:,16},'Distribution','Binomial');

abcdyhr_mdl_01_y2_etg = fitglm(table2array(P_mean(patientOrderFinal_hair_2y_idx,:)),abcdyhr01mod_final_01_hair_2y{:,9},'Distribution','Binomial');
abcdyhr_mdl_01_y2_cotinine = fitglm(table2array(P_mean(patientOrderFinal_hair_2y_idx,:)),abcdyhr01mod_final_01_hair_2y{:,10},'Distribution','Binomial');
abcdyhr_mdl_01_y2_sum = fitglm(table2array(P_mean(patientOrderFinal_hair_2y_idx,:)),abcdyhr01mod_final_01_hair_2y{:,16},'Distribution','Binomial');

abcdyhr_mdl_01_y3_etg = fitglm(table2array(P_mean(patientOrderFinal_hair_3y_idx,:)),abcdyhr01mod_final_01_hair_3y{:,9},'Distribution','Binomial');
abcdyhr_mdl_01_y3_cotinine = fitglm(table2array(P_mean(patientOrderFinal_hair_3y_idx,:)),abcdyhr01mod_final_01_hair_3y{:,10},'Distribution','Binomial');
abcdyhr_mdl_01_y3_sum = fitglm(table2array(P_mean(patientOrderFinal_hair_3y_idx,:)),abcdyhr01mod_final_01_hair_3y{:,16},'Distribution','Binomial');

[abcdyhr_mdl_01_y1_etg_Ftest(2,1),abcdyhr_mdl_01_y1_etg_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y1_etg);
[abcdyhr_mdl_01_y1_cotinine_Ftest(2,1),abcdyhr_mdl_01_y1_cotinine_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y1_cotinine);
[abcdyhr_mdl_01_y1_sum_Ftest(2,1),abcdyhr_mdl_01_y1_sum_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y1_sum);

[abcdyhr_mdl_01_y2_etg_Ftest(2,1),abcdyhr_mdl_01_y2_etg_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y2_etg);
[abcdyhr_mdl_01_y2_cotinine_Ftest(2,1),abcdyhr_mdl_01_y2_cotinine_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y2_cotinine);
[abcdyhr_mdl_01_y2_sum_Ftest(2,1),abcdyhr_mdl_01_y2_sum_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y2_sum);

[abcdyhr_mdl_01_y3_etg_Ftest(2,1),abcdyhr_mdl_01_y3_etg_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y3_etg);
[abcdyhr_mdl_01_y3_cotinine_Ftest(2,1),abcdyhr_mdl_01_y3_cotinine_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y3_cotinine);
[abcdyhr_mdl_01_y3_sum_Ftest(2,1),abcdyhr_mdl_01_y3_sum_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y3_sum);


% year 1
dummy_var = [abcdyhr_mdl_01_y1_etg_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y1_etg.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y1_etg.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B2','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H2','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y1_cotinine_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y1_cotinine.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y1_cotinine.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B10','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H10','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y1_sum_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y1_sum.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y1_sum.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B18','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H18','Sheet','Substances');

% year 2
dummy_var = [abcdyhr_mdl_01_y2_etg_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y2_etg.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y2_etg.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B26','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H26','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y2_cotinine_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y2_cotinine.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y2_cotinine.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B34','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H34','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y2_sum_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y2_sum.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y2_sum.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B42','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H42','Sheet','Substances');

% year 3
dummy_var = [abcdyhr_mdl_01_y3_etg_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y3_etg.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y3_etg.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B50','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H50','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y3_cotinine_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y3_cotinine.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y3_cotinine.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B58','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H58','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y3_sum_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y3_sum.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y3_sum.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B66','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H66','Sheet','Substances');


%%%%%% NORMALIZED DATA

% Grades - LinRes
grades_mdl_y2_norm = fitlm(normalize(table2array(P_mean(patientOrderFinal_saag_2y_idx,:))),normalize(abcdsaag01mod_final_saag_2y{:,6}));
grades_mdl_y3_norm = fitlm(normalize(table2array(P_mean(patientOrderFinal_saag_3y_idx,:))),normalize(abcdsaag01mod_final_saag_3y{:,6}));

[grades_mdl_y2_norm_Ftest(2,1),grades_mdl_y2_norm_Ftest(1,1)] = coefTest(grades_mdl_y2_norm);
[grades_mdl_y3_norm_Ftest(2,1),grades_mdl_y3_norm_Ftest(1,1)] = coefTest(grades_mdl_y3_norm);

dummy_var = [grades_mdl_y2_norm_Ftest; cell2mat(struct2cell(grades_mdl_y2_norm.Rsquared))];
writematrix(grades_mdl_y2_norm.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B18','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H18','Sheet','Grades');

dummy_var = [grades_mdl_y3_norm_Ftest; cell2mat(struct2cell(grades_mdl_y3_norm.Rsquared))];
writematrix(grades_mdl_y3_norm.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','B26','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm.xlsx'],'Range','H26','Sheet','Grades');



%% Voxel-based analyses pub_noAnom_new_common - Regression USING 6 PCS - Requires Matlab 2019b or later

clear; clc

% You first have to have run the first part of a simple CCA model with the
% same train-test sets as our main model. That way you can use the SVD
% results while respecting the same train-test sets.
svd_framework = 'C:\Users\Konstantinos Tsirlis\Dropbox\My Desktop\Mourao - Miranda\ABCD results\abcdAll_noA\framework\cca_holdout1-0.20_common_PCA_CCA_9_correct_outmat';
cd([svd_framework '\load\svd']);

to_load = dir('*_svd*');
load(to_load(1).name,'RY'); RYtest = RY;
load(to_load(2).name,'RY'); RYtrain = RY;


basedir = 'C:\Users\Konstantinos Tsirlis\Dropbox\My Desktop\Mourao - Miranda\ABCD results\abcdAll_noA\framework\spls_fair20-5_fair_modeA_pub_noAnom_new_common_bmi\res\';
load([basedir(1:end-4) 'outmat_1.mat'],'oteid')

RY = zeros(length(oteid),620);
RY(oteid,:) = RYtest;
RY(~oteid,:) = RYtrain;

RY = (RY-mean(RY))./std(RY);
P_mean = table();
for i = 1:6 % change according to the number of significant modes
P_mean.(['PC' num2str(i)]) = RY(:,i);
end

cd(basedir)
load('C:\Users\Konstantinos Tsirlis\Dropbox\My Desktop\Mourao - Miranda\ABCD Preprocessing\Behavioral Stuff\Data_and_Code\Predict Future Outcomes\future_outcomes_checkpoint1_new_common.mat')


%%%%%% RAW DATA

% Grades - LinRes
grades_mdl_y2 = fitlm(table2array(P_mean(patientOrderFinal_saag_2y_idx,:)),abcdsaag01mod_final_saag_2y{:,6});
grades_mdl_y3 = fitlm(table2array(P_mean(patientOrderFinal_saag_3y_idx,:)),abcdsaag01mod_final_saag_3y{:,6});

[grades_mdl_y2_Ftest(2,1),grades_mdl_y2_Ftest(1,1)] = coefTest(grades_mdl_y2);
[grades_mdl_y3_Ftest(2,1),grades_mdl_y3_Ftest(1,1)] = coefTest(grades_mdl_y3);

dummy_var = [grades_mdl_y2_Ftest; cell2mat(struct2cell(grades_mdl_y2.Rsquared))];
writematrix(grades_mdl_y2.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B2','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H2','Sheet','Grades');

dummy_var = [grades_mdl_y3_Ftest; cell2mat(struct2cell(grades_mdl_y3.Rsquared))];
writematrix(grades_mdl_y3.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B10','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H10','Sheet','Grades');


% Suicidality - Binomial LogRes
suicidality_mdl_01_y1 = fitglm(table2array(P_mean(patientOrderFinal_suicidality_1y_idx,:)),suicidality01mod_final_01_suicidality_1y{:,10},'Distribution','Binomial');
suicidality_mdl_01_y2 = fitglm(table2array(P_mean(patientOrderFinal_suicidality_2y_idx,:)),suicidality01mod_final_01_suicidality_2y{:,10},'Distribution','Binomial');
[suicidality_mdl_01_y1_Ftest(2,1),suicidality_mdl_01_y1_Ftest(1,1)] = coefTest(suicidality_mdl_01_y1);
[suicidality_mdl_01_y2_Ftest(2,1),suicidality_mdl_01_y2_Ftest(1,1)] = coefTest(suicidality_mdl_01_y2);

dummy_var = [suicidality_mdl_01_y1_Ftest; cell2mat(struct2cell(suicidality_mdl_01_y1.Rsquared))]; dummy_var(5:end)=[];
writematrix(suicidality_mdl_01_y1.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B2','Sheet','Suicidality');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H2','Sheet','Suicidality');

dummy_var = [suicidality_mdl_01_y2_Ftest; cell2mat(struct2cell(suicidality_mdl_01_y2.Rsquared))]; dummy_var(5:end)=[];
writematrix(suicidality_mdl_01_y2.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B10','Sheet','Suicidality');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H10','Sheet','Suicidality');


% Hair - Binomial LogRes

abcdyhr_mdl_01_y1_etg = fitglm(table2array(P_mean(patientOrderFinal_hair_1y_idx,:)),abcdyhr01mod_final_01_hair_1y{:,9},'Distribution','Binomial');
abcdyhr_mdl_01_y1_cotinine = fitglm(table2array(P_mean(patientOrderFinal_hair_1y_idx,:)),abcdyhr01mod_final_01_hair_1y{:,10},'Distribution','Binomial');
abcdyhr_mdl_01_y1_sum = fitglm(table2array(P_mean(patientOrderFinal_hair_1y_idx,:)),abcdyhr01mod_final_01_hair_1y{:,16},'Distribution','Binomial');

abcdyhr_mdl_01_y2_etg = fitglm(table2array(P_mean(patientOrderFinal_hair_2y_idx,:)),abcdyhr01mod_final_01_hair_2y{:,9},'Distribution','Binomial');
abcdyhr_mdl_01_y2_cotinine = fitglm(table2array(P_mean(patientOrderFinal_hair_2y_idx,:)),abcdyhr01mod_final_01_hair_2y{:,10},'Distribution','Binomial');
abcdyhr_mdl_01_y2_sum = fitglm(table2array(P_mean(patientOrderFinal_hair_2y_idx,:)),abcdyhr01mod_final_01_hair_2y{:,16},'Distribution','Binomial');

abcdyhr_mdl_01_y3_etg = fitglm(table2array(P_mean(patientOrderFinal_hair_3y_idx,:)),abcdyhr01mod_final_01_hair_3y{:,9},'Distribution','Binomial');
abcdyhr_mdl_01_y3_cotinine = fitglm(table2array(P_mean(patientOrderFinal_hair_3y_idx,:)),abcdyhr01mod_final_01_hair_3y{:,10},'Distribution','Binomial');
abcdyhr_mdl_01_y3_sum = fitglm(table2array(P_mean(patientOrderFinal_hair_3y_idx,:)),abcdyhr01mod_final_01_hair_3y{:,16},'Distribution','Binomial');

[abcdyhr_mdl_01_y1_etg_Ftest(2,1),abcdyhr_mdl_01_y1_etg_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y1_etg);
[abcdyhr_mdl_01_y1_cotinine_Ftest(2,1),abcdyhr_mdl_01_y1_cotinine_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y1_cotinine);
[abcdyhr_mdl_01_y1_sum_Ftest(2,1),abcdyhr_mdl_01_y1_sum_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y1_sum);

[abcdyhr_mdl_01_y2_etg_Ftest(2,1),abcdyhr_mdl_01_y2_etg_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y2_etg);
[abcdyhr_mdl_01_y2_cotinine_Ftest(2,1),abcdyhr_mdl_01_y2_cotinine_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y2_cotinine);
[abcdyhr_mdl_01_y2_sum_Ftest(2,1),abcdyhr_mdl_01_y2_sum_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y2_sum);

[abcdyhr_mdl_01_y3_etg_Ftest(2,1),abcdyhr_mdl_01_y3_etg_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y3_etg);
[abcdyhr_mdl_01_y3_cotinine_Ftest(2,1),abcdyhr_mdl_01_y3_cotinine_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y3_cotinine);
[abcdyhr_mdl_01_y3_sum_Ftest(2,1),abcdyhr_mdl_01_y3_sum_Ftest(1,1)] = coefTest(abcdyhr_mdl_01_y3_sum);


% year 1
dummy_var = [abcdyhr_mdl_01_y1_etg_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y1_etg.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y1_etg.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B2','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H2','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y1_cotinine_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y1_cotinine.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y1_cotinine.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B10','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H10','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y1_sum_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y1_sum.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y1_sum.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B18','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H18','Sheet','Substances');

% year 2
dummy_var = [abcdyhr_mdl_01_y2_etg_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y2_etg.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y2_etg.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B26','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H26','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y2_cotinine_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y2_cotinine.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y2_cotinine.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B34','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H34','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y2_sum_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y2_sum.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y2_sum.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B42','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H42','Sheet','Substances');

% year 3
dummy_var = [abcdyhr_mdl_01_y3_etg_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y3_etg.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y3_etg.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B50','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H50','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y3_cotinine_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y3_cotinine.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y3_cotinine.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B58','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H58','Sheet','Substances');

dummy_var = [abcdyhr_mdl_01_y3_sum_Ftest; cell2mat(struct2cell(abcdyhr_mdl_01_y3_sum.Rsquared))]; dummy_var(5:end)=[];
writematrix(abcdyhr_mdl_01_y3_sum.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B66','Sheet','Substances');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H66','Sheet','Substances');


%%%%%% NORMALIZED DATA

% Grades - LinRes
grades_mdl_y2_norm = fitlm(normalize(table2array(P_mean(patientOrderFinal_saag_2y_idx,:))),normalize(abcdsaag01mod_final_saag_2y{:,6}));
grades_mdl_y3_norm = fitlm(normalize(table2array(P_mean(patientOrderFinal_saag_3y_idx,:))),normalize(abcdsaag01mod_final_saag_3y{:,6}));

[grades_mdl_y2_norm_Ftest(2,1),grades_mdl_y2_norm_Ftest(1,1)] = coefTest(grades_mdl_y2_norm);
[grades_mdl_y3_norm_Ftest(2,1),grades_mdl_y3_norm_Ftest(1,1)] = coefTest(grades_mdl_y3_norm);

dummy_var = [grades_mdl_y2_norm_Ftest; cell2mat(struct2cell(grades_mdl_y2_norm.Rsquared))];
writematrix(grades_mdl_y2_norm.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B18','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H18','Sheet','Grades');

dummy_var = [grades_mdl_y3_norm_Ftest; cell2mat(struct2cell(grades_mdl_y3_norm.Rsquared))];
writematrix(grades_mdl_y3_norm.Coefficients{:,:},[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','B26','Sheet','Grades');
writematrix(dummy_var,[basedir 'future_outcomes_common_bmi_norm_PC.xlsx'],'Range','H26','Sheet','Grades');