在R中复制SAS GLIMMIX对比语句

Question

我正在与一位使用SAS的合作者一起工作，而我更熟悉R。我的合作者已经创建了一个GLIMMIX过程类型III测试，它具有自定义假设测试的固定效果(使用对比语句)。我已经找到了一种方法来重新创建固定效果的类型III测试，使用：

car::Anova(model, type = 3, test = "F"

我的问题是，有没有办法复制对比声明？具体来说，如果我的模型是：

glm(Response ~ Group, data = dataset, family = "binomial")

我得到了如上所述的组变量固定效应的III检验，有没有办法在组变量的水平之间进行测试？

SAS中的对比语句输出给出了一个任意标签，Num DF，Den DF，F value和Pr >F。我在R中寻找这种类型的输出，以便将每个实验组与控制组进行对比。

Answer 1

多亏了Ben Bolker的建议，我找到了可以工作的代码！

model <- glm(Response ~ Group, data = dataset, family = binomial)

emm.object <- emmeans(model, ~ Group)
joint_tests(emm.object) #Produces similar to type III test, but check
contrastlist <- list ("Control vs. Group 1" = c(1, -1, 0, 0), "Control vs. Group 2" = c(1, 0, -1, 0), "Control vs. Group 3" = c(1, 0, 0, -1))
tcontrast <- data.frame(contrast(emm.object, contrastlist))

F-values <- tcontrast$t.ratio^2

注意:有时候，对比会给你一个z.ratio而不是一个t.ratio。不确定这是否与我们的目的不同。F值就是https://stats.stackexchange.com/questions/55236/prove-f-test-is-equal-to-t-test-squared的t.ratio平方。

我找到了一个emmeans对比度的例子，上面提到的就是F值的t.ratio的平方，但我再也找不到合适的网页了。