сравнение 5 групп по качественному признаку |
Здравствуйте, гость ( Вход | Регистрация )
сравнение 5 групп по качественному признаку |
7.04.2018 - 18:41
Сообщение
#16
|
|
Группа: Пользователи Сообщений: 902 Регистрация: 23.08.2010 Пользователь №: 22694 |
Автор пакета "Exact" и здесь то же повторяет: https://stats.stackexchange.com/questions/1...hloo-or-barnard Так что это, в отличие от x13, не ошибка какая-то, но я пока не понимаю что это значит Поразбираться можно, но очень всё это ресурсозатратно. Видно проще хроники написать, чем нормальный учебник: https://errorstatistics.files.wordpress.com...ston-mielke.pdf Я ворчал только по поводу формулы. Вот оригинальная статья Boschloo (1970), но требуется регистрация: https://onlinelibrary.wiley.com/doi/pdf/10.....1970.tb00104.x Что касается "Хроник...", то - да, книженция любопытная, мне она попадалась в Google.books с очень малым числом страниц для просмотра. Сообщение отредактировал 100$ - 7.04.2018 - 18:43 |
|
7.04.2018 - 22:33
Сообщение
#17
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
|
|
7.04.2018 - 22:56
Сообщение
#18
|
|
Группа: Пользователи Сообщений: 902 Регистрация: 23.08.2010 Пользователь №: 22694 |
Да, спасибо, хотя Nokh уже угостил ссылкой в #15. Правда, жизнь и впрямь слишком коротка, чтобы осиливать подобные талмуды. По существу вопроса попутно могу сообщить, что в случае биномиального распределения при p1=p2=p Boschloo сконструирован так, что точной верхней гранью рассматриваемого критериального множества как раз и является результат точного метода Фишера (ТМФ). |
|
5.06.2018 - 11:45
Сообщение
#19
|
|
Группа: Пользователи Сообщений: 14 Регистрация: 12.03.2017 Пользователь №: 29479 |
Пока в R такого не делал, но подскажу. Когда я последний раз интересовался этим 3-5 лет назад, в R точным методом можно было обсчитать только таблицы 2 х 2, а пакета coin вроде не было. Я точным методом считал давно в StatXact, а последнее время довольствуюсь Монте-Карло в бесплатном PAST (https://folk.uio.no/ohammer/past/). Когда освоите, можете поделиться функциями здесь; в руководстве coin видел какой-то пример, но не было времени вникать. В пакете coin есть возможность провести тесты. Я надеюсь кому-нибудь пригодится. library(coin) tmp <- matrix( c(73, 23, 36, 86), nrow = 2, byrow = TRUE, dimnames = list( "Outcome" = c("Outcome 1", "Outcome 2"), "Factor" = c("Factor 1", "Factor 2"))) tmp <- as.table(tmp) # Permutation test, Monte Carlo test tmp <- chisq_test(tmp, distribution = approximate(B = 10000)) pvalue(tmp) midpvalue(tmp) # Exact permutation test tmp <- chisq_test(tmp, distribution = "exact") pvalue(tmp) midpvalue(tmp) |
|