diff --git a/analysis/matlab/variations/boxa_a2_d0_10/learning_curve.png b/analysis/matlab/variations/boxa_a2_d0_10/learning_curve.png
index 3cf72c5..2b5c991 100644
Binary files a/analysis/matlab/variations/boxa_a2_d0_10/learning_curve.png and b/analysis/matlab/variations/boxa_a2_d0_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/boxa_a2_d0_10/learning_curve_rate.png b/analysis/matlab/variations/boxa_a2_d0_10/learning_curve_rate.png
index f3e0092..5f7ea1e 100644
Binary files a/analysis/matlab/variations/boxa_a2_d0_10/learning_curve_rate.png and b/analysis/matlab/variations/boxa_a2_d0_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/boxa_a2_d0_10/plotcurve.m b/analysis/matlab/variations/boxa_a2_d0_10/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/boxa_a2_d0_10/plotcurve.m
+++ b/analysis/matlab/variations/boxa_a2_d0_10/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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/boxa_a2_d0_5/learning_curve.png b/analysis/matlab/variations/boxa_a2_d0_5/learning_curve.png
index 5cb589f..34d34af 100644
Binary files a/analysis/matlab/variations/boxa_a2_d0_5/learning_curve.png and b/analysis/matlab/variations/boxa_a2_d0_5/learning_curve.png differ
diff --git a/analysis/matlab/variations/boxa_a2_d0_5/learning_curve_rate.png b/analysis/matlab/variations/boxa_a2_d0_5/learning_curve_rate.png
index 4aa6fca..876c2dd 100644
Binary files a/analysis/matlab/variations/boxa_a2_d0_5/learning_curve_rate.png and b/analysis/matlab/variations/boxa_a2_d0_5/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/boxa_a2_d0_5/plotcurve.m b/analysis/matlab/variations/boxa_a2_d0_5/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/boxa_a2_d0_5/plotcurve.m
+++ b/analysis/matlab/variations/boxa_a2_d0_5/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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/boxa_a2_d6_10/learning_curve.png b/analysis/matlab/variations/boxa_a2_d6_10/learning_curve.png
index fcbb813..4a60ff2 100644
Binary files a/analysis/matlab/variations/boxa_a2_d6_10/learning_curve.png and b/analysis/matlab/variations/boxa_a2_d6_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/boxa_a2_d6_10/learning_curve_rate.png b/analysis/matlab/variations/boxa_a2_d6_10/learning_curve_rate.png
index c1c8d4a..67f3ac5 100644
Binary files a/analysis/matlab/variations/boxa_a2_d6_10/learning_curve_rate.png and b/analysis/matlab/variations/boxa_a2_d6_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/boxa_a2_d6_10/plotcurve.m b/analysis/matlab/variations/boxa_a2_d6_10/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/boxa_a2_d6_10/plotcurve.m
+++ b/analysis/matlab/variations/boxa_a2_d6_10/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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/boxa_b2_d0_10/learning_curve.png b/analysis/matlab/variations/boxa_b2_d0_10/learning_curve.png
index 99834fa..b58ea4e 100644
Binary files a/analysis/matlab/variations/boxa_b2_d0_10/learning_curve.png and b/analysis/matlab/variations/boxa_b2_d0_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/boxa_b2_d0_10/learning_curve_rate.png b/analysis/matlab/variations/boxa_b2_d0_10/learning_curve_rate.png
index 34e00b2..597e785 100644
Binary files a/analysis/matlab/variations/boxa_b2_d0_10/learning_curve_rate.png and b/analysis/matlab/variations/boxa_b2_d0_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/boxa_b2_d0_10/plotcurve.m b/analysis/matlab/variations/boxa_b2_d0_10/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/boxa_b2_d0_10/plotcurve.m
+++ b/analysis/matlab/variations/boxa_b2_d0_10/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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/boxa_b2_d0_5/learning_curve.png b/analysis/matlab/variations/boxa_b2_d0_5/learning_curve.png
index 3ace26d..8843432 100644
Binary files a/analysis/matlab/variations/boxa_b2_d0_5/learning_curve.png and b/analysis/matlab/variations/boxa_b2_d0_5/learning_curve.png differ
diff --git a/analysis/matlab/variations/boxa_b2_d0_5/learning_curve_rate.png b/analysis/matlab/variations/boxa_b2_d0_5/learning_curve_rate.png
index cde4eaa..e899e80 100644
Binary files a/analysis/matlab/variations/boxa_b2_d0_5/learning_curve_rate.png and b/analysis/matlab/variations/boxa_b2_d0_5/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/boxa_b2_d0_5/plotcurve.m b/analysis/matlab/variations/boxa_b2_d0_5/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/boxa_b2_d0_5/plotcurve.m
+++ b/analysis/matlab/variations/boxa_b2_d0_5/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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/boxa_b2_d6_10/learning_curve.png b/analysis/matlab/variations/boxa_b2_d6_10/learning_curve.png
index 766bfc9..b263c69 100644
Binary files a/analysis/matlab/variations/boxa_b2_d6_10/learning_curve.png and b/analysis/matlab/variations/boxa_b2_d6_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/boxa_b2_d6_10/learning_curve_rate.png b/analysis/matlab/variations/boxa_b2_d6_10/learning_curve_rate.png
index c3e955c..de91998 100644
Binary files a/analysis/matlab/variations/boxa_b2_d6_10/learning_curve_rate.png and b/analysis/matlab/variations/boxa_b2_d6_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/boxa_b2_d6_10/plotcurve.m b/analysis/matlab/variations/boxa_b2_d6_10/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/boxa_b2_d6_10/plotcurve.m
+++ b/analysis/matlab/variations/boxa_b2_d6_10/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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/matched_current_d0_3/learning_curve.png b/analysis/matlab/variations/matched_current_d0_3/learning_curve.png
index 04c2af3..2ebcc47 100644
Binary files a/analysis/matlab/variations/matched_current_d0_3/learning_curve.png and b/analysis/matlab/variations/matched_current_d0_3/learning_curve.png differ
diff --git a/analysis/matlab/variations/matched_current_d0_3/learning_curve_rate.png b/analysis/matlab/variations/matched_current_d0_3/learning_curve_rate.png
index 81a1872..a960ff8 100644
Binary files a/analysis/matlab/variations/matched_current_d0_3/learning_curve_rate.png and b/analysis/matlab/variations/matched_current_d0_3/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/matched_current_d0_3/plotcurve.m b/analysis/matlab/variations/matched_current_d0_3/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/matched_current_d0_3/plotcurve.m
+++ b/analysis/matlab/variations/matched_current_d0_3/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 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_10/learning_curve.png b/analysis/matlab/variations/naive_a2_d0_10/learning_curve.png
index 456dcfa..fc020c1 100644
Binary files a/analysis/matlab/variations/naive_a2_d0_10/learning_curve.png and b/analysis/matlab/variations/naive_a2_d0_10/learning_curve.png differ
diff --git a/analysis/matlab/variations/naive_a2_d0_10/learning_curve_rate.png b/analysis/matlab/variations/naive_a2_d0_10/learning_curve_rate.png
index a094a92..973cc57 100644
Binary files a/analysis/matlab/variations/naive_a2_d0_10/learning_curve_rate.png and b/analysis/matlab/variations/naive_a2_d0_10/learning_curve_rate.png differ
diff --git a/analysis/matlab/variations/naive_a2_d0_10/plotcurve.m b/analysis/matlab/variations/naive_a2_d0_10/plotcurve.m
index 193a44b..8fb7f44 100644
--- a/analysis/matlab/variations/naive_a2_d0_10/plotcurve.m
+++ b/analysis/matlab/variations/naive_a2_d0_10/plotcurve.m
@@ -38,8 +38,9 @@ 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, 'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
-xtickangle(90);                  % vertical x-axis tick labels
+set(gca, 'XTick', xd, 'XTickLabelRotation', 90, ...   % vertical x-axis tick labels
+    'FontName', 'Arial', 'FontSize', 13, 'LineWidth', 1.5, 'Box', 'off');
+drawnow;                         % force layout so the rotation is applied headless
 exportgraphics(fig, outFile, 'Resolution', 150);
 close(fig);
 fprintf('%s: wrote %s (anodal N=%d, control N=%d)\n', vname, outFile, na, nc);