proj_erythritol_4a-code/firstlevel/ery_4a_firstlevel_s3_diagnose_model.m

%% ery_4a_firstlevel_s3_diagnose_model.m
%
% This script will run diagnostics on first level models and publish them
% as an html report using the CANlab function scn_spm_design_check, and
% save variance inflation factors as a .mat file
%
% DOCUMENTATION
% help scn_spm_design_check
% 
% USAGE
% Script should be called from LaBGAScore_firstlevel_s1_options_design_struct.m, it
% is not for standalone use
%
%__________________________________________________________________________
%
% authors: Lukas Van Oudenhove
% date:   March, 2022
%
%__________________________________________________________________________
% @(#)% LaBGAScore_firstlevel_s3_diagnose_model.m         v1.0        
% last modified: 2022/03/22


%% VIFS AND DESIGN

vifs = scn_spm_design_check(subjfirstdir,'events_only','vif_thresh',LaBGAS_options.mandatory.vif_thresh);
drawnow, snapnow
save('vifs','vifs');
close all

%% FIRST LEVEL MAPS

if LaBGAS_options.display.plotmontages

    load(fullfile(subjfirstdir,'SPM.mat'));

    tmapnames = dir(fullfile(subjfirstdir,'spmT_*.nii'));
    tmapspaths = cell(1,size(tmapnames,1));
    tmapsobj = cell(1,size(tmapnames,1));
    montages = cell(1,size(tmapnames,1));

    [~,maskname,maskext] = fileparts(LaBGAS_options.display.mask);
    mask = [maskname,maskext];

    fprintf('\nShowing results at p < %s %s, k = %s, mask = %s\n',num2str(LaBGAS_options.display.input_threshold),LaBGAS_options.display.thresh_type,num2str(LaBGAS_options.display.k),mask); 

    for tmap = 1:size(tmapnames,1)
        tmapspaths{tmap} = fullfile(tmapnames(tmap).folder,tmapnames(tmap).name);
        tmapsobj{tmap} = statistic_image('image_names',tmapspaths{tmap},'type','T','dfe',SPM.xX.erdf);
        tmapsobj{tmap} = threshold(tmapsobj{tmap},LaBGAS_options.display.input_threshold,LaBGAS_options.display.thresh_type,'k',LaBGAS_options.display.k,'mask',LaBGAS_options.display.mask);
        create_figure('fmridisplay');
        wh_montage = 5;
        axis off
        figtitle = DSGN.contrastnames{tmap};
        o3 = canlab_results_fmridisplay([],'outline','linewidth',0.5,'montagetype','compact','overlay','mni_icbm152_t1_tal_nlin_sym_09a_brainonly.img');
        o3 = addblobs(o3,tmapsobj{tmap},'splitcolor',{[.1 .8 .8] [.1 .1 .8] [.9 .4 0] [1 1 0]});
        [o3,title_handle] = title_montage(o3,wh_montage,figtitle);
        set(title_handle,'FontSize',18);
        fighan = activate_figures(o3);
        f3 = fighan{1};
        f3.Tag = figtitle;
        f3.WindowState = 'maximized';
        drawnow,snapnow
        close(f3)
        clear figtitle o3 title_handle fighan f3
    end
    
end