Dear R users,

I have a dataset with two ordered variables, tr_x1 and tr_y1. A crosstable of them can bee seen below.

     tr_x1
tr_y1   -1    0    1
   -1  629  100  629
   0  1396 4353 1443
   1   668  126  655

It is clear that if tr_x1 is 0, it has an effect on tr_y1. A chi-square statistic is clearly showing this with a low p-value.
Is there a regression-based method you would offer? I tried polr from MASS package but without finding a significant coefficient, because the columns for tr_x1 and tr_y1 are similar.

Thank you for your help!

daniel

	[[alternative HTML version deleted]]

