Как провести точный попарный тест Фишера? - PullRequest
0 голосов
/ 01 февраля 2019

Я пытаюсь сделать pairwise.fisher.exact тест, но у меня есть некоторые проблемы с ним.К сожалению, я не могу найти какие-либо решения в Интернете.Таким образом, я надеюсь, что кто-то может мне помочь.Это мои данные:

dput(test)
structure(list(yes = c(2L, 2L, 5L, 5L, 1L, 0L, 0L, 12L, 4L, 7L, 
4L, 1L, 6L, 25L, 29L, 66L, 18L, 77L, 53L, 19L, 13L, 15L, 7L, 
20L, 58L, 34L, 38L, 39L, 11L, 45L, 33L, 50L, 48L, 42L, 29L), 
    no = c(898L, 1438L, 1435L, 1435L, 1439L, 1440L, 1440L, 1428L, 
    1436L, 1433L, 1436L, 1439L, 1434L, 1415L, 1351L, 1374L, 1422L, 
    1363L, 1387L, 1421L, 1427L, 1425L, 1433L, 1420L, 1382L, 1406L, 
    1402L, 1401L, 1429L, 1395L, 1407L, 1390L, 1392L, 1398L, 907L
    ), sum = c(900L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 
    1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1380L, 1440L, 
    1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 
    1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 1440L, 
    936L)), class = "data.frame", row.names = c(NA, -35L))

Сначала я сделал fisher.test, чтобы увидеть, есть ли разница между "да" и "нет".fisher.test(test,simulate.p.value=TRUE) R дает мне существенную разницу, поэтому я хочу знать, какая разница между «да» и «нет».Я написал:

pairwise.fisher.test(c(2,2,5,5,1,0,0,12,4,7,4,1,6,25,29,66,18,77,53,19,13,15,7,20,58,34,38,39,11,45,33,50,48,42,29),c(900,1440,1440,1440,1440,1440,1440,1440,1440,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,936),p.adjust.method = "bonferroni") 

R показывает мне разницу между «да» и «нет», но только для первых 30 значений.Там написано:

"достигла getOption (" max.print ") - пропущено 5 строк"

Таким образом, последние 5 строк больше не отображаются.Я не знаю, как решить эту проблему.

1 Ответ

0 голосов
/ 02 февраля 2019

результат <- pairwise.fisher.test (c (2,2,5,5,1,0,0,12,4,7,4,1,6,25,29,66,18,77,53,19,13,15,7,20,58,34,38,39,11,45,33,50,48,42,29), c (900,1440,1440,1440,1440,1440,1440, 1440,1440,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400,1400, 1400,1400,936), p.adjust.method = "bonferroni") </p>

Просмотр (результат [["p.value"]])

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...