Home > marsbar > @mardo_2 > ui_report.m

ui_report

PURPOSE ^

method for SPM2 design reporting

SYNOPSIS ^

function varargout = ui_report(D, varargin)

DESCRIPTION ^

 method for SPM2 design reporting

 Copied with minor edits from:  
 @(#)spm_DesRep.m    2.31 Andrew Holmes 03/03/28
 
 See that file for detailed commentary

 $Id$

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function varargout = ui_report(D, varargin)
0002 % method for SPM2 design reporting
0003 %
0004 % Copied with minor edits from:
0005 % @(#)spm_DesRep.m    2.31 Andrew Holmes 03/03/28
0006 %
0007 % See that file for detailed commentary
0008 %
0009 % $Id$
0010 
0011 %-Format arguments
0012 %-----------------------------------------------------------------------
0013 if nargin < 2
0014   action = 'desrepui'; 
0015 else
0016   action = varargin{1};
0017 end
0018 
0019 %-Generic CallBack code
0020 %-----------------------------------------------------------------------
0021 cb = 'tmp = get(findobj(''Tag'', ''DesRepUI''),''UserData''); ';
0022 
0023 % simplify access to design
0024 SPM = des_struct(D);
0025 
0026 % Add empty fields where necessary
0027 try
0028   SPM.xC;
0029 catch
0030   SPM.xC = {};
0031 end
0032 try
0033   SPM.xsDes;
0034 catch
0035   SPM.xsDes = [];
0036 end
0037 
0038 switch lower(action)
0039 
0040 %=======================================================================
0041 case 'desrepui'                                    %-Design reporting UI
0042 %=======================================================================
0043 % h = ui_report(D, 'DesRepUI')
0044 % h = ui_report(D, 'DesRepUI',SPM)
0045 
0046 %-Add a scaled design matrix to the design data structure
0047 %-----------------------------------------------------------------------
0048 if ~isfield(SPM.xX,'nKX')
0049     SPM.xX.nKX = spm_DesMtx('Sca',SPM.xX.X,SPM.xX.name);
0050 end
0051 
0052 % put back into design object
0053 D = des_struct(D, SPM);
0054 
0055 %-Draw menu
0056 %=======================================================================
0057 
0058 %-Get Interactive window and delete any previous DesRepUI menu
0059 %-----------------------------------------------------------------------
0060 Finter = spm_figure('GetWin','Interactive');
0061 delete(findobj(get(Finter,'Children'),'flat','Tag','DesRepUI'))
0062 
0063 %-Draw top level menu
0064 %-----------------------------------------------------------------------
0065 hC      = uimenu(Finter,'Label','Design',...
0066          'Separator','on',...        
0067          'Tag','DesRepUI',...
0068          'UserData',D,...
0069          'HandleVisibility','on');
0070 
0071 %-DesMtx
0072 %-----------------------------------------------------------------------
0073 hDesMtx = uimenu(hC,'Label','Design Matrix','Accelerator','D',...
0074         'CallBack',[cb,...
0075         'ui_report(tmp, ''DesMtx'')'],...
0076         'UserData',hC,...
0077         'HandleVisibility','off');
0078 
0079 %-Design matrix orthogonality
0080 %-----------------------------------------------------------------------
0081 h = uimenu(hC,'Label','Design orthogonality','Accelerator','O',...
0082         'CallBack',[cb,...
0083         'ui_report(tmp, ''DesOrth'')'],...
0084         'UserData',hC,...
0085         'HandleVisibility','off');
0086 
0087 %-Explore design
0088 %-----------------------------------------------------------------------
0089 hExplore = uimenu(hC,'Label','Explore','HandleVisibility','off');
0090 
0091 switch modality(D)
0092 case 'pet'
0093     hFnF = uimenu(hExplore,'Label','Files and factors','Accelerator','F',...
0094         'CallBack',[cb,...
0095         'ui_report(tmp, ''Files&Factors'')'],...
0096         'UserData',hC,...
0097         'HandleVisibility','off');
0098     hCovs = uimenu(hExplore,'Label','Covariates','Accelerator','C',...
0099         'CallBack',[cb,...
0100         'ui_report(tmp, ''Covs'')'],...
0101         'UserData',hC,...
0102         'HandleVisibility','off');
0103     if isempty(SPM.xC), set(hCovs,'Enable','off'), end
0104 case 'fmri'
0105     for j = 1:length(SPM.Sess)
0106         h = uimenu(hExplore,'Label',sprintf('Session %.0f ',j),...
0107             'HandleVisibility','off');
0108         for k = 1:length(SPM.Sess(j).Fc)
0109             uimenu(h,'Label',SPM.Sess(j).Fc(k).name,...
0110                  'CallBack',[cb,...
0111             sprintf('ui_report_fmri(tmp,%d,%d);',j,k)],...
0112                  'UserData',hC,...
0113                  'HandleVisibility','off')
0114         end
0115     end
0116 end
0117 
0118 
0119 %-Clear, Quit, Help
0120 %-----------------------------------------------------------------------
0121 uimenu(hC,'Label','Clear','Accelerator','L','Separator','on',...
0122     'CallBack','spm_results_ui(''Clear'')',...
0123     'HandleVisibility','off');
0124 uimenu(hC,'Label','Help','Separator','on',...
0125     'CallBack','spm_help(''spm_DesRep'')',...
0126     'HandleVisibility','off');
0127 
0128 %-Pop open 'Interactive' window
0129 %-----------------------------------------------------------------------
0130 figure(Finter)
0131 
0132 %-Return handle of menu
0133 %-----------------------------------------------------------------------
0134 varargout = {hC};
0135 
0136 
0137 %=======================================================================
0138 case 'files&factors'                         %-Summarise files & factors
0139 %=======================================================================
0140 % ui_report(D, 'Files&Factors',fnames,I,xC,sF,xs)
0141 fnames  = image_names(D);
0142 if isempty(fnames)
0143   fnames = cell(size(SPM.xX.X, 1), 1);
0144 end
0145 
0146 I       = SPM.xX.I;
0147 xC      = SPM.xC;
0148 sF      = SPM.xX.sF;
0149 xs      = SPM.xsDes;  %-Structure of description strings
0150 
0151 [fnames,CPath] = spm_str_manip(fnames,'c');    %-extract common path component
0152 nScan          = size(I,1);            %-#images
0153 bL             = any(diff(I,1),1);         %-Multiple factor levels?
0154 
0155 %-Get graphics window & window scaling
0156 Fgraph = spm_figure('GetWin','Graphics');
0157 spm_results_ui('Clear',Fgraph,0)
0158 FS     = spm('FontSizes');
0159 
0160 %-Display header information
0161 %-----------------------------------------------------------------------
0162 hTax = axes('Position',[0.03,0.85,0.94,0.1],...
0163     'DefaultTextFontSize',FS(9),...
0164     'XLim',[0,1],'YLim',[0,1],...
0165     'Visible','off');
0166 
0167 text(0.5,1,'Statistical analysis: Image files & covariates...',...
0168     'Fontsize',FS(14),'Fontweight','Bold',...
0169     'HorizontalAlignment','center')
0170 
0171 dx1 = 0.05;
0172 dx2 = 0.08;
0173 
0174 x = 0; text(x+.02,.1,'image #','Rotation',90)
0175 if bL(4), x=x+dx1; text(x+.01,.1,sF{4},'Rotation',90), end
0176 if bL(3), x=x+dx1; text(x+.01,.1,sF{3},'Rotation',90), end
0177 if bL(2), x=x+dx1; text(x+.01,.1,sF{2},'Rotation',90), end
0178 if bL(1), x=x+dx1; text(x+.01,.1,sF{1},'Rotation',90), end
0179 
0180 for j = 1:length(xC)
0181     n = size(xC(j).rc,2);
0182     if n>1, tmp=xC(j).cname; else, tmp={xC(j).rcname}; end
0183     for k=1:n
0184         x=x+dx2;
0185         text(x,.1,tmp{k},'Rotation',90,'Interpreter','TeX')
0186     end
0187 end
0188 
0189 x=x+dx2;
0190 text(x,0.65,'Base directory:','FontWeight','Bold')
0191 text(x,0.5,CPath,'FontSize',FS(8))
0192 text(x,0.2,'filename tails...')
0193 
0194 line('XData',[0 1],'YData',[0 0],'LineWidth',3,'Color','r')
0195 
0196 %-Tabulate file & covariate information
0197 %-----------------------------------------------------------------------
0198 hAx = axes('Position',[0.03,0.05,0.94,0.8],...
0199     'DefaultTextFontSize',FS(8),...
0200     'Units','points',...
0201     'Visible','off');
0202 AxPos = get(hAx,'Position'); set(hAx,'YLim',[0,AxPos(4)])
0203 
0204 dy = FS(9); y0 = floor(AxPos(4)) -dy; y  = y0;
0205 
0206 for i = 1:nScan
0207 
0208     %-Scan indices
0209     x = 0; text(x,y,sprintf('%03d',i))
0210     if bL(4), x=x+dx1; text(x,y,sprintf('%02d',I(i,4))), end
0211     if bL(3), x=x+dx1; text(x,y,sprintf('%02d',I(i,3))), end
0212     if bL(2), x=x+dx1; text(x,y,sprintf('%02d',I(i,2))), end
0213     if bL(1), x=x+dx1; text(x,y,sprintf('%02d',I(i,1))), end
0214 
0215     %-Covariates
0216     for j = 1:length(xC)
0217         for k=1:size(xC(j).rc,2)
0218             x=x+dx2;
0219             text(x,y,sprintf('%6g',xC(j).rc(i,k)),...
0220                 'HorizontalAlignment','Center')
0221         end
0222     end
0223 
0224     %-Filename tail(s) - could be multivariate
0225     x=x+dx2;
0226     text(x,y,fnames{i})
0227     y=y-dy;
0228 
0229     %-Paginate if necessary
0230     if y<dy
0231         text(0.5,0,sprintf('Page %d',spm_figure('#page')),...
0232             'FontSize',FS(8),'FontAngle','italic')
0233         spm_figure('NewPage',[hAx;get(hAx,'Children')])
0234         hAx = axes('Units','points','Position',AxPos,...
0235             'DefaultTextFontSize',FS(8),'YLim',[0,AxPos(4)],...
0236             'Visible','off');
0237         y = y0;
0238         text(y,0,'continued...','FontAngle','Italic')
0239     end
0240 end
0241 
0242 line('XData',[0 1],'YData',[y y],'LineWidth',3,'Color','r')
0243 
0244 
0245 %-Display description strings
0246 % (At bottom of current page - hope there's enough room!)
0247 %-----------------------------------------------------------------------
0248 if ~isempty(xs)
0249     y = y - 2*dy;
0250     for sf = fieldnames(xs)'
0251         text(0.3,y,[strrep(sf{1},'_',' '),' :'],...
0252             'HorizontalAlignment','Right','FontWeight','Bold',...
0253             'FontSize',FS(9))
0254         s = getfield(xs,sf{1});
0255         if ~iscellstr(s), s={s}; end
0256         for i=1:prod(size(s))
0257             text(0.31,y,s{i},'FontSize',FS(9))
0258             y=y-dy;
0259         end
0260     end
0261 end
0262 
0263 %-Register last page if paginated
0264 if spm_figure('#page')>1
0265     text(0.5,0,sprintf('Page %d/%d',spm_figure('#page')*[1,1]),...
0266         'FontSize',FS(8),'FontAngle','italic')
0267     spm_figure('NewPage',[hAx;get(hAx,'Children')])
0268 end
0269 
0270 %-Pop up the Graphics window
0271 %-----------------------------------------------------------------------
0272 figure(Fgraph)
0273 
0274 
0275 
0276 %=======================================================================
0277 case {'desmtx','desorth'} %-Display design matrix / design orthogonality
0278 %=======================================================================
0279 % ui_report(D, 'DesMtx',xX,fnames,xs)
0280 % ui_report(D, 'DesOrth',xX)
0281 
0282 xX      = SPM.xX;
0283 fnames  = image_names(D);
0284 xs      = SPM.xsDes;  %-Structure of description strings
0285 
0286 desmtx = strcmp(lower(varargin{1}),'desmtx');
0287 
0288 
0289 %-Locate DesMtx (X), scaled DesMtx (nX) & get parameter names (Xnames)
0290 %-----------------------------------------------------------------------
0291 if isfield(xX,'xKXs') & ...
0292         ~isempty(xX.xKXs) & isstruct(xX.xKXs)
0293     iX = 1;
0294     [nScan,nPar] = size(xX.xKXs.X);
0295 elseif isfield(xX,'X') & ~isempty(xX.X)
0296     iX = 0;
0297     [nScan,nPar] = size(xX.X);
0298 else
0299     error('Can''t find DesMtx in this structure!')
0300 end
0301 
0302 if isfield(xX,'nKX') & ~isempty(xX.nKX)
0303     inX = 1; else, inX = 0; end
0304 
0305 if isfield(xX,'name') & ~isempty(xX.name)
0306     Xnames = xX.name; else, Xnames = {}; end
0307 
0308 
0309 %-Compute design orthogonality matrix if DesOrth
0310 %-----------------------------------------------------------------------
0311 if ~desmtx
0312     if iX
0313     tmp  = sqrt(sum(xX.xKXs.X.^2));
0314     O    = xX.xKXs.X'*xX.xKXs.X./kron(tmp',tmp);
0315         tmp  = sum(xX.xKXs.X);
0316     else
0317     tmp  = sqrt(sum(xX.X.^2));
0318     O    = xX.X'*xX.X./kron(tmp',tmp);
0319         tmp  = sum(xX.X);
0320     end
0321     tmp = abs(tmp)<eps*1e5;
0322     bC  = kron(tmp',tmp);
0323 end
0324 
0325 
0326 %-Display
0327 %=======================================================================
0328 
0329 %-Get graphics window & FontSizes
0330 %-----------------------------------------------------------------------
0331 Fgraph = spm_figure('GetWin','Graphics');
0332 spm_results_ui('Clear',Fgraph,0)
0333 FS = spm('FontSizes');
0334 
0335 
0336 %-Title
0337 %-----------------------------------------------------------------------
0338 hTax = axes('Position',[0.03,0,0.94,1],...
0339     'DefaultTextFontSize',FS(9),...
0340     'XLim',[0,1],'YLim',[0,1],...
0341     'Visible','off');
0342 
0343 str='Statistical analysis: Design'; if ~desmtx, str=[str,' orthogonality']; end
0344 text(0.5,0.95,str,'Fontsize',FS(14),'Fontweight','Bold',...
0345     'HorizontalAlignment','center')
0346 
0347 line('Parent',hTax,...
0348     'XData',[0.3 0.7],'YData',[0.92 0.92],'LineWidth',3,'Color','r')
0349 
0350 
0351 %-Display design matrix
0352 %-----------------------------------------------------------------------
0353 hDesMtx = axes('Position',[.07 .4 .6 .4]);
0354 if inX        %-Got a scaled DesMtx
0355     hDesMtxIm = image((xX.nKX + 1)*32);
0356 elseif iX    %-No scaled DesMtx, DesMtx in .xKXs structure
0357     hDesMtxIm = image((spm_DesMtx('sca',xX.xKXs.X,Xnames) + 1)*32);
0358 else        %-No scaled DesMtx, no .xKXs, DesMtx in .X
0359     hDesMtxIm = image((spm_DesMtx('sca',xX.X,     Xnames) + 1)*32);
0360 end
0361 
0362 STick = ui_report(D, 'ScanTick',nScan,32);
0363 PTick = ui_report(D, 'ScanTick',nPar,32);
0364 
0365 set(hDesMtx,'TickDir','out',...
0366     'XTick',PTick,'XTickLabel','',...
0367     'YTick',STick,'YTickLabel','')
0368 if desmtx
0369     xlabel('parameters'), ylabel('images')
0370 else
0371     set(get(hDesMtx,'Xlabel'),...
0372         'Position',get(get(hDesMtx,'Ylabel'),'Position'),...
0373         'Rotation',90')
0374     xlabel('design matrix')
0375 end
0376 
0377 %-Parameter names
0378 if ~isempty(Xnames)
0379     axes('Position',[.07 .8 .6 .1],'Visible','off',...
0380         'DefaultTextFontSize',FS(8),'DefaultTextInterpreter','TeX',...
0381         'XLim',[0,nPar]+0.5)
0382     for i=PTick, text(i,.05,Xnames{i},'Rotation',90), end
0383 end
0384 
0385 %-Filenames
0386 % ( Show at most 32, showing every 2nd/3rd/4th/... as necessary to pair )
0387 % ( down to <32 items. Always show last item so #images is indicated.   )
0388 if desmtx & ~isempty(fnames)
0389     axes('Position',[.68 .4 .3 .4],'Visible','off',...
0390         'DefaultTextFontSize',FS(8),...
0391         'YLim',[0,nScan]+0.5,'YDir','Reverse')
0392     for i = STick
0393         try
0394             str  = fnames(i,:);
0395         catch
0396             str  = fnames{i};
0397         end
0398         text(0,i,spm_str_manip(str,'Ca35'));
0399     end
0400 end
0401 
0402 %-Setup callbacks to allow interrogation of design matrix
0403 %-----------------------------------------------------------------------
0404 if iX,     set(hDesMtxIm,'UserData',...
0405     struct('X',xX.xKXs.X,'Xnames',{Xnames},'fnames',{fnames}))
0406 else,     set(hDesMtxIm,'UserData',...
0407     struct('X',xX.X,     'Xnames',{Xnames},'fnames',{fnames}))
0408 end
0409 set(hDesMtxIm,'ButtonDownFcn',[cb 'ui_report(tmp, ''SurfDesMtx_CB'')'])
0410 
0411 
0412 if desmtx
0413     %-Parameter estimability/uniqueness
0414     %---------------------------------------------------------------
0415     hPEstAx   = axes('Position',[.07 .315 .6 .025],...
0416             'DefaultTextInterpreter','TeX');
0417     if iX,    est = spm_SpUtil('IsCon',xX.xKXs);
0418     else,    est = spm_SpUtil('IsCon',xX.X); end
0419     hParEstIm = image((est+1)*32);
0420     set(hPEstAx,...
0421         'XLim',[0,nPar]+.5,'XTick',[1:nPar-1]+.5,'XTickLabel','',...
0422         'YLim',[0,1]+.5,'YDir','reverse','YTick',[],...
0423         'Box','on','TickDir','in','XGrid','on','GridLineStyle','-');
0424     xlabel('parameter estimability')
0425     text((nPar+0.5 + nPar/30),1,...
0426         '(gray \rightarrow \beta not uniquely specified)',...
0427         'Interpreter','TeX','FontSize',FS(8))
0428     set(hParEstIm,'UserData',struct('est',est,'Xnames',{Xnames}))
0429     set(hParEstIm,'ButtonDownFcn',[cb 'ui_report(tmp, ''SurfEstIm_CB'')'])
0430 else
0431     %-Design orthogonality
0432     %---------------------------------------------------------------
0433     hDesO   = axes('Position',[.07 .18 .6 .2]);
0434     tmp = 1-abs(O); tmp(logical(tril(ones(nPar),-1))) = 1;
0435     hDesOIm = image(tmp*64);
0436     
0437     set(hDesO,'Box','off','TickDir','out',...
0438         'XaxisLocation','top','XTick',PTick,'XTickLabel','',...
0439         'YaxisLocation','right','YTick',PTick,'YTickLabel','',...
0440         'YDir','reverse')
0441     tmp = [1,1]'*[[0:nPar]+0.5];
0442     line('Xdata',tmp(1:end-1)','Ydata',tmp(2:end)')
0443 
0444     xlabel('design orthogonality')
0445     set(get(hDesO,'Xlabel'),'Position',[0.5,nPar,0],...
0446         'HorizontalAlignment','left',...
0447         'VerticalAlignment','top')
0448     set(hDesOIm,...
0449         'UserData',struct('O',O,'bC',bC,'Xnames',{Xnames}),...
0450         'ButtonDownFcn',[cb 'ui_report(tmp, ''SurfDesO_CB'')'])
0451 
0452     if ~isempty(Xnames)
0453         axes('Position',[.69 .18 0.01 .2],'Visible','off',...
0454             'DefaultTextFontSize',FS(10),...
0455             'DefaultTextInterpreter','TeX',...
0456             'YDir','reverse','YLim',[0,nPar]+0.5)
0457         for i=PTick
0458             text(0,i,Xnames{i},'HorizontalAlignment','left')
0459         end
0460     end
0461 
0462 end
0463 
0464 %-Design descriptions
0465 %-----------------------------------------------------------------------
0466 if desmtx
0467     str = 'Design description...';
0468     line('Parent',hTax,...
0469         'XData',[0.3 0.7],'YData',[0.28 0.28],'LineWidth',3,'Color','r')
0470     hAx = axes('Position',[0.03,0.05,0.94,0.22],'Visible','off');
0471 else
0472     str = '';
0473     line('Parent',hTax,...
0474         'XData',[0.3 0.7],'YData',[0.14 0.14],'LineWidth',3,'Color','r')
0475     hAx = axes('Position',[0.03,0.05,0.94,0.08],'Visible','off');
0476     xs = struct('Measure',    ['abs. value of cosine of angle between ',...
0477                  'columns of design matrix'],...
0478             'Scale',    {{    'black - colinear (cos=+1/-1)';...
0479                     'white - orthogonal (cos=0)';...
0480                     'gray  - not orthogonal or colinear'}});
0481 end
0482 
0483 
0484 if ~isempty(xs)
0485     set(hAx,'Units','points');
0486     AxPos = get(hAx,'Position');
0487     set(hAx,'YLim',[0,AxPos(4)])
0488     
0489     dy = FS(9); y0 = floor(AxPos(4)) -dy; y = y0;
0490 
0491     text(0.3,y,str,...
0492         'HorizontalAlignment','Center',...
0493         'FontWeight','Bold','FontSize',FS(11))
0494     y=y-2*dy;
0495     
0496     for sf = fieldnames(xs)'
0497         text(0.3,y,[strrep(sf{1},'_',' '),' :'],...
0498             'HorizontalAlignment','Right','FontWeight','Bold',...
0499             'FontSize',FS(9))
0500         s = getfield(xs,sf{1});
0501         if ~iscellstr(s), s={s}; end
0502         for i=1:prod(size(s))
0503             text(0.31,y,s{i},'FontSize',FS(9))
0504             y=y-dy;
0505         end
0506     end
0507 end
0508 
0509 %-Pop up the Graphics window
0510 %-----------------------------------------------------------------------
0511 figure(Fgraph)
0512 
0513 
0514 %=======================================================================
0515 case 'covs'                %-Plot and describe covariates (one per page)
0516 %=======================================================================
0517 % ui_report(D, 'Covs',xX,xC)
0518 
0519 xX = SPM.xX;
0520 xC = SPM.xC;
0521 
0522 if ~length(xC), spm('alert!','No covariates!',mfilename), return, end
0523 
0524 %-Get graphics window & window scaling
0525 Fgraph = spm_figure('GetWin','Graphics');
0526 spm_results_ui('Clear',Fgraph,0)
0527 FS = spm('FontSizes');
0528 
0529 %-Title
0530 %-----------------------------------------------------------------------
0531 hTax = axes('Position',[0.03,0,0.94,1],...
0532     'DefaultTextFontSize',FS(9),...
0533     'XLim',[0,1],'YLim',[0,1],...
0534     'Visible','off');
0535 
0536 text(0.5,0.95,'Statistical analysis: Covariates',...
0537     'Fontsize',FS(14),'Fontweight','Bold',...
0538     'HorizontalAlignment','center')
0539 
0540 text(0.5,0.82,'(covariates plotted over transposed design matrix)',...
0541     'FontSize',FS(8),'HorizontalAlignment','center')
0542 
0543 line('XData',[0.3 0.7],'YData',[0.92 0.92],'LineWidth',3,'Color','r')
0544 line('XData',[0.3 0.7],'YData',[0.44 0.44],'LineWidth',3,'Color','r')
0545 
0546 
0547 %-Design matrix (as underlay for plots) and parameter names
0548 %-----------------------------------------------------------------------
0549 [nScan,nPar]   = size(xX.X);
0550 if isfield(xX,'name') & ~isempty(xX.name)
0551     Xnames = xX.name; else, Xnames = {}; end
0552 
0553 %-Design matrix
0554 hDesMtx = axes('Position',[.1 .5 .7 .3]);
0555 if isfield(xX,'nKX') & ~isempty(xX.nKX)
0556     image(xX.nKX'*32+32)
0557 elseif isfield(xX,'xKXs') & ~isempty(xX.xKXs)
0558     image(spm_DesMtx('sca',xX.xKXs.X,Xnames)*32+32)
0559 else
0560     image(spm_DesMtx('sca',xX.X,Xnames)*32+32)
0561 end
0562 set(hDesMtx,'Visible','off')
0563 
0564 %-Parameter names
0565 hParAx = axes('Position',[.8 .5 .2 .3],'Visible','off',...
0566     'DefaultTextFontSize',FS(8),'DefaultTextInterpreter','TeX',...
0567     'YLim',[0.5,nPar+0.5],'YDir','Reverse');
0568 hPNames = zeros(nPar,1);
0569 for i = 1:nPar, hPNames(i) = text(.05,i,Xnames{i}); end
0570 
0571 
0572 %-Covariates - one page each
0573 %-----------------------------------------------------------------------
0574 for i = 1:length(xC)
0575 
0576     %-Title
0577     %---------------------------------------------------------------
0578     hSTitle = text(0.5,0.87,sprintf('%d : %s',i,xC(i).rcname),...
0579             'Parent',hTax,...
0580             'HorizontalAlignment','center',...
0581             'FontSize',FS(13),'FontWeight','Bold');
0582 
0583     %-Plot
0584     %---------------------------------------------------------------
0585     hAx = axes('Position',[.1 .5 .7 .3],...
0586             'TickDir','out','Box','off','Color','none',...
0587             'NextPlot','add',...
0588             'XLim',[0,nScan]+0.5);
0589     plot(xC(i).rc,'LineWidth',2)
0590     if nScan<48, plot(xC(i).rc,'.k','MarkerSize',20); end
0591     xlabel('image #')
0592     ylabel('covariate value')
0593 
0594 
0595     %-Descriptions
0596     %---------------------------------------------------------------
0597     hDAx = axes('Position',[0.03,0.1,0.94,0.30],'Visible','off');
0598     
0599     set(hDAx,'Units','points');
0600     tmp = get(hDAx,'Position');
0601     set(hDAx,'YLim',[0,tmp(4)])
0602     
0603     dy = FS(9); y0 = floor(tmp(4)) -dy; y = y0;
0604 
0605     %-Description strings from xC(i).descrip
0606     text(0.3,y,'Details :',...
0607         'HorizontalAlignment','Right',...
0608         'FontWeight','Bold','FontSize',FS(9))
0609     s = xC(i).descrip;
0610     if ~iscellstr(s), s={s}; end
0611     for j=1:prod(size(s))
0612         text(0.31,y,s{j},'FontSize',FS(9))
0613         y=y-dy;
0614     end
0615     y=y-dy;
0616 
0617     %-Key (if block of covariates entered)
0618     %---------------------------------------------------------------
0619     if size(xC(i).rc,2)>1
0620         ColorOrder = get(hAx,'ColorOrder');
0621         text(0.3,y,'Key :',...
0622             'HorizontalAlignment','Right',...
0623             'FontWeight','Bold','FontSize',FS(9))
0624         for j = 1:size(xC(i).rc,2)
0625             color = ColorOrder(mod(j-1,size(ColorOrder,1))+1,:);
0626             if size(xC(i).rc,2)==length(xC(i).cname)
0627                 str = xC(i).cname{j};
0628             else
0629                 str = sprintf('column %d',j);
0630             end
0631             text(0.31,y,str,'FontSize',FS(9),...
0632                 'Color',color)
0633             text(0.5,xC(i).rc(1,j),[str,' \rightarrow'],...
0634                 'Parent',hAx,...
0635                 'FontSize',FS(8),'FontWeight','Bold',...
0636                 'HorizontalAlignment','Right',...
0637                 'Interpreter','TeX',...
0638                 'Color',color)
0639             y=y-dy;
0640         end
0641         y=y-dy;
0642     end
0643 
0644 
0645     %-Associated parameters
0646     %---------------------------------------------------------------
0647     text(0.3,y,'Design matrix columns :',...
0648         'HorizontalAlignment','Right',...
0649         'FontWeight','Bold','FontSize',FS(9))
0650     if isempty(xC(i).cols)
0651         text(0.31,y,'(none)','FontSize',FS(9))
0652     else
0653         for j = xC(i).cols
0654             text(0.31,y,sprintf('%d : %s',j,Xnames{j}),...
0655                 'FontSize',FS(9),'Interpreter','TeX')
0656             y=y-dy;
0657         end
0658     end
0659     y=y-dy;
0660 
0661 
0662     %-Highlight parameter names
0663     %---------------------------------------------------------------
0664     hCurPNames = hPNames(xC(i).cols);
0665     set(hCurPNames,'Color','r','FontWeight','Bold','FontSize',FS(8))
0666 
0667 
0668     %-Paginate (if more than one covariate)
0669     %---------------------------------------------------------------
0670     if length(xC)>1
0671         spm_figure('NewPage',[hSTitle; hAx; get(hAx,'Children');...
0672             hCurPNames; hDAx; get(hDAx,'Children')]);
0673     end
0674 
0675 end
0676 
0677 %-Pop up the Graphics window
0678 %-----------------------------------------------------------------------
0679 figure(Fgraph)
0680 
0681 
0682 %=======================================================================
0683 case 'scantick'
0684 %=======================================================================
0685 % ui_report(D, 'ScanTick',nScan,lim)
0686 % ( Show at most 32, showing every 2nd/3rd/4th/... as necessary to pair )
0687 % ( down to <32 items. Always show last item so #images is indicated.    )
0688 if nargin<3, lim=32; else, lim=varargin{3}; end
0689 if nargin<2, error('insufficient arguments'), end
0690 nScan = varargin{2};
0691 
0692 p = max(1,ceil(nScan/lim));
0693 s = 1:p:nScan; s(end)=nScan;
0694 
0695 varargout = {s,lim};
0696 
0697 
0698 %=======================================================================
0699 case {'surfdesmtx_cb','surfdesmtxmo_cb','surfdesmtxup_cb'} %-Surf DesMtx
0700 %=======================================================================
0701 % ui_report(D, 'SurfDesMtx_CB')
0702 % ui_report(D, 'SurfDesMtxMo_CB')
0703 % ui_report(D, 'SurfDesMtxUp_CB')
0704 
0705 h    = get(gca,'Xlabel');
0706 
0707 if strcmp(lower(varargin{1}),'surfdesmtxup_cb')
0708     UD = get(h,'UserData');
0709     set(h,'String',UD.String,'Interpreter',UD.Interpreter,...
0710         'UserData',UD.UserData)
0711     set(gcbf,'WindowButtonMotionFcn','','WindowButtonUpFcn','')
0712     return
0713 end
0714 
0715 
0716 if strcmp(lower(varargin{1}),'surfdesmtx_cb')
0717     UD = struct(    'String',    get(h,'String'),...
0718             'Interpreter',    get(h,'Interpreter'),...
0719             'UserData',    get(h,'UserData'));
0720     set(h,'UserData',UD)
0721     set(gcbf,'WindowButtonMotionFcn',[cb 'ui_report(tmp, ''SurfDesMtxMo_CB'')'],...
0722          'WindowButtonUpFcn',    [cb 'ui_report(tmp, ''SurfDesMtxUp_CB'')'])
0723 end
0724 
0725 mm  = [get(gca,'YLim')',get(gca,'XLim')']+[.5,.5;-.5,-.5];
0726 ij  = get(gca,'CurrentPoint');
0727 ij  = round(min(max(ij(1,[2,1]),mm(1,:)),mm(2,:)));
0728 
0729 istr = 'none';
0730 switch get(gcbf,'SelectionType')
0731 case 'normal'
0732     try, str = sprintf('X(%d,%d) = %g',ij(1),ij(2),...
0733         subsref(get(gco,'UserData'),...
0734         struct('type',{'.','()'},'subs',{'X',{ij(1),ij(2)}})));
0735     catch, str='(no cached design matrix to surf)'; end
0736 case 'extend'
0737     try, str = sprintf('Image %d: %s',ij(1),...
0738         spm_str_manip(...
0739         subsref(get(gco,'UserData'),...
0740         struct('type',{'.','()'},...
0741             'subs',{'fnames',{ij(1),':'}})),'Ca40'));
0742     catch, str='(no cached image filenames to surf)'; end
0743 case 'alt'
0744     try, str = sprintf('Parameter %d: %s',ij(2),...
0745         subsref(get(gco,'UserData'),...
0746         struct('type',{'.','{}'},'subs',{'Xnames',{ij(2)}})));
0747         istr = 'tex';
0748     catch, str='(no cached parameter names to surf)'; end
0749 case 'open'
0750     try,    assignin('base','ans',subsref(get(gco,'UserData'),...
0751             struct('type',{'.'},'subs',{'X'})))
0752         evalin('base','ans')
0753     catch,    fprintf('%s GUI: can''t find design matrix\n',mfilename)
0754     end
0755     return
0756 end
0757 
0758 set(h,'String',str,'Interpreter',istr)
0759 
0760 
0761 %=======================================================================
0762 case {'surfestim_cb','surfestimmo_cb','surfestimup_cb'}  %-Surf ParEstIm
0763 %=======================================================================
0764 % ui_report(D, 'SurfEstIm_CB')
0765 % ui_report(D, 'SurfEstImMo_CB')
0766 % ui_report(D, 'SurfEstImUp_CB')
0767 
0768 h    = get(gca,'Xlabel');
0769 
0770 if strcmp(lower(varargin{1}),'surfestimup_cb')
0771     UD = get(h,'UserData');
0772     set(h,'String',UD.String,'Interpreter',UD.Interpreter,...
0773         'UserData',UD.UserData)
0774     set(gcbf,'WindowButtonMotionFcn','','WindowButtonUpFcn','')
0775     return
0776 end
0777 
0778 if strcmp(lower(varargin{1}),'surfestim_cb')
0779     UD = struct(    'String',    get(h,'String'),...
0780             'Interpreter',    get(h,'Interpreter'),...
0781             'UserData',    get(h,'UserData'));
0782     set(h,'UserData',UD)
0783     set(gcbf,'WindowButtonMotionFcn',[cb 'ui_report(tmp, ''SurfEstImMo_CB'')'],...
0784          'WindowButtonUpFcn',    [cb 'ui_report(tmp, ''SurfEstImUp_CB'')'])
0785 end
0786 
0787 mm  = [get(gca,'XLim')]+[.5,-.5];
0788 i   = get(gca,'CurrentPoint');
0789 i   = round(min(max(i(1,1),mm(1)),mm(2)));
0790 
0791 istr = 'none';
0792 switch get(gcbf,'SelectionType')
0793 case 'normal'
0794     try, tmp = {' (not unique)',' (unique)'};
0795     str = sprintf('Parameter %d : %s%s',...
0796         i,...
0797         subsref(get(gco,'UserData'),...
0798             struct('type',{'.','{}'},'subs',{'Xnames',{i}})),...
0799         tmp{subsref(get(gco,'UserData'),...
0800             struct('type',{'.','()'},'subs',{'est',{i}}))+1});
0801         istr = 'tex';
0802     catch, str='(no cached data to surf)'; end
0803 case {'extend','alt'}
0804     return
0805 case 'open'
0806     try,    UD = get(gco,'UserData');
0807         assignin('base','ans',...
0808             subsref(get(gco,'UserData'),...
0809                 struct('type',{'.'},'subs',{'est'})))
0810         evalin('base','ans')
0811     catch,    fprintf('%s GUI: can''t find design orthogonality\n',mfilename)
0812     end
0813     return
0814 end
0815 
0816 set(h,'String',str,'Interpreter',istr)
0817 
0818 
0819 
0820 %=======================================================================
0821 case {'surfdeso_cb','surfdesomo_cb','surfdesoup_cb'}    %-Surf DesOrthIm
0822 %=======================================================================
0823 % ui_report(D, 'SurfDesO_CB')
0824 % ui_report(D, 'SurfDesOMo_CB')
0825 % ui_report(D, 'SurfDesOUp_CB')
0826 
0827 h    = get(gca,'Xlabel');
0828 
0829 if strcmp(lower(varargin{1}),'surfdesoup_cb')
0830     UD = get(h,'UserData');
0831     set(h,'String',UD.String,'Interpreter',UD.Interpreter,...
0832         'UserData',UD.UserData)
0833     set(gcbf,'WindowButtonMotionFcn','','WindowButtonUpFcn','')
0834     return
0835 end
0836 
0837 if strcmp(lower(varargin{1}),'surfdeso_cb')
0838     UD = struct(    'String',    get(h,'String'),...
0839             'Interpreter',    get(h,'Interpreter'),...
0840             'UserData',    get(h,'UserData'));
0841     set(h,'UserData',UD)
0842     set(gcbf,'WindowButtonMotionFcn',[cb 'ui_report(tmp, ''SurfDesOMo_CB'')'],...
0843          'WindowButtonUpFcn',    [cb 'ui_report(tmp, ''SurfDesOUp_CB'')'])
0844 end
0845 
0846 mm  = [get(gca,'YLim')',get(gca,'XLim')']+[.5,.5;-.5,-.5];
0847 ij  = get(gca,'CurrentPoint');
0848 ij  = round(min(max(ij(1,[2,1]),mm(1,:)),mm(2,:)));
0849 if ij(1)>ij(2), return, end
0850 
0851 istr = 'none';
0852 switch get(gcbf,'SelectionType')
0853 case 'normal'
0854     try
0855         UD = get(gco,'UserData');
0856         if abs(abs(UD.O(ij(1),ij(2)))-1) < eps*1e1
0857              str = '{\bf colinear}';
0858         elseif abs(UD.O(ij(1),ij(2))) < eps*1e1
0859             str = '{\bf orthogonal}';
0860         else
0861             str = '{\bf not orthogonal}';
0862         end
0863         if ~diff(ij), str=[str,' {\it(same column)}']; end
0864         if UD.bC(ij(1),ij(2)), tmp=' ={\it r}'; else, tmp=''; end
0865         str = {    sprintf('{\\bf %s} (col %d) & {\\bf %s} (col %d)',...
0866                 UD.Xnames{ij(1)},ij(1),...
0867                 UD.Xnames{ij(2)},ij(2)),...
0868             sprintf('cos(\\theta)%s = %1.2f',...
0869                 tmp,UD.O(ij(1),ij(2))),...
0870             ['\rightarrow ',str]};
0871         istr = 'tex';
0872     catch, str='(no cached data to surf)'; end
0873 case {'extend','alt'}
0874     return
0875 case 'open'
0876     try,    UD = get(gco,'UserData');
0877         assignin('base','ans',UD.O)
0878         evalin('base','ans')
0879     catch,    fprintf('%s GUI: can''t find design orthogonality\n',mfilename)
0880     end
0881     return
0882 end
0883 
0884 set(h,'String',str,'Interpreter',istr)
0885 
0886 
0887 %=======================================================================
0888 case {'surfcon_cb','surfconmo_cb','surfconup_cb'}        %-Surf Contrast
0889 %=======================================================================
0890 % ui_report(D, 'SurfCon_CB')
0891 % ui_report(D, 'SurfConOMo_CB')
0892 % ui_report(D, 'SurfConOUp_CB')
0893 
0894 cUD = get(gco,'UserData');
0895 if ~isstruct(cUD) | ~isfield(cUD,'h')
0896     warning('contrast GUI objects setup incorrectly'), return
0897 end
0898 h    = cUD.h;
0899 
0900 if strcmp(lower(varargin{1}),'surfconup_cb')
0901     UD = get(h,'UserData');
0902     set(h,'String',UD.String,'Interpreter',UD.Interpreter,...
0903         'UserData',UD.UserData)
0904     set(gcbf,'WindowButtonMotionFcn','','WindowButtonUpFcn','')
0905     return
0906 end
0907 
0908 if strcmp(lower(varargin{1}),'surfcon_cb')
0909     UD = struct(    'String',    get(h,'String'),...
0910             'Interpreter',    get(h,'Interpreter'),...
0911             'UserData',    get(h,'UserData'));
0912     set(h,'UserData',UD)
0913     set(gcbf,'WindowButtonMotionFcn',[cb 'ui_report(tmp, ''SurfConMo_CB'')'],...
0914          'WindowButtonUpFcn',    [cb 'ui_report(tmp, ''SurfConUp_CB'')'])
0915 end
0916 
0917 mm  = [get(gca,'YLim')',get(gca,'XLim')']+[.5,.5;-.5,-.5];
0918 ij  = get(gca,'CurrentPoint');
0919 ij  = round(min(max(ij(1,[2,1]),mm(1,:)),mm(2,:)));
0920 
0921 istr = 'none';
0922 switch get(gcbf,'SelectionType')
0923 case 'normal'
0924     try
0925         if cUD.i>0, str = sprintf('%d',cUD.i); else, str = ''; end
0926         switch get(gco,'Type')
0927         case 'image'
0928             str = sprintf('%s\\{F\\}: {\\bf%s} (%d,%d) = %.2f',...
0929                 str,cUD.xCon.name,ij(2),ij(1),...
0930                 cUD.xCon.c(ij(2),ij(1)));
0931         case 'patch'
0932             str = sprintf('%s\\{T\\}: {\\bf%s} (%d) = %.2f',...
0933                 str,cUD.xCon.name,ij(2),...
0934                 cUD.xCon.c(ij(2)));
0935         otherwise, error('unexpected object type')
0936         end
0937         istr = 'TeX';
0938     catch, str='(no cached data to surf)'; end
0939 case {'alt','extend'}
0940     return
0941 case 'open'
0942     try,    assignin('base','ans',cUD.xCon.c')
0943         evalin('base','ans')
0944     catch,    fprintf('%s GUI: can''t find contrast\n',mfilename)
0945     end
0946     return
0947 end
0948 
0949 set(h,'String',str,'Interpreter',istr)
0950 
0951 
0952 %=======================================================================
0953 otherwise                                        %-Unknown action string
0954 %=======================================================================
0955 error(['Unknown action string: ',varargin{1}])
0956 
0957 
0958 
0959 %=======================================================================
0960 end

Generated on Thu 17-Mar-2016 01:14:46 by m2html © 2005