diff --git a/analysis/matlab/variations/naive_a2_d0_13/learning_curve.png b/analysis/matlab/variations/naive_a2_d0_13/learning_curve.png
index 10ec60a..5e03d21 100644
Binary files a/analysis/matlab/variations/naive_a2_d0_13/learning_curve.png and b/analysis/matlab/variations/naive_a2_d0_13/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_a2_d0_13/learning_curve_rate.png b/analysis/matlab/variations/naive_a2_d0_13/learning_curve_rate.png
index 40e1c33..e0da25e 100644
Binary files a/analysis/matlab/variations/naive_a2_d0_13/learning_curve_rate.png and b/analysis/matlab/variations/naive_a2_d0_13/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_a2_d0_13/plotcurve.m b/analysis/matlab/variations/naive_a2_d0_13/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_a2_d0_13/plotcurve.m
+++ b/analysis/matlab/variations/naive_a2_d0_13/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_a2_d0_5/learning_curve.png b/analysis/matlab/variations/naive_a2_d0_5/learning_curve.png
index a133800..a9113a7 100644
Binary files a/analysis/matlab/variations/naive_a2_d0_5/learning_curve.png and b/analysis/matlab/variations/naive_a2_d0_5/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_a2_d0_5/learning_curve_rate.png b/analysis/matlab/variations/naive_a2_d0_5/learning_curve_rate.png
index 80ceb5c..8acbb21 100644
Binary files a/analysis/matlab/variations/naive_a2_d0_5/learning_curve_rate.png and b/analysis/matlab/variations/naive_a2_d0_5/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_a2_d0_5/plotcurve.m b/analysis/matlab/variations/naive_a2_d0_5/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_a2_d0_5/plotcurve.m
+++ b/analysis/matlab/variations/naive_a2_d0_5/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_a2_d6_10/learning_curve.png b/analysis/matlab/variations/naive_a2_d6_10/learning_curve.png
index f272551..5131881 100644
Binary files a/analysis/matlab/variations/naive_a2_d6_10/learning_curve.png and b/analysis/matlab/variations/naive_a2_d6_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_a2_d6_10/learning_curve_rate.png b/analysis/matlab/variations/naive_a2_d6_10/learning_curve_rate.png
index 445e3d2..3c2fc37 100644
Binary files a/analysis/matlab/variations/naive_a2_d6_10/learning_curve_rate.png and b/analysis/matlab/variations/naive_a2_d6_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_a2_d6_10/plotcurve.m b/analysis/matlab/variations/naive_a2_d6_10/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_a2_d6_10/plotcurve.m
+++ b/analysis/matlab/variations/naive_a2_d6_10/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_a2_d6_13/learning_curve.png b/analysis/matlab/variations/naive_a2_d6_13/learning_curve.png
index 69ba7ce..4850011 100644
Binary files a/analysis/matlab/variations/naive_a2_d6_13/learning_curve.png and b/analysis/matlab/variations/naive_a2_d6_13/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_a2_d6_13/learning_curve_rate.png b/analysis/matlab/variations/naive_a2_d6_13/learning_curve_rate.png
index 80456ee..1e7519e 100644
Binary files a/analysis/matlab/variations/naive_a2_d6_13/learning_curve_rate.png and b/analysis/matlab/variations/naive_a2_d6_13/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_a2_d6_13/plotcurve.m b/analysis/matlab/variations/naive_a2_d6_13/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_a2_d6_13/plotcurve.m
+++ b/analysis/matlab/variations/naive_a2_d6_13/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_boxa_d0_10/learning_curve.png b/analysis/matlab/variations/naive_boxa_d0_10/learning_curve.png
index 1c1abc5..fec8036 100644
Binary files a/analysis/matlab/variations/naive_boxa_d0_10/learning_curve.png and b/analysis/matlab/variations/naive_boxa_d0_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d0_10/learning_curve_rate.png b/analysis/matlab/variations/naive_boxa_d0_10/learning_curve_rate.png
index 69697ac..8d9a0c2 100644
Binary files a/analysis/matlab/variations/naive_boxa_d0_10/learning_curve_rate.png and b/analysis/matlab/variations/naive_boxa_d0_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d0_10/plotcurve.m b/analysis/matlab/variations/naive_boxa_d0_10/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_boxa_d0_10/plotcurve.m
+++ b/analysis/matlab/variations/naive_boxa_d0_10/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_boxa_d0_13/learning_curve.png b/analysis/matlab/variations/naive_boxa_d0_13/learning_curve.png
index 7ee39e3..fcf0362 100644
Binary files a/analysis/matlab/variations/naive_boxa_d0_13/learning_curve.png and b/analysis/matlab/variations/naive_boxa_d0_13/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d0_13/learning_curve_rate.png b/analysis/matlab/variations/naive_boxa_d0_13/learning_curve_rate.png
index 662ec6b..9d14c87 100644
Binary files a/analysis/matlab/variations/naive_boxa_d0_13/learning_curve_rate.png and b/analysis/matlab/variations/naive_boxa_d0_13/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d0_13/plotcurve.m b/analysis/matlab/variations/naive_boxa_d0_13/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_boxa_d0_13/plotcurve.m
+++ b/analysis/matlab/variations/naive_boxa_d0_13/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_boxa_d0_5/learning_curve.png b/analysis/matlab/variations/naive_boxa_d0_5/learning_curve.png
index b9bb2ec..b7e73a4 100644
Binary files a/analysis/matlab/variations/naive_boxa_d0_5/learning_curve.png and b/analysis/matlab/variations/naive_boxa_d0_5/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d0_5/learning_curve_rate.png b/analysis/matlab/variations/naive_boxa_d0_5/learning_curve_rate.png
index b3bc2b3..49917f8 100644
Binary files a/analysis/matlab/variations/naive_boxa_d0_5/learning_curve_rate.png and b/analysis/matlab/variations/naive_boxa_d0_5/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d0_5/plotcurve.m b/analysis/matlab/variations/naive_boxa_d0_5/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_boxa_d0_5/plotcurve.m
+++ b/analysis/matlab/variations/naive_boxa_d0_5/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);
diff --git a/analysis/matlab/variations/naive_boxa_d6_10/learning_curve.png b/analysis/matlab/variations/naive_boxa_d6_10/learning_curve.png
index f36a51b..184c44b 100644
Binary files a/analysis/matlab/variations/naive_boxa_d6_10/learning_curve.png and b/analysis/matlab/variations/naive_boxa_d6_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d6_10/learning_curve_rate.png b/analysis/matlab/variations/naive_boxa_d6_10/learning_curve_rate.png
index 04109bb..e18e8ba 100644
Binary files a/analysis/matlab/variations/naive_boxa_d6_10/learning_curve_rate.png and b/analysis/matlab/variations/naive_boxa_d6_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_boxa_d6_10/plotcurve.m b/analysis/matlab/variations/naive_boxa_d6_10/plotcurve.m
index 8fb7f44..8399de2 100644
--- a/analysis/matlab/variations/naive_boxa_d6_10/plotcurve.m
+++ b/analysis/matlab/variations/naive_boxa_d6_10/plotcurve.m
@@ -38,9 +38,22 @@ legend([e1 e2], {sprintf('anodal, N = %d', na), sprintf('control, N = %d', nc)},
     'Location', 'northwest', 'Box', 'off');
 xlabel('training day'); ylabel(ylab);
 title(vname, 'Interpreter', 'none');
-set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+% Vertical x-axis tick labels. XTickLabelRotation / xtickangle are NOT honored
+% by exportgraphics on a headless invisible figure, so draw the labels manually
+% as rotated text objects (which always render).
+set(gca, 'XTick', xd, 'XTickLabel', {}, 'Position', [0.13 0.20 0.80 0.70], ...
     'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-drawnow;                         % force layout so the rotation is applied headless
+drawnow;
+yl = ylim(gca);
+yb = yl(1) - 0.015 * (yl(2) - yl(1));
+for k = 1:numel(xd)
+    text(xd(k), yb, num2str(xd(k)), 'Parent', gca, 'Rotation', 90, ...
+        'HorizontalAlignment', 'right', 'VerticalAlignment', 'middle', ...
+        'FontName', 'Arial', 'FontSize', 13, 'Clipping', 'off');
+end
+xlh = get(gca, 'XLabel');            % push the axis label below the vertical ticks
+set(xlh, 'Units', 'normalized', 'Position', [0.5 -0.15 0]);
+drawnow;
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);