Извлечение необработанных p-значений из функции glm glht (вместо скорректированных Tukey p-значений) - PullRequest
0 голосов
/ 29 октября 2018

Мне дали приведенный ниже код и попросили извлечь необработанные p-значения, а не скорректированные по Тьюки значения (так как мы будем корректировать множественные сравнения с использованием Homes-Bonferroni на более позднем этапе), но я не уверен, что заменить "Tukey" на (я новичок в использовании R) .....

        res=glht(x, linfct=mcp(Letter="Tukey")
        out=summary(res)
        out

1 Ответ

0 голосов
/ 07 ноября 2018

Я нашел ответ. Для всех, кто заинтересован ... Опция "Tukey" для функции glht в пакете multcomp фактически не использует коррекцию Tukey, она просто устанавливает все парные сравнения. Это не делает p-значения; для этого вам нужно summary.glht. Чтобы получить необработанные значения p, вы используете test = Adjust («none»).

res=glht(x, linfct=mcp(Letter="Tukey")
   out=summary(res, test = adjusted("none"))
   out
...