我有一个我想进行的分析,我已经计划了对比(不是事后比较!)我想在治疗组之间进行。治疗组变量有 (k =) 4 个水平。我计划总共进行 3 次不同的比较,因此,如果我理解正确的话,我不需要对计算的 p 值进行任何调整,因为比较是 k-1。
我想在 R 中使用multcomp
orlsmeans
包来做到这一点。我的问题是:有谁知道是否可以在不对置信区间(和 p 值)进行任何调整的情况下进行这种计划比较?据我从我看过的小插图和我看过的例子中可以看出,该summary.glht()
功能会作为默认值进行调整,我不清楚哪个选项会撤消此操作。
如果有人需要一个可重现的例子,他们可以使用我在http://www.ats.ucla.edu/stat/r/faq/testing_contrasts.htm上找到的这个例子:
library(multcomp)
hsb <- read.csv("http://www.ats.ucla.edu/stat/data/hsb2.csv")
m1 <- lm(read ~ socst + factor(ses) * factor(female), data = hsb)
summary(m1)
K <- matrix(c(0, 0, 1, -1, 0, 0, 0), 1)
t <- glht(m1, linfct = K)
summary(t)