numrical report

imapLMMmcc.m

@@ -29,7 +29,6 @@
 StatMap.Pmask0 = StatMap.Pmask;
 % clear old fields
 StatMap.Pmask  = [];
-StatMap.report = [];
 mask   = isnan(LMMmap.MSE)==0;
 tbl    = LMMmap.Variables;
 nonnan = find(mask);

imapLMMreport.m

@@ -0,0 +1,67 @@
+function imapLMMreport(StatMap_c,varargin)
+if nargin>1
+    % a 4d matrix with (Npredictor,beta[95%CI],xSize,ySize)
+    betamaps  = varargin{1};
+    plotcondi = 1;
+else
+    plotcondi = 0;
+end
+%%
+clc
+label = StatMap_c.label;
+Pmask = StatMap_c.Pmask;
+Fmap  = StatMap_c.map;
+dof   = StatMap_c.df;
+if isfield(StatMap_c,'beta')
+    betaall   = StatMap_c.beta;
+    betaCIall = StatMap_c.betaCI;
+end
+for ilabel = 1:length(label)
+    Pmasktmp = squeeze(Pmask(ilabel,:,:));
+    if sum(Pmasktmp(:))>0
+        disp(['Displaying result for ',label{ilabel}])
+        [maskbw,cluster]=bwlabel(Pmasktmp);
+        figure;imagesc(maskbw);axis equal off;colorbar
+        for ic = 1:cluster
+            disp(['  Numrical report for cluster ' num2str(ic)])
+            cluster_sel = maskbw == ic;
+            Fvalueall = Fmap(ilabel,cluster_sel);
+
+            disp(['    The local maximum within the cluster: F(',num2str(dof(ilabel,1)),...
+                ',',num2str(dof(ilabel,2)),') = ',num2str(max(Fvalueall))])
+            if isfield(StatMap_c,'beta')
+                maxloc  = Fvalueall==max(Fvalueall);
+                beta    = betaall(ilabel,cluster_sel);
+                betaCI  = squeeze(betaCIall(ilabel,:,cluster_sel));
+                disp(['    The beta contrast at the local maximum is ', num2str(beta(maxloc)),...
+                    ' [',num2str(betaCI(1,maxloc)),', ',num2str(betaCI(2,maxloc)),']'])
+            end
+            if plotcondi==1
+                for ic2 = 1:size(betamaps,1)
+                    beta1=squeeze(betamaps(ic2,:,cluster_sel));
+                    maxloc  = Fvalueall==max(Fvalueall);
+                    disp(['    Condition ',num2str(ic2),': ',num2str(beta1(1,maxloc)),...
+                        ' [',num2str(beta1(2,maxloc)),', ',num2str(beta1(3,maxloc)),']'])
+                end
+            end
+
+            disp(['    The local minimum within the cluster: F(',num2str(dof(ilabel,1)),...
+                ',',num2str(dof(ilabel,2)),') = ',num2str(min(Fvalueall))])
+            if isfield(StatMap_c,'beta')
+                minloc  = Fvalueall==min(Fvalueall);
+                beta    = betaall(ilabel,cluster_sel);
+                betaCI  = squeeze(betaCIall(ilabel,:,cluster_sel));
+                disp(['    The beta contrast at the local minimum is ', num2str(beta(minloc)),...
+                    ' [',num2str(betaCI(1,minloc)),', ',num2str(betaCI(2,minloc)),']'])
+            end
+            if plotcondi==1
+                for ic2 = 1:size(betamaps,1)
+                    beta1=squeeze(betamaps(ic2,:,cluster_sel));
+                    minloc  = Fvalueall==min(Fvalueall);
+                    disp(['    Condition ',num2str(ic2),': ',num2str(beta1(1,minloc)),...
+                        ' [',num2str(beta1(2,minloc)),', ',num2str(beta1(3,minloc)),']'])
+                end
+            end
+        end
+    end
+end