Оформление результатов стат. анализа, какие данные вводить в таблицу? |
Здравствуйте, гость ( Вход | Регистрация )
Оформление результатов стат. анализа, какие данные вводить в таблицу? |
12.03.2018 - 15:57
Сообщение
#1
|
|
Группа: Пользователи Сообщений: 1 Регистрация: 9.03.2018 Пользователь №: 31077 |
Добрый день.
Прошу помочь разобраться в таком вопросе. При анализе значимости с помощью непараметрического критерия Манна-Уитни. Какие результаты нужно представлять в таблице. Медиану и интерквартильный размах или среднее и ошибку среднего. По механизму самого анализа вроде медиана и интерквартильный размах. С уважением... |
|
30.05.2018 - 13:52
Сообщение
#2
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
https://stat.ethz.ch/pipermail/r-help/2012-...ary/303977.html
Тут все ответы, Разность медиан (median(A)-median(B)) и ДИ к разности медиан - раньше я использовала package pairwiceCI, а сейчас пакeт iNZightVIT (R должен быть на машине), который рисует график после мультика, сопровождающего bootstrap распределение медианы. Пришпиленную картинку я сделала по тестовому примеру, который дан в ссылке. Вы можете также взять этот же пример и проверить результат своего кода. ДИ median(A-B ) дает wilcox.test: > I=wilcox.test(A,B,conf.lev=0.95,conf.int=TRUE,exact=F,correct=T) > inf=I$conf.int[1] > sup=I$conf.int[2] > C<-(I) > C Тут все однозначно и сходится с ручным расчетом по методу Altman, который также описан в русскоязычной статье ИНТЕРВАЛЬНАЯ ОЦЕНКА МЕДИАНЫ И ЕЕ АВТОМАТИЗАЦИЯ |
|
5.06.2018 - 11:35
Сообщение
#3
|
|
Группа: Пользователи Сообщений: 14 Регистрация: 12.03.2017 Пользователь №: 29479 |
https://stat.ethz.ch/pipermail/r-help/2012-...ary/303977.html Тут все ответы, Разность медиан (median(A)-median(B)) и ДИ к разности медиан - раньше я использовала package pairwiceCI, а сейчас пакeт iNZightVIT (R должен быть на машине), который рисует график после мультика, сопровождающего bootstrap распределение медианы. Пришпиленную картинку я сделала по тестовому примеру, который дан в ссылке. Вы можете также взять этот же пример и проверить результат своего кода. ДИ median(A-B ) дает wilcox.test: > I=wilcox.test(A,B,conf.lev=0.95,conf.int=TRUE,exact=F,correct=T) > inf=I$conf.int[1] > sup=I$conf.int[2] > C<-(I) > C Тут все однозначно и сходится с ручным расчетом по методу Altman, который также описан в русскоязычной статье ИНТЕРВАЛЬНАЯ ОЦЕНКА МЕДИАНЫ И ЕЕ АВТОМАТИЗАЦИЯ Спасибо за ответ. Мой код выдал A =c(619, 600, 490, 1076, 654, 955, 563, 955, 827, 873, 1253) B =c(346, 507, 598, 228, 576, 338, 1153, 354, 560, 517, 381) A <- apply(matrix(sample(A, rep=TRUE, 5000*length(A)), nrow=5000), 1, median) B <- apply(matrix(sample(B, rep=TRUE, 5000*length(B)), nrow=5000), 1, median) quantile(A-B, c(.025, 0.975)) 2.5% 97.5% 59 601 Вроде не сильно отличается. |
|