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

Доброй ночи.
Если кому-то интересно. Написал Харреллу, сказал, что это ошибка http://stats.stackexchange.com/q/191063.
Чтобы с Github'a загрузить фикс, использовать следующий код:

Код

require(rms)
getRs('contrast.s', grepo='rms', dir='R', put='source')

[p2004r]

Доброй ночи.
Если кому-то интересно. Написал Харреллу, сказал, что это ошибка http://stats.stackexchange.com/q/191063.
Чтобы с Github'a загрузить фикс, использовать следующий код:

Код

require(rms)
getRs('contrast.s', grepo='rms', dir='R', put='source')

И какую разницу показал расчет с помощью поправленного кода?