从glm-glht函数中提取原始p值(而不是Tukey调整的p值)



我得到了下面的代码,并被要求提取原始p值,而不是Tukey调整后的值(因为我们将在稍后阶段使用Homes Bonferroni进行多次比较调整(,但我不确定用什么来代替"Tukey"(我是使用R的新手(。。。。。

res=glht(x, linfct=mcp(Letter="Tukey")
out=summary(res)
out
我找到了答案。对于其他感兴趣的人。。。multcomp包中glht函数的"Tukey"选项实际上并没有使用Tukey校正,它只是设置了所有的成对比较。它不做p值;为此,您需要summary.glht。要获得原始p值,请使用test=adjusted("none"(。
res=glht(x, linfct=mcp(Letter="Tukey")
out=summary(res, test = adjusted("none"))
out

最新更新