我已在我拥有的数据库中使用了lme()
功能。之后,我使用了glht()
来查看均值的多重比较。 glht()
的摘要在线性假设中给出了4列,其中最后一列(Pr(> | z |))是我目前感兴趣的那一列。
因为我在10个不同的小组中处理20个参数,4种不同的治疗方法,所以我们计划创建一个热图"在summary(glht())
的最后一列中获得的 p 值。
我最初的计划是将这些值复制并粘贴到Excel中,然后我就可以制作这个热图。但是,R中的summary()
将列分隔为空格,这对发送到Excel没有帮助。
我有2个问题,然后:(1)是否可以以这种方式更改summary()
所有值都以制表符分隔? (2)有没有人有更好/更快的方法来准备这张热图?对于最后一个问题,我已经考虑用R中的 p 值填充矩阵作为启动器。但是,我不确定如何从summary()
函数调用此Pr(> | z |)列。
非常感谢您的帮助!