hmrG_stats (897c409d) · Commits · github_fork / Homer3

FuncRegistry/UserFunctions/hmrG_t_paired_baseline_HRF.m

+12 −12

Original line number	Diff line number	Diff line
		@@ -28,8 +28,6 @@ nSubj = length(yAvgSubjs);

		for iBlk = 1:nDataBlks



		for iSubj = 1:nSubj

		yAvg = yAvgSubjs{iSubj}(iBlk).GetDataTimeSeries('reshape');
		@@ -47,8 +45,8 @@ for iBlk = 1:nDataBlks
		end
		end

		baseline_yAvg(iSubj,:,:,:) = mean(yAvg(round(fq(tHRFrange(1) + abs(min(tHRF)))):round(fq(tHRFrange(2) + abs(min(tHRF)))),:,:,:),1);
		mean_yAvg(iSubj,:,:,:) = mean(yAvg(1:round(fq*abs(min(tHRF))),:,:,:),1);
		baseline_yAvg(iSubj,:,:,:) = mean(yAvg(1:round(fq*abs(min(tHRF))),:,:,:),1);
		mean_yAvg(iSubj,:,:,:) = mean(yAvg(round(fq(tHRFrange(1) + abs(min(tHRF)))):round(fq(tHRFrange(2) + abs(min(tHRF)))),:,:,:),1);

		end

		@@ -66,8 +64,6 @@ for iBlk = 1:nDataBlks
		end
		end
		end
		end


		% output
		hmrstatsG_base_cond.pval = pval;
		@@ -75,4 +71,8 @@ hmrstatsG_base_cond.hval = hval;
		hmrstatsG_base_cond.cval = cval;
		hmrstatsG_base_cond.tstats = tstats;
		hmrstatsG_base_cond.ml = ml;
		hmrstatsG_base_cond.mean_yAvg = mean_yAvg;

		end

+8 −3

Original line number	Diff line number	Diff line
		@@ -68,7 +68,7 @@ if ~isempty(nCond)
		end
		end
		end
		end


		% output
		hmrstatsG_contrast.pval = pval;
		@@ -77,3 +77,8 @@ hmrstatsG_contrast.cval = cval;
		hmrstatsG_contrast.tstats = tstats;
		hmrstatsG_contrast.ml = ml;
		hmrstatsG_contrast.c_vector = c_vector;
		hmrstatsG_contrast.mean_yAvg = mean_yAvg;

		end