Краскел-Уоллис как регрессия, Все про 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]

Добрый вечер.
1. Однофакторный анализ (y~x). Соответственно можно думать о (y~x1*x2) и тп.
2. http://r.789695.n4.nabble.com/kruskal-wall...-td4288008.html - "The Kruskal-Wallis test is a special case of the proportional odds ordinal logistic model. You can get any contrast you want by testing regression coefficients. In a couple of weeks the rms package's contrast function will allow for individual confidence intervals of effects that together have a 0.05 type I error, by using the multcomp package (called automatically from contrast.rms)".