 Краскел-Уоллис как регрессия, Все про R |
Здравствуйте, гость ( Вход | Регистрация )
 27.12.2015 - 11:28
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 107 Регистрация: 27.12.2015 Пользователь №: 27815  |
Всем добрый день.
Возник вопрос касательно порядковых данных и применения одно-двух-повторного анализа данных. В 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" предполагается)? |
 |
 
|
 |
|
 |
Ответов
|
23.01.2016 - 02:25
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 107 Регистрация: 27.12.2015 Пользователь №: 27815 |
Доброй ночи.
Если кому-то интересно. Написал Харреллу, сказал, что это ошибка http://stats.stackexchange.com/q/191063. Чтобы с Github'a загрузить фикс, использовать следующий код: Код require(rms)
getRs('contrast.s', grepo='rms', dir='R', put='source') |
|
|
|
|
|
|
23.01.2016 - 09:29
Сообщение
#3
|
|
 Группа: Пользователи Сообщений: 1091 Регистрация: 26.08.2010 Пользователь №: 22699 |
Цитата(comisora @ 23.01.2016 - 02:25)  Доброй ночи. Если кому-то интересно. Написал Харреллу, сказал, что это ошибка http://stats.stackexchange.com/q/191063. Чтобы с Github'a загрузить фикс, использовать следующий код: Код require(rms) getRs('contrast.s', grepo='rms', dir='R', put='source') И какую разницу показал расчет с помощью поправленного кода?  |
|
|
|
|
|
|
23.01.2016 - 14:11
Сообщение
#4
|
|
|
Группа: Пользователи Сообщений: 107 Регистрация: 27.12.2015 Пользователь №: 27815 |
Цитата(p2004r @ 23.01.2016 - 10:29) И какую разницу показал расчет с помощью поправленного кода? Добрый день. Разницы между orm и lrm нет. В частных случаях результаты функций совпадают с результами wilcox.test() и dunn.test(). |
|
|
|
|
|
Сообщений в этой теме
comisora Краскел-Уоллис как регрессия 27.12.2015 - 11:28
p2004r Цитата(comisora @ 27.12.2015 - 11:28... 27.12.2015 - 21:27 comisora Добрый вечер.
1. Однофакторный анализ (y~x). Соотв... 27.12.2015 - 21:54
p2004r Цитата(comisora @ 27.12.2015 - 21:54... 27.12.2015 - 22:55 comisora kruskal.test(y~x1*x2) не работает, а lrm(y~x1*x2) ... 27.12.2015 - 23:21 p2004r Цитата(comisora @ 27.12.2015 - 23:21... 28.12.2015 - 23:59 comisora Цитата(p2004r @ 28.12.2015 - 23:59) ... 29.12.2015 - 01:28 p2004r Цитата(comisora @ 29.12.2015 - 01:28... 29.12.2015 - 16:01 comisora Добрый день.
Ничего не мешает запустить. Выдаёт хи... 11.01.2016 - 19:08 nokh Добавлю свои 5 копеек.
Во-первых, dunn.test - это ... 11.01.2016 - 22:33 comisora Цитата(nokh @ 11.01.2016 - 22:33) Во... 12.01.2016 - 11:22 |