ComputeIQConnections.m
b90a4357cb97master

ComputeIQConnections.m
View Options

			function [n_NIQ,IDX_NIQ1,IDX_NIQ2] = ComputeIQConnections(ISFC_MW,ISFC_RS,ISFC_PRS,Group_Label,Group_Label_pureRS,IQ,Motion_regressor,alpha)

			n_CTR = sum(Group_Label==1);
			n_ASD = sum(Group_Label==2);

			n_ROI = 270;
			n_conn = 36315;

			% Computes significant excursions separately for TD and ASD subjects, as in
			% our former paper (Bolton et al. 2018)
			tmpRS_TD = ISFC_RS(Group_Label==1);
			tmprRS2_TD = ISFC_PRS(Group_Label_pureRS==1);
			tmpMW_TD = ISFC_MW(Group_Label==1);

			SignMat_TD = FindExcursions(tmpRS_TD,tmpMW_TD,alpha,'Subject-Marginalized-Added-RS',tmprRS2_TD);

			tmpRS_ASD = ISFC_RS(Group_Label==2);
			tmprRS2_ASD = ISFC_PRS(Group_Label_pureRS==2);
			tmpMW_ASD = ISFC_MW(Group_Label==2);

			SignMat_ASD = FindExcursions(tmpRS_ASD,tmpMW_ASD,alpha,'Subject-Marginalized-Added-RS',tmprRS2_ASD);

			IQ_ASD = nanmean(IQ(Group_Label == 2));
			IQ_TD = nanmean(IQ(Group_Label == 1));

			IQ_avg = nan(n_CTR+n_ASD,1);
			IQ_avg(Group_Label == 1) = IQ_TD;
			IQ_avg(Group_Label == 2) = IQ_ASD;

			n_excursions = nan(n_conn,n_CTR+n_ASD);
			n_excursions(:,Group_Label==1) = sum(abs(SignMat_TD),2);
			n_excursions(:,Group_Label==2) = sum(abs(SignMat_ASD),2);


			EXC_mot = nan(size(n_excursions'));
			for c = 1:n_conn
			[~,EXC_mot(:,c)] = y_regress_ss(n_excursions(c,:)',Motion_regressor');
			EXC_mot(:,c) = EXC_mot(:,c) + mean(n_excursions(c,:));
			end
			EXC_mot = EXC_mot';

			IDX_NIQ1 = FindIQRelatedConns(n_excursions,IQ_avg,Group_Label,IQ);
			IDX_NIQ2 = FindIQRelatedConns(EXC_mot,IQ_avg,Group_Label,IQ);

			n_NIQ(1) = sum(IDX_NIQ1);
			n_NIQ(2) = sum(IDX_NIQ2);
			end