Catanova for two-way contingency tables with ordinal variables using orthogonal polynomials