The two chisq.test p values differ when the contingency table is transposed! (PR#3486)

I think it works to use

           STATISTIC <- sum(sort((x - E) ^ 2 / E, decreasing = TRUE))

instead: this starts by summing the big values, and hence if at all
slightly 'underestimates' the real value, which is fine for the
comparisons.