clear; clc; sbj = readtable('301_subjs_corrected_stats_forAI.xlsx'); %sbj = sbj(:,2:end); sbj_ag = sbj(:,end-1:end); %age and gender columns sbj_b = sbj(:,2:end-2); %brain data only sbj_b = table2array(sbj_b); %convert table to matrix for division sbj_adj=sbj_b./sum(sbj_b,2); %adjust values by dividing by total brain volume sbj_adj = array2table(sbj_adj); %convert back to table (adjusted brain volumes) sbj_ag.AgeRound = fix(sbj_ag.Age./10)*10; %round age, make 53 and 59 = 50 sbj_ag.AgeGroup2(sbj_ag.Age<=25)='Y'; %young sbj_ag.AgeGroup2(sbj_ag.Age>25)='A'; %adult sbj_ag.AgeGroup3(sbj_ag.Age<=25)='Y'; %young sbj_ag.AgeGroup3((sbj_ag.Age>25) & (sbj_ag.Age<60))='A'; %adult sbj_ag.AgeGroup3(sbj_ag.Age>=60)='O'; %old sbj_table = [sbj_adj,sbj_ag]; %brain data, age, gender, ageRound, ageGroup %sbj_table_cc = sbj_table(:,[123,124,125,126,127,128,280,281,282,283,284]); %gcc,bcc,scc,age,gender,ageRound,ageGroup %sbj_table_cc.Properties.VariableNames=["GCC_L","GCC_R","BCC_L","BCC_R","SCC_L","SCC_R","AGE","GENDER","AGE_ROUND","AGE_GROUP2","AGE_GROUP3"]; sbj_table_cc = sbj_table(:,[123,124,125,126,127,128,281,283,284]); %gcc,bcc,scc,age,gender,ageRound,ageGroup2, ageGroup3 sbj_table_cc.Properties.VariableNames=["GCC_L","GCC_R","BCC_L","BCC_R","SCC_L","SCC_R","GENDER","AGE_GROUP2","AGE_GROUP3"]; %sbj_a.Properties.VariableNames = ["SFG_L","SFG_R","SFG_PFC_L","SFG_PFC_R","SFG_pole_L",... %"SFG_pole_R","MFG_L","MFG_R","MFG_DPFC_L","MFG_DPFC_R","IFG_opercularis_L","IFG_opercularis_R",... %"IFG_orbitalis_L","IFG_orbitalis_R","IFG_triangularis_L","IFG_triangularis_R","LFOG_L","LFOG_R",... %"MFOG_L","MFOG_R","RG_L","RG_R","PoCG_L","PoCG_R","PrCG_L","PrCG_R","SPG_L","SPG_R",... %"SMG_L","SMG_R","AG_L","AG_R","PrCu_L","PrCu_R","STG_L","STG_R","STG_L_pole","STG_R_pole",... %"MTG_L","MTG_R","MTG_L_pole","MTG_R_pole","ITG_L","ITG_R","PHG_L","PHG_R","ENT_L","ENT_R",... %"FuG_L","FuG_R","SOG_L","SOG_R","MOG_L","MOG_R","IOG_L","IOG_R","Cu_L","Cu_R","LG_L",... %"LG_R","rostral_ACC_L","rostral_ACC_R","subcallosal_ACC_L","subcallosal_ACC_R","subgenual_ACC_L",... %"subgenual_ACC_R","dorsal_ACC_L","dorsal_ACC_R","PCC_L","PCC_R","Insula_L","Insula_R","Amyg_L","Amyg_R",... %"Hippo_L","Hippo_R Caud_L Caud_R Put_L Put_R GP_L GP_R Thalamus_L Thalamus_R HypoThalamus_L HypoThalamus_R BasalForebrain_L BasalForebrain_R NucAccumbens_L NucAccumbens_R RedNc_L RedNc_R Snigra_L Snigra_R CerebellumGM_R CerebellumGM_L CP_L CP_R Midbrain_L Midbrain_R CST_L CST_R SCP_L SCP_R MCP_L MCP_R PCT_L PCT_R ICP_L ICP_R ML_L ML_R Pons_L Pons_R Medulla_L Medulla_R ACR_L ACR_R SCR_L SCR_R PCR_L PCR_R %GCC_L GCC_R BCC_L BCC_R SCC_L SCC_R PVWl_L PVWl_R ALIC_L ALIC_R PLIC_L PLIC_R RLIC_L RLIC_R EC_L EC_R CGC_L CGC_R CGH_L CGH_R Fx/ST_L Fx/ST_R Fx_L Fx_R IFO_L IFO_R PTR_L PTR_R SS_L SS_R SFO_L SFO_R SLF_L SLF_R Cl_L Cl_R BasalForebrain_L BasalForebrain_R BasalForebrain_L BasalForebrain_R BasalForebrain_L BasalForebrain_R BasalForebrain_L BasalForebrain_R Mammillary_R Mammillary_L LV_Frontal_L LV_body_L LV_atrium_L LV_Occipital_L LV_Inferior_L LV_Frontal_R LV_body_R LV_atrium_R LV_Occipital_R LV_Inferior_R III_and_IV_ventricle PVWa_L PVWa_R PVWp_L PVWp_R SFWM_L SFWM_R SFWM_PFC_L SFWM_PFC_R SFWM_pole_L SFWM_pole_R MFWM_L MFWM_R MFWM_DPFC_L MFWM_DPFC_R IFWM_opercularis_L IFWM_opercularis_R IFWM_orbitalis_L IFWM_orbitalis_R IFWM_triangularis_L IFWM_triangularis_R LFOWM_L LFOWM_R MFOWM_L MFOWM_R RGWM_L RGWM_R PoCWM_L PoCWM_R PrCWM_L PrCWM_R SPWM_L SPWM_R SMWM_L SMWM_R AGWM_L AGWM_R PrCuWM_L PrCuWM_R STWM_L STWM_R STWM_L_pole STWM_R_pole MTWM_L MTWM_R MTWM_L_pole MTWM_R_pole ITWM_L ITWM_R FuWM_L FuWM_R SOWM_L SOWM_R MOWM_L MOWM_R IOWM_L IOWM_R CuWM_L CuWM_R LWM_L LWM_R rostralWM_ACC_L rostralWM_ACC_R subcallosalWM_ACC_L subcallosalWM_ACC_R subgenualWM_ACC_L subgenualWM_ACC_R dorsalWM_ACC_L dorsalWM_ACC_R PCCWM_L PCCWM_R CerebellumWM_R CerebellumWM_L MCP_cb_L MCP_cb_R ICP_pons_L ICP_pons_R FrontSul_L FrontSul_R CentralSul_L CentralSul_R SylFrontSul_L SylFrontSul_R SylTempSul_L SylTempSul_R SylParieSul_L SylParieSul_R ParietSul_L ParietSul_R CinguSul_L CinguSul_R OcciptSul_L OcciptSul_R TempSul_L TempSul_R Caudate_tail_L Fimbria_L Caudate_tail_R Fimbria_R Chroid_LVetc_L Chroid_LVetc_R Age Gender %sbj_a = removevars(sbj_adj,{'Gender'}); %sbj_g = [sbj_adj,sbj_ag(:,2)]; %brain data + gender %sbj_a2=sbj_a; %create a variable to round age %sbj_a2.Age = roundn(sbj_a2.Age,1); %round age, make 53 = 50, make 59 = 60 %sbj_a2.Age = fix(sbj_a2.Age./10)*10; %round age, make 53 and 59 = 50 %sbj_a3 = sbj_a; %young<25, adult<60 old<80 %sbj_a3 = sbj_a3; %sbj_a3.AgeG = rand(301,1); %sbj_a3.AgeG(sbj_a3.Age<=25)='Y'; %sbj_a3.AgeG((sbj_a3.Age>25) & (sbj_a3.Age<60))='A'; %sbj_a3.AgeG(sbj_a3.Age>=60)='O'; %sbj_a3.Age=[]; %figure(1); hist(sbj_a2{:,end}); %sbj_a_26 = sbj_a2( (sbj_a2{:,end} >= 20) & (sbj_a2{:,end} < 80), : ); %remove certain ages %figure(2); hist(sbj_a_26{:,end}); %GCC 124 125 %BCC 126 127 %SCC 128 129 %sbj_a_gcc= sbj_a2(:,[123,124,280]); %gcc, age %sbj_a_bcc= sbj_a2(:,[125,126,280]); %bcc, age %sbj_a_scc= sbj_a2(:,[127,128,280]); %scc, age % sbj_a3_gcc= sbj_a3(:,[123,124,280]); %gcc, age % sbj_a3_bcc= sbj_a3(:,[125,126,280]); %bcc, age % sbj_a3_scc= sbj_a3(:,[127,128,280]); %scc, age % sbj_a3_cc = sbj_a3(:,[123,124,125,126,127,128,280]); %gcc,bcc,scc,age % sbj_a3_cc.Properties.VariableNames=["GCC_L","GCC_R","BCC_L","BCC_R","SCC_L","SCC_R","AGE"]; % sbj_g_gcc= sbj_g (:,[123,124,280]); %gcc, gender % sbj_g_bcc= sbj_g (:,[125,126,280]); %bcc, gender % sbj_g_scc= sbj_g (:,[127,128,280]); %scc, gender % sbj_g_cc = sbj_g (:,[123,124,125,126,127,128,280]); %gcc,bcc,scc,gender %sbj_g_cc.Properties.VariableNames=["GCC_L","GCC_R","BCC_L","BCC_R","SCC_L","SCC_R","GENDER"]; %sbj_ag_cc = [sb_bj(:,[123,124,125,126,127,128,280]); %sbj_ag_cc = [sbj_a3_cc , sbj_g_cc(:,7)]; %scatter(sbj_a2{:,280},sbj_a2{:,124}) %[corr(sbj_a2{:,123},sbj_a2{:,280}), corr(sbj_a2{:,124},sbj_a2{:,280}), corr(sbj_a2{:,125},sbj_a2{:,280}),... %corr(sbj_a2{:,126},sbj_a2{:,280}), corr(sbj_a2{:,127},sbj_a2{:,280}), corr(sbj_a2{:,128},sbj_a2{:,280})] %[corr(sbj_g{:,123},sbj_g{:,280}), corr(sbj_g{:,124},sbj_g{:,280}), corr(sbj_g{:,125},sbj_g{:,280}),... %corr(sbj_g{:,126},sbj_g{:,280}), corr(sbj_g{:,127},sbj_g{:,280}), corr(sbj_g{:,128},sbj_g{:,280})] %%%%%%%%%%%%%% STAT %%%%%%%%%%%%%%%%%%%%%%%%%%% %plot(sbj_table_cc.GCC_L, sbj_table_cc.GENDER);