interim solution: applying pop_rmbase() to grand averages

Author: Buch

code/analyses-erp/OUTPUT.txt

@@ -43,31 +43,31 @@
 fprintf('=== VISIBLE-TARGET ANALYSES ===\n\n');
 
 %% ERN - test if detectable
-% social-FE p = 0.044 → DETECTABLE
+% social-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_112_social_vis_FE', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_FE', 'Soc_Corr'}, ...
     'visible_soc_FE_ERN_detection', [0 100], ern_cluster);
 
-% social-NFE p = 0.275 → NOT DETECTABLE
+% social-NFE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_113_social_vis_NFE', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_NFE', 'Soc_Corr'}, ...
     'visible_soc_NFE_ERN_detection', [0 100], ern_cluster);
 
-% nonsocial-FE p = 0.591 → NOT DETECTABLE
+% nonsocial-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_212_nonsoc_vis_FE', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_FE', 'Nonsoc_Corr'}, ...
     'visible_nonsoc_FE_ERN_detection', [0 100], ern_cluster);
 
-% nonsocial-NFE p = 0.865 → NOT DETECTABLE
+% nonsocial-NFE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_213_nonsoc_vis_NFE', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_NFE', 'Nonsoc_Corr'}, ...
     'visible_nonsoc_NFE_ERN_detection', [0 100], ern_cluster);
 
-%% delta-ERN ANOVA (error type × social condition) → no interaction, no error type effect, social condition effect p = .005
+%% delta-ERN ANOVA (error type × social condition)
 design = table({'FE'; 'FE'; 'NFE'; 'NFE'}, ...
     {'social'; 'nonsocial'; 'social'; 'nonsocial'}, ...
     'VariableNames', {'ErrorType', 'SocialCondition'});
@@ -79,31 +79,31 @@
     'ranova', design);
 
 %% Pe - test if detectable
-% social-FE p = 0.028 → DETECTABLE
+% social-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_112_social_vis_FE', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_FE', 'Soc_Corr'}, ...
     'visible_soc_FE_Pe_detection', [200 500], pe_cluster);
 
-% social-NFE p = 0.069 → NOT DETECTABLE
+% social-NFE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_113_social_vis_NFE', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_NFE', 'Soc_Corr'}, ...
     'visible_soc_NFE_Pe_detection', [200 500], pe_cluster);
 
-% nonsocial-FE p = 0.068 → NOT DETECTABLE
+% nonsocial-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_212_nonsoc_vis_FE', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_FE', 'Nonsoc_Corr'}, ...
     'visible_nonsoc_FE_Pe_detection', [200 500], pe_cluster);
 
-% nonsocial-NFE p = 0.169 → NOT DETECTABLE
+% nonsocial-NFE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_213_nonsoc_vis_NFE', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_NFE', 'Nonsoc_Corr'}, ...
     'visible_nonsoc_NFE_Pe_detection', [200 500], pe_cluster);
 
-%% delta-Pe ANOVA (error type × social condition) → all p > .05
+%% delta-Pe ANOVA (error type × social condition)
 testing_meanAmplitude(char(diffWaves), ...
     {'diffWave_soc_vis_FE', 'diffWave_nonsoc_vis_FE', ...
      'diffWave_soc_vis_NFE', 'diffWave_nonsoc_vis_NFE'}, ...
@@ -115,26 +115,26 @@
 fprintf('\n--- Collapsed Errors Analysis ---\n');
 
 %% ERN - collapsed errors
-% social collapsed errors p = 0.151 → NOT DETECTABLE
+% social collapsed errors
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_110_social_vis_error', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_Error_Collapsed', 'Soc_Corr'}, ...
     'visible_soc_collapsed_ERN_detection', [0 100], ern_cluster);
 
-% nonsocial collapsed errors p = 0.769 → NOT DETECTABLE
+% nonsocial collapsed errors
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_210_nonsoc_vis_error', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_Error_Collapsed', 'Nonsoc_Corr'}, ...
     'visible_nonsoc_collapsed_ERN_detection', [0 100], ern_cluster);
 
