Matrices, common in statistics because of applications of linear algebra, have long maintained the convention of i,j indexes starting with one.

It's not. Zero based is wise for offsets (CS) and fields where you have a zero element, and a few corner cases like factorial, but by and large any time you're enumerating an order over a nonempty set you start with one in math, because, max(indices) = |set|