Hi Branko,
The differences in the results shared by you, specifically in the Sum of squares regression are significant. I am assuming that you might be using some kind of distribution with fixed parameters but varied samples for the two cases. It can be the reason of different sum of squares value but not very significant impact in the results.
I have tried to reproduce the difference at my end, but I am getting exact results for both the releases. I am attaching the results for reference.
SBPerror = [1.2, 0.8, 1.5, 1.0, 1.3, 0.9, 1.1, 1.4, 1.6, 1.2];
SBPK = [1, 2, 3, 1, 2, 3, 1, 2, 3, 1];
PP = [1, 1, 1, 2, 2, 2, 3, 3, 3, 3];
MP = [1, 2, 2, 1, 1, 2, 2, 1, 2, 2];
IABP_dP = [1, 2, 1, 2, 1, 2, 1, 2, 1, 2];
CP_PeakSlopeOW = [1, 1, 2, 2, 1, 1, 2, 2, 1, 1];
% Perform 5-way ANOVA
P = anovan(SBPerror, {SBPK, PP, MP, IABP_dP, CP_PeakSlopeOW}, 'continuous', [1 2 3 4 5], 'varnames', strvcat('SBPK', 'PP', 'MP', 'IABP_dP', 'CP_PeakSlopeOW'));
The results for R2020b
The results for R2023b
Both the results are exactly same.
If you think your data is an exception and the behaviour is for exactly same dataset, please share the exact dataset.
