Files
Experiments DB Dev 963fd889b2 feat(matlab): windowed paper-LME replication for mergeNaive (d0_10, d0_13)
Make tdcs_paper_lme window-aware (windowKey 'full'|'d0_10'|'d0_13',
default 'full' -> unchanged filename). Windowed runs give the two anchor
arms equal day coverage, so the stim:day interaction is not confounded by
the full-range coverage imbalance; the report now states whether coverage
is equal and adds an INTERPRETATION block. Add switch cases
paper_mergeNaive_d0_10 / _d0_13 and include them in run_all.

Finding: unlike mergeA2 (fair-window interaction n.s.), the pooled-naive
control keeps the interaction significant in the fair d0_13 window
(F(1)=7.07, p=0.0088, +1.80 [+0.46,+3.14]; equal on Day 1 p=0.20),
reproducing the paper without the coverage confound; d0_10 is borderline
(p=0.051).

Tests: tLme asserts equal coverage / no warning / one full summary for
both windows and a significant positive d0_13 interaction; bad windowKey
errors. Suite 41/41.

Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
2026-07-20 17:49:09 -04:00

23 lines
1.1 KiB
Matlab

%RUN_ALL Run all tDCS scenarios and write results/<scenario>.txt for each.
% Usage: matlab -batch "run_all"
% Covers the 6 GLMM scenarios and their 6 linear "days x tDCS" (lme_*) variants.
scenarios = {'unmerged_full', 'unmerged_d0_10', ...
'mergeA2_full', 'mergeA2_d0_10', ...
'mergeB2_full', 'mergeB2_d0_10', ...
'lme_unmerged_full', 'lme_unmerged_d0_10', ...
'lme_mergeA2_full', 'lme_mergeA2_d0_10', ...
'lme_mergeB2_full', 'lme_mergeB2_d0_10', ...
'lme_unmerged_d0_13', 'lme_mergeA2_d0_13', 'lme_mergeB2_d0_13', ...
'phase_unmerged', 'phase_mergeA2', 'phase_mergeB2', ...
'paper_unmerged', 'paper_mergeA2', 'paper_mergeB2', ...
'mergeNaive_full', 'paper_mergeNaive', 'phase_mergeNaive', ...
'paper_mergeNaive_d0_10', 'paper_mergeNaive_d0_13'};
for i = 1:numel(scenarios)
fprintf('\n\n### Running scenario: %s ###\n', scenarios{i});
tdcs_glm(scenarios{i});
end
fprintf('\n\nAll %d scenarios complete. See analysis/matlab/results/*.txt.\n', numel(scenarios));