-% delta-ERN t-test (social vs nonsocial) → significant: p = .008 
+% delta-ERN t-test (social vs nonsocial)
 testing_meanAmplitude(char(diffWaves), ...
     {'diffWave_soc_vis_error', 'diffWave_nonsoc_vis_error'}, ...
     {'Soc_Error_Collapsed', 'Nonsoc_Error_Collapsed'}, ...
     'visible_collapsed_deltaERN_comparison', [0 100], ern_cluster);
 
 %% Pe - collapsed errors
-% social collapsed errors p = 0.013 → DETECTABLE
+% social collapsed errors
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_110_social_vis_error', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_Error_Collapsed', 'Soc_Corr'}, ...
@@ -158,25 +158,25 @@
 fprintf('\n=== INVISIBLE-TARGET ANALYSES ===\n\n');
 
 %% ERN - test if detectable
-% social-FE p = 0.991 → NOT DETECTABLE
+% social-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_102_social_invis_FE', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_Invis_FE', 'Soc_Vis_Corr'}, ...
     'invisible_soc_FE_ERN_detection', [0 100], ern_cluster);
 
-% social-NFG p = 0.082 → NOT DETECTABLE
+% social-NFG
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_104_social_invis_NFG', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_Invis_NFG', 'Soc_Vis_Corr'}, ...
     'invisible_soc_NFG_ERN_detection', [0 100], ern_cluster);
 
-% nonsocial-FE p = 0.015 → DETECTABLE (but reversed direction!)
+% nonsocial-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_202_nonsoc_invis_FE', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_Invis_FE', 'Nonsoc_Vis_Corr'}, ...
     'invisible_nonsoc_FE_ERN_detection', [0 100], ern_cluster);
 
-% nonsocial-NFG p = 0.102 → NOT DETECTABLE
+% nonsocial-NFG
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_204_nonsoc_invis_NFG', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_Invis_NFG', 'Nonsoc_Vis_Corr'}, ...
@@ -194,25 +194,25 @@
     'ranova', design_invis);
 
 %% Pe - test if detectable
-% social-FE p = 0.420 → NOT DETECTABLE
+% social-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_102_social_invis_FE', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_Invis_FE', 'Soc_Vis_Corr'}, ...
     'invisible_soc_FE_Pe_detection', [200 500], pe_cluster);
 
-% social-NFG p = 0.519 → NOT DETECTABLE
+% social-NFG
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_104_social_invis_NFG', 'grandAVG_111_social_vis_corr'}, ...
     {'Soc_Invis_NFG', 'Soc_Vis_Corr'}, ...
     'invisible_soc_NFG_Pe_detection', [200 500], pe_cluster);
 
-% nonsocial-FE p = 0.019 → DETECTABLE
+% nonsocial-FE
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_202_nonsoc_invis_FE', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_Invis_FE', 'Nonsoc_Vis_Corr'}, ...
     'invisible_nonsoc_FE_Pe_detection', [200 500], pe_cluster);
 
