Appropriate anova function in MATLAB to check PDL::Stats results?

[etj]

As can be seen at https://github.com/PDLPorters/PDL-Stats/blob/474f0f56dc8cdbe52036a9cf82b7df1058fa7339/t/stats_glm.t#L596-L607, I am trying to fix PDL::Stats's anova functions. You can see in the linked section I have added what I think is the correct MATLAB (in fact Octave) code to replicate the PDL::Stats code, but I am not at all sure.

Can a statistics wizard help a poor stats neophyte out? Especially by creating the right Octave code and checking the tested-for values starting line 623 to tell me what's correct?

[etj]

A note from the future: as guided by the knowledgeable zmughal, I stopped trying to use MATLAB/Octave, and switched to using R, which is certainly the go-to solution for free statistics software. This is the R equivalent of the test in the given PDL::Stats test file:

library(data.table)
library(rstatix)
tdata <- data.frame(
 stringsAsFactors = FALSE,
 dv = c(3.0,2.0,1.0,5.0,2.0,1.0,5.0,3.0,1.0,4.0,1.0,2.0,3.0,5.0,5.0,3.
+0,4.0,2.0,1.0,5.0,4.0,3.0,2.0,2.0),
 id = c(0L,1L,2L,3L,0L,1L,2L,3L,0L,1L,2L,3L,0L,1L,2L,3L,0L,1L,2L,3L,0L
+,1L,2L,3L),
 w = c(0L,0L,0L,0L,0L,0L,0L,0L,1L,1L,1L,1L,1L,1L,1L,1L,2L,2L,2L,2L,2L,
+2L,2L,2L),
 b = c(0L,0L,0L,0L,1L,1L,1L,1L,0L,0L,0L,0L,1L,1L,1L,1L,0L,0L,0L,0L,1L,
+1L,1L,1L)
)
as.data.table(tdata)
tdata <- tdata %>% convert_as_factor(id, w, b)
as.data.table(tdata)
res.aov <- anova_test(
  data = tdata, dv = dv, wid = id,
  within = c(w), between = c(b), detailed = TRUE
  )
res.aov
get_anova_table(res.aov, correction = "none")

And after the code change, PDL::Stats now gives the same results as R! If you want to run the above yourself without "installing" any "software", you can enter it in this extremely handy service: https://rdrr.io/snippets/.

[etj]

And with much further effort and maths-figuring, and some experimenting, the rest of the problem is now fixed in:
• https://github.com/PDLPorters/PDL-Stats/commit/fc4609139d32a669566a46429fb93b62c6abf209
• https://github.com/PDLPorters/PDL-Stats/commit/7ddc45c643e0b210ce500a0dea8d08625814c9d7
• https://github.com/PDLPorters/PDL-Stats/commit/cee6320e224b66e668c38d3198a58ba8e1dce388
As noted on https://rt.cpan.org/Ticket/Display.html?id=97925, as soon as I've put together a demo of all the quite neat stats-plotting stuff (doing statistical analysis is really aided by looking at plots), a new version will go up on CPAN.