Dear colleagues with particular interests in r software
I have completed an ordinal regression model using two different r functions (polr - MASS package and clm - ordinal package). The output was measured using a 5 level Likert scale (0-4). Predictors are numerical as well as categorical. N=722 obs. Using both functions, results for coefficients values and p value are very similar, except for intercepts. In one case (polr) intercepts are significately different from zero. In the other case (clm) 3 out of 4 intercepts values are not different from zero. Difference come obvisouly from the calculation of standard error for intercepts (intercept values are the same in both of the models). Is anyone had a similar experience ? How to decide wich function is the most appropriate ?
Many thanks in advance
Dr Jeoffrey DEHEZ