-% nonsocial-NFG p = 0.015 → DETECTABLE 
+% nonsocial-NFG
 testing_meanAmplitude(char(grandAver), ...
     {'grandAVG_204_nonsoc_invis_NFG', 'grandAVG_211_nonsoc_vis_corr'}, ...
     {'Nonsoc_Invis_NFG', 'Nonsoc_Vis_Corr'}, ...
@@ -248,13 +248,6 @@
 end
 fprintf('Moved %d invisible-target statistics files\n', length(invisFiles));
 
-%% generate html report
-fprintf('\nGenerating HTML report...\n');
-publish('batch_erp_analyses.m', ...
-    'format', 'html', ...
-    'outputDir', reportsDir, ...
-    'evalCode', false);
-
 %% create summary file
 summaryFile = fullfile(outputDir, 'analysis_summary.txt');
 fid = fopen(summaryFile, 'w');
@@ -277,4 +270,12 @@
 fprintf('\n========================================\n');
 fprintf('Analysis Complete!\n');
 fprintf('Output saved to:\n%s\n', outputDir);
-fprintf('========================================\n');
+fprintf('========================================\n');
+
+
+% %% generate html report
+% fprintf('\nGenerating HTML report...\n');
+% publish('batch_erp_analyses.m', ...
+%     'format', 'html', ...
+%     'outputDir', reportsDir, ...
+%     'evalCode', true);

code/analyses-erp/batch_erp_analyses.asv

@@ -34,13 +34,13 @@ visual_comperp({'grandAVG_111_social_vis_corr', 'grandAVG_211_nonsoc_vis_corr',
 
 % ERN
 pimpfigure('ERP', [-150 400 -3.5 2.5], {'1' '1:' '3' '3:' '5' '5:'}, ...
-           {'correct social', 'correct nonsocial', 'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial (N = 23)'}, ...
            [0 100], [vis_cor; vis_cor; vis_FE; vis_FE; vis_NFE; vis_NFE]);
 title('Ne/ERN - visible-target - grand averages');
 
 % Pe
 pimpfigure('ERP', [-100 600 -2 3], {'1' '1:' '3' '3:' '5' '5:'}, ...
-           {'correct social', 'correct nonsocial', 'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial (N = 23)'}, ...
            [200 500], [vis_cor; vis_cor; vis_FE; vis_FE; vis_NFE; vis_NFE]);
 title('Pe - visible-target - grand averages');
 
@@ -51,15 +51,15 @@ visual_comperp({'diffWave_soc_vis_FE', 'diffWave_nonsoc_vis_FE', ...
 
 % ERN
 pimpfigure('ERP', [-150 400 -2.5 2.5], {'1' '1:' '3' '3:'}, ...
-           {'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial'}, ...
            [0 100], [vis_FE; vis_FE; vis_NFE; vis_NFE]);
-title('∆ERN - visible-target - difference waves');
+title('∆ERN - visible-target');
 
 % Pe
 pimpfigure('ERP', [-100 600 -3 2], {'1' '1:' '3' '3:'}, ...
            {'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
            [200 500], [vis_FE; vis_FE; vis_NFE; vis_NFE]);
-title('∆Pe - visible-target - difference waves');
+title('∆Pe - visible-target');
 
 
 %% additional tests - visible-target - collapsed errors (FE + NFE combined)
@@ -70,13 +70,13 @@ visual_comperp({'grandAVG_111_social_vis_corr', 'grandAVG_211_nonsoc_vis_corr',
 
 % ERN
 pimpfigure('ERP', [-150 400 -3.2 2], {'1' '1:' '3' '3:'}, ...
-           {'correct social', 'correct nonsocial', 'error social', 'error nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [0 100], [vis_cor; vis_cor; vis_err; vis_err]);
 title('Ne/ERN - visible-target - collapsed across error type');
 
 % Pe
 pimpfigure('ERP', [-100 600 -1.5 2.75], {'1' '1:' '3' '3:'}, ...
-           {'correct social', 'correct nonsocial', 'error social', 'error nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [200 500], [vis_cor; vis_cor; vis_err; vis_err]);
 title('Pe - visible-target - collapsed across error type');
 
@@ -86,13 +86,13 @@ visual_comperp({'diffWave_soc_vis_error', 'diffWave_nonsoc_vis_error'}, ...
 
 % ERN
 pimpfigure('ERP', [-150 400 -2 2], {'1' '1:'}, ...
-           {'error social', 'error nonsocial'}, ...
+           {'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [0 100], [vis_err; vis_err]);
 title('∆ERN - visible-target - collapsed across error type');
 
 % Pe
 pimpfigure('ERP', [-100 600 -3.1 2], {'1' '1:'}, ...
-           {'error social', 'error nonsocial'}, ...
+           {'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [200 500], [vis_err; vis_err]);
 title('∆Pe - visible-target - collapsed across error type');
 
@@ -109,13 +109,13 @@ visual_comperp({'grandAVG_111_social_vis_corr', 'grandAVG_211_nonsoc_vis_corr',
                 char(grandAver));
 
 % ERN
-pimpfigure('ERP', [-150 400 -3.5 1.5], {'1' '1:' '3' '3:' '5' '5:'}, ...
+pimpfigure('ERP', [-150 400 -2.5 1.5], {'1' '1:' '3' '3:' '5' '5:'}, ...
            {'vis-correct social', 'vis-correct nonsocial', 'FE social', 'FE nonsocial', 'NFG social', 'NFG nonsocial'}, ...
            [0 100], [vis_cor; vis_cor; invis_FE; invis_FE; invis_NFG; invis_NFG]);
 title('Ne/ERN - invisible-target - grand averages');
 
 % Pe
-pimpfigure('ERP', [-100 600 -2 3], {'1' '1:' '3' '3:' '5' '5:'}, ...
+pimpfigure('ERP', [-100 600 -1 2.75], {'1' '1:' '3' '3:' '5' '5:'}, ...
            {'vis-correct social', 'vis-correct nonsocial', 'FE social', 'FE nonsocial', 'NFG social', 'NFG nonsocial'}, ...
            [200 500], [vis_cor; vis_cor; invis_FE; invis_FE; invis_NFG; invis_NFG]);
 title('Pe - invisible-target - grand averages');
@@ -124,13 +124,46 @@ title('Pe - invisible-target - grand averages');
 visual_comperp({'diffWave_soc_invis_FE', 'diffWave_nonsoc_invis_FE', ...
                 'diffWave_soc_invis_NFG', 'diffWave_nonsoc_invis_NFG'}, ...
                 char(diffWaves));
+
 % ERN
-pimpfigure('ERP', [-150 400 -3.5 2.5], {'1' '1:' '3' '3:'}, ...
+pimpfigure('ERP', [-150 400 -1.2 2], {'1' '1:' '3' '3:'}, ...
            {'FE social', 'FE nonsocial', 'NFG social', 'NFG nonsocial'}, ...
            [0 100], [invis_FE; invis_FE; invis_NFG; invis_NFG]);
-title('∆ERN - invisible-target - difference waves');
+title('∆ERN - invisible-target');
+
 % Pe
-pimpfigure('ERP', [-100 600 -2 3], {'1' '1:' '3' '3:'}, ...
+pimpfigure('ERP', [-100 600 -2.6 1.5], {'1' '1:' '3' '3:'}, ...
            {'FE social', 'FE nonsocial', 'NFG social', 'NFG nonsocial'}, ...
            [200 500], [invis_FE; invis_FE; invis_NFG; invis_NFG]);
-title('∆Pe - invisible-target - difference waves');
+title('∆Pe - invisible-target');
+
+
+
+
+
+%% MASTER LEGEND
+figure('Position', [100 100 400 300]);
+hold on;
+
+% visible conditions
+plot(NaN, NaN, '-', 'Color', vis_cor/255, 'LineWidth', 4); 
+plot(NaN, NaN, ':', 'Color', vis_cor/255, 'LineWidth', 4);
+plot(NaN, NaN, '-', 'Color', vis_FE/255, 'LineWidth', 4);
+plot(NaN, NaN, ':', 'Color', vis_FE/255, 'LineWidth', 4);
+plot(NaN, NaN, '-', 'Color', vis_NFE/255, 'LineWidth', 4);
+plot(NaN, NaN, ':', 'Color', vis_NFE/255, 'LineWidth', 4);
+plot(NaN, NaN, '-', 'Color', vis_err/255, 'LineWidth', 4);
+plot(NaN, NaN, ':', 'Color', vis_err/255, 'LineWidth', 4);
+
+% invisible conditions
+plot(NaN, NaN, '-', 'Color', invis_NFG/255, 'LineWidth', 4);
+plot(NaN, NaN, ':', 'Color', invis_NFG/255, 'LineWidth', 4);
+
+legend({'Correct (nonsocial)', 'Correct (social)', ...
+        'Flanker Error (nonsocial)', 'Flanker Error (social)', ...
+        'Non-Flanker Error (nonsocial)', 'Non-Flanker Error (social)', ...
+        'Error collapsed (nonsocial)', 'Error collapsed (social)', ...
+        'Non-Flanker Guess (nonsocial)', 'Non-Flanker Guess (social)'}, ...
+        'Location', 'bestoutside', 'FontSize', 12);
+title('Legend - Line Types & Colors');
+axis off;

code/analyses-erp/batch_erp_analyses.m

@@ -34,13 +34,13 @@
 
 % ERN
 pimpfigure('ERP', [-150 400 -3.5 2.5], {'1' '1:' '3' '3:' '5' '5:'}, ...
-           {'correct social', 'correct nonsocial', 'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial (N = 23)'}, ...
            [0 100], [vis_cor; vis_cor; vis_FE; vis_FE; vis_NFE; vis_NFE]);
 title('Ne/ERN - visible-target - grand averages');
 
 % Pe
 pimpfigure('ERP', [-100 600 -2 3], {'1' '1:' '3' '3:' '5' '5:'}, ...
-           {'correct social', 'correct nonsocial', 'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial (N = 23)'}, ...
            [200 500], [vis_cor; vis_cor; vis_FE; vis_FE; vis_NFE; vis_NFE]);
 title('Pe - visible-target - grand averages');
 
@@ -51,13 +51,13 @@
 
 % ERN
 pimpfigure('ERP', [-150 400 -2.5 2.5], {'1' '1:' '3' '3:'}, ...
-           {'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial (N = 23)'}, ...
            [0 100], [vis_FE; vis_FE; vis_NFE; vis_NFE]);
 title('∆ERN - visible-target');
 
 % Pe
 pimpfigure('ERP', [-100 600 -3 2], {'1' '1:' '3' '3:'}, ...
-           {'FE social', 'FE nonsocial', 'NFE social', 'NFE nonsocial'}, ...
+           {'FE social (N = 25)', 'FE nonsocial (N = 26)', 'NFE social (N = 22)', 'NFE nonsocial (N = 23)'}, ...
            [200 500], [vis_FE; vis_FE; vis_NFE; vis_NFE]);
 title('∆Pe - visible-target');
 
@@ -70,13 +70,13 @@
 
 % ERN
 pimpfigure('ERP', [-150 400 -3.2 2], {'1' '1:' '3' '3:'}, ...
-           {'correct social', 'correct nonsocial', 'error social', 'error nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [0 100], [vis_cor; vis_cor; vis_err; vis_err]);
 title('Ne/ERN - visible-target - collapsed across error type');
 
 % Pe
 pimpfigure('ERP', [-100 600 -1.5 2.75], {'1' '1:' '3' '3:'}, ...
-           {'correct social', 'correct nonsocial', 'error social', 'error nonsocial'}, ...
+           {'correct social (N = 29)', 'correct nonsocial (N = 29)', 'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [200 500], [vis_cor; vis_cor; vis_err; vis_err]);
 title('Pe - visible-target - collapsed across error type');
 
@@ -86,13 +86,13 @@
 
 % ERN
 pimpfigure('ERP', [-150 400 -2 2], {'1' '1:'}, ...
-           {'error social', 'error nonsocial'}, ...
+           {'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [0 100], [vis_err; vis_err]);
 title('∆ERN - visible-target - collapsed across error type');
 
 % Pe
 pimpfigure('ERP', [-100 600 -3.1 2], {'1' '1:'}, ...
-           {'error social', 'error nonsocial'}, ...
+           {'error social (N = 26)', 'error nonsocial (N = 26)'}, ...
            [200 500], [vis_err; vis_err]);
 title('∆Pe - visible-target - collapsed across error type');