Краскел-Уоллис как регрессия, Все про R Опции

[comisora]

Всем добрый день.
Возник вопрос касательно порядковых данных и применения одно-двух-повторного анализа данных. В R есть пакет "dunn.test" с помощью которого можно провести одномерный вариант с попарными сравнениями (dunn.test(x,g, "BY") - обычный KW. Есть пакет rms с функцией orm, частным случаем которой должен быть KW для одного вектора с порядковыми данными с разделением по подгруппам. Результат функции он выдает в виде регрессии с коэффициентами, однако возникло желание привести результат к "привычному" виду. Возникла мысль использовать конструкцию contrasts(Glm(M), list(g="a"), list(g="b")) потом "a" и "c" и "b" и "c". Результат схож (но не идентичен) с результатом dunn.test, однако уверенности в том, что так можно делать нет.
1. Насколько правомерно так тестировать контрасты для данных, измеренных в порядковой шкале?
2. Можно ли таким способом тестировать двух-трех-многовходовые порядковые данные для вычленения парных различий (использование поправки "BY" предполагается)?

[comisora]

kruskal.test(y~x1*x2) не работает, а lrm(y~x1*x2) и orm(y~x1*x2) работают. Для одномерного случая результаты обеих функций неплохо бьются с результатами dunn.test. У lrm контрасты аналогичны попарным сравнениям, а у orm нет, хотя они обе выдают одно и тоже. Хочется разобраться ЧЯДНТ.