2

按照说明,我想出了以下内容:

library(papaja)
library(datasets)
library(stats)
library(multcomp)


bla <- glm('MPG.highway ~ DriveTrain*Origin',data = Cars93)
contrast.matrix <- rbind('main' = c(0,1,0,0,0.5,0), 'int' = c(0,0,0,0,1,0))
blac <-  glht(bla, linfct = contrast.matrix)
apa_print.glht(blac,test = multcomp::adjusted())
apa_print.summary.glht(summary(blac,test=adjusted(type="bonferroni")))

这不起作用:“不能对不存在的列进行子集化。x 列p.value不存在。”

有更好的想法吗?或者由于该软件包是实验性的,它根本不起作用?关于如何将 glht 的输出打印到由 knitr 生成的 htlm 文档的提示也被接受。

4

2 回答 2

2

apa_print.glht()笔记的文档

这些方法没有经过适当的测试,应该被认为是实验性的。

似乎有一些上游更改导致此方法失败。对此感到抱歉,但现在应该已修复。请尝试安装最新的开发版本papaja

remotes::install_github("crsh/papaja@devel")

apa_print(blac,test = multcomp::adjusted())
于 2020-07-16T14:39:52.713 回答
-1

安装最后一个开发后。papaja 的版本,并使用“emmeans”对象运行 apa_print,在我的情况下也会发生相同的错误:

<错误/vctrs_error_subscript_oob>

错误:不能对不存在的列进行子集化。x 列p.value不存在。

于 2020-08-03T10:32:28.033 回答