按照说明,我想出了以下内容:
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 文档的提示也被接受。