Я нашел ответ. Для всех, кто заинтересован ...
Опция